ぎょーむ日誌 2002-03-22

2002 年 03 月 22 日 (金)

• 0655 起床． 6 時間半寝たんだけど 脳内ばてが解消しないかんぢがする． やはり呪われ言語 C++ のせいで疲れてるのかなあ．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー．
• 0750 自宅発． 晴れ． 特急に乗って 0801 京急平和島発． 車中で 「来週はどういう発表すりゃあいいんだろう」 などと今ごろになって考え出す． 普通に乗り換えて 0828 京急上大岡発． 0845 研究所着．
• 朝からいきなり毎木シミュレイターを直す作業． ふたつの森林， 小川 (冷温帯林) ⇔ Pasoh (熱帯林) をいったりきたりしてるうちに， すっかりややこしくなってしまった 計算系統を修繕するため．
• で， 1100 までかかってこの直しばかりやっていた． たしかに数時間前よりは単純化されたけど， まだまだ複雑さをとどめている． しかし生態学会大会前のシミュレイター大改装は， おそらくここまで．
• 試験運転 …… 最大化対数尤度 (MLL) が 4.857e+3? なんかみょーに上がってるな． パラメーターは少し変えたが． ざっと点検してもヘンなところ見つからなかったんで再試行．
• えーと， 現在の小川毎木シミュレイターの場合， データがないところを乱数でごまかす， たとえば樹冠の「カタチ」とか． という方針を採用しているので仮想森林作り直すごとに 推定量がびみょーに異なる ……

[乱数樹冠] これのカタチが毎回変わる． 開空度 (らしきもの) 計算→パラメーター推定に 大きく影響してるみたいだ．

• 現在のコード・パラメーターセットで計算させた場合， 6ha 約 5000 本の個体の開空度を計算するのに約 5 分 (それを観測年 = {1987, 1989, 1991} で計算するので 15 分)， 24 樹種のパラメーター推定に約 1 分で， 合計 16 分ぐらいかかる (BookPC, PentiumIII 733MHz)． メモリの使用量はおよそ 160MB …… まぁ， これらの数値は「樹冠のカタチ」の複雑さを決めるパラメーターなんかに 強く依存してはいるんだけど．
• しかしメモリ使用量はまだ減らせるんではないかな． 今後の課題のひとつとして， データベイス管理システム (DBMS) との連動を考えなくては．
• などと書いてたら推定試行おわった． 最大化対数尤度 4.877e+3． ますます高い (つまりモデルによくあっている)． 改装したときに混入していたバグを直してしまったのだろうか? ともあれ， MLL のばらつきがでかくなる， とわかったんでもう一回再試行してみるか． 午前もちょうど終ったことだし．

[仮想森林内の固定観測点] タテ軸が「明るさ」 ヨコ軸が順位． しかし実測観測値と比べると ちょっと明るすぎる． 上の図のよーに， マガマガしい樹冠を生成してなお …… 実際はどうなってるんだ?

• 計算命じておいてから 1210 トイレでジャージ上下に着替えて 独房群内の人目につきにくい経路を選んで脱柵． いつのまにか曇天． 富岡総合公園 前峰までゆっくり走って登り， 頂上でストレッチング． 反対側に降りて主峰に登る． 谷間の公園はサクラ満開で子供をつれたお母さん多数． 南東展望台をまわってから引き返して 1245 独房群帰着．
• 最大化対数尤度 4.842e+3． やはりばらつき大きいな． 昨日以前よりは高い値ばかりが出ているけれど． 昼飯食ってるあいだにもう一度再計算．
• 最大化対数尤度 4.874e+3． ばらつきでかい． 樹冠の形状に敏感といってよいのだろう …… さて， どうしたもんかね． まぁ， 開空度計算系統とパラメーター推定系統は完全に分離してるから， ある「ある日たまたま得られた明るさ観測値セット」のもとでの さまざまな最尤推定， というふうに示せばいいんだけど．
• 何もとくに思いつかないので， とりあえずここまでのファイルいくつかをこちらに保存．
  • 樹種名一覧: tree_name.txt
  • 森林設定: FOREST.txt
  • 樹木設定: tree_common.txt
  • 確率論的な直径成長モデル: growth_function.cc
  • 推定結果の一例: growth.parameter
• BookPC がヒマそうにしてるんで， 同じ設定でしつこく再計算を命じてみる． 最大化対数尤度 4.866e+3， 4.863e+3， 4.854e+3， 4.867e+3, 4.862e+3, 4.852e+3, 4.850e+3, 4.867e+3, 4.874e+3, …… MLL の分布は漸近的にχ二乗分布に収束， ってことらしいんだけど．
• で， こうやってどんどん計算させてるあいだ私は何してたかというと …… ぎょ−むさぼって引っ越し準備． 東京電力 のペイジで電気代清算申込． しかしこの登録フォームのできが悪くて なかなかいらいらさせられる． しかも確認のメイルを送ってこない． 東京水道局 ， これも同じようなフォ−ムで． こっちはすぐに確認メイル来た． 東京ガス …… あ， ガスに関しては 「お客様番号」なる数列書いた領収書用意してなかった． これは最悪の場合， 電話でどうこうせんといかんのかしらん? 東京ガスの場合 19 時まで受けつけてくれるらしいけど (0570-002211)．
• クレジットカ−ドはネット上で登録情報変更できるようになってたんで， とりあえずそのアカウントだけとって後日に変更することに． 銀行での住所変更手続きはあとから郵送でのんびりやるか． 郵便局に転送依頼 …… ネット上ではダメか． いやはや．
• あてもなく BookPC には計算を継続させつつ， 人間のほうは来週の発表の構成を考えてみる． 当日 (3/28) の プログラム を見るに， S9 データマイニングと同時に 植物生態学系のシンポジウムもいくつか開催されるし， その時間は招待講演で Stuart Chapin 御大しゃべってるから …… やはり聴衆には植物集団の動態まにあなヒトはあまりいない， と考えるべきだろう． そもそもこの S9， 発表者からして動植半々だし．
• データほりおこしシンポジウムでは 20-25 分ほどしゃべってよろしい， ということらしいんだけど …… 先日に弥生講堂でやったみたいなハナシやってたら， ほんとに上っ面の部分だけで終ってしまうな． 小川毎木モデリングはむしろ例題のよーに取り扱い， このモデリング過程を軸とする …… 何か統計学・データ解析法バナシというふうに再構成しないといけない．
• ありゃ， WindowManager (WM) の設定ちょっと変えたら WM がとんでしまった． うーむ， 謎． いま使ってるのは WindowMaker-0.80.0-0vl7． WindowMaker の問題か， あるいは私のかなりいいかげんなカスタマイズが原因か． まぁ， しばらく放置．
• 独房区画内闇ルーター BookPC で最尤推定計算を継続しつつ， 闇ホスト ThinkPad240Z に闇 login して 発表内容 HTML メモを作ってみる． なかなかはかどらない．
• 1815 研究所発． 雨か． 傘無し． そのまま強行突破して 1825 歯科医院着． すぐに診療開始． 今日は上の列の歯石とり． すぐに終る． 歯医者先生， 私の左上の親不知が相当にヘンなかたちで生えてるから まだ虫歯になってないけど抜いてはどうかとの提案を． いや私は今月末で京浜流刑地から仮釈放になって 北辺更正労働地区に移されるんで， と申し上げると， じゃあそちらで治療続けてくださいとのこと． たしかにヘンな歯は取ってしまいたいんだけど， 抜歯か …… あれ， 抜いてしばらくは痛いんだよねえ． 本日の診療費 500 円．
• 普通に乗って 1838 京急杉田発． 快特に乗り換えて 1849 京急上大岡発． 普通に乗り換えて 1910 京急川崎発． 1930 帰宅． 体重 70.4kg． 最尤やせ．
• どうも夜に仕事ができない．
• 今日の食卓
  • 朝 (0710): 米 0.6 合． ブリ・ダイコン煮物．
  • 昼 (1255): 弁当． 米 0.8 合． 朝と同じ． ただしダイコンのみ．
  • 晩 (1950): 米 0.6 合． 朝と同じ．