ぎょーむ日誌 2000-02-10

2000 年 02 月 10 日 (木)

• 0550 起床． 0650 ごろ明るくなる． 晴れ． 10-20 cm ほど雪が積もったか．
• 朝から昨日のぎょーむ日誌を書いてみたり， 会議のメモに手を入れたりする．
• バグ取り続行． 早起きしたせいか， 現実逃避のためか， 10 時から 11時までまた寝てしまう．
• あー， ばかばかしい． バグだバグだと思ってた逆行列計算関数には 実はまったくバグがなかったことが判明した． 計算試験用に私が与えた行列がたまたま Gauss-Jordan 法では処理できない 行列だったのである． すなわち， Micro$oft ふうに述べるならば， GJ 法の「バグではなく仕様」ということになる． まあ， 最尤法の計算では行列の要素がすべて整数になり かつそれが「解けない」組み合わせである 確率は限りなくゼロに近いだろうから 大丈夫だろう …… 抜本的には LU 分解を 使うようにしなければならない． いかなる場合でも必ず解けるようにする関数を 作るのは結構めんどーだろう． よって後回しだ．
• 逆行列問題を解決しても， うまく計算できない． 今度は数値偏微分行列の問題か． 数値微分の数値微分を二重に計算しなければ ならないところがあるので， テンポラリーな Parameter 型変数を使う アルゴリズムに変更してみる …… うーん， まだうまくいかないなぁ． Parameter Class まわりの細かいバグを つぶしていく．
• 昼飯を食いながら， なんとなくラジオをつけてみる． 国会中継やってる． たぶん 「統計データーを出すと具体性がでて良い」 というような 「質問マニュアル」があるんだろう． 数字がやたらと質問の中に出てくるんだが …… なかなか悪どいぞ． 「君にも出せる有意差」(cf. 粕谷ピンク本) というほど高尚なものではない． 「日本の『食糧自給率』は 60% だから， 輸入が止まると国民の 4 割が飢える」 …… まぁ， この程度だ． 論理的な議論というより プロパガンダの表明合戦なんだろうから， これでも良いのか．
• 昼飯食ってもバグ取りは続く． この最尤法のプログラムは C++ で書いてるんだけど， こいつは結構ややこしい決まり事がいろいろとある． たとえば文法的に正しい書き方がわからなくなったとする． どうしたらいいか． 私の場合， ごく短い実験用コードを書いてみる． そういう実験用コードを生成するためのプログラムを あらかじめ準備しているほどだ． あいまいな点を調べられるようなプログラムを作り， コンパイルして実際に動かしてみる． 迂遠なようだけど， 確信もてる書き方を 意外と効率よく見つけ出すことできる．
• 1630 自宅発． 1654 JR 琴似駅発． 1710 研究室着．
• 1715 より隔週の怪しげセミナー． 担当は清野さん． アロメトリーの話のようなので 今回はまともなのかと思いきや， やはりキワモノ的であった． 著者のよくわからん見解によれば Corner の法則 (ついに何がホントのCorner の法則かわからんかった) は「当年生シュートの太さと葉面積は正比例」 ということになったり (パイプモデルみたいですな)， 主成分分析と共分散分析を混合してみたり， どうも恣意的なんじゃないかと疑われるよーな 種間比較やってみたり …… ま， やりたい放題なわけです． 反パイプモデル闘士の牧さんが 「こんなんパイプモデルとは関係ありませんよ」 とすさんでいました．
• Ngraph の魔術師・田代さんによると (参照： かとーさんの掲示板 )， Ngraph でも塗りつぶし棒グラフが書けるようです．
• 最尤法プログラミング続行． 完成がどんどん遅れていくにつれ， この仕事をとりまく情勢が変化していく． いまや元祖発注者である浦口さんは 私の連日の納期延長の言い訳など もはや全くアテにせずに修論どんどん書き進めるし， 来週に研究の中間報告する牧さんは 「パラメーターの推定値がなきゃ， そういう発表にします」と 早くもリスクヘッジを考え出すし， 進んでいない「しれとこモデリング」発注者の 西村さんには 「いやー，いま作ってる推定プログラムができれば， 針葉樹の成長モデルなんてすぐにできますよ」 と 未完成製品に面倒を押しつけるし ……
• くそう面倒だ． しかし過程を追ってバグをつぶしていくしかない． すでに最尤化の計算が走り始める状態には なっているんだが， パラメーターセットが 明後日の方向にすっとんで行く という問題にわずらわされている． どうも今回の新機軸である数値微分のところに 問題ありそうだ．
• 昨日フォントを入れ換えた牧さんの <蒼龍>． 新しく導入された 「丸文字フォント」を設定した Netscape Communicator で この「ぎょーむ日誌」など表示してみて はしゃいでおられる． 「女子高生が書いた日記みたーい」
• 2330 研究室発． 2350 JR 札幌駅発． 0020 帰宅．
• 終り無きバグ取りは続く．
• 今日の食卓
  • 朝 (0700): 昨夜の残りのダイコン・豆腐・キャベツの 味噌汁．
  • 昼 (1330): 朝の味噌汁と飯の残りで雑炊を作ってみる．
  • 晩 (0050): ダイコンの雑炊． 帰宅してみるとおかずはないのに， 冷やご飯だけがちょっと残っている． 昆布でだしをひき， 拍子木に切った大根を煮る． 味噌と冷やご飯を投入して 雑炊とする．