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 6.8: Parameter estimates (bootstrapped standard errors) <strong>of</strong> the four component<br />

restricted covariance model<br />

Component<br />

θ<br />

1<br />

θ<br />

2<br />

θ<br />

3<br />

θ<br />

12<br />

θ<br />

13<br />

θ<br />

23<br />

p<br />

j<br />

1 6.4384 0.0152 8.5477 0.0000 0.0000 2.4015 0.0143<br />

(0.3673) (0.0318) (0.3180) (0.0000) (0.0000) (0.1516)<br />

2 0.8485 0.1696 13.5921 0.0000 0.0000 0.0000 0.1213<br />

(0.0277) (0.0070) (0.0778) (0.0000) (0.0000) (0.0000)<br />

3 1.9083 0.4127 2.8167 0.0000 0.0000 0.0000 0.3575<br />

(0.0376)<br />

4 0.8075<br />

(0.0249)<br />

(0.0055)<br />

0.1545<br />

(0.0045)<br />

(0.0285)<br />

0.0819<br />

(0.0049)<br />

(0.0000)<br />

0.0000<br />

(0.0000)<br />

(0.0000)<br />

0.0000<br />

(0.0000)<br />

(0.0000)<br />

0.0000<br />

(0.0000)<br />

0.5069<br />

For the restricted covariance model it is also observed that the components <strong>of</strong> the model<br />

with small mixing proportions have the large standard errors.<br />

6.4.2 Results <strong>for</strong> the different <strong>multivariate</strong> Poisson <strong>hidden</strong> Markov <strong>models</strong><br />

Similar to the <strong>multivariate</strong> Poisson finite mixture <strong>models</strong>, <strong>for</strong> the <strong>multivariate</strong> Poisson<br />

<strong>hidden</strong> Markov <strong>models</strong> all three <strong>models</strong>, that is, the local independence model, the<br />

common covariance model, and the model with restricted covariance structure were<br />

fitted sequentially <strong>for</strong> 1 to 7 components ( k =1,…,7). Furthermore, in order to overcome<br />

the well-known drawback <strong>of</strong> the EM algorithm, i.e. the dependence on the initial<br />

starting values <strong>for</strong> the model parameters, 10 different sets <strong>of</strong> starting values were chosen<br />

at random. In fact, the transition probabilities ( P ij<br />

) were uni<strong>for</strong>m random numbers with<br />

constraint ∑<br />

m<br />

P ij<br />

j=<br />

1<br />

= 1 , 1 ≤ i ≤ m . The λ ’s were generated from a uni<strong>for</strong>m distribution<br />

123

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

Saved successfully!

Ooh no, something went wrong!