ぎょーむ日誌 2004-11-01

2004 年 11 月 01 日 (月)

• 0750 起床． 朝飯． コーヒー． 0850 自宅発． 曇． 0900 研究室着．
• 52 時間ほどかけて苫小牧直径成長推定計算 (新データ) が終了． 予想したより計算時間が短かった， と思ったら …… ミズナラ (QC) 計算がトバされてる． データをみると， 個体 2120 の 2002 年の DBH だけが欠測になってるのが原因， と特定できた． とりあえず， そういう個体は除外して再計算をやらせてみると同時に， 苫小牧ボスにお尋ねメイルだす．
• じつは上の「計算とばし」は wiki 上に結果をアップロードしてる 最中に気づいたんだよねえ …… このあたりの作業， 半自動化になってるわけだがそのへんの関連プログラムの使いかた， すっかり失念していた． いろいろぢたばたしたり， とばしに気づいたり， お尋ねメイルとか出したりしてると， 時刻はすでに 1030． いやはや．
• ミズナラ推定再計算， 所用時間は 6-7 時間， というところか． つまり今日の夕方終了だな．
• またぢりぢりと鹿児島準備．
• いつものごとく， 苫小牧ボスから修正策のご指示いただいたので， それにあわせて計算やりなおし． 時刻は 1107．
• 1305-1335 北大構内走． 曇天． 昼飯．
• 昼飯後もぢりぢりと準備． 作図プログラムといったゆーがなコトやってる 時間も失われてしまったので， ドロー系ソフトウェア Tgif (これって Unix 界では 20 年ぐらい使われてんぢゃなかろーか …… と思ったら最初の version は 1990 年だった) なんかで書いて EPS 出力， LaTeX にとりこみ．

• 図はひととーりそろってきたか?
• ばててきたんで 1920 研究室発． 1930 帰宅． 体重 73.4kg． 晩飯．
• 晩飯後も作業．
• 2140 ミズナラ再計算終了． 予想より計算時間ながく， 10.5 時間というところ． つまり， 新データのもとで 7 樹種推定 (アオダモ & ミズキが却下され， ウダイカンバが追加された)， 合計約 62 時間かかる， ということか．
• さっそく Perl & R に作らせた報告と図を Wiki 上に更新 (全樹種そろわんと更新できないわけで)． また作図プログラム・感度解析プログラムの 微修正に苦闘する．

[ミズナラ 6 月の成長]

樹木個体数を増やすと

• 気象変量に振り回される年変動
• 非単調なサイズ依存性

両方を「個体差」の影響こみで推定できている， とゆーかんぢで．

• すると， 苫小牧ボスよりさっそくメイルいただき， ミズナラはこれで良さそうだけど， オオモミジに怪しい一個体あり， とのこと． 送っていただいた「外れ値チェック図」 はなかなか巧妙で， たしかにアブなそう．
• オオモミジは個体数おおくないんで， この影響は大きいかも， ということで …… また計算やりなおすことに (オオモミジだけ)． まあ， 明朝までには終わるか．
• これとは別に， モデル選択 (変数選択) の手順はこれでいいのか， という懸念もある． 自作関数でこのあたり処理してるわけだが …… まあ， 高速化を実現したあとでないと， 改善策を検討できないな．
• もう 2400 すぎですよ． 夜中にバグとりやったりしたんで， 現在アタマの「電源の切れる準備」中 (時間かかる) ……
• [今日の運動]
  • 北大構内走 1305-1335． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0820): 米麦 0.8 合． キャベツ・ブナシメジ・豆腐のシチュー．
  • 昼 (1350): 研究室お茶部屋． 米麦 0.8 合． キャベツ・ブナシメジ・豆腐のシチュー．
  • 晩 (2030): 米麦 0.8 合． キャベツ・ブナシメジ・豆腐のシチュー．