ぎょーむ日誌 2005-07-25

2005 年 07 月 25 日 (月)

• 0750 起床． 朝飯． コーヒー． 0930 自宅発． 晴． 0945 研究室着．
• アカマツ原稿修正に着手．
• すると D1 井田君がやってきてウコンウツギの訪花データ解析こんさる． ふーむ， 場所と日づけの二重 nested か …… まぁ， 実際の数値計算においては， どんな多変量確率分布も「離散化」すりゃー 一次元わーるどに持ちこめるけどね …… と書いてみて， ようやくにして， われながらボケてると気づいた．
  • 日づけ 1
  • 場所 A
  • 場所 B
  • 日づけ 2
  • 場所 A
  • 場所 B
  と階層的になってる (つまり nest してる) 場合には， 「一次元化」などできやしない (実際には日づけ 2 では場所 C と D なんだが)． 兇悪というしかないちょー難問だ …… ということで， 昨日かんがえてたよーな「尤度の期待値の対数」 vs 「対数尤度の期待値」 といった問題を検討する必要とかもでるわけですよ (いや， そんなの考えてもたいして改善されないか?)．
• しかしここまで来てしまうと， もういっそのことべいづで， という心境にもなりますな．
• またアカマツ原稿の修正．
• 昼飯．
• D1 井田君のウコンウツギのハナシのつづき …… 午前中， こいつは「二重混合モデル」か， と危惧してたんだけど意外な方面から モデルを簡単化する材料が提供されてしまった． 観測の交絡 (confounding) である． 井田君の調査方法では観測日が変わると観測対象の樹木個体がかわってしまう． そして観測対象となる樹木個体は一日に 1-2 本だ． これは何を意味するのか? すなわち 「樹木個体の random effects」 と 「日づけの random effects」 の区別がつかん， ということである．
• 一般に「交絡」は実験・観測において忌避されるべきもの， とされるわけだが …… われわれはこの「区別がつかん」という実験計画法的欠陥を逆に 悪用 してみせることにした． つまり混合モデルにいれるべき random effects はただひとつで OK だと正当化可能になり， それは 「樹木個体・日づけの random effects」 と呼ばれるものになる． これで計算が劇的に簡単になる． ははは
• で， あとは詳細をツメるばかり．
  • 株訪問回数の統計モデル: これは混合ポアソン分布でうまく説明できそうだ． しらべてみると， 各反復での株訪問回数の平均と分散はほぼ比例関係にあるんで． これは井田君担当． R の glmmML とか使えば計算できることなんで．
  • 株内の連続訪花数: これは平均の二乗と分散が比例関係にある， とわかった …… これは何を意味するのかといえば， 回数は幾何分布にしたがうことになる． 幾何分布 とは何かを考えてみれば， この確率分布が 株内のハチの行動を説明しうる簡単で合理的なモデル になってることがわかる． さらに， 「連続訪花回数はゼロ」 というのはありえないので幾何分布から回数ゼロをぬいた 条件つき確率分布を考えないといけない． さらに random effects も必要である． ということで， 「混合」「条件つき」幾何分布の最尤推定となる． しょうがないんでこれは久保担当．
  いやー， 井田君これを進化学会大会とやらで発表するつもりらしくて， ですね． まにあうのか?
• 1600 また別の「発表せとぎわ」なこんさる． こちらはより深刻で IMC9 でエゾシカ発表をする フィールド科学センター M2 の伊藤君である． といってもこの問題は (残り時間を考えると) すでに統計モデリングは終了してるんだが …… 何やら「割り算」思考からのがれられんよーで． ということで， 先日の負の二項分布モデルの解説のやりなおし，
• おーぷん大学雑用． みなさんちょっとした図を作るのになぜ PoorPoint なんぞ使うのか． そしてそれをそのまま当方に送りつけるのか， OpenOffice.org 2.0beta の解毒能力にだって限界はあるんだけど (とくに Mac 版のおひすファイルは苦しい) …… といぶかしみつつ， とりあえず ネット版 研究紹介を更新していく．
• 1825 研究室発． 1840 帰宅． 通販で買ってしまった自宅作業用の椅子がとどく． 組み立て作業．
• 作業修了後いきなりばててしまった． 夜中におきて， おーぷん大学雑用のつづき． 晩飯．
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0840): 米麦 0.7 合． キャベツ・モヤシ・ショウガ・ホタテの炒めもの．
  • 昼 (1310): 研究室お茶部屋． 米麦 0.5 合． キャベツ・モヤシ・ショウガ・ホタテの炒めもの．
  • 晩 (2420): 米麦 0.7 合． キャベツ・モヤシ・ショウガ・ホタテの炒めもの．