ぎょーむ日誌 2004-06-18

2004 年 06 月 18 日 (金)

• 0800 起床． 朝飯． コーヒー． 0900 自宅発． 曇． 0910 研究室着．
• 昨晩やらせていた推定計算は 7 時間ほどかけて 0318 終了， か． 計算のやりかたにもう少し工夫が必要か?
• で， 強制的に dbh.log 項に依存する， と仮定した計算だと． 推定値は

  	scale.log
          -1.02724964
  	const            dbh.log        dbh.log.sq
  	-6.06808777      5.04546525     -0.71576079
  	dbh.log:ps.prev  mb6 
  	-0.00322468      -1.05195074 
  	
  	temp.center=16.0/slope=0.3
  	ppfd: factor=0.02
  	vpd: factor=0.20
  	temp.mb: center=20.0/slope=1.2
  

  といった奇妙なものになるか． 定式化の段階からよく検討せんといかんな． 最大化対数尤度が -629 あたり， と尤度的には悪くはないんだが (降水量モデルとの比較において)．
• 念のため， 新しくした気象フィルターについての検算やってみるけど， まぬけミスは無し． やはり推定計算が「ヘン」なのだ． えーと， dbh.log をふつーの項としてあつかった場合は dbh.log:mb6 (+), ps.prev(-) とゆーような組み合わせが最良だったかな． というかどういう気象フィルターパラメーターにおいても， だいたいにおいてこの組み合わせだった． モデル選択関数は大丈夫なのかしらん．
• とりあえずの策がないので， vpd ぬきな気象値で計算させてみるか? 今回は dbh.log についてはおまかせにして， と． サイズ依存項 dbh.log を「必修」にしてしまうと， やはりちょっと推定結果がぎこちないような気がする．
• いやいやいっそのこと光合成・代謝の両面について， 気温だけが影響するという極端モデルはどうだろう? もはや苫小牧ボスの意向無視しまくりではあるんだけど， 私としてもよーやくにして 気象データなるものの複雑さ・面倒さが理解できつつあるんで． 気温だけなら 2 時間弱で片づきそうだな．
• まーたまた， 釧路大会を阻止するために わーどファイルを抱えて飛びこんでくる自爆テロが発生． なぜそこまでして Bill Gates 師への忠誠を誓うのか …… 彼が天上の楽園にいたる階悌を準備してくれているとでも? あるいは何か弱みをおさえられてしまっているとか． ともあれ， 放置だ．
• 気温だけに依存モデル， 1 時間で計算結果でた． ふーむ， 新しく導入した単調増加型の気温フィルター， 文字どーり完全にウラ目に出てるね (最良モデルに関しては)． vpd うんぬんの問題以前に， 気温フィルターがよろしくないかも， と特定できた． さーて， どうしたもんかいのう． 時刻は 1135．
• じゃ， 気温 → 光合成なモデルだけやってみる． これは 7 分ほどで片づく． やっぱぜんぜんダメか． フィルターのパラメーター範囲を変更してやりなおし …… ふーむ， やはり光合成に関しては， 「気温の高すぎる日」 は捨てなきゃいけないんじゃないの． 環境が統制された実験室で測定される光合成速度と， イヤらしい交絡だらけの 野外のそれのあいだになンか違いがあるような気がするな． 気温のキキかたについて． あるいは瞬間速度と長期にわたる同化量の違いとかに．
• 気温なるものが (これまた日平均気温なんだが) ある一日の光合成量を規定しており， さらに「ある年の夏によく光合成してるなら翌年の成長はよい」 ことが正しいと仮定する． するとこのような気温フィルターが最良， ということになる． つまり 18°C あたりがベストで 気温が低すぎたり高すぎたりするのはダメ， ということだ．

• もちろん， 上の仮定が成立しない場合もあるだろう． たとえば， 「ある年の夏によく光合成してるなら翌年の成長はよい」 …… これは一見もっともらしいけどホントなのか? そうでない場合もありうる． 例をでっちあげるなら …… ある年の夏によく光合成できたなら， その年の秋の繁殖に使ってしまうんでは? すると翌年の成長が悪くても不思議ではないよな． このよーに 「今年の光合成量が多いと翌年は成長がよくない」 なるひねくれきったハナシが正しいとすると， 気温は高けりゃ高いほどよい， という推定結果が最良モデルとなる． さー， どっちだ? 私にはわからんねえ ……
• 考えてもしょうがないので， 以上のような事態を報告にまとめて， 苫小牧ボスの裁可をあおぐことにする． 関数型が違うから単純には比較できんのよ (めんどくさくなったらそのまま比較するけれど)．
• 1425-1505 北大構内走． 曇 …… のはずだったんだけど走ってるうちに晴れてきた． 昼飯．
• まーた， 釧路大会がらみの自爆テロ …… 口頭発表をもうしこんでおきながら いまごろになって参加を取り消すハタ迷惑なすっとぼけが． たぶん最初からこれをネラって申し込みをやったにちがいない． プログラム確定後のキャンセル処理過程はまだ作ってないのに． えーい， 爆散屍体は放置だ． それにしても， 釧路大会ってのはどーしてこれほどまでに憎まれてしまってるんだろね?
• さて， 北大構内走前からやらせておいた推定計算の結果だが …… 代謝速度への気温の影響についても， じつに巧妙としかいいようのない最尤推定値を発見してくれますねえ． 成長のとても良い年を選んで 気象値が劣悪になる気象フィルターを器用にも構成しやがったよ． 針の穴を通すよーな， といいますか． 神の企みに悪意はない のかもしれないケド， 冗談はきつすぎやしませんかね?
• 代謝を表現する気象フィルターの探索レンジを限定． なぜ生態学のデータ解析で Bayes 統計学が必要とされつつあるのか， つくづく身をもって理解できたよーな気分．
• またキャンセリング自爆テロか． いや， テロ予告だな． 一般講演とりけし …… こいつも狙ってやがったな．
• 推定計算ぢりぢりとやりなおしてみたり．
• しかしチカラつきたので， 10 時間ぐらいかかる計算を命じておいて 1940 研究室発． 2020 帰宅． 体重 73.8kg．
• 釧路大会の遅延要旨爆弾突入テロ． どうやら環境省とゆーお役所では 上司に代わって下っぱが爆発物を製作して 自爆してみせるのが礼儀にかなった作法らしい． たぶん上司のほうは公開講演会壇上で自爆するんだろう． 私のまわりはどんどん死屍累々になってんだけど， ほーちだ．
• 今日はテロ厄日か …… 金はらったけど参加登録してませんどうにかしろテロ． さしづめ人質はとったから言うこと聞けやこら， というところか． 放置．
• [今日の運動]
  • 北大構内走 1425-1505． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0820): 米麦 0.7 合． タマネギ・トマト・ピーマン・コマツナ・エノキダケ・ホタテ の炒めもの．
  • 昼 (1520): 研究室お茶部屋． 食パン．
  • 晩 (2130): 米麦 1.0 合． ニラ・モヤシの炒めもの．