ぎょーむ日誌 2007-06-23

2007 年 06 月 23 日 (土)

• 0900 起床． コーヒー． 朝飯． 怠業．
• 昼飯． 洗濯． 怠業． アリばてか?
• 1515 自宅発． 曇． 1530 研究室着．
• また岩倉さん修論の石狩浜アリ世界の深淵をさまよう． ここは慎重に行動すべき作業フェイズなので， 時間かかるのはしょうがない．
• 岩倉さん修論読みつつ， 今までの階層ベイズモデリング反芻してるうちに， 「推定結果の説明を簡単にするため」 その他の目的のため， もう一個よけーなパラメーターを追加したほうがよいかも とゆー気がしてきた …… むろん， こんなことをすると計算時間がそのぶんよけーに必要になり， しかも他のパラメーターの収束わるくする危険性もあるんだけど．
• まあ石狩浜右往左往してるついでの副業ということで， と休眠状態の A801 院生部屋の Dell desktop 機をたたきおこし， 最新の formica.yessensis セットを rsync で送りこんで BUGS プログラムと R2WinBUGS なドライヴァーを改造． 再計算を命じる．
• MCMC 計算まち時間 (1000 秒ちょい) にも岩倉さん修論の解読は進捗していく …… methods に関しては小うるさそうな大統領ちぇっくが入ってないか， あるいは方法論節の重要性に関して ……
• で， 再計算で得られた事後分布だけど， やはり予想どーり攻撃確率に関するパラメーターの 収束わるくなった． この「あちこち交絡しまくり実験計画」のすばらしさのあらわれ， とでも申しましょうか ……
• そうだなあ， 「策」としては …… 無情報事前分布の「無情報ぶり」をやや低下させてみるか? これは fixed effects の「過小評価」につながるので， それほど unfair な策ではない． むしろ (たとえば「ゆーい差」決戦主義なヒトたちから見れば) 自分で自分のクビをしめてるようなものだろう． しかしながらおそらく， 事前分布ニラみつつ観測データに「あわせて」いく MCMC 計算によって生成される事後分布というものは， 私の期待としては， いま少し ``robust'' なものなのである．
• また 1000 秒ほど岩倉さん修論に専念． Methods section の書きなおしはあとまわしだ． Results section の解読にとりかかる．
• ネラいどーり全パラメーターで無事に収束した． そして予想どおりの推定結果というか， 「30 m 離れたところで拉致られたアリ」も またヨソものとしてちくちく攻撃されてる (ただし supercolony 内のアリは同種内平和ぼけなので 攻撃確率 0.01 - 0.1 ぐらい)， と判明した (print(ant.bugs, digits.summary = 2))． よしよし．

• そして R & LaTeX の自動作図・作表のすばらしさは， かかる再計算時の図表つくりなおしが make いっぱつで終了することである …… ただし， 実際には推定結果が変わると図のレイアウトに影響する場合があり， 作図プログラムの微調整とかは必要になる．
• ついでに， これら図表でもちいられている術語のいくつかがあまりよろしくない， と岩倉さん修論よんで気づいたので修正しておく．
• しかしまあ， 今回の階層ベイズモデルで推定してみせた攻撃確率なるものも， 自分で計算しといてなんだけど， ワケわからんものだなぁ． 1 contact における攻撃確率とやらが 0.05 だったとしよう． これは「低い」のか「高い」のか? 拉致られアリとしては各回の遭遇において 無傷のうちに脱出できる見こみは低くない …… しかしある nest の防衛圏から脱出するまでには多数の地元アリに 遭遇してしまうわけで， そうだとすると噛まれたりせずに防衛圏から脱出できる確率は かなり小さくなってしまう (ましてやうっかり地元巣に入ってしまったら ……)， と．
• 「拉致られアリに対する遭遇ごとの攻撃確率 0.05 でした」 といった「一回ごとの確率の低さ」によって 「すーぱーころにー」なるものは特徴づけられているんだろうけど …… 識別能力が高めの地元アリが 1-2 割いるだけでその巣の総合的な 他者排除能力ってのは高くなるような気がするんだけどなあ．
• よくわからぬまま撤退． 2115 研究室発． 2130 帰宅． 体重 67.6kg． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0940): イソップベイカリーのクロワッサン．
  • 昼 (1300): スパゲッティー． レトルトパウチドのカルボナーラソース．
  • 晩 (2330): 米麦 0.7 合． レタス・ネギ・卵炒飯． トマト．