ぎょーむ日誌 2005-11-05
2005 年 11 月 05 日 (土)
-
0740 起床．
朝飯．
コーヒー．
-
屋久島照葉樹葉っぱ問題にとりくむ．
統計モデルの修繕，
というアレ．
昨日
「1 年めに脱落した確率○，
2 年めに脱落した確率○，
3 年めに脱落した確率○，
……」
といった場合わけ尤度計算をくみこんだ統計モデルをつくり，
とか書いたけどこれはまちがいで，
もっと簡単だとわかった
(昨晩のうちに)．
-
つまり齡 a -1 から a のあいだの生存確率を
p(a)
とすると，
齡 a まで生き残る確率は Π p(a) ．
すると齡 a までのどこかの時点で死んでしまう確率は
その余事象だから単純に 1 - Π p(a) ．
つまり二項分布にすぎない．
-
しかしながら
p(a)
は logit でもなければ probit でもないので
R
の glm()
は使えない．
optim()
を駆使して最尤推定する必要あり．
-
で，
ぢたばたしてるうちになんとか推定計算・作図ができるようになった．
うーむ．
しかし何かずれてるような
……
-
昼飯をはさんで計算チェックなど．
問題なさそう．
1400 自宅発．
晴．
1415 研究室着．
-
上の図でなんだかずれてるみたいだ，
となってた原因わかった．
(モデルが良いかどうかは別として)
推定計算 (推定結果で線で表示) 自体には単純なまちがいはなかった．
まずかったのは，
観測値 (点々表示) は年枝ごとにぷーるした腐れ割算値である．
-
いやはや，
またまた割算もしくは平均値の呪い
とでもいうべき状況だ．
なぜこういうでたらめ統計量がうまくいかないのかというと，
個体差
が強すぎるためである．
下の図のごとし．
-
つまり同一樹種で同一の場所 (明るい・暗い) であっても，
枝によっては早々に葉っぱを脱落させるやつもいれば，
いつまでたっても葉っぱを落とさないのもいる，
ということ．
-
ここで使っている統計モデルは
まだそういう個体差を考慮した推定になってないので，
いずれ改訂せねばならぬ．
混合モデル化もしくはべいづ推定化．
これはなかなかたいへんな計算になるのであとまわしだ．
今回はたぶんべいづ推定なモデルにしたてあげていかないと
ダメなんぢゃないかな
……
とゆーか，
たぶん避けられない．
-
今回の対象樹種，
全 22 種に関する自動作図 & 自動作表関数を完成させる．
なお，
いまやってのるは樹種ごとに独立の推定．
はい．
とりあえず．
-
とりあえず
「ここまで計算できました」
報告を牛原さんに．
どっちが院生でどっちがせんせーだかわからなくなってきた．
-
推定計算各部の細かい考えかたを説明すると時間かかる
……
が，
ここは理解してもらわねばならぬところ．
大学院生や研究者にはふたとおりあり，
わからないときにわかりませんというタイプと
わかってるふりをするタイプ．
牛原さんが前者なので助かる．
-
お茶部屋で説明だの議論だのしてるわけだが
……
ここは「食堂」でもあるわけで．
こちらが尤度だの環境要因のうまい定式化うんぬんと検討してるかたわらで，
院生たちが旺盛な食欲を発揮しつつ
「いやー，
あたしって植物生態と食べものぐらいしか話題ないんですよねー」
といったおしゃべりに熱中してるのは，
なかなかに猥雑でよいのかもしれない．
-
夜になって屋久島相談はひとまず終了．
-
院生部屋 A802 の電力事情が悪化している
(A801 も?)．
つぎたし工法で構築された
(さらにメンテナンスがへっぽこな)
A 棟 7-8F は北海道の建築物とは思えぬほど保温性が悪い．
しかも寒さに敏感なヒトが多いので冬は局所暖房のため
電気ヒーターなどが使われる．
これが停電の原因となる．
しかもいつのまにか大電力消耗するデスクトップ機が増えてしまって
(ええ，私が予算消化のため買ってしまったマシンもそれにすごく貢献してます)，
ますます停電おこりやすい状況になってしまった．
-
計算機の分散でもはかるか?
A802 在住で電力大食い計算機つかう女性大学院生に，
試みに
「電力的に余裕のある A803B (かとーオフィス) にひっこします?
私のうしろの大きな机が空いてるんだけど」
とお尋ねしてみたところ，
ただちに「絶対にイヤです」と強く断られてしまったので，
この案はぼつ．
-
となると，
あとは A802 に隣接する部屋 (かとーオフィスもしくはお茶部屋)
から電線ひっぱるしかないのかな．
昨年，
電流増大案とやらが出たんだけど，
いつのまにやらうやむやになってしまったし．
-
2000 研究室発．
2030 帰宅．
2225 自宅発北大構内走．
2310 帰宅．
体重 72.2kg．
うんどうぶそくのせいか，
重くなった．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
米麦 0.4 合．
ハクサイ・タマネギ・エノキダケ・大豆のシチュー．
- 昼 (1220):
蕎麦．
- 晩 (2340):
スパゲッティー．
トマト・ピーマン・ニンニク・バジル・サンマのソース．
ナメコ・海藻スープ．
