ぎょーむ日誌 2005-04-30

2005 年 04 月 30 日 (土)

• 1120 起床． いかん， かんぜんに周期がずれてる． 朝飯． コーヒー．
• 怠業 …… とゆーか， また昨日のカクレ問題 にとらっぷされてしまう． 私のアタマの中の (外部からは観測できない) カクレ状態がそうさせているのだろう．

• かかる Hidden Markov Model と観測データのつじつまあわせの方法， つまり最尤推定法の尤度方程式をどういうふうに構成すればよいのか， よーやくにして考えついた． この問題においては， MCMC 法 によって生成される カクレ状態を sampling しておくだけでは不完全で， 各場所における occurence probability P_i も「記録」せんといかん， とゆーことだ …… これをもちいて「対数尤度の期待値」を計算して， そいつを最大化するパラメーターセットを探索すればよろしい (すなわち EM アルゴリズム )， そういうことだ．
• これにて一件落着． 実用的によい推定結果が得られるかどうかはともかく， 概念的にはこの手順でもって長年の懸案 (たんに時間かかりすぎただけなんだが) が矛盾なく解決した， ということで． そして現用の計算機でそれほど手間どらずに計算可能である …… たぶん， 例によって生態学ではない他分野 (空間統計学とか computer science 系など) では， 誰かが同じような内容を発表ずみである可能性は高そうだけど．
• とゆーか， 上の図は (なぜか hidden layer 内で相互作用のある) ニューラルネットそのもの， と見えてしまうかも． 視覚情報処理とか． たしかにこういうふうに図化してみると， Hidden Markov Model (HMM) と Neuralnet Model のあいだに関連ありそう， というのがわかるな． ふーむ
• さてさて， これをちょっとヒネれば， 苫小牧 leafminer 問題にもただちに応用できるし， 大阪大会でさんざんみせられた いろいろな動物・植物の分布推定問題のたぐいも ちょっとはマシなものになるだろう．
• むろんこのあたりを Bayesian inference にもっていくことは容易である …… が， パラメーターの事後分布に関する簡単な Gibbs sampler を構成できないように見えるんで， Rejection sampling とか Metropolis-Hastings algorithm でやることになり， かえってめんどう感あるかも．
• いっぽうで パラメーターではなく各場所の状態の「事後分布」に関しては， Gibbs sampler が簡単なものになっている …… したがって， 上のようにおてがる MCMC 法駆動な EM algorithm が使えるわけだ． しかしながら， 先日あれこれと検討したように， この方式でパラメーター最大化をやるのはそれなりに手間がかかり， 計算量としてはべいぢあんなやり口のそれに近づいてしまうのかもしれない．
• ひとくぎりついたんで， 1745 自宅発北大構内走． 曇天． やや寒い． 1815 帰宅． 体重 73.8kg． 昨日と比較すると 1.2kg 減 …… まあ， 90 分間も走れば多少は脱水するだろう， というよーな． いずれにせよ， 日単位で生じるのは水の出はいりだけだ．
• 本日 90 分間も走ったのは， 体調・天気のぐあいがよかったのと， 明日はうんどう休養日になる可能性があるため． 予定としては， 明朝は「生活周期ずれ」解消のために 0700 ごろにおきて， (家でじっとしてるとまた寝てしまうだろうから) 大学の研究室にいってうだうだしてみる …… というような． 睡眠ぶそくでふらふらしてる状態では， 北大構内走とか無理かもしれんわけで．
• 買いもの． 晩飯．
• またしょうこりもなく寒天でざーとを作ってしまった． 昨日のはすでに食ってしまったんで． 買いものにいったら乾燥ブルーベリーが賞味期限ぎれ接近で安売りされてたので (75g 150 円ぐらい)， これを材料に． 棒寒天， あらかじめ一時間ぐらい水につけておく． これに上記ブルーベリーと料理酒を加えてぐらぐらと煮たてる． 寒天をよくとかし， 酒精をトバすため． なぜか当家では料理酒だけはすこし 良さげなもの を使ってたりする． 角形の容器がないのでボウルにながしこむ． 水浴であらねつとってから冷蔵庫へ． 一時間もするとかたまった (容積 700cc ぐらい)． 全体が濃い紫色でおもしろい． 味もまあまあ …… 砂糖をほとんど加えてないので， いつものごとくウス味なんだけど． ふと思いついて， プレインヨーグルトであえてみると， これはなかなか良いくみあわせとわかった．
• やはり真夜中すぎても眠くならん．
• [今日の運動]
  • 北大構内走 90 分間．
• [今日の食卓]
  • 朝 (1150): ヨーグルト．
  • 昼 : 食ってない …… ことにするか．
  • 晩 (2030): コンブ納豆あえ． キュウリ． スパゲッティー． キャベツ・ニンジン・セロリ・ブナシメジ・豆腐・タコのカレー．