02.09.2014 Views

multivariate poisson hidden markov models for analysis of spatial ...

multivariate poisson hidden markov models for analysis of spatial ...

multivariate poisson hidden markov models for analysis of spatial ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Table 7.6 Bacterial counts by 3 samplers in 50 sterile locations<br />

X<br />

1<br />

X<br />

2<br />

X<br />

3<br />

X<br />

1<br />

X<br />

2<br />

X<br />

3<br />

X1<br />

X<br />

2<br />

X<br />

3<br />

X1<br />

X<br />

2<br />

X<br />

3<br />

X<br />

1<br />

X<br />

2<br />

X<br />

3<br />

1 2 11 3 6 6 3 8 2 7 10 5 22 9 6<br />

8 6 0 3 9 14 1 1 30 2 2 8 5 2 4<br />

2 13 5 4 2 25 4 5 15 3 15 3 2 0 6<br />

2 8 1 9 7 3 7 6 3 1 8 2 2 1 1<br />

5 6 5 5 4 8 8 10 4 4 6 0 4 6 4<br />

14 1 7 4 4 7 3 2 10 8 7 3 4 9 2<br />

3 9 2 7 3 2 6 8 5 6 6 6 8 4 6<br />

7 6 8 1 14 6 2 3 10 4 14 7 3 10 6<br />

3 4 12 2 13 0 1 7 3 3 3 14 4 7 10<br />

1 9 7 14 9 5 2 9 12 6 8 3 2 4 6<br />

The observed covariance matrix and the observed correlation matrix are given below:<br />

⎡15.0714<br />

⎢<br />

0.2755<br />

⎢<br />

⎢⎣<br />

− 3.6939<br />

0.2755<br />

13.6428<br />

− 7.7347<br />

− 3.6939⎤<br />

− 7.7347<br />

⎥<br />

⎥<br />

32.6122 ⎥⎦<br />

⎡1<br />

⎢<br />

⎢<br />

⎢⎣<br />

0.0192<br />

1<br />

− 0.1666⎤<br />

− 0.3667<br />

⎥<br />

⎥<br />

.<br />

1 ⎥⎦<br />

Since the correlations were not very high, only the independence covariance structure is<br />

considered and the loglikelihood and the AIC values are recorded. According to the<br />

AIC criterion (Table 7.7), the model with the highest AIC value is selected<br />

(loglikelihood <strong>of</strong> –382.2335). This model, the seven-component local independence<br />

<strong>multivariate</strong> Poisson finite mixture model, gives a reasonable fit with respect to the<br />

loglikelihood.<br />

153

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!