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)