ぎょーむ日誌 2000-02-25

2000 年 02 月 25 日 (金)

• 0730 起床． 雪． ちょードカ雪． ラジオによると一晩で 33cm の積雪， JR の列車は 9本も運休． 運行状態はたぶん無茶苦茶なんだろう． やはり寒波が来てたそーだ． のんびりすることに．
• のんびりしすぎて 0900-1200 机で寝てしまう． ここしばらく早起きだったためか? うーむ， 間抜けすぎるぞ．
• 列車も相変わらずまともに動いてないようなので， 自宅に閉じ籠って作業することに． 今日は北大は入試の日らしいんで， 近づくとまた予備校に雇われた マヌケ兄ちゃんに入学案内を押しつけられそうだし．
• 札幌市内でも雪のふりかたに差があるらしく， 私が住んでる西区とその北側の手稲区で たくさん積もったそーだ． ああ， つまり山の近くということですね．
• さてプログラム改良の方向だが …… ともあれ曲がりなりにも「種間比較型」の 尤度比検定はできるようになったので， 他のタイプの尤度比検定ができるように 最尤推定部を拡張したい． つまりまた面倒な入力部の改良ですね． あー， たいへんだ．
• いままで， 浦口さんの種間比較問題を前提に設計していたところを， もっと一般的な形式の問題に対応できるように 書き換えていく． すると ソースファイルがどんどんばらばらにされていく． ああ， 大丈夫なのか．
• いまのところ， こういうふうに分割してる．
  • function.cc : 尤度関数を定義する関数． 基本的に， このファイルを取り換えるだけで さまざまな種類の問題に対応できる．
  • main.cc : メインルーチンと 引数処理， そして なぜかデーター読み込み関数など．
  • matrix.cc : 行列と逆行列をセットにしたクラス Matrix の実装部．
  • matrix.h : Matrix の宣言部．
  • ml.h : この最尤推定プログラムの 全モジュールで使われる include 文と グローバル定数．
  • newton.cc : Newton-Raphson 法の過程を 統括している関数が書かれている．
  • output.cc : 出力ストリームの関数をまとめている． って， いまのところ Parameter しかないが．
  • parameter.cc : Newton-Raphson 法で使うのに 便利な関数をもつクラス Parameter の実装部．
  • parameter.h : Parameter の宣言部．
• C++ の STL (standard template library) を使ってて 腹立たしいのは， ちょっと書き間違えただけで コンパイルエラーが何画面ぶんもだーっと 出力されてしまうことだ． いやはや． Perl でプログラムを作っていると 「自分はなんて機知にとんだ書き方ができるんだ」 という錯覚にひたれるが， C++ でコーディングしてると 「あ， おれって世界最悪のプログラマーかも」 とぜつぼーしたくなる．
• いままでパラメーターの集まりは vector<Parameter> Plist という 単なる可変長配列だったんだけど， 新しく Class Plist を定義してみる． その結果， たとえば いままでこう書いていたのが

  	( ( *plist )[ i ] ).Add( deltap * ( ( *plist )[ i ] ).Dfactor() ); 

  こうなってしまいました．

  	plist->PP( i )->Add( deltap * plist->PP( i )->Dfactor() ); 

  …… うーん， あまり良くなったカンジがしないので， なにもかも Parameter に押しつけることに． こういうふうに． やっぱダメダメかな．

  	plist->PP( i )->NewtonChange( deltap ); 

• 最初からきちんと設計すべきだったのか? それとも， 一度作ったものにこだわらず 一から作り直すべきだったのか? Perl か何かでプロトタイプを作っておいてから， 必要な部分だけを C++ 化すべきだったのか? …… どうやっても面倒そう．
• 1630 ごろ ひととおりコンパイルエラーが消えたけど， 最尤化計算を途中で投げてしまうようになった． なんでなんで? 皿洗いしてみたり， コーヒーをいれたりして落ち着こうとする．
• 昨日， ostream の scientific 形式では 有効桁数が変えられんのかとか ぎゃあぎゃあ騒いでた問題は 私の勘違いであること判明． setprecision(int) をストリームに流し込むだけで 簡単に変更できる．
• ラジオのニュースによると， この一日は札幌 2 月における観測史上で 最もすごい積雪だったそうだ．
• ようやく思い通りに動きはじめる． しかし今度はデーターに 「パラメーター減らしたモデル」 あてはめるのが結構むずかしい， と発見する． 考えてみれば当り前のことか …… などと納得したつもりにならないで， さらに懐疑にとりつかれてコードを見直すと …… ほーら， 間違いが見つかった． パラメーターの数を変えたぐらいで， 収束速度が桁違いに遅くなったりするもんか． 莫迦．
• 天才ならざる者は 何事も信用しない 強い猜疑心を持たねばならない， か． いやはや．
• よーやく stepwise 型の尤度比検定プログラムが 完成した． 最後の難関はχ二乗分布の積分． 自由度が奇数で 尤度比が極端に大きい場合 計算がでたらめになることがある， という部分の補正． 時刻は 0100．
• あとは尤度比検定の結果を受け取って描画する， という部分だが …… まぁ， これはしばらくは放っておいてよさそうな気が． ともかく今日の仕事はおしまい．
• いろいろと金の請求が． 部屋代 43000 円． びんぼー人にはぜーたくすぎるか． もっと安いところ探すか． ガス代 3751 円． 水道・下水道代 4032 円． ま， このあたりは毎月変わらない． 電気代もこれぐらいだったかな． 最悪なのは電話代． 7339 円． 当家の電話回線に人語が通ることはほぼ皆無で， これぜーんぶ 計算機と計算機の間のお話のために 支払われるのである． はあ． テレホーダイ入るかなぁ． 少なくとも時間外には 接続したくなくなるかもしれない． ともかく現状のままでは いかん． と思いつつも， 今日もしばしば PPP してました． ネット依存症だなー．
• 今日の食卓
  • 朝 (0800): 味噌汁の残りと ホーレンソウとネギの炒めもの．
  • 昼 (1210): 朝の ホーレンソウとネギの炒めもの を弁当にしておいたのだが， またまたそれを家で食うことに なってしまった． とほほ． 何やってんでしょうね， 私．
  • 晩 (1930): 味噌汁の残りに タマネギ・ネギ・シイタケを 追加してみる．