ぎょーむ日誌 2000-02-23

2000 年 02 月 23 日 (水)

• 0950 起床． うう，眠い． 晴れ．
• 尤度比検定のプログラミング， 完成まであと少しだ， と盲信してぢりぢりと葡匐前身 ……
• 最尤推定プログラムにパラメーター読ませて 尤度だけ計算させるオプションをつけるためだけに， 入力処理の部分を大幅に作り直していく．
• 尤度比の確率分布を計算する部分以外は， 尤度比検定のプログラムがほぼ完成しつつある． これの主要コンポーネントは 今月の中ごろにようやく稼働した 最尤推定プログラムである．
• 結局， この部品の入力部の改善に時間の大半を取られた． 要するにこの 数値偏微分→逆行列→Newton-Raphson 法の 計算モジュール複合体の動作を 制御するくだりをいじっていただけなのである． そんなに複雑なことやらせたかったわけではない． 「てきとーなパラメーター初期値と データーセットのファイルの場所を教えてやるから， 最尤推定しろ」 「あ， データーセットは標準入力から渡すかもしれん」 「それから尤度比検定にも使いたいんで， 帰無仮説のパラメーターセット受け取ったら 尤度だけ計算しろ」 …… たったこれだけのことにすぎないのだが．
• なぜ時間をくっているのか?
  • 私の設計能力が低い： まぁ， 何と言いますか， いろいろなレヴェルで 設計がズサンなんです． 尤度比検定のこと考えて 最尤推定プログラム作ってない， 引数処理やファイル入力用の 汎用クラスをあらかじめ準備してなくて プロジェクトごとに作り直してる， などなど．
  • C++ は型チェックがうるさい： うるさいというより いつも愛用してる Perl みたいに てきとーに取り計らってくれない …… コマンドラインから さまざまな型の引数与えるときに めんどーなわけです． このあたりは継承やテンプレイトで吸収したり， 引数用に 「Perl スカラー変数みたいなクラス」 作るべきなんですが …… うーん， PipeTree では ある程度それを実現してたんだけど， まだまだ試作段階なんだよねえ．
• さーてと， あとは尤度比の確率分布の計算だが …… 尤度比は極限でχ二乗分布で漸近する …… ってのを， 最初は使ってみるかな． 安直ではあるけれど． 自由度はパラメーター数 - 1 で良かったのかしらん …… おっと， 学校に行かねばならぬ時間だ．
• 1430 自宅発． 1445 地下鉄琴似駅発． 1500 大通り公園近くの銀行で ThinkPad560E の代金を 振り込み． 1520 研究室着．
• 1530 PRA 輪読会． 今日は第 2 章で担当は鍋島さん． 1830 ごろまで続く．
• 久しぶりにホーエルの「入門数理統計学」に 目を通してみる． 手持ちの統計学の本の中では 最も平易に尤度比検定について説明してあった． なるほど帰無仮説のパラメーター数を 自由度とするχ二乗分布に漸近するわけね． たしかに言われてみれば， そういう状況について調べてるわけだし．
• またプログラムのデーター構造を 書き直すところが出現． とほほ．
• 2105 研究室発． 2130 JR 札幌駅発． 晴れている． それゆえにすごく冷える． 2220 帰宅．
• 帰宅してすぐに料理に取り組む． 少し休んだりすると， ジャガイモなんかの皮をむいたりする気力が 失せるんで ……
• プログラミング続けようとしたものの， 眠くなりつつある．
• 今日の食卓
  • 朝 (1030): ホーレンソウ・ネギ・エノキダケ・卵 の炒めもの． あ， それからワカメ (1 パック 190 円) に酢醤油かけたもの． おかずが二品も． なんと豪勢な …… って昨日の昼・晩と 同じモノ食っております． はい．
  • 昼 (1320): 朝と同じ． 学校で食べようと弁当にしておいたのに， 仕事にハマってしまい， 自宅で弁当を食べるというマヌケなことに ……
  • 晩 (2330): ジャガイモ・ニンジン・シイタケ・豆腐 の味噌汁． だしは煮干と昆布でひく． やはり煮干だしはうまい． 昼の残りのワカメもついたので， これで一汁一菜 というやつを実現してしまった．