ぎょーむ日誌 2004-06-21

2004 年 06 月 21 日 (月)

• 0740 起床． 朝飯． コーヒー． 0850 自宅発． 曇． 0900 研究室着．
• 苫小牧直径成長モデリング …… 今朝がた得られた計算結果はなんだか奇妙だぞ? またモデリング wiki 上に報告まとめつつ 「どっかに計算まちがいありそーな」 と考えていたら …… しばらくしてから， かなり阿呆阿呆な失敗をやったと気づいた． うむ． これは釧路大会雑用後遺症にちがいない．
• えーとですね， 直径成長の混合モデル (樹木ごとの個体差を考慮したモデル) の最尤推定計算やらせるときに， 数値積分を含む尤度方程式の評価を 拡張モジュール でやらせてんですよ． 時間かかるから呪われ言語 C を使って． で， このときに double 型の値の絶対値が必要なんだけど， これを <math.h> の abs() でやらせちゃったんだよね …… これは「整数型」の絶対値を返すやつだ． 正しくは fabs()， ですね． 計算やりなおし．
• で， これまでの計算結果とかの整理にちょっと手間どってみたり．
• M2 平林さんからいつものごとく難しい計算問題をだされて， 数時間これにトラップされてしまう． そして解決できてない． 期待値が given である多項分布の推定， というやつ …… 三とーりの多項分布を考えんといかんのだけど， そのうち一つだけがとてもややこしい． 尤度方程式を作ることができれば R の optim() で強引に解けるんだけど， それすらできない． いまのところ．
• 1355 北大構内走発． 今日も曇天ときおりぱらぱら雨． 1430 戻る． 昼飯． 宣教師殿から office24 札幌駅前通店 なる 24 時間営業の印刷屋の存在を教えてもらう． たしかに半光沢紙にしただけで， 見てくれがかわるな． ただし公費つかえんとのこと．
• 本日も釧路大会を阻止するべくさまざまな自爆メイルテロが． 最小限の処理だけして放置．
• 苫小牧直径成長モデリング wiki ペイジ整備をススめる．
• fabs() を使ってやりなおした計算結果， これで問題なさそうに見えるんだが …… 昨日の計算結果は再現されんなぁ． そのあたりを調べる計算をもういっぺんだけやらせる． これでまた 7.5 時間かかったりするわけだが．
• 2040 研究室発． 2050 帰宅． 体重 73.6kg． なぜか減らない． 晩飯．
• [今日の運動]
  • 北大構内走 1355-1430． ストレッチング．
• [今日の食卓]
  • 朝 (0830): 米麦 0.7 合． キャベツ・ニンジン・ブナシメジ・煮干の味噌汁． トマト．
  • 昼 (1455): 研究室お茶部屋． 食パン 2 枚．
  • 晩 (2120): 米麦 1.0 合． キャベツ・ニンジン・ブナシメジ・煮干の味噌汁． トマト．