ぎょーむ日誌 2003-04-04

2003 年 04 月 04 日 (金)

• 0830 起床． どうもびみょーに不調ぎみだな． 夜中とか目がさめて寝られなくなったりするし． 朝飯． コーヒー． 0900 自宅発． 曇． 0910 研究室着．
• オショロコマ のメタ個体群動態で終わった一日．
• 午前中は農学部 前川研 の D3 小泉さんからサカナ・観測データ・モデルの説明をうける． 支流を local population と見立てて Hanski ふう メタ個体群動態モデル を作ろうとするも， パラメーター推定がうまくいってないようで， 長さ数十キロメートルこえる流域のメタ個体群全体があっさりと 数年で絶滅しちまう， と．
• で， Hanski グループ謹製の パラメーター推定・シミュレイター統合プログラム に観測データ与えて憶測された数値を詳しく調べてみると …… なんつーか， そっこーで全滅するのも当り前， というシロモノで．
• まずはここから改善しないといけない， しかしながら， この北欧製プログラムはブラックボックスだ …… というより表示から憶測するに怪しい計算やってるような気がする． ということで， データもらって自分で推定計算やりなおしてみる．
• Hanski 教授があたかもイヤがらせのごとく， 性質のよくない関数型ばかり選んでくれているので， 収束計算はなかなか難しい …… とはいえ， 極東某国にも迎撃の用意あり． 頑強さだけが売り物の Nelder-Mead 最大化 Perl モジュール， そしてこの手の呪われきった最尤推定の戦訓から得られた工夫あれこれ．
• 豪雨の中で転がりながらの組みうち， というかんぢのぢたばた試行錯誤の結果 …… 夕方ごろにまともそうな推定値が得られた． それを小泉さんにメイルで連絡． 先方がチェックしてるあいだに北大構内走にでかける． 路上の雪はほとんど消え， 日長も長くなったので夕方でも走れる．
• もどってみると， どーも良くない結果のようで …… さきほど送ったメイルを見直してみると， 推定値は間違ってないけど， 計算法を説明するところで単位系を間違えてたことに気づく． うかつ．
• 修正依頼を出してから帰宅． 1950 研究室発． 2010 帰宅． 体重 72.0kg． 晩飯．
• お， こんどはうまく計算できたようで． これにてひとまずは一件落着． まぁ， まだメタ個体群シミュレイション結果をどう解析するか， だとか， Hanski プログラムのどのあたりに問題あったか， というよーな問題あれこれが残存しているわけだが．
• 土日はニセコのスキー実習にでかける． 雪のシーズンの終わりになって， ようやく初めてスキーがはける …… ともあれ楽しいことです． しかし一年ぶりだから何も進歩してなくて， あいかーらずヘタなままとゆーのがなんとも．
• [今日の運動]
  • 北大構内走 1740-1815． ストレッチング．
• [今日の食卓]
  • 朝 (0840): 米麦 0.7 合． コマツナ・マイタケ・豆腐の炒めもの．
  • 昼 : 研究室お茶部屋． なんとなく食ってない． ミカンとか．
  • 晩 (2200): 米麦 0.7 合． コマツナ・マイタケ・豆腐の雑炊．