ぎょーむ日誌 2002-08-30
2002 年 08 月 30 日 (金)
- 0830 起床．
どうもいまいち．
シャワー．
0840 自宅発．
曇．
0850 研究室着．
朝飯．
コーヒー．
- 昨日のぎょーむ日誌をほとんど書いてなかったんで，
でれでれと書く．
今日の午前中は口頭試問あるんだけど，
下っぱは出なくてよいでしょうという理由で怠業する予定．
その準備だけは昨日のうちにあれこれと．
- 0955 やはり雪野さんがその件に関してかとーオフィスに来られたんだけど，
あーだこーだ言って難を逃れる．
- 1030 から 1230 まで某あやしげ財団に金をせびる紙切れを準備する．
金額としは 20-30 万円だし使えるのは来年度からだけど，
9 月中に採否が決まるんで今後の見通しがたちやすいだろう，
というようなかんぢで．
- ひとまず片づいたので昼飯．
PD 石井さんの研究費に余剰が生じつつある公算大，
との情報に接したので，
無一文な私としてはそれを横取りしようという計画を考える．
- 午前中に書いたあやしげ作文には研究科長 (ま，学部長ごときもの)
のハンコが必要らしいんで，
雪野さんのお手をわずらわせて取ってきていただく．
封筒の宛名がきで漢字が手書きでは書けなくて苦労する．
- かとー先生と相談しつつ，
石井さんに無心する本の一覧というのを作っていく．
- 1500-1830，
お茶部屋で下野さんと (言うところの)
logistic 回帰の問答．
いろいろと自分で調べた下野さんのココロの中にある
logistic モデル像の説明とそれへの質問でナゾの部分を解明していこう
とゆー試み．
「わからない」という状況からも逃げ出さない下野さんの粘り，
かとー先生・石井さん・谷さんたちの
適切な合いの手によって当方はたいへんに勉強になった．
そう，
おそらく私がもっともトクをしている．
- で，色々なことがわかってきたんだが
……
たとえば統計学ソフトウェアに「ロジスティック回帰」をやらせると，
尤度比検定による帰無仮説との比較があるのは当然として，
個々の推定値にもイヤらしい p 値が計算されている．
これは何の計算しているのか?
……
皆で調べてみると Wald 検定だということがわかった．
- これは何らかの方法で推定値の誤差を評価してやり，
それで標準化した最尤推定値は
これまたイヤらしい標準正規分布に漸近する
(帰無仮説のもとでは)，
という性質を使ってるんだね．
ということで，
ある仮定のもとで個々のパラメーターが N(0,1)
からの標本になってそうかどうかを見ている．
他のパラメーターとは完全に独立であると仮定してるわけね．
- さらに面白いのは，
この Wald 検定は
尤度方程式全体を評価している尤度比検定と漸近的に一致
……
ふーむ．
- S 言語つかった S-Plus だの R だのでは，
対比行列 (contrast matrix)
という概念が重要になってくる
……
というところまでわかったんだけど，
この対比行列なるものの全容は解明しきれなかった．
今後の課題とゆーことで．
- SAS 系ソフトウェアでは飽和モデル (saturated model)
との尤度比が計算される
……
これは，
言わば「これ以上パラメーター数を増やせない」
というところまでパラメーターを増やしたモデルで，
ここからどこまで減じていけるかを見るわけだ．
- ……
というような問答で午後が終ってしまった．
統計学問答ばて．
1850 研究室発．
1900 帰宅．
体重 72.4kg．
うーむ．
- ばて．
- [今日の素読]
-
NHK ラジオ「ビジネス英会話」
Lesson 22 ``Volunteer Families'' (8/26-8/31)．
- [今日の運動]
- [今日の食卓]
- 朝 (0900):
研究室お茶部屋．
パン．
コーヒー．
- 昼 (1240):
弁当．
研究室お茶部屋．
米麦 0.7 合．
- 晩 (2130):
パン．