ぎょーむ日誌 2002-12-25
2002 年 12 月 25 日 (水)
- 0800 起床．
朝飯．
コーヒー．
0850 自宅発．
晴．
0900 研究室着．
- ちかごろは朝にぎょーむ日誌を書いていることが多いような気がする．
よいのか悪いのか．
- 明日は苫小牧の leaf miner 寄生蜂の推定問題をやらねばならんので，
今日はまた昨日の複数確率分布推定問題のつづき．
こっちの推定問題発注者がとりあえず使える段階まで持っていく，
ということで．
- たたみこみ数値積分の改良．
これが重要だとわかったんで．
数値計算本などながめる．
なるほど，
これまたいろいろなヴァリアントがあるようで．
- それ以前にまずは，
Gamma 分布 (自然が作ったと仮定しているゆらぎ)
のまわりで積分するというわれながらまぬけな
計算方針はやめて Gauss 分布 (人間の測定誤差)
のまわりでの積分に切りかえる．
人間の測定誤差のほうが相対的に小さいと仮定してるんで．
もしそうでなければ，
二種類のばらつきを分離して推定などできやしない
……
この変更だけで計算速度が 5 倍ほど向上した．
- というような変更と入出力まわりの改善だけで午前が終わった．
1230 北大構内走発．
晴．
寒い．
1300 研究室戻る．
昼飯を食っていたらいきなり窓の外がホワイトアウト．
突発的な降雪である．
- 午後も数値積分の改良．
- たしかに拡張 Simpson 則を使ったほうが結果はマシになる．
劇的な改良，
というほどでもないけど．
# ホントの値
p0 -3.0e+00
p1 5.0e-02
beta 1.0e-01
sigma 1.0e-01
# Trapezoidal rule
estimated_p0 -3.253e+00
estimated_p1 5.131e-02
estimated_beta 7.255e-02
estimated_sigma 1.084e-01
log_likelihood 6.192070e+01
# Extended Simpson's rule
estimated_p0 -3.227e+00
estimated_p1 5.201e-02
estimated_beta 7.527e-02
estimated_sigma 1.060e-01
log_likelihood 6.236328e+01
とはいえ計算の手間を増やさずに
推定結果を改善できるので，
拡張 Simpson 則を採用する．
まあ，
積分の区切りを大きくしても結果がそれほど
悪くはならない，
というところに御利益ありそうだけど．
- とゆーことで，
とりあえず発注者の浦口さんに連絡して，
この推定問題はたなあげ．
すでに 1630 か．
数値計算ばて．
- うだうだとしてからメイル書きなど．
2020 研究室発．
北 12 生協で買いものしてから
2030 帰宅．
体重 73.4kg．
- 晩飯食ってからさっさと寝てしまった．
ここしばらく寝不足ぎみだったので．
- [今日の素読]
- [今日の運動]
-
北大構内走 1230-1300．
ストレッチング．
-
腹筋運動 30 × 3 回．
- [今日の食卓]
- 朝 (0840):
研究室お茶部屋．
米麦 0.5 合．
ハクサイ・キャベツ・ナス・
タマネギ・シイタケ・マイタケ・豆腐
の野菜のカレー．
- 昼 (1320):
弁当．
研究室お茶部屋．
米麦 0.8 合．
ハクサイ・キャベツ・ナス・
タマネギ・シイタケ・マイタケ・豆腐
の野菜のカレー．
- 晩 (2150):
米麦 0.8 合．
ハクサイ・キャベツ・ナス・
タマネギ・シイタケ・マイタケ・豆腐
の野菜のカレー．