ぎょーむ日誌 2009-05-02

2009 年 05 月 02 日 (土)

• 0800 起床． 朝飯． コーヒー． 怠業．
• 1100 自宅発． 晴． 1115 研究室着．
• ちょっとまた Web Site 生成 Perl すくりぷとまわりの改良に トラップされてしまった ……
• 昼飯調達のため， 北大構内をぷらぷら． よーやく当地の植物も開葉・開花しはじめた． 観光客多数．
• 歩いてるうちになぜかしら， ウミガメ モデルについて検討してしまった …… このモデルでは MCMC 計算 (ちょー) 高速化のため， あちこちに正規乱数が入っている． これで計算は格段に速くなるのだが， 推定された時系列あれこれはけっこうがたがた …… アンサンブル平均なんかはそうでもないけど．
• で， 歩いてるうちに思いついたのは， 「MCMC の各ステップでアンサンブル平均のようなものをとればいいのでは?」 といったもの． これなら計算速度はそんなに変わらないだろう ……
• で， 研究室にもどって， BUGS code をこのようにしてみたのだが ……

  for (i in 1:N.sample) {
  	x[LV[i]] ~ dnorm(mean.x[i], tau[2])
  	mean.x[i] <- x[LV[i] - 1] + r[i]
  	r[i] <- mean(rs[LV[i],])
  	for (j in 1:N.rs) {
  		rs[LV[i], j] ~ dnorm(rs[LV[i] - 1, j], tau[3])
  	}
  }
  

  やはり， そんなにムシのいいハナシはない， ともうしますか …… もくろみどーり， 計算時間は 300 秒ほどで終了したので (つまり以前とほぼ同じ), しめしめと思ったのだが …… 推定結果のほうも以前と同じ， つまり， 下のように結果として得られた事後分布も 以前 と同じでほとんど変わらなかった． いやはや．

• おとなしく， R 作図わざ本の翻訳原稿みなおし作業のつづき．
• A802 院生部屋のプリンター， CANON LBP5400 の K トナーぎれ． 研究室雑用 Wiki によると， 3 月に交換したばかりなのだが? …… そして予備トナーカートリッジがないので， 連休あけまで運転停止．
• R 作図わざ本， 今日も 11 ペイジほどススんだ． 1940 研究室発． 2000 帰宅． 晩飯の準備． ちょっとうんどう． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 背筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回． ストレッチ．
• [今日の食卓]
  • 朝 (0930): イソップベイカリーのクロワッサン．
  • 昼 (1400): スパゲッティー． リンゴ．
  • 晩 (2130): 米麦 0.8 合． ダイコン・ニンジン・豆腐・鶏肉のスープ． キャベツ・ミズナ・ネギ・ハムのサラダ． ヨーグルト．