ぎょーむ日誌 2001-11-30

2001 年 11 月 30 日 (金)

• 0650 起床． どうも改良中のモデルが気になってよく眠れなかった．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー．
• 0740 自宅発． 曇り． モデリングのこと考えてたので 間違えて 0747 羽田行き急行に乗ってしまった． 本線・羽田線の分岐である 京急蒲田で降りて後続の特急を待つ． 特急に乗って 0755 京急蒲田発． 普通に乗り換えて 0818 京急上大岡発． 0835 研究所着．
• う， もう午前終了か． ここまで， 観測された a 齢茎のサイズ分布と a+1 齢のそれを比較しながら， a → a+1 の間に生じるサイズ依存的な死亡率の関数型を データから推定する方法を考えてた． なかなか難しい． 単純に切断してしまえ， というユーワクあるんだけどね． しかし茎のサイズ分布って対数正規分布みたいなんだよね． つまり分散が x に依存するやつだ．
• いやまてよ …… 対数正規分布にこじつけた場合， x ではなく log(x) の世界では期待値と分散が完全に独立なのか． このへんに小細工の余地あるか?
• ひとつのアイデアとしては， a → a+1 となるときに， 分散が小さくなるように (どれぐらい圧縮すべきかは観測データから憶測されるはずなんだが) 分布の下側が切断する (これはかなりがさつな手法だ)． このやりかたは工夫すれば分布のカタチに依存しないな． たぶん． これでいくか? ドロくさいのはしょうがないとして， こういうモデリングで「おのぞみの結果」とやら一 式を算出できるかどうか， もうちょっとデータを見直したほうがよさそうだな．
• 昼飯食ってると 私の作った mlfitting というパラメーター推定ソフトウェアの製造物責任を問うメイルが …… 北大の浦口さんからのものである． 近ごろは， ベータ分布を使ったパラメーター推定に取り組んでおられる．
• さっそくプラグインモジュールのコード (これは呪われ C++ で書かないといけない， という点でかなり使用者非友好的である) やパラメーターファイル， データなど拝見してみる．
• 午後はこの浦口コードを少しばかり手直しなどして， 1600 ごろにようやくそれらしい推定結果が得られるようになった． 問題点は …… 一週間ほど前に私が出したいんちきメイルが原因だったのである． いろいろと． ベータ分布推定用の プラグインモジュールそれ自体は かなりよくできている．
• 私が見落としていた問題点は …… たとえば， ベータ分布の平均値は常に (0, 1) でなければならないので， これが定数ではなく何かの関数である場合には， 関数型に注意を払う必要がある．
• あるいは， ベータ分布 z^(a - 1) (1 - z)^(b - 1) / Beta(a, b) において
  • (平均) = a / (a + b)
  • (分散) = a b / ((a + b)^2 (a + b + 1))
  となってるんだけど， 計算してみればわかるように， つねに (平均) > (分散) が成立している (ベータ分布のパラメーター a, b はどちらも正である)． このあたりも気をつけないといけない． うーん， 勉強になる．
• とりあえずベータ分布を使った推定は私にとって初めてなので， つい嬉しがって推定結果の図を出してみたり ……

[初めてのベータ分布推定] 浦口さん が苦労して開発されたベータ分布用の mlfitting プラグインモジュールである plugin_function.so を横取りして (極悪ですなぁ)， 計算させてみた結果． ただしデータはまだ調整中のもの． 縦軸の値がベータ分布に従う． 曲線はその分布の平均で (1 - exp(-a * x)) * exp(-b * x) という関数型である． この (a^, b^) を最尤推定してみた (もひとつ分散を決めるパラメータも含めて)．

何度もいうように， このベータ分布推定用モジュールのコードは浦口さんが 工夫したものであり， 私はそれを借用して試験運転してるだけなのである．
• それにしても mlfitting の プラグインモジュール開発の面倒さはどうにかならんものかな …… これはまったく解決策を思いつかない． この部分だけ外部のスクリプト言語で記述してよいようにする? ただでさえ遅い計算速度がさらに耐えがたいものになりそうだけど． いやいや， Perl とかで尤度関数を書けるようになったからといって， 値がヘンな領域に飛んだりするのを防ぐことできないか．
• 1750 研究所発． 普通に乗って 1800 京急杉田発． 快特に乗り換えて 1810 京急上大岡発． 普通に乗り換えて 1830 京急川崎発． 1900 帰宅．
• 夕方から茎の死亡率関数の推定方法を考えているんだが， なかなかまとまらない． 「Variance が減衰するよーに」 パラメーターを選ぶのは安直そうなんだけど (しかしこの推定値をまぢめに計算するとなると約 100 万項 の非線形方程式を解かないといかんな)， しかしあまりにも露骨に現象論的すぎる．
• ふーん， PerlによるExcelファイルの読み取り／書き込み か …… Spreadsheet::ParseExcel というような解毒 Perl モジュールが 開発されつつあるらしい． とはいえ， 私はしばらく使わんだろうな． うけとったおひすファイルが研究に必要なものであれば， wvware とかでテキストに変換してから …… あれれ? wvWare とか xlHTML のサイトがいつのまにか無くなってる? 謎． れどもんどの圧力か．
• 今日の食卓
  • 朝 (0710): 米 0.6 合． タマネギ・ニンニク茎・エリンギの炒めもの．
  • 昼 (1220): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1930): 米 0.6 合． 朝の残りで雑炊．