KuboLog 2007-07-08

苦情・お叱りは，たいへんお手数かけて恐縮ですが，久保 (kubo@ees.hokudai.ac.jp) までお知らせください．

2007 年 07 月 08 日 (日)

0900 起床．コーヒー．朝飯． 1055 自宅発．快晴． 1115 研究室着．

岩倉さんからアリ拉致実験についていろいろとメイルで教えていただいた．「たぶん，こうだろう」という憶測で作ったモデルだけど，「ほぼ正しい」とわかった．

[行動連鎖の仮定]

まず antennation する・しないがあって，次に攻撃する・しない，と仮定しているわけだが …… 今さらながら，この仮定が正しいかどうかは元データをあたればよい，と気づいた．で grep してみると ……

$ lv allinone.txt | grep " , a" 1 a a "b , a" a i jozankei jozankei30m o

…… また定山渓アリ の凶行か．全 1650 接触のうち Jozankei vs. Jozankei の一ヶ所で「攻撃 → antennation」があるな．この場合は地元アリが (30 m 離れた地点で拉致された) 拉致られアリに (瞬間的な判断で) 噛みついて，さらに antennation しているのである．「また定山渓」というのは，やはり全接触のうち一ヶ所だけ「拉致られアリが地元アリ (Hoshioki) に逆襲」という event があり，その「逆襲した拉致られアリ」が定山渓出身なのである (ぎょーむ日誌 2007-06-11)．

まあ，これはまたあとでデータ変換からやりなおして，階層ベイズモデルの MCMC 計算もやりなおしておくことにしよう．どうせここ一ヶ所かえても結果はほとんど変わらないし (なにしろ Jozankei vs. Jozankei 状況における event だし)， R によって 作図・作表は自動化 (投稿可能しあがり水準) されてるので何べんで推定計算やりなおしもどうってことない．

といった返信メイルかいて昼飯．

で，昼飯中に推定の再計算は終了 …… したと思ったけど，やはり「30 m 離れたところで拉致」「Jozankei vs. Jozankei」の交絡は依然として「スキあらば収束悪化させる」要因だとわかった．

で，事前分布を工夫して，いわゆる「切片」こと定数項パラメーターではない fixed effects パラメーターは無情報ではなく超事前分布に分散を制約される (IID) 事前分布をもつパラメーター，としてみる …… これで収束はうまくいき，以前と同じ結果になった (事後分布 summary table)．そして全ての図表の自動再生成，と．

といった試行錯誤の計算まち時間にも多巣アリ総説論文をちらちらと読みつつ．なんか，アルゼンチンアリに関しては (正しいのかどうかよくわからんけど) 遺伝的多様性でハナシをまとめようとしてるね．だとしたら，遺伝的なんちゃらでは場所による挙動のちがい説明できそうにないエゾアカヤマアリはおもしろそう，ということなのか?

1650 研究室発． 1710 帰宅．晩飯の準備．晩飯．

1925 より，札幌駅すぐ北の小映画館蠍座で「フラガール」みた．前回みた「インディアン」をいかにもな「男の子」映画だとすると，こちらは由緒正しい「女の子 (集団)」映画，とするべきか．好奇心旺盛で姿勢ただしく柔軟性と実行力にとみ理性的だけど情理をわきまえ思いやりのある「女の子」たち (中には高齢の「女の子」も含まれる) が果敢に問題にとりくみ，悲しいことがあって泣き，みんなで助けあい，みんなどんどんキレイになっていって，問題を解決して嬉しくて泣く映画 …… まあ，今回もまた正調なるエンターテイメントです．男どもはひきたて役というか道化役で，立ちションして水かけられたり自分でハサんでしまったり，ハワイ樹木ふぇちになってヤシに抱きついたり，硬直した議論ふりまわして周囲を困惑させてみたり，せいぜい借金とり相手にツルハシふりまわして借用書を食ってしまう，ぐらいの役まわりで．
(後記) 書き忘れてたけど，この映画みると平八郎提督の建白による北海道ハワイ化計画って実現するためにはそうとうの気合いが必要だとしみじみとわかった．

[今日の運動]

うんどう休養日

[今日の食卓]

朝 (0940): イソップベイカリーの巨大クロワッサン．
昼 (1330): 研究室お茶部屋．米麦 0.5 合．フリーズドライ春雨スープ．
晩 (1840): 米麦 0.6 合．麻婆豆腐，ためしてガッテンふう．キャベツ・ニンジン・タマネギ・キュウリ・レタスのサラダ．

ぎょーむ日誌 2007-07-08

2007 年 07 月 08 日 (日)