File talk:Em old faithful.gif

The's a easier version, you don't have to do image conversion manually:

library(animation)

saveMovie({

 ani.options(autobrowse=FALSE)
 ## plot initial contours
 iter <- 1
 plot.em(theta)
 ## run EM and plot
 for (iter in 2:30){
   T <- E.step(theta)
   theta <- M.step(T)
   plot.em(theta)
   ani.pause()
 }

}, movie.name="em_faithful.gif", outdir=getwd())

--Gong-Yi Liao 21:07, 19 January 2012 (UTC)