ぎょーむ日誌 2005-07-09
2005 年 07 月 09 日 (土)
-
0810 起床．
朝飯．
コーヒー．
-
怠業，
とゆーか，
昨日のぎょーむ日誌の数式画像どもをてぎわよく生成するため
tex2png.pl
(自作 Perl スクリプト)
の改良．
-
日射はそれほど強くないんだけど
……
北大構内走にでる気にならず．
お，
すこし雲ってきたかな，
と
1430 自宅発北大構内走．
とおもったらまた晴れてきた．
1520 帰宅．
体重 73.8kg．
-
北大構内走中，
また統計学問題を考える．
昨日のつづき．
生物の「いる・いない」問題を検討するときに，
「場所 A とその近傍の場所 B での植物 X のいる・いないは独立ではない
(相関がある)」
といったたぐいの問題 (他にも同所的な群集の構成問題とか)
を考えるときにマルコフ連鎖モンテカルロ
(MCMC)
法を応用した解析が今後はどんどん使われるようになるだろう．
-
で，
昨日検討したよーに，
それには
-
MPLE 補正 (Geyer & Thompson,
尾形たち)
-
EM アルゴリズムを使った最尤推定 (MCEM)
-
ベイズ推定
といった手法が使えそうだ．
現在，
生態学まわりでしょぼしょぼと使われてるのは 1. だけで，
これすらほとんど知られていないだろう．
基本的には MPLE を計算し，
それでもって MCMC を生成，
さらにその結果を使って MPLE を「補正」するとう手法である．
ちょい難解なので，
あまり広まらないんではないかな．
-
とはいえ，
2. と 3. もそれほど簡単ではないかもしれん．
直射日光さけつつ北大構内を走りながらそのあたりを考える．
まず，
これらの方法で MCMC の使いかたがちょっと違うように思える．
くだんの Gibbs 分布
(昨日しめしたよーに)
で「いる・いない」問題を考え，
やはり昨日の Gibbs sampler を使うところまでは同じなんだけど，
ここから先が違ってくるよーな．
-
(ここで考えてる MCMC 援用型の)
EM アルゴリズム
は (全部かぞえてしまうと) 莫大な個数になる「可能な配置」に関する期待値を，
Gibbs 分布からのランダムサンプル
(条件つきランダムサンプルとでもいうべきか?)
の標本平均というか期待値に置き換えてしまう，
という計算方法だと理解している．
つまり配置あるいはパターンが random effects であるよーな
混合モデルを考えていて，
まぢになって積分するのは苦しいので MCMC で代用している，
と．
-
いっぽうで，
べいづな推定
では
β
だの
γ
だのといったパラメーターの確率分布 (事前分布) を分母・分子にかけて，
さらに分母ではこれをパラメーターに関して積分して「つぶして」
(つまり分母を完全に定数にしてしまって)，
分子だけながめてその Gibbs 分布からの sampling やって，
パラメーターに関する確率分布 (事後分布) を生成させる，
という手ぐちだ．
-
そうですなぁ
……
将来は (とゆーかすでに現在?)
WinBUGS
みたいなソフトウェアが発展して，
「なんだかわかんないけど，
べーづ推定してみました」
という方向性なんだろうなぁ．
まー，
それはそれで楽しいかも
……
近ごろちょくちょく見かける
「なんだかわかんないけど，
glm() してみました」
だって当方の想像力をいろいろと刺激してくれているしね．
-
1705 自宅発．
曇．
いまごろになって曇，
か．
1720 研究室着．
-
しかし研究室にいるとなぜか研究ススまず．
ということで A 棟 8F 闇ネットの久保製の劣悪配線の部分的な修理に着手．
これは私が作ったへっぽこケイブルが原因で，
10M なら OK だけど 100M では通信できんという状態になっている．
ケイブルの作りかた
を復習し
(まんなか 4 本のうち 2 本を「交叉」させるところがぽいんと)，
とりあえずかとーオフィス内のそれを調べる
……
が意外にも，
このあたりはまともなケイブルばかりで (かとーさん作?)
一部の配線を組み換えるだけに．
[LAN ケイブルつくりなおし]
かしめ工具と先日購入したケイブルテスター．
-
お茶部屋は劣悪ケイブルばかりだったので全部つくりなおし．
闇ルーターまわりも，
意外にも OK か
(というか以前にこのあたり苦労させられた)．
A801 室もたぶん大丈夫だろう．
A804 セミナー室のケイブル点検にいったら小菅せんせーの「隠し」 PowerBook
が無人の室内で稼働していて不気味だった．
本人によると A801 からアクセスすることで隠匿気分を楽しむ，
とのこと．
-
お茶部屋で院生とデータ構造雑談
……
してるはずが，
広島風お好み焼きにもやしをいれるかどーかとか，
オタフクソースの万能性とか，
東日本における塩味のつよさとか．
2200 研究室発．
2230 帰宅．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0850):
クロワッサン．
- 昼 (1245):
蕎麦．
- 晩 (2255):
米麦 0.8 合．
コマツナ・ネギ・ショウガ・ニンニク・牛レバの炒めもの．