ぎょーむ日誌 2009-01-30

2009 年 01 月 30 日 (金)

• 0830 起床． 4 時ごろにいったん起きて， 2 時間ぐらい眠れなかったのだが， また寝てしまって起きたらこの時刻だった …… うーむ， 不調パターンがつづきますな． 朝飯． コーヒー． 0930 自宅発． 晴． 0945 研究室着．
• さっさと公開講座 (私の分担部分はいわば「ゐんばぐす Why & How」) の準備しなければならないのだが …… ぢゅぎょー準備ばかりでイヤになってしまったのか， たまには勉強しなくてはとかつぶやきつつ， 鍋嶋・石井解説文 「樹高成長の制限とそのメカニズム」 を勉強． これはよくまとまっていますな． そして樹高成長なるものは依然としてナゾばかり， というかんぢだ．
• 昼飯．
• …… いかん， イヤな雑用メイルかきに着手しなくては， と考えていたのだけど， なかなか着手できなくて (とうめんの仕事にも雑用にも関係しない) プログラミング言語の最新事情あれこれについて調べてしまった． アタマぼけぎみですな．
• とりあえず B 棟 1F のネットワーク室に現状を偵察にいく．

[ネットワーク室といっても]

こういう巨大な「スイッチ」(手前右) のハコが鎮座してるだけで ……

[動いてない server 機なんかが]

ぞろぞろあるだけ． じつは机の下の一台だけがまともに動いている．

[このへん片づければ]

A 棟の server 機器を疎開させることができるかも． なんだか北大内の地球環境のネットワークに関するはりがみがあったんだけど， A 棟と大計が光ケイブルで直結してるのか， 情報教育棟から B/C 棟と分岐しているのか， よくわからないかきかただった．

• とりあえず情報基盤センターに状況・実現したいこと・要望したいこと を最小限にまとめたメイルを送ってみる． どうなることやら． とにかく交渉あるのみ， だ．
• ついでにもう一件ネット雑用． 研究室の (あまり使われていない) 無線 LAN アンテナの調査 …… うーむ， とくに問題ないように思うけど．
• 統数研公開講座ねたというか架空な例題を考えてみる …… これが思いつかなくて， 準備がとどこおっているんだよね．
• 講義のーとにのってるやつは簡単すぎるような気がするし， かといって複雑な例題は教える側にとっては自滅的というほかないし ……
• そうそう， 昨日 car.normal() 解説文を改訂したんだけど， MCMC とは何かといった説明がわかりにくい， と指摘された． もともとは

  尤度が良くなる、 つまり観測データへのあてはまりが良くなるように パラメーターを変化させていくという目的は ベイズモデルの MCMC 法でも最尤推定法と同じである。 しかしながら、 MCMC 法では「尤度が良くなる方向」にばかり進むのではなく、 何か特定の値を探索するのではなく尤度が良いときの 「モデルの状態」にむしろ興味がある、 というところが異なっている。 この手法を適用することで、 ベイズ統計モデルで定義された事後分布そのものを直接に推定はできないが、 事後分布からの多数のランダムサンプルが得られる (以下、サンプリングとよぶ)。 これによってわれわれは事後分布の性質を知ることになる。

  といったかんぢで， MCMC 計算の中での パラメーターの値の 「動き」みたいなのに着目していたんだけど， たしかに上の説明は「何が何やら?」いまいちなかんぢだ． そこで以下のように改訂してみた．

  (MCMC 計算は) 最尤推定法と同じく尤度 (統計モデルの観測データへのあてはまり) にもとづくものである。 数値的な最尤推定法では試行錯誤によって 尤度最大となるパラメーターの値を推定する。 ベイズモデルの MCMC による推定でも パラメーターの値を試行錯誤で変化させる。 しかしながら、 MCMC 法ではこの試行錯誤 (以下、サンプリングとよぶ) で得られたパラメーターの値のセットが 事後分布からの無作為抽出標本となるような、 巧妙な「値の変化のさせかた」のルールが適用されている。 このようにしてサンプリングされた値のセットによって、 われわれはパラメーターの事後分布の性質 (事後平均値など) を推定できる。

  より抽象的な記述になっているわけだが， まあ中途半端に具体的な「動き」なんかを書くより， 「どういった値が得られるのか」 に焦点をあててみたわけだが …… どうでしょうかね?
• てなことを書いてたら， まさにこの 昨日提出した原稿ファイルに関して， めんどうな指示が …… この作業フェイズさえ終えてしまえば， ほぼすべて完了のはず， とがんばってみる．
• しかし最後までは到達せず撤退． 1920 研究室発． 買いもの． 1950 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • うー ……
• [今日の食卓]
  • 朝 (0900): 米麦 0.8 合． 鶏レバのショウガ煮． キムチ．
  • 昼 (1310): 研究室． 米麦 0.8 合． 鶏レバのショウガ煮． コーンスープ．
  • 晩 (2200): 米麦 0.8 合． ダイコン・ネギ・鶏肉のスープ． 鶏レバのショウガ煮． ヨーグルト． リンゴ．