ぎょーむ日誌 2002-12-22
2002 年 12 月 22 日 (日)
- 0920 起床．
5 時間ぐらいしか寝てない．
朝飯．
コーヒー．
- 1600 ごろまでひたすら数値計算プログラムの改良．
ばてたので昼飯．
昼飯食うと猛烈に眠くなって
……
意識が回復すると 2120．
やれやれ．
また生活周期がずれそう．
- うーむ，
今日いろいろやってみたんだが，
汎用最尤推定計算プログラム
mlfitting
はけっこう使えるやつなんだが
……
今回のややこしー推定には歯がたたぬようだ．
-
mlfitting
の要点は多次元の微分方程式群を
-
数値微分→勾配行列
-
勾配行列の逆行列計算
-
Newton-Raphson 法
とゆーようなよくある数値計算技法の
コンビネイション技でしんぷるに解いてしまうことにある
(加えて最大化したい尤度関数を
plugin_function.so
として動的にプラグインできるのもウリ)．
- 今日はさらにここに「最適な『歩幅』の自動探索」
なるものを加えたんだが
……
やはりそれでも解けなかった．
そんなに複雑な確率分布ではないはずなんだがなぁ
(下の図のようなガンマ分布と正規分布のたたみこみ)．
- かなり逡巡してから
……
まったく
別の計算方式である修正 Powell 法の
プログラムを作ってみることにする．
とりあえずのプロトタイピングとて
Perl モジュールとして書いてみることに．
いやはや，
数値計算ものぐるい状態かも．
- 修正 Powell 法モジュールつくっていく
……
コード量的には 2/3 ぐらい書けたか?
時刻はすでに 2800 すぎ．
しかしアタマは稼働状態のままなので，
ぬるい風呂に入ってでれでれしてから寝る．
- 今日は部屋から一歩も出ないで終わった．
ひたすらに数値計算な一日．
- [今日の食卓]
- 朝 (1030):
米麦 0.8 合．
ハクサイ・ネギ・シイタケ・卵の炒めもの．
- 昼 (1630):
パンケイキ．
コーヒー．
- 晩 (2230):
スパゲッティー．
タマネギのカレー．
