ぎょーむ日誌 2004-05-02

2004 年 05 月 02 日 (日)

• 0930 起床． 朝飯． コーヒー．
• 怠業．
• おそい昼飯くってから 1620 自宅発． 曇． 1630 研究室着．
• 苫小牧直径成長モデル． 昨晩ふらふらと歩きながら考えついた例の着眼， つまり負の二項分布を導出するところで置換積分が必要， にもとづいて推定プログラムを書き直してみる …… うん? あいかわらずヘンだぞ．
• で， よくよく調べてみると， じつはこの置換積分がぜんぜん簡単ではない， ということがわかってしまいましたよよん． 置換できん， というより分解不可能． 積分したいパラメーターが p で積分したくないパラメーターが q だとすると， どうしても exp(pq) というような部分がでてくるぢゃありませんか． せ き ぶ ん で き ん． えーい， どうしろってゆーわけ?!
• そりゃ， まあ， 数値積分はできますよ． しかし R でそんなことやったら， どうなることやら． ああ， せっかくの (Poisson)*(Gamma) = (負の二項分布) 作戦が …… だれか助けてー
• しょうがないので

  integrate(
  	function(c) (
  		dpois(growth, lambda = c * mu.nonconst)
  		* dgamma(c, shape = shape, scale = scale)
  	),
  	0, Inf, stop.on.error = FALSE
  )
  

  てなかんぢで， 数値積分を含むカタチに書き換えてみる …… 計算時間は 4 倍以上を要するようになってしまった． どひゃー
• しかも， みょーに簡単なモデルばかりが選択されるし． 定式化のまわりに， まだ何か間違いが残存してるってことか?
• このあたりの高速化にとりくむ． とりあえず dnbinom() 版の 2 倍程度の計算時間， というところまでは回復した． ふう． しかしあいかわらず単純化されたモデルばかりが選ばれる． 気象の影響より「個体差」のほうが重要， ということなのか? うう， それは困るんだけどなぁ．
• 計算ばてしてきたので 2130 研究室発． 2150 帰宅． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0940): 「北欧」バゲット．
  • 昼 (1500): 米麦 0.7 合． チンゲンサイ・卵の炒飯．
  • 晩 (2240): 米麦 1.0 合． チンゲンサイ・タマネギ・エノキダケ・ホタテの炒めもの．