二階堂愛 (Itoshi NIKAIDO) が綴る Hacking is beliveing な日常_ [R]Rでグラフをアニメーションさせてみたよ R News 2008年2号 (PDF) や useR!2008での面白そうな発表をいくつか上げてみたで紹介されていたR の animation library を使ってみたよ。 まずは、非線形モデル y = mesor + a * cos(2*pi*(t-acrophase)/P) にN(0,0.5)のガウシアンノイズを加えたデータを用意し、これを測定されたデータと考える。そのデータに対して非線形回帰するために、Nelder-Mead法で目的関数sum( (y-yhat)^2 )を最大化する。いわゆる cosinor analysis ですね。その計算過程を各ステップごとにプロットし、animation libraryを使ってアニメ化する