ぎょーむ日誌 2005-02-05
2005 年 02 月 05 日 (土)
-
0820 起床.
外はかなり雪がふってる.
朝飯.
コーヒー.
洗濯.
提出用ぱいぷ樹木原稿関連ファイルなあれこれ処理.
-
1040 自宅発.
雪.
1100 研究室着.
怠業.
昼飯.
北大生協書籍部まで雪中散歩.
-
統計学勉強.
多変量正規分布 (multivariate normal distribution)
の条件つき分布について.
他の変量が決まってるという条件のもとで,
ある変量の確率分布は
(共分散行列で決定される)
ある定数分散の正規分布となる.
うーむ
……
これは私のへっぽこというかアテにならぬ思い込み
(分散も増減すると思ってた
--- 平均は当然ながら与えられた条件によって増減するわけだが)
がまた間違ってたわけで.
いやはや.
-
蛇足ながら
R
で逆行列を計算する関数は
solve()
(これはちょっといやはやな名前ですな).
あるいは MASS package には ginv()
なる
ムーア・ペンローズ型一般化逆行列を使った
計算関数もある.
-
ということで,
Markov Chain Monte Carlo (MCMC) 法において,
多変量正規乱数の Gibbs sampler はとても簡単で
計算負荷の少ないものとなる,
ということで.
-
そして,
空間相関ある GLMM みたいな状況で使う場合には,
線形予測子全体が多変量正規分布になる
……
と考えたほうが,
やはり私には理解しやすい.
理解しやすい,
というかそういうふうに計算しないと,
ダメだ.
Random effects の平均がゼロにならない,
というのはなんだかヘンなかんぢではあるけれど.
いや,
集団全体ではゼロになってんのかな.
-
ともあれ,
このような技法でもって苫小牧寄生蜂問題だの,
林冠動態ギャップ拡大問題
(こちらは私にとっては 10 年来の難問だった)
だのといった愉快な面倒は解明されていくのでありましょう
……
とゆーよーなハナシを「空き時間の埋め草」
として
進化植物研究会
(つい先ほど手直しして公開)
でやらんといかんかな,
と思っていたんだけど,
西村さんの尽力で大原研のかたがしゃべってくれそうな
(現時点ではまだ演題などいれてない)
……
ということで,
道内の生態学研究者たちが
MCMC 洗礼をあびるのは先のばしされた.
もとをただせば小菅君が小菅人脈
(これは旭川の陸自第二師団内部から
故郷たる横須賀の近所の総理実家その他にまでつながってるらしい)
を
「こんなことのために貴重な人脈を使うのは」
と出し惜しみしたのが,
今回の混乱のそもそもの発端といってよいだろう.
そうにちがいない.
-
1900 研究室発.
雪.
1910 帰宅.
晩飯.
早めに寝てしまったものの,
夜中に目覚めてしまった.
あいかわらずの降雪.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0850):
米麦 0.5 合.
トマト・豆腐のスープ.
- 昼 (1230):
研究室お茶部屋.
米麦 0.8 合.
おでんの残り
……
には固形物が含まれてはいなかったのだが.
- 晩 (2000):
米麦 1.0 合.
トマト・豆腐のスープ.