ぎょーむ日誌 2004-04-25
2004 年 04 月 25 日 (日)
-
1100 起床．
ほぼ回復．
朝飯．
コーヒー．
-
1610 自宅発．
曇．
1620 研究室着．
-
釧路大会小雑用．
一般講演 & 参加申し込み数が増えてきた．
4/30 (金) しめきりなもんで．
で，
それにともなって「訂正依頼」ってのも増えてきた．
登録者が自分で自分の間違いを訂正できるシステムにしとくべきだった．
しかしもはや改造する気にもならんので，
不名誉なる「手動訂正」に終始しそうだ．
そう，
これは手間の問題ではなく名誉の問題だ．
-
苫小牧直径成長モデリング，
ボスたる鍋嶋さんの的確なるご指示により，
一年の中のある月の成長量にだけ焦点をあてた推定計算に切りかえてみる．
パラメーター数はめでたく 12 → 8 となり，
このおかげで
R
の
optim()
による負の二項分布の対数尤度最大化収束計算が
……
おお，
1 秒未満で片づくようになった．
いいねえ．
呪われ言語プログラミングに手をそめなくても良いかもしれない．
-
蛇足ながら補足しておくと，
R で
family = 負の二項分布
とする一般化線形モデルの推定計算は
glm.nb(MASS)
でできる．
なぜ私が迂遠なる optim(base)
なんぞ使ってるか，
と言いますと，
苫小牧直径成長モデリングは
glm.nb(MASS)
では片づかない計算もふくんでいるからだ．
-
で，
そのひとつは parameter mapper
(これは読みにくい R コードで書けた)
と連動する model selector とか．
-
……
で，
うだうだと効率のワルく作業して，
optim()
な汎用 model selector の骨組みを書いてみる．
簡単な試験運転，
パス．
時刻は 2150．
-
2210 研究室発．
2220 帰宅．
晩飯．
-
[今日の運動]
-
腹筋運動 30 ×
3 回．
腕立ふせ 3 ×
3 回．
-
[今日の食卓]
- 朝 (1130):
米麦 0.5 合．
チンゲンサイ・ニンジン・タマネギ・エノキダケのカレー．
- 昼 (1500):
米麦 0.5 合．
チンゲンサイ・ニンジン・タマネギ・エノキダケのカレー．
- 晩 (2240):
米麦 1.0 合．
チンゲンサイ・ニンジン・タマネギ・エノキダケのカレー．