ぎょーむ日誌 2009-02-23

2009 年 02 月 23 日 (月)

• あ， これは昨晩の雪の新千歳空港だ．
• 0500 起床． 東京西麻布． 3 時間弱ねてしまった． 本日 1000 から公開講座の準備． はい． まだできてません．
• ぢりぢりと準備のつづき ……

• なんとかできたよーに思えたので， 0900 宿発． 雨． 0910 南麻布の 統数研 着． 2F 講堂の公開講座 マルコフ連鎖モンテカルロ法の基礎と実践 会場へ． 公開講座を支えてくださってるスタッフの皆さんにご挨拶．
• 伊庭さん はまだみえていない …… 開始時間直前にあらわれ， 教室前方で PC の準備． そして立ったまま朝食のおにぎりとお茶を． なんだかスゴい．
• 受講者は 80 名ぐらい． 北海道からくるとちょっと異様に思えたのは， マスク着用なヒトが多いとゆーあたり． スギ花粉がすでにとびちってるようで．
• 公開講座開始． 1000-1100 伊庭さん「アルゴリズムとしての MCMC」． 内容的には 岩波「統計科学のフロンティア12 計算統計 II マルコフ連鎖モンテカルロ法の基礎」 の p.27 あたりまで． Java によるイジングモデルのデモとか． 最後にベイズと MCMC の関係について． MCMC は「統計学の範囲に限っても階層ベイズへの応用だけではない」 わけだが， 「階層ベイズ + MCMC のセットが当たったのは事実 (+ WinBUGS の 3 点セット)」． Metropolis-Hastings 法は， ちょっと Gibbs sampler よりな Metropolis， といったハナシも．
• 1100-1230 久保「ベイジアンモデリングと MCMC」． 生態学ねたではない架空例題で階層ベイズモデルを説明しようとしたのだが …… いまいちだったか? 90 分のうち 75 分あたりで終わってしまったのだが， その後に質問がいろいろと出た． いきなり「ベイズでモデル選択するには」 といったちょー難問だったので， いささかの迷いもなく伊庭さんに回答をお願いする． 伊庭さんの回答の結論もやはり「難しいですね」 ということなのだが， そこにいたるまでにいろいろとありがたいハナシが聞けるのである．

[久保説明風景]

受講者内の知人のかたがとってくださいました． どうもありがとうございます． うーむ， いかにもアヤしい講師だ …… しゃべってる本人にはちょっと「へろへろ」疲労感あり， コトバつまりぎみ． 辞任したさきの財務大臣 (この時事ねたは一年後には忘却されて意味不明だろう) のキモチが少しわかった． 私はカゼ薬も酒ものんでないけど．

• 昼休み． 昼飯くいつつ午後の投影資料のチェック．
• 午後の部． 1330-1430 久保「R と WinBUGS の使いかた」 ということで午前にやった「打率の推定」 階層ベイズモデルを推定計算させるための ひたすらコマかいハナシ． 午前・午後の久保のハナシの投影資料・ファイルは こちら で公開してます (配布資料はかなり修正する必要があるのでまだアップロードしてない)．
• これまた 60 分もち時間のうち 45 分ぐらいで終わったのだけど， 質疑応答があって， 伊庭さんからのまとめコメントで無事に終了． 質問時など， スタッフのかたがマイクの支援などやってくださった．
• 1430-1600 伊庭さん「アルゴリズムとしての MCMC」． 内容的には 岩波「統計科学のフロンティア12 計算統計 II マルコフ連鎖モンテカルロ法の基礎」 の p.28-98 あたりの範囲のものが多かった． ベイズの補足， マルコフ連鎖の収束定理 (エルゴート性)， レプリカ交換法とその背景， MCMC のベイズじゃない応用 (魔方陣の探索とか) ……
• 16 時すぎ， 公開講座ぶじに終了． 本日きいた伊庭さんコメントいくつか …… 例によって久保の聞き書きなので， まちがってたりする部分もあります．
  • 統数研ではベイズや MCMC に関する長い研究の歴史があって， そもそも赤池さん自身「AIC の次の研究」 としてベイズにとりくんでいた． AIC によるモデル選択は離散的というか， 「モデル A かモデル B か，どちらか一方」 といった「かくかく」したものなので， ベイズによってそのあたりなめらかにならないか， といった方向性で．
  • このように「モデル選択はやらなくてもよいように」 といったことがベイズを研究する理由のひとつだった． しかしモデル選択のようなことができれば便利な場合もある． モデル選択については「ベイズ統計と統計物理」 でもふれている． Bayes factor が使えるかも， といったハナシはあるのだが， 事前分布に無情報事前分布があると ……
  • はずれ値は t 分布などスソの長い確率分布で ベイズモデル化できるけれど， その MCMC 計算は難しいものになりがち． 変化点の検出だの， グループ数を指定しないグループわけだのも同様．
  • しかし予備検定とかそういうことをやらずに， まとめてひとつの統計モデルであつかいたい， というのがベイズな方向性のひとつ．
  そういえば， 伊庭さんの講義の中で， まだ MCMC わざが統計物理なヒトたちだけが使うものだった時代， その計算技法はマンツーマンで伝授されていくものだったんだけど， こんにちのように MCMC 使いが急増してしまってだいじょうぶかなぁ， といったハナシも．
• 終了後， 統数研-ドイツ大使館前の坂をくだった スープ屋 さんで伊庭さんとスープ (なぜ?) とコーヒー． ふーむ， 岩波書店ってそんなに「攻撃的」なところだったんですか ……
• ふと， 南麻布の統数研にくるのはこれが最後かも， と思ったり． 今年 移転 するそーで．
• 1738 日比谷線広尾発． 1742 恵比寿着． 1748 JR 山手線に乗って発． 羽田空港まで 550 円． 1757 品川着． 1801 京急にのりかえて発 …… のはずが， 乗りそこなってしまった． 8 両編成なのにホームの一番前方で待ってたため． そして人ごみなので移動できず． 1805 久里浜いき快特にのって発． 1815 蒲田着． 1822 京急羽田線に乗って発． 1832 羽田空港着．
• 1915 ADO の B767-300 に搭乗． 15 分おくれで 1930 taxing 開始． なかなか離陸しなかった． 機内で晩飯， あとはひたすら休む． 2110 新千歳空港着． 2130 空港発． 2210 JR 札幌着． 気温 -9°C． 2223 帰宅． ぐったり．
• 帰宅したら不幸のメイルが． A 棟 7-8F 内の PC 一台がウィルスのたぐいにやられて暴走， あちこちの IP アドレス 445 番ポートに対する無差別大量アクセスやってるので， 北大ネットから切り離されてしまった． やれやれ， 明日の朝一番はこれに対処しないといけないよーで ……
• [今日の運動]
  • ばてました
• [今日の食卓]
  • 朝 (0700): こんびに握り飯．
  • 昼 (1240): 統数研の公開講座会場で こんびに握り飯．
  • 晩 (2010): AIRDO 機内で羽田空港サンドイッチ． 機内さーびすのゴボウスープ．