ぎょーむ日誌 2002-12-24

2002 年 12 月 24 日 (火)

• 0820 起床． 朝飯． コーヒー． 0855 自宅発． 曇． 0905 研究室着．
• ばたばたとメイル書きなど．
• Perl モジュール化した Nelder-Mead 式 Polytope 法を 最尤推定の道具とするべく調整を続ける． まずは簡単にとける推定問題 …… 先日やってた quasi-likelihood 推定問題で Polytope 法の試験を行う． なかなか悪くない．
• おっと， 最尤推定計算ばかりにかまけていないで 生態学会大会の要旨をさっさと提出しなくては． なにしろ 明日 必着だからなぁ． 共著者・甲山さんに検分してもらって発送の準備． 1130 研究室発． 1140 郵便局着． 速達料 390 円． ぎりぎりまで提出しなかった罰金か． しかもここ北日本から関東の辺境たるつくばまで， 郵便が到達するまで何日ようするか， じつはわかったもんではない， という …… ついでに生態学会会費支払う． 11200 円． 高い．
• 1220 北大構内走発． この冬初の雪上転倒をやってしまった． むこうからよたよたと走ってくる自動車をかわそうとして コケた． 圧雪の上で重心急移動をやってはいけない． 雪の上は危険だ． 1255 研究室もどる． 昼飯．
• ぎりぎりになって速達にて送信してしまう小心な私などより よほど大物である浦口さんの 講演要旨の「下書き」を拝読している自分を発見した．
• さて， いよいよ

  (観測された値) = (ある確率分布の標本) + (測定誤差)

  というような複数の確率分布の結合した状況 (というかスケイルパラメーターが複数ある状況) での推定問題にとりくむ． これは関数のたたみこみ， つまり数値積分を含むので， そのあたりをまず Perl コード化しないといけない．
• 再掲すると， こういう確率分布を考えているわけで．

[Gamma + Gauss 分布] Gamma 分布のばらつきを決める beta と Gauss 分布のばらつきを決める sigma を同時に決めないといけない (左の図では sigma = 0.1)． もちろん Gamma の平均も． Gauss 分布は人間の測定誤差をあらわすので 平均はゼロと仮定している (ホントにゼロなんだろうか)．

• 苫小牧から第三次攻撃通告の電文． えーい． 借金でクビがまわらなくなる， とはこのことか．
• という状況であるにもかかわらず， お茶部屋で D1 平尾君と統計学雑談を長々とやってしまった …… いやはや．
• 2030 よーやく長年の懸案であった， 複数確率分布モデルの最尤推定問題に解決への第一歩となる 計算結果が得られた． 下の図のタテ軸のばらつきがくだんの joint distribution (上の図) から得られたランダムサンプルになってるわけで．

[それなりに良い推定結果の例] 青線が推定で緑線がホントの平均値． 標本数 152 コ． まぁ， よくがんばりましたとゆーことで． まだまだ検討・改良が必要ではあるけれど． 推定値 ホントの値 p0 -3.14 -3.0 p1 0.048 0.06 beta 0.083 0.1 sigma 0.115 0.1

• Perl + Polytope 法というおよそ数値計算むけではない環境なんだけど， いまのところは「遅い」という気がしない． この問題だと 100-300 反復ぐらいで収束してるようだ． 時間にして 1-3 分程度 (PentiumIII 600MHz の ThinkPad X20 で)． とはいえ， これは標本数が 152 コ， とそんなに多くないから， この程度で片づいているんだろうなぁ．
• ただしこの収束過程を見てると Polytope 法が 「無駄計算ばかりやってやがる」 というのがひしひしと実感される． まぁ， それゆえの安心感みたいなモノはあるんだが．
• ひとまずは良い結果が得られて， すっかりうかれてしまったので， biometry というメイリングリストによけーなメイルを送ってしまった …… 不幸中の幸いでこれは biometry あてとはなっていなくて， 粕谷さんあての私信となったわけだが (で，その返信にてまたいろいろと 親切なるご教示をいただいてしまった)．
• 腹もへったし帰るとしますか． 2130 研究室発． 北 12 生協で買いものなどして 2150 帰宅． 体重 73.2kg． 晩飯．
• 「おお， この最尤推定の結果の善し悪しは， たたみこみの数値積分の精度にかなり依存してしまっているのか!」 というような計算あそびを続けてしまったので， また寝る時間が遅くなってしまった．
• [今日の素読]
  • 数値計算じゃんきーな日々が続くと， かかる敵国語研究がススまぬ， というような ……
• [今日の運動]
  • 北大構内走 1220-1255． ストレッチング．
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (0840): 米麦 0.8 合． キャベツ・タマネギ・マイタケの炒めもの．
  • 昼 (1310): 弁当． 研究室お茶部屋． 米麦 0.8 合． キャベツ・タマネギ・マイタケの炒めもの．
  • 晩 (2220): 米麦 0.8 合． 残りもの野菜のカレー …… ハクサイ・キャベツ・ナス・ タマネギ・シイタケ・マイタケ・豆腐． なんなのだろうか．