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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

5.3.3.1 Properties <strong>of</strong> the EM algorithm<br />

The EM algorithm is certainly one <strong>of</strong> the most accepted algorithms to estimate finite<br />

mixture <strong>models</strong> due to some <strong>of</strong> its attractive properties listed below:<br />

• The most important advantage <strong>of</strong> the EM algorithm is surely its convergence<br />

towards the optimum parameter values. This means that, given the recent<br />

mixture model parameters, a single EM iteration provides new parameter<br />

estimates, which are proven to increase the loglikelihood <strong>of</strong> the model<br />

(Dempster et al., 1977; McLachlan et al., 1997). The convergence <strong>of</strong> the EM<br />

algorithm is proven by Meilijson (1989) and Wu (1983).<br />

• The EM algorithm ensured that the estimated parameters are within the required<br />

range (admissible range). This means that, <strong>for</strong> example <strong>for</strong> the Poisson<br />

distribution, the parameter values are zero or positive and cannot take negative<br />

values.<br />

• The EM algorithm is fairly easy to program.<br />

However, apart from these appealing properties <strong>of</strong> the EM algorithm, some limitations<br />

have been identified as well:<br />

• The setback with the EM estimation is that the procedure may converge to a<br />

local but not a global optimum (McLachlan et al., 1988; Titterington et al.,<br />

1985). It is generally accepted that the best way to avoid a local solution is to<br />

use multiple sets <strong>of</strong> starting values <strong>for</strong> the EM algorithm and to observe the<br />

evolution <strong>of</strong> final likelihood <strong>for</strong> the different restarts <strong>of</strong> the EM algorithm.<br />

84

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

Saved successfully!

Ooh no, something went wrong!