ぎょーむ日誌 2000-12-21
2000 年 12 月 21 日 (木)
- 0700 起床.
5 時間ぐらい寝た.
- 朝飯・弁当の準備.
朝飯.
コーヒー.
シャワー.
- 0800 自宅発.
今日も曇天.
ローソンで電気代支払う.
2027 円.
0819 東京モノレイル流通センター発.
浜松町の気温 8 ℃.
0840 研究所着.
- 以前から考え続けている最尤推定プログラム
mlfitting の機能拡張について検討する.
これは推定に関する機能を強化するのではなく,
推定した結果を図なんかにする補助的なものを
考えている.
-
森林総研
の新山さんから毎木データベイス化作業の
委員会でアイデア出してくれない,
という依頼.
呪われ E*cel でバグ混入な森林総研データには
苦闘の記憶ばかりあるので,
うけたまわる.
ここはリレイショナルデータベイス (RDB) 化,
だろうな.
勉強しておかなくては.
- 昨晩,
九大山岳部部
メイリングリストに
現役リーダーが遭難対策に関する
阿呆な見解を投稿したんだけど,
それに対して顧問の中溝先生が
珍しく厳しい返信.
冬山遭難対策か.
そうだ,
そろそろ装備の点検やらんといかんな.
- Pasoh の光環境データないんから
ボルネオのそれで代用しようか,
という甲山さんからの提案.
いつもながら曲芸じみてますなぁ.
とはいえ,
熊本県立大の
山田さん
(面識ないが)
が感光紙法で取ったという光データと
精密毎木調査の結果というのはちょっと見てみたいな.
- 小川開空度データのほうは ofm-ml 上で
情報やりとり続く.
柴田さんの尽力で話が進みつつある.
- メイルのやりとりで中断されてしまった
mlfitting 機能拡張の検討であるが
……
うーん,
--plot mode というのを作ろうとしてるんだけど,
例によって入出力をどうしようかと迷う.
- 答えが出ないまま昼飯の時間.
- 昼飯食ってから郵便の投函に出る.
昨晩書いた共立出版への住所変更届である.
天気,
晴れ.
今日は冬至だけど気温は高い.
路上は昼飯のため出た人々で混雑している.
浜松町は基本的におやぢな街である.
しかしながら,
昼間は若い女性がたくさんいる.
朝のモノレイルから道路を俯瞰するとわかるんだけど,
都心部における女性雇用力には驚愕させられる
(とはいえ,
ホントに必要な需要なんだろうかとも思う).
となりのコンヴィニエンスストア「ポプラ」
で 80 円切手購入.
ちょっとあたりを見回すとポスト発見.
なぜかは知らないが
この周辺は郵便ポストの
密度が異常に高い.
- う.
昼からもメイルのやりとりが続く.
プログラム改良の考えがまとまらん.
- ふと,
われにかえるとネット上で
二次曲面・楕円放物面について
一生懸命しらべている.
えーと,
何でこんなもん調べてるんだろ
……
思い出せん.
- かなり時間が立ってから,
3 次元的な内挿を楕円放物面でやるとしたら,
いったい何点のデータがあれば
曲面の形が決まるのか調べていたのであった.
地形データ補間作業で GMT を廃絶したい,
とういアレね.
ようやく思い出した.
ボケてるなぁ.
あ,
でも考えてみればこういう種類の内挿点ならば,
無理に given な点を通らなくても
いいわけだから
……
- おっと,
考えてみればこの問題自体,
直面してるそれとは何の関係もなかった.
いやはや.
- 逃避から戻って来て,
mlfitting のコードと向き合う.
一から全部書き直したほうがマシだなぁ
……
思いつつ,
--plot オプションを指定すると
flag が立ち output_plot 関数を呼ぶように
作りかえる.
- C++ の場合,
どうも徹頭徹尾クラスという観念で固めていかないと,
どうにもプログラムコードが
わかりにくくなるような気がする.
以前に書いたように複雑な経緯を有する
mlfitting は多分に FORTRAN 的・BASIC 的な
気質が残るコードに
Parameter や Matrix といったクラスたちが
混入していて見ばえが良くないのである.
- 1730 output_plot 関数の中身はまだ書いてないけれど,
これまでの版とのズレを最小に抑えつつ
Plugin 関数のつけかえといった難問は解決できたので
帰宅準備.
今日は早めに帰って,
平和島駅近くの床屋に寄ってうっとおしくなってきた
髪を切る予定なのである
……
と思ったら,
北日本某大学の院生よりメイル.
内容はむろん mlfitting で
パラメーターの最尤推定がうまくできません,
というものである.
できの悪いプログラムを作る愚か者は
その愚かしさに見合った罰を受けねばならない.
- さてそのデータファイルとパラメーターファイルを
当方の mlfitting にかけると
……
たしかにあさっての方向に 3 つのパラメーターが
トんでいきますな.
- で,
ちょっと試行錯誤するとうまく収束する
初期値セットが得られた.
このあたりは
思考による解決というより
経験に基づく直観によるのかもしれない.
なんとならば,
最尤法ふくむ非線形パラメーター推定の世界は
基本的に不合理で,
そう呪われているからである.
- メイル返信.
で,
帰る準備.
今日は某地域生態系学の研究室は忘年会が
始まっている時刻なんで,
今晩は追い撃ちがかかる公算は低いためである.
- 1815 研究所発.
1834 東京モノレイル浜松町発.
なーんか腹減ったな.
1940 帰宅.
- 晩飯食ってからメイル取ると,
ちょーど研究所出た時刻に発信された追い撃ち一通.
すばやい.
えーと
……
あのパラメーター初期値セットでは
うまく収束しませんでした ?
うーむ,
謎.
- 謎が深すぎて 2 時間ほど寝てしまった.
目をさましてココアなど飲む.
- ああ,
時差が.
ま,
謎解決は明日でいいか,
と妥協的な態度.
東京流刑地では
時差拡大が気になるんで夜に仕事できん.
はななだ不都合.
- 今日の食卓
- 朝 (0720):
米 0.3 合.
豆腐・ワカメ・切り干しダイコンの雑炊.
だしはサバぶし.
味噌味.
- 昼 (1215):
弁当.
米 0.8 合.
ホウレンソウおひたし.
- 晩 (2010):
米 0.7 合.
雑炊.
昆布・切り干しダイコン・高野豆腐.
なぜか血迷ってカレー味.
残念ながら,
あまりうまくない.