ぎょーむ日誌 2004-05-06

2004 年 05 月 06 日 (木)

• 0900 起床． う， やはり周期がずれてる． 朝飯． コーヒー． 0930 自宅発． 晴． 0940 研究室着．
• さて， 昨晩からやらせている苫小牧直径成長モデルの推定計算だが …… 10 時間 40 分で気象値 972 セットか． やはり Dell 機は計算おそい? で， 結果をみると当年の代謝速度みたいな値がマイナスに影響してる． これは気温と降水量の背反する関係を反映しているんだろうな． いまさらながら， 気温・降水量の影響をそれぞれ分離して推定するのは難しそうだ．
• 釧路大会雑用にとりくんでみる．
• 進捗しないうちに， 都立大の加賀屋さんから， offset つき Poisson モデルと logit モデルの関係についての スルどい質問うけたので， あわてて Alan Agresti カテゴリカルデータ解析本など 調べなおす． いや， ちょっと気にはなっていたもんで． よりどころは 「対数線形モデルと logit モデルの等価性」 といったところなんだけど …… offset 項が小さいときに (極端には log(1) = 0)， びみょーにずれが生じるんではなかろうか． 0 < λ << 1 かつ Prob[x > 1] = O(λ2)? ふーむ
• 雪野さんがとある申請書類の Word ファイルをおひす98 for Mac で扱おうとして苦労されている． たしかに起動して 10 秒以内に停止するなぁ． 何かよろしくない機種依存文字のたぐいが入ってるんでは．
• 1300-1335 北大構内走． 晴． 昼飯． 院生のみなさん， そろそろ調査がたいへんになってきたようで．
• お茶部屋 iMac に MacOS X 10.3 を再インストール． 修理から帰ってきた時点では 10.2 にされていたんで． ついでに X11 や Fink もインストール．
• 釧路大会大雑用は進捗せず小雑用の対応におわれる．
• イヤになったんで， 苫小牧直径成長モデルの推定プログラムで遊ぶ． てきとーにモデルを単純化して計算させてみると， なンかこの複雑な関係の理解が少しマシになるよーな …… 想像だにしなかったのだけど， ひょっとすると単純に 「雨が少なかった年」 というのが 「良い前年」なのか? うーむ， これは検討してみなくては．
• で， いろいろな状況で推定プログラムを走らせたところ， なかなか興味ぶかいパターンが …… もしかしたら， ですよ， 観測塔上にて測定された変動はげしい 光量子密度量 (PPFD) の観測値の大きさなんぞよりも 「その夏，雨がどれだけ降らなかったか」 のほうが， 翌年のイタヤカエデの成長速度を知る上では 信頼性の高い情報源でありうる， というような ……
• まったくもって現象論的なハナシというべきなんだろーけど， 確率論的モデルのあてはまりは格段に異なる． 人間から見ると 「今年はあまり降らんなぁ」 という年であっても， 苫小牧の樹木たちにとっては これがもしすでに十分なる水供給であったとするなら …… いやー， じつにトリッキーでじつに面白い．
• この見とおしは正しいのか正しくないのか． 役に立つのか立たないのか． 少なくとも， 数日前から気になってた推定プログラムの 「奇妙な推定値」 の謎は完全に解明された．
• というようなことを苫小牧ボスに御報告する作文にて 本日のぎょーむはほぼ終ってしまった． おっと， 帰宅する前に 気象庁サイト から 3287 日ぶんの 24 時間観測データをダウンロードせよ， という秘密指令を自動処理してくれる Perl スクリプトをしかけておかなくては． これってかなり「くらっきんぐ」に近いよーな．
• 2050 研究室発． えーい， 釧路大会大雑用を放置するのは危険だな． 2100 帰宅． 晩飯．
• [今日の運動]
  • 北大構内走 1300-1335． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0910): 米麦 0.8 合． ハクサイ・コンニャク・サケあらの味噌汁．
  • 昼 (1350): 弁当． 研究室お茶部屋． 米麦 0.8 合． ハクサイ・コンニャク・サケあらの味噌汁．
  • 晩 (2200): 米麦 1.0 合． ハクサイ・コンニャク・サケあらの味噌汁．