ぎょーむ日誌 2009-05-04
2009 年 05 月 04 日 (月)
-
0840 起床．
ねむい．
朝飯．
コーヒー．
1000 自宅発．
晴．
1015 研究室着．
-
さて，
昨晩やらせておいたウミガメ時系列データ解析の
「脱のいず」改良版モデルの計算だが
……
6900 秒ほどで終了．
なにしろ 5 万 MCMC step とかだからなぁ
……
ちなみにさきの盛岡大会
ポスター発表
の計算は 2000 MCMC step でしかもそのうち半分は burn-in
といったありさまで．
つまり，
「のいず」まみれな推定結果
(じつはこれは見てくれが悪いだけで，
言えることは「脱のいず」版とほぼ同じ)
でよければ数分で計算終了，
となるわけで．
-
で，
時間のかかった MCMC 計算だけど，
部分的に収束がいまいち．
そこで，
「脱のいず」
のための微小のいずの SD を
sqrt(1 / 1.0E+5)
から
sqrt(1 / 5.0E+4)
に変更して再計算を命じてみる．
また 2 時間ばかりがんばってください．
-
R 作図わざ本の翻訳原稿の点検作業
……
があまりススまぬうちにいったん帰宅して昼飯．
北海道もいよいよ春，
とゆーかんぢで．
-
昼飯からもどると長閑なるウミガメ計算 6800 秒が終了してたわけだが
……
やはり 5 万ではびみょーに足りない．
ということで，
9 万 MCMC step (予定計算時間は 3 時間弱)
に変更してまた計算させる．
-
また R 作図わざ本の作業．
-
1510 第 5 章終了．
残り (原書で) 90 ペイジ弱．
もうちょい続ける．
-
一時間ほど作業．
R 作図わざ本の翻訳原稿チェック，
本日は 13 ペイジほどススんだ．
-
「ちょっとのあいだ放置」
(うわー，二週間も経過!)
してた苫小牧シュート論文の作業のつづき．
まづは，
使っていると不幸きわまりない気分になる MS わーど上において，
ファイルの統合作業とか．
-
ウミガメモデル改良版計算，
予想よりもやや長い計算時間，
11800 秒 (3.2 時間)
を費やして終了した．
9 万 MCMC step の 2000 MCMC step とばし，
といった私としてはやったことがないほどのちょー長時間サンプリング．
さすがにぴしっと収束している
(BUGS code，
事後分布表)．
-
r[t]
の推定なんかも，
以前のモデルのもの
とはちょっと「かんぢ」がちがうなぁ
……
-
しかしまあ，
面白いのは最近 10 年間の
r[t]
経路幾何平均の分布なんかは，
以前とはあまり変わらないわけで．
-
「予測」できることも，
以前とはあまり変わらない．
つまりまあ，
この観測データに対して，
どこかに random walk 的なものが入った統計モデルで
抽出しうる情報なんてそんなには多くない，
ということだろう
……
モデル中「海岸ノイズ」とよばれてる部分がすごく大きいわけだしね．
ウミガメ上陸数の予測不可能性，
ともうしますか．
-
とりあえずウミガメ再計算の試行錯誤，
結果あれこれをきちんと保存して，
作業はしばらく休止．
-
苫小牧シュート論文，
あまりススまなかった
……
1900 研究室発．
1915 帰宅．
2050 自宅発北大構内走．
2120 帰宅．
ちょっとうんどう．
体重 74.0 kg
……
うーむ，
体重変動が大きい．
晩飯の準備．
晩飯．
-
[今日の運動]
-
腹筋運動 30 ×
3 回．
背筋運動 30 ×
3 回．
腕立ふせ 10 ×
3 回．
北大構内走 30 分間．
ストレッチ．
-
[今日の食卓]
- 朝 (0920):
食パン．
キャベツ・ミズナ・ネギ・ハムのサラダ．
- 昼 (1240):
スパゲッティー．
- 晩 (2210):
米麦 0.8 合．
ダイコン・ホッケの煮物．
キャベツ・ミズナ・ネギ・ハムのサラダ．
ヨーグルト．
