ぎょーむ日誌 2001-02-28

2001 年 02 月 28 日 (水)

• 0700 起床． 5 時間ぐらい寝た． 体力的には少し回復．
• 朝飯・弁当の準備． 朝飯． シャワー． ひげそり． コーヒー．
• 0743 自宅発． 曇． 0802 東京モノレイル流通センター発． 浜松町の気温 6 ℃． 0820 研究所着．
• 呪われて動かない ADSL 高速接続． (昨晩も試行したけどダメだった)． 東京めたりっく通信 のペイジを見てると， サーヴィスに関しては support@metallic.co.jp だけどモデムについては help@metallic.co.jp に， などとややこしいことがさりげなく書かれている． support のほうにしかメイル送ってなかったな． まるで返事が来ないんで help にもメイル送ってみる． たぶん人数が足りなくて破綻してるんじゃないのかな．
• もひとつメイル送る． これは次の次ぐらいの仕事の布石である．
• さて， ここの報告会とやらで 話す内容は森林三次元モデルの正当化みたいなことである． 地球表面のどろどろねばねばした流体にしか 興味持ち得ない病気に取り付かれた人々で構成されている やくざ組織に， 何故あーきてくちゃーぐるーぷなどというものがあるのか， という売り込みですかね． あー， やだやだ．
• で， 一昨日あたりからいろいろと手を入れているのが 三次元モデルというやつである． 森林動態モデルと言えば， 現在の主流は無次元か一次元． そこで， あまり意味ないんだが (ゼロ次元だろうが 26 次元だろうが 何でも勝手に好きなようにやってくれれば 良いんだけど)， 今回は三次元モデルと一次元モデルの比較， ということを行う． 残念ながら他に手がないためである．
• その前の準備段階として mlfitting を呼び出す estimator.pl というスクリプトを改造して， と．
• あ， また北大ネットワークが通信途絶． 毎朝毎朝のように不調になってるな …… というのは大間違い． こっちの建物のルーターが死んでいる． メイルも受信できない状況だ． さっきまではつながっていたのになぁ． やれやれ．
• 腐れルーターは死んでも， 私の机の上に構築されている 闇ネット・闇ルーターは健在である． 当然ながら． しかし， これはちょっと嬉しい．
• 自分でもまったく忘れていた機能なんだけど， 最尤推定プログラム mlfitting に --estimator (integer) という コマンドラインオプションをつけると 推定すべきパラメーター数を実行時に指定できる． もちろんこれは最尤帰無仮説を計算するために つけられている機能に他ならない．
• make growth_null で最尤帰無仮説をまとめて計算させてみる． この帰無仮説は光分布が林内で一様， と仮定してる無次元モデルに該当する． これと比較すると， 三次元光分布モデルを使った計算では， たしかに全体的に尤度は上がっている． 昨晩ちょっと確認したとーりに． さて， それでは一次元モデルはどのあたりに落ち着くのか ? 今日はこれを計算しなければならない．
• 他にありそうなモデルとして …… 甲山さんが自分の趣味に合致している 「準三次元モデル」 みたいなやつも比較しろ， と言いそうだけど， 報告会とやらではやらない． 聞いてるほうはもはや何が何やら区別つかんだろうから． そういう計算機能はつけるけどね．
• いやはや． 領域長の安岡さんの命令で 金曜日に「報告会の練習」とか称して， このあたりの計算結果を 喋らされることになってしまった． 修論発表会じゃあるまいし． 大学のせんせいにはこういうヒトがいるんだよねえ …… とはいえ， たいへん残念なことに， 私だけそれをボイコットするほどの 勇気も持ち合わせていない． 金曜日は要旨提出もやらんといかんのになぁ． さて， 間に合うのか． 立腹ぱわーで間に合わせるか．
• Bad クラスというのを作る． 一次元モデル計算用である． 名前に他意はなく basal area distribution の略である． ふだんはこういう略称はクラス名には使わないんだけど (あとで見たときわからないから)， これはなぜか気に入ってしまったのである．
• ネットワークは 1055 に回復．
• Bad クラスを組み込むと コンパイラーが 呪われ STL 関連エラーを大量に出す． 正確にはリンカー， か． なんだろ， これは …… そうか STL で降順 sort するには

	sort( container.begin(), container.end(), greater<Type>() )

とするんだった． 最後の () を忘れていた．
• 1145 ようやくにして一次元 (つまり垂直) BA 分布クラス Bad を組み込んだ小川シミュレイターの 試験運転完了． みょーに手こずったな． いやいややってるからか ?
• 縦深メイル防御陣より警報． うわっ， 文字化けメイル …… とおもったらメイラー Mew はちゃんと解読してくれた． 何を考えているのか本文が base64 で encode されていたのである． 東京めたりっく通信からようやく返信． みょーに丁寧なおわび． ADSL モデムが呪われてしまっている原因は不明． 自力で解決できないなら担当のモノを派遣 …… 指示された web ペイジ見ると， 平日のみ作業してるとのこと． つまり ADSL 接続のために 有給休暇をとれってことですか． くそう．
• ちょー高速常時接続のためには涙を飲んで 3/12 (月曜日) に予約をいれてみる．
• とりあえず， Bad を使ったパラメーター推定作れるように 推定系を作り直して行く． まずは静的シミュレイターから 値を出し入れするところ． おっ， 小川ではふたつの指標の間に 多少は相関あるかな ?

[さて ?]
横軸が開空度
らしき量で，
縦軸は一次元
の合計胸高断
面積．

しかし縦軸はサイズとほぼ同義なんだよね (上下が逆だけど)．
• パラメーター推定部に手をいれる前に， 昼飯にしようかな．
• 相変わらず， fj.comp.lang.c や fj.soc.law なんかでは， 人間にとって狂気とはごくごく身近なものだと 感じさせる投稿が続いている． もちろん真面目に読んでるわけでなく， 「あ，今日もやってるな」 と ざっとながめてるだけだが (量が尋常ではないんですよ)．
• さっさと午後の仕事に取り組む． ちょっと試しにいくつかの樹種でパラメーター推定計算． なーんか収束までえらく時間のかかるやつがいるなぁ． どうも胸高断面積の数値だけが 他のものより大きなケタであることが原因なんだろうか．
• ばからしー …… と思われるかもしれないけど， これは経験的にそうなのである． 試しにそこを意味もなく 1/100 にしてみると …… ほーら， すぐに収束するようになった． ええっと， これは mlfitting 内部の どういう計算に影響されてたんだっけ．
• 1335 さて， これで
  • 無次元モデル (定数とサイズだけ)
  • 一次元モデル (定数とサイズと上側胸高断面積)
  • 三次元モデル (定数とサイズと開空度らしきもの)
  の推定結果が出そろった． 当然ながら， 下にいくほど尤度 (モデルによって説明されている度合) は高い． というのも， 上で書いているように 一次元的な Bad というのは サイズの単調な関数に過ぎないからだ． それにたいして 「開空度らしきもの」 はいま少し複雑で 「同じサイズであっても異なる値」 を取りうる． これは「うまく」計算してやれば， 樹木の成長量の説明に寄与したりするわけだ．
• 三次元モデルのほうが良いのは 半ば当然なのだけど， 次の問題はこれが「十分に良いものであるか」 という問題だ． モデルごとに 全樹木個体わたって 対数尤度の和を計算してみるか．
• ちょっと昨日の「ぎょーむ日誌」など直したりしていて 計算がおわったのが 1430． うーん． どうだろ．

無次元モデルに比べて一次元モデルの 尤度は 100 以上上昇している． パラメーター数は 24 増えただけだから， AIC 的に評価するなら， これを使うことは 「よりよいモデル選択」となる． 尤度比検定でもよゆーで無次元モデルは 棄却できるだろう． まぁ， 三次元モデルは一次元モデルと同じパラメーター数なのに さらによい， ということで頑張ってるわけですかね．
• 正体不明の三次元モデルは まだまだ改良の工夫の余地ありそうだな．
• 昨日考えた計算結果はこれで出た． 面白みがないけど， これで一応発表にはなるか． どうせ 15 分だし． あとで時間あったら結果を もうちょっといじってみるか．
• とりあえず ``あーきてくちゃーグループ長'' 甲山さんにあーきてくちゃーな計算結果らしきものが 出ました， というメイルを送る．
• 金曜日にやらされる発表「練習」の準備のため， またしょうこりもなく HTML スライドを作る． 半分は先月の使いまわしになるだろう． この作業は闇ネットを通じて ThinkPad 内で行うので， 手すきになった BookPC は ナゾの三次元モデルの性質を解明させるべく， ごりごりと計算を続行させてみる．
• 三次元モデルの尤度はさらに上がる (上の表はその数値を使って修正済)． その調整に使ったパラメーターは， そう， 多く見積もっても 5 個かな．
  • Voxel サイズ (1.0, 1.0, 1.0) (m)
  • 樹冠サイズ係数 0.5
  • 樹冠密度係数 3.0e-2
  ちなみにこのようなセットで樹冠を生成して 三次元光分布を計算した場合， それは開空度の「観測値」とはまるっきり 合致してない のである． 私にとっては 少しばかり衝撃的である．

[でたらめモデル]
しかしこれが樹木の生長速度をもっ
ともよく説明しているのである．あ
えて議論するなら，観測値は地表 1m
に設置された観測点のものである
(縦軸もそのようにして計算した)．
しかしパラメーター推定するときは
そんな低いところにモデル観測点を
おかない (各樹木のてっぺんにつけ
る)．ということは ……

さらに樹冠モデルはこうなってしまっている．

[きのこの森]
…… この樹冠がもっともよく ……

• 1755 研究所発． 1814 東京モノレイル浜松町発． 1850 帰宅．
• 晩飯食ったあとに少し濃い目に淹れた紅茶を飲んでいても 眠りこんでしまいそう． うう． いかんいかん．
• 眠気覚ましもかねて食糧調達に出る． 米がきれかけている． いつものごとく， 一番安い米 10Kg 2990 円買ってしまう． これで一月もつはずだから， 一日約 100 円か．
• 何をやっても脳内ホワイトアウト状態なので， 早々に寝ることにする． 予定の計算結果も得られたことだし， よしとするか． 2200 すぎに就寝
• 今日の食卓
  • 朝 (0715): 米 0.6 合． ホウンレンソウ茹でて 麺つゆとキムチの素で味付けしたもの． エリンギをオリーブオイルで炒めて 塩胡椒で味付けしたもの．
  • 昼 (1230): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1910): 米 0.6 合． 朝の残りにニラを加えて野菜雑炊にする．