ぎょーむ日誌 2004-02-29

2004 年 02 月 29 日 (日)

• 0850 起床． 朝飯． コーヒー．
• 朝から昨日の計算の続きにしつこくとりくむ． 時間もあまり残されていないのに …… ThinkPad X31 (Pentium M 1.6GHz) は入院中なんで， ThinkPad 240Z (Cerelon 0.5GHz) で計算させる． アカマツ stem 数 150 ぐらい (つまり 1 個体) 500-1000 iterations で収束して 0.5-1 分ぐらい． R プログラミングの下手さが計算速度に悪影響あるかも， と思うと気分わるい．
• 重量成長モデルにサイズ依存性を組み込むと， 推定結果がマシになると確認できた． このあたりはモデリングわざが要求されるところだ． ただし， 推定そのものの不安定さは (これまた昨晩よりはマシになったとはいえ) まだ解消されていない． とくにサイズ依存的な死亡の推定が難しい． 他のアカマツ個体でどうなるか， をまだ確認してないけど， 場合によっては死亡モデルは抜いていいのかもしれん …… いや， やっぱりダメかな．

[なかなか面妖なる]

サイズ依存性を示していて， その不思議さを楽しめる …… けど， 推定がうまくいってないところはコホート全体の重量成長の点からは， 重要ではないんで， まあ， このまま計算しちまうというテもあるわけ．

• 死亡に関するパラメーターが初期値依存なのは， このサブモデルがどうでもいいからだ． 複合分布の推定への寄与が少ない． いっぽうで， 成長モデルに関しては， そちらにはあまり影響される， だいたい同じ値に収束してるようだ． むろん完全に独立してるわけではないけれど．
• 1310 自宅発北大構内走． 晴． 1355 帰宅． 体重 73.2kg． 昼飯．
• 1600 自宅発． 曇． 寒い． 1610 研究室着．
• またアカマツ stem のサイズ分布を推定する問題にとりくむ．
• 小菅君が 進化植物学研究会 関連でうろうろしている …… どうも講座からの申込者が少ない， と思ったら連絡がいってなかったようで． このあたりは Trendy セミナー の ML で通知されてたんだけど， 今の修士課程院生は (M2 大澤君をのぞいて) 誰も登録してないんだよね (さらにその後に M1 平林さんも学外アカウントから trendy-ml を取っていた， と判明)． いやはや， と思いつつ M1 & M2 まとめて強制登録．
• 一個体ごとの計算はだいたいうまくいくようになってきたんで， R プログラミング Tips 大全 読みつつプログラムの改善にとりくむ …… このときバックアップをとらずにいきなり改良に着手してしまったわけだが …… たいへん幸運なことに， データ構造の変更のワナにはまることなく， かつ数倍の高速化が実現した． しかし， じつはまだまだ効率が悪い， というおそるべきことも理解できてしまった． これってかなり面倒なプログラミング言語ですよ．
• いろいろと計算したいんだけど， 時間切れ． 帰ろうとすると， お茶部屋で小菅君と大澤君が Trendy セミナーの今後についてうだうだと話している …… 大学院生たちを集団としてみた場合， その「やる気」変動はけっこう大きいように思うんで， こういうのは大学院生まかせにだけはできんのかもしれんな． 講座サーヴァー管理を院生に押しつけるとかえってたいへんであるように (日本の大学の多くの研究室では院生がサーヴァー管理奴隷やらされてんだが)． 2220 研究室発． 2240 帰宅．
• [今日の運動]
  • 北大構内走 1310-1355． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (1000): 米麦 0.7 合． 納豆． キャベツ・ニンジン・豆腐・コンブの味噌汁．
  • 昼 (1420): 米麦 0.7 合． 納豆． キャベツ・ニンジン・豆腐・コンブの味噌汁．
  • 晩 (2340): うどん． キャベツ・ニンジン・豆腐・コンブの味噌汁．