ぎょーむ日誌 2000-10-05

2000 年 10 月 05 日 (木)

• 0730 起床． まぁ， 6 時間ぐらいまとめて寝れたんで， かなりすっきりはしている． やはり晩飯食ったあとに眠気を解除する 工夫が重要だったのではなかろーか．
• 朝飯・弁当の準備． 朝飯． シャワー． コーヒー飲む．
• 0805 自宅発． どんよりとした曇天． ゴミ出し． 0825 東京モノレイル流通センター発． 浜松町の気温 20 ℃． 0845 研究所着．
• 小川モデリングのパラメーター推定系の 第一段階としてアロメトリーなるものを 推定せねばならぬ． これは 「直径 100 cm の樹木の高さはどれぐらいか?」 ということを憶測するためである． こーいう計算のためには 観測データーと 太さ-高さ関係式を用いて パラメーターを推定しなければならない．
• 自作最尤推定コードである mlfitting を使えば， たしかに樹種ごとに一見もっともらしい 数値は得られる．
• ああ， しかしこの関係式である 拡張相対成長式とかいうのが 生物学的にインチキくさいものであれば， それをデーターに非線形あてはめるという 手法も統計学的に疑わしいシロモノである． コードを書いていて 怒りがふつふつとわいてくる． いつかはこのあたりを 何もかも無茶苦茶に打ち砕いてやる， というキモチになる．
• とはいえ， なかなか良い代替案がないんだよねぇ． 太さ成長と高さ成長を同時に推定する， とか． しかし高さ成長が推定できん． うーん …… くそう， しばらくは旧来の方法使うしかないのか． 100 年前と本質的に なんら変わることのない ……
• 小川 24 樹種を扱う tree_name.pm という Perl モジュールをちょっと改良． 樹種一覧を返す機能をつけた． 一行スクリプトで無理矢理呼び出す場合の 呪文はこうだ．

	perl -I$OFR_HOME/lib -Mtree_name -e \
		'print TreeName->new->NameList( "Abbreviation" )'
	(結果) ACM ACP ACR ACS ACT BTG CLB COC COK CRC CRJ CRL
	       CRT CSC FGC FGJ FRL MLM OSJ PRV QRM QRS STO ZZZ

ZZZ は「その他おおぜい」に割り振られた略号．
• 太さ-高さデーターを mlfitting に渡す一歩手前まできて 一休み． 昼飯の弁当食いつつ netnews 読む．
• 近ごろ fj.* ではやるもの (自分的に)．
  • fj.lang.c: 乱数の生成． めるせんぬついすたー法． 周期は 2^19937 -1 !!! (<cstdlib> の rand は 2^16 ぐらいらしい)． 慶應大学の松本さんのペイジ に解説とコードあります．
  • fj.comp.lang.pascal: Pascal は コンパイル速い言語仕様かどーか， というよくわからん議論．
  • fj.comp.lang.perl: X 上で 選択した文字列の取得． kono 先生の手筋を 一行スクリプトで書くと perl -MTk -e 'print MainWindow->new->SelectionGet()' ． 例によって Perl/Tk である．
  • fj.comp.misc: Apple 株大暴落 というスレッドから 上の Pascal 議論は 派生しているのである． わけわからん．
  • fj.comp.mobile: JR の どの車両なら内部で 電力窃盗できるか， という情報交換．
  • fj.life.health: タマネギが 高血圧に良いとか， 水にさらしちゃ効果ありません とか．
  • fj.os.linux: Linux マシンを クラスター化するには ディスクレスが良いとのこと．
  • fj.rec.aerospace: 飛行機が 墜落したときの 安全姿勢について． 足を椅子にあげるのは 意味ないかも． しかし御鷹巣山的な hard landing だと足切りを まぬがれるかも．
  • fj.rec.games.roguelike: なんと Java 化された rogue である ``jrogue''．
  • fj.sci.math, fj.sci.philosophy: あいかわらず アチラの世界に ……
• まーた， 連絡手ちがいかい． しくじるのはそちらの勝手だけど わあわあ騒ぐのはやめてくれ． おじーちゃん (うるさくないほう) は その横でよく熟睡できるなぁ． 私は修行が足りないのかしらん．
• あるいは Mac の使いかたがわからなくなって 毎日毎日 SE に電話する， とか． 「世界一ィィィ」 なる 虚しいスローガンのもとで すぱこんの建造進めてる科学技術庁な ヒトたちなのに ソフトウェアのマニュアルが自分では読めん ってことはないよな． お役人がそれでは， IT あいちーと騒ぎまくってる そーり大臣にブン殴られますよ．
• ええーい， こういう下らぬ雑念は頭から締め出して， mlfitting のドライヴァー作りに 集中しようとする．
• わりとあっさりできた． ま， mlfitting はすでに一定の汎用性を付与されてるからな． そうかこの日のために， dh_function.cc を plugin モジュール化したり， いろいろと苦労してたんだね．
• あり ? なんかパラメーター推定がヘン． 気分を落ち着けるために となりのコンヴィニエンスストアー「ポプラ」で 野菜ジュース買って飲む．
• なーんだ， タテヨコ逆じゃん． ばか． dh_function.cc を書き直す． make plugins． そう， plugin だけコンパイルやり直せばよいのである． mlf_driver.pl …… ほーら， ちゃんとした結果になった．
• うーん， しかし全 24 樹種のうち何種かで うまく推定できてないな． ACT (ヒナウチワカエデ) とか CLB (リョウブ) とか QRS (コナラ) とか． 拡張相対成長式なるものが破綻しているんだろーな．
• mlfitting に「食わせる」 初期値ファイルのパラメーター値を変えたり， 樹種によってその初期値ファイルを変更するように mlf_driver.pl を書き換えたり．
• ofr1999 で推定に用いていた polytope 法は あてはまりがよくない場合でも 無理矢理に曲線をコジつけていたんだが， 局所的に解を探索することになる 最尤法ではそのような無茶はちょっとやりにくい． 尤度を最大化するような polytope 法コードを作るべきか? という誘惑にかられる．
• Mac 使いのDennis さんが 「おススメのノート PC は ?」 という質問． 後で答えることにして ネット上を調べる． 外国人に勧められるノート PC は 日本では入手が難しい． 英語版 Window$ を正式にサポートしてるところが 少ないためである． たしか Window$ はキーアサインを交換できないので， キーボードも US 配列にしておかなくては …… という条件を満たすのは以下の三社かな ?
  • Toshiba Direct
  • Gateway2000
  • Fujitsu
• ひととおり調べあげてから， パラメーター推定に戻る． 時刻はすでに 1900 すぎ．
• 「モバイルギア」なる小さな計算機で PocketBSD + jvim を使うという うらやましい環境を構築している 京大の長谷川さんから， US 鍵盤の CE マシンならこんなにありますぜ， という悪魔の誘惑のごときメイル． ここでそれに乗ってしまっては 金も時間もどんどん失われるので， それらのマシンは現在は入手困難であるとか 電池が 3 時間しかもたないからダメですとか 理由をつけて手を出さないことにする．
• 最尤法の強みで， 推定に用いる確率分布を取り換える． これも dh_function.cc を直すだけでよい． 観測値が各サイズにおいて正規分布ではなく (正規分布なワケないだろ! と正規分布を憎むワタシ)， 対数正規分布であったなら ? …… あ， うまくいった． 問題になっている樹種だけでなく， 他の樹種のパラメーターをより短いステップ数で 収束するようになった．
• えーと， 対数正規分布である， という仮定の正当化は ……
  • 樹高とは正の値である と仮定している． 正規分布仮定だと マイナス値の樹高も考えている ことになる．
  • 低い樹高は正確に測定される一方で 高い樹高は不正確になりがちである と仮定している． 正規分布仮定だと 誤差はどの樹高でも一定である．
  • 樹高は過小推定より過大推定されがちである と仮定している． 正規分布仮定だと これらは同程度である．
  一見どれももっともらしいが …… うさんくさいですな．
• そもそも「樹高」などという定義の明確ではない， 生物学的意味の明らかではない， ノイズの混じりやすい， しかも 測定困難な指標を 我々はいつまで使い続けるのだろうか． いまから 100 年後の生態学者もまた樹高を 測定しているんだろうか． あやしげなアロメトリー式にデーターを こじつけているんだろうか．
• ま， 今日はここまでにしておこう． 2000 研究所発． 浜松町の本屋 dan で少しうろうろ． 2025 東京モノレイル浜松町発． 腹減ったなぁ． 2104 帰宅．
• …… 晩飯食ったら寝てしまった． おそらく 2200-0100 の間か． とほほ． これはこのあと眠れなくなるパターンだ．

[深夜に空腹になったんで]
小麦粉 100g ・ベイキング
パウダー・塩少々でかよう
なものを焼いてかじる．う
う．貧乏人は麦を食え．麦
ちがいか．

• 上の写真のようなものをかじりつつ， これからどうなるんだろうね， と思う． なんとなく明日の朝飯の準備．
• 今日の食卓
  • 朝 (0750): 米 0.6 合． 味噌汁の残り． ピーマン・ニンニクの茎・マッシュルームの 炒めもの． カレールウひとかけらで味つけ．
  • 昼 (1215): 弁当． 米 0.7 合． ピーマン・ニンニクの茎・マッシュルームの 炒めもの．
  • 晩 (1930): 米 0.5 合． 朝の残り．