ぎょーむ日誌 2005-09-27
2005 年 09 月 27 日 (火)
-
0740 起床．
(書類上の) 夏休み第二日目．
朝飯．
コーヒー．
洗濯．
ネット雑用などかたづけて
1010 自宅発．
晴．
1025 研究室着．
-
なんだかどたばたしつつ，
一般化推定方程式 (GEE; generalized estimating equations)
の勉強のつづき．
-
やっぱりこれって正真正銘の quasi-likelihood 法だ
……
分散の指定のしかたがいささかコっているだけである．
-
で，
その分散を指定する working correlation matrix
の意味なんだが
……
自分で簡単な例を作ってるうちにだいたいわかった．
しかし，
これって時系列とかそういうハナシではないように思う．
というのも対称行列になってるんで，
これだと「過去」だけでなく「未来」からも影響されることになるよ．
まあ，
非対称行列にすればそれで時系列の効果を表現できると思うけど．
-
昼飯．
亀山さんから大雪山ツガザクラ調査に関する驚きのハナシをいろいろとうかがう．
じつは今まで
(私などには)
エゾツガザクラ
として聞かされてた樹種はぢつはほとんどが
コエゾツガザクラ
である，
とか．
で，
雑種だろう疑惑が濃厚なんだけどコエゾは「品種」あつかい，
とか．
うーむ
……
-
まあ，
ツガザクラ雑種形成シミュレイションはもうちょい放置することにして，
と
(いささか落ちつかなくなってきたが)．
-
昼飯後，
一時間ほど高田さんと来年度新入生の雑談につきあう．
-
また GEE 勉強つづき．
1600 すぎになって，
よーやくこれがどういうふうに推定計算されて，
応用されてるのか理解できるようになってきた．
これまた (他の quasi likelihood な方法たちと同じく)
Bayesian の普及とともに衰亡していきそうな一族に含まれてるな．
こういう ``practical'' な方法を考えだしたヒトたちは
ホントにすごいと思うけど．
-
アタマの中の整理と称して 1615 北大構内散歩発．
曇天．
いまにも雨ふりそう．
北大生協にいって VMware アカデミック版の購入申請書かき．
うろうろ歩きつつ自分自身に GEE とは，
と説明してみる．
-
あ，
書き忘れてたけど，
これは Calgary 大学在の
石井さんから出された宿題に回答するための検討．
で，
おもに参考にした文献は
(これまた石井さんに reading assignment
として指示された)
Horton, NJ and Lipsitz, SR. 1999.
Review of Software to Fit Generalized Estimating Equation Regression Models.
The American Statistician 53: 160-169.
また，
計算法については Yudi Pawitan 尤度本
``In All Likelihood''
の Chapter 14 が参考になる．
ネット上の教材だと
Højsgaard
の
GEE 解説
とかかなあ．
-
ともかくひととーりわかったんで，
宿題の回答メイルかき．
結論: GEE つかってもご利益ないんで (これまでどーり)
一般化線形混合モデル (GLMM) でいいんぢゃないでしょうか．
-
1940 研究室発．
1955 帰宅．
2135 自宅発北大構内走．
2230 帰宅．
体重 73.4kg．
晩飯．
-
「(書類上の) 夏休み」第二日目終了．
また夏休み
うそ日記
を書く．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0840):
米麦 0.6 合．
ニンニク茎・ショウガ・鶏ハツ炊きこみ飯．
マイタケ・ワカメ・豆腐のスープ．
- 昼 (1240):
研究室お茶部屋．
米麦 0.6 合．
ニンニク茎・ショウガ・鶏ハツ炊きこみ飯．
海藻スープ．
- 晩 (2330):
米麦 0.7 合．
ニンニク茎・ショウガ・鶏ハツ炊きこみ飯．
キャベツ・ピーマン・豆腐のカレースープ．