ぎょーむ日誌 2006-04-13
2006 年 04 月 13 日 (木)
-
0830 起床．
朝飯．
コーヒー．
洗濯．
1040 自宅発．
曇．
1055 研究室着．
-
本日からかとー室長が筑波大がよいなので，
A803B はさしずめ「かとー留守オフィス」とでも呼ばれるべき状態になった．
かとーさん荷物はまだいろいろとあるわけで
……
しかしすばやい大学院生がかとー高級椅子を「借りて」いってしまった．
-
苫小牧樹木直径モデル論文原稿にふたたび取り組んでみる．
起動状態よろしくない．
-
昼飯．
-
1328 地震．
また
震度 1
ぐらいなんだけど，
A 棟 8F はぐらぐらする．
-
論文原稿よみ．
ぢりぢりと鍋嶋わーるどにはいりこむ．
自分で書いた部分にさしかかると，
いつものことながらアタマいたくなる．
-
この推定計算のわくぐみでは「樹種間比較」はかなり難しいわけだが
……
ふと「全樹種まったく同じパラメーター」
という再単純モデルなら簡単に計算できることに気づく．
たぶんすごくあてはまりの悪いんだろうけど
……
というのも，
かなり樹種差があるのに「差がない」と仮定したモデルなんで．
-
A801 室「苫小牧ボス」こと鍋嶋さんに拝謁し，
かかる計算の勅許をいただく．
Dell desktop 機 (P4 2.4GHz) で 25 分弱で計算終了．
やはり「樹種べつ」モデルの合計 AIC = 5470 にくらべて
かなりよろしくないモデルだ (AIC = 5535)．
-
しかしまあ，
いちおうこれで「樹種間比較」最初の一歩にはふみこんでいることになる．
これ以上のまぢめな計算は難しい．
それでも論文原稿修正の，
多少の役にはたつかもしれぬ．
-
M1 大学院生が突然かとー留守オフィスにやってくる．
「昨日のぎょーむ日誌に関する質問あるのですが」
「はい ……」
(いきなりなにか花粉流動まにあっくな質問だろうか……)
「以前から疑問に思ってたのですが」
「はい ……」
(どきどき……)
「『今日の食卓』の『電磁波酒蒸し』って何ですか，あれは?」
かなりずっこけつつ電子れんぢすなわち microwave oven
でつくる酒蒸しであります，
と回答さしあげる．
-
しかしその後はあまり進捗もなく撤退．
1910 研究室発．
1930 帰宅．
晩飯．
-
昨日は敵前逃亡したこんさるメイル書き．
二通だけなんだが
……
やはり時間かかる．
-
さて後藤さんと昨日検討していたハナシの一部をこのへんに書いてみる．
Burczyk 的モデルとはこういったかんぢで
(自殖ナシの場合):
が遺伝子型
h
の母親から
g
という種子ができる条件つき確率，
それが系外花粉による場合
と父親
j
の花粉による場合
;
が (あれこれとアヤしまれてる)
「系外からの花粉である確率」，
そして (系内限定の) 父親間競争は
と多項 logistic 的になっている．
は母親
i
における父親
j
の「競争能力」で樹間距離依存性
(これが dispersal kernel できまる)
とか父親のサイズとか開花状況の関数である．
-
で，
後藤さんたちの新モデルは
が 花粉 dispersal kernel
のカタチと母親樹の plot 内位置で決まるんぢゃないの，
という変更をいれたもの
……
私はその確率計算についてはあまり正しいと思っていない．
-
で，
こういう系外・系内分割のやりかたがそもそもおかしいのでは，
という意見もありえるわけで．
どこからきた花粉であろうと同じように柱頭の上で競争させればいいはず，
ということで
……
たとえば代替案としては，
とでもしてみるか?
は系外花粉の競争能力で，
ここをどう定式化するかがとてもむずかしいんだけど
(つまり自分でもどうしたらいいかわかってないんだけど)，
-
花粉 dispersal kernel のカタチ
-
母親樹の plot 内位置
-
系外樹木密度とサイズ分布
みたいなものどもに依存している，
と．
-
さらに花粉親樹木あるいは相性問題を random effects でくみこむ，
とか，
は系外遺伝子頻度に依存しているんだけどこれを
EM アルゴリズムとか Bayes 推定的に構成しながら計算する，
とか，
このあたりのモデリングはまだまだ考えるべきことがあると思うんだよね．
嗚呼，
難問山積
……
-
ありがたいことに
富田君
が重要な難問に関する研究を一歩ススめてくれた
……
Burczyk は
「ぶぁーちっく」と読むのが適当
ではなかろーか，
という考究である．
なるほど，
「いんちきドイツ語ふう読み」
とでもいうべき久保説「ぶるくつぃく」
よりはもっともらしいような気がする
……
たしかに
Wikipedia でも
``cz''は「ちゃ」音，
とあるよな．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0920):
米麦 0.6 合．
海藻スープ．
サバ電磁波酒蒸し．
- 昼 (1310):
研究室お茶部屋．
米麦 0.6 合．
レトルトパウチドのトマトソース．
- 晩 (2030):
米麦 0.7 合．
タマネギ・豆腐・タラすりみカレー．
