ぎょーむ日誌 2001-05-08

2001 年 05 月 08 日 (火)

• 0600 起床． 昨晩は 2400 に「消灯」して， 予想したより早めに眠ることができた． おかげで今朝はいつもより眠くない．
• 朝飯・弁当の準備． 弁当の菜を昨晩作っておかなかったので， 今朝になって野菜を切り炒める． 時間はある． 朝飯． コーヒー． シャワー． ひげそり．
• 0720 自宅発． 曇． 0725 京急平和島着． ID カードを忘れてきたのに気づく． しょうがないんで取りに帰る． 時間はある． 0730 自宅を再出発． 0737 京急平和島発． 普通に乗り換えて 0801 京急上大岡発． 0820 研究所着．

[今朝の忘れ物] これが無いと独房群に入れない． 将来は建物内移動すらできなくなるらしい． 明日からホルダーは独房内におき， カードだけ財布に入れることにする．

• いま下請けぎょーむでやってる Pasoh 解析問題． 推定された値はまったく信用できないにせよ， パラメーター推定までのワクぐみは 自動化が終了している． で， その先がまっくらな状態だ (下請けから見ると)． さらにセットで同時に推定される その値は独立ではなく 「こっちが増えるとこっちが減る」 というような釣り合い関係にある． これもよくわからぬところだ．
• 今朝， 歩きながらそのあたりをつらつらと考えてみる． ああいうふうに推定されてしまうのは， かぎられた領域に分布する点の集まりに， これまた「形」がずいぶんと制限されている 確率分布をあてはめよう， としてる最尤推定計算の結果だなぁ． たぶん． であるならば， この特性を逆用すれば， 推定値セットなんかをどう調べるんだよう， という そもそもの問題にも多少の光明はあるわけだ． なるほど．
• メイル書きなど． Pasoh 発注者の甲山さんがやはり 「99% サイズ」 の計算注文だして来たんで， すぐに計算して送り返す． 任意の頻度サイズが計算できるように Perl スクリプトを作っておいたのである (当然ながら)． そして言うまでもなく 95%サイズに 比べて 99% での推定値は信頼区間のハバがでかい． はてさて， どうするつもりか．
• 1100 までまた ssh 問題に取り組む． $HOME/.ssh/authorized_keys を設定することで (たんに host 側 で ssh-keygen した $HOME/.ssh/identity.pub を追加するだけ)， RSA 認証できることはわかった． しかしパスワード無しの (あるいは最初だけパスフレイズ入力の) 接続にはいたらない． どーも， これは机上闇ルーターの BookPC 表・闇ふたつの 顔を持つことに関係ありそう． うう， めんどくさい． たぶん， 自宅闇ネットならもう少し簡単だろう． ということにして， 今日の ssh 研究はヤメ． ふう．
• パラメーター推定のときにパラメーター数を 調節する問題． 一年間 (って換算もでたらめなものだが) にどれだけ樹木が太るか， という成長モデルにいくつのパラメーターが必要か， という問題である． いまのところ， 定数・サイズ依存・明るさ依存 の 3 つがある． そしてこの順番に怪しさが増大する．
• で， 科学の基本として怪しいパラメーターの数は減じたい， ということで， パラメーター数を増減させたモデルを比較するわけだ． この場合だとパラメーター数は 1 から 3 まで 変えられる．
• さて私の自作最尤推定コードは --estimator コマンドラインオプションつけることで， 推定値の数を変えることができる． ということで， 単に 1,2,3 パラメーターモデルそれぞれの 推定値を得るのか簡単である． ただし， 計算時間がいままでの約 3 倍に増大するので イライラさせられるけど．
• 推定計算の試験運転が終わらないまま昼飯． うーむ， 問題は最大化対数尤度が出たあとの処理か． AIC か尤度比検定か． しかしまだその段階でもないような気がする． 推定値なんて， 樹冠モデルや成長モデルを変えただけで ころころ転じていくもんだしね． だいたいモデル選択したところで， その先の処理は示されていないのである．
• でれでれと昼飯食いつつメイルなど書いたりしてると， 昼休みが終わりそうになる． しかし雨が降ってきたので， 今日は外に出られないな． ありがたいことに今日は昼飯食っても眠気がきざさない．
• 2 パラメーターモデル (定数とサイズ依存) と 3 パラメーターモデル (加えて明るさ依存) の 比較をしてみる． とはいえ， 「明るさ」 というのが信用できんので試行段階なんだが …… 471 樹種の中で (明るさ無視モデルと明るさ依存モデルの間での) 対数尤度の差が 1 未満であるものが 230 樹種ほどだった (むろんもの数は樹冠モデル・成長モデルに依存して ころころ変わってしまう)．
• AIC によるモデル選択の規準からすると， これら対数尤度差が 1 未満のものは 「明るさ無視モデルのほうがマシ」となる． というのも， 比較してるモデル間の パラメーター数の差が 1 であるからだ．
• いっぽうで尤度比検定だとどうなるか． 尤度比検定では 対立仮説と帰無仮説の尤度の比が χ² 分布 (自由度は「パラメーター数の差」) に漸近する， というリクツを使う． 自由度 1 の χ² 分布 を右から積分していって 0.05 となる点は …… 3.84 か． これの log をとると (つまり尤度比→対数尤度差と変換)， だいたい 1.3 ぐらいかな． ということで， AIC の場合と同じように， 230 ぷらす何種かにおいて 帰無仮説 (「明るさ無視モデル」) が棄却できない．
• …… ということで， この「モデル選択」の部分は自動化はまだだけど， 処理の方法に関するメドはたった． もし自動化が必要なら Perl スクリプト 一個書けばよいだけだ．
• 夕方まで発注者の甲山さんとのメイルやりとり続く …… なるほど． 全体的には強引に押し切るワケね． 面白くない． されど当方も代案なし． ということで， 残された問題は細かいつじつまあわせばかり， か．
• 1810 研究所発． 雨． 1829 京急杉田発． 特急に乗り換えて 1850 京急上大岡発． 1920 帰宅．
• 晩飯食ったらまた 1 時間ほどオチてしまった．
• 昨日書いてたコロリョフ最後のロケット <N-1> の ペイジ 見つけた． なんかすごいぞ ……
• 生活に「消灯時間」導入しようと思ってるんだけど …… 今晩はうまくいかなかった． 「眠れる」という確信ないと 電気を消せないのである． 消せば眠れる， とは思うんだけど ……
• 今日の食卓
  • 朝 (0630): 米 0.6 合． ミブナ・ニンジン・シイタケの野菜炒め． ミブナはなぜか埼玉県産．
  • 昼 (1220): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1930): 米 0.7 合． 朝と同じ． 干しシイタケ・乾燥ワカメの味噌汁．