EM algorithm

Initialize EM algorithm

Compute EM trajectory

Visualize EM solution

Initial model

EMAllPlot[initialModel, gdata, pr] ;

[Graphics:../HTMLFiles/index_16.gif]

Final model

EMAllPlot[em // Last, gdata, pr] ;

[Graphics:../HTMLFiles/index_18.gif]

Plot log-likelihood of data given the model as a function of EM iteration

PlotLogLikelihood[em, xd] ;

[Graphics:../HTMLFiles/index_20.gif]

Save EM data

Animation

EMAllPlot[#, gdata, pr] & /@ em2 ;

[Graphics:../HTMLFiles/index_30.gif]


Created by Mathematica  (September 8, 2003)