ぎょーむ日誌 2010-02-(01-10)

2010 年 02 月 01 日 (月)

• 0800 起床． ねむい． 朝飯． コーヒー． 0850 自宅発． 晴． 0910 研究室着．
• 午前中は小山さんとやちパワー解析のつづき． 先日の結果の説明． 「悪い時期」 を特定できるような改善． 再計算 …… などなど． 無事に終了．

• 露崎さんからゐんどーづ 7 情報をいろいろ教えてもらっていたので， 小山さん PC 上における R の挙動が少し改善された …… 先月 20 日 では install.packages() はできるけど， update.packages() しようとすると Program_Files 内の上書き禁止がどーのこーので うまく動作しなかった．
• 解決策は R アイコンを右クリックして「管理者」として R を起動すること． これは Vista, 7 以降の Windows で変更されたユーザー権限が原因でそうなる．
• ちょっとめんどくさいのは library(R2WinBUGS) を使って R から WinBUGS を呼びだして使いたい場合もこの「管理者」起動が必要とされる， といったあたり． WinBUGS が何やら Program_Files 内のファイルを書き換えたりするんだよね．
• …… といったあたりの情報も ちょっと追記 してみる．
• 北部生協で昼飯． その後， 地環研までてけてけと歩いていき， 「過半数代表候補者選挙」 なるものをすませてくる． これは棄権すると， 再選挙の危険性がたかまるめんどうなものなので， 昼飯後の散歩と考えてでかけてみた次第．
• さーて， いよいよ岩波本作業にもどろう (一週間なにもできなかった) …… と思ったのだけど， あれこれとまたわりこんできて， ですね．
• 今月は二回も東京出張があるのだけど， 二番目について旅の手配を何もしてなかったことに気づいたので， そのあたりをうだうだと． なお， 一番目の出張は次の日曜日に出発しなければならない．
• 空間統計モデリングの問題にはまる． 空間統計もののおそろしいところは， その問題の検討がおもしろすぎることだ． おそるべきイキおいで時間が消耗されていく ……
• 空間統計ほど楽しくはないが， ROC だの AUC だのの検討． こんなアヤしげなものどもに何か意味あるんだろうか …… と考えてみると， 少しはおもしろいところもある問題かも， と思えてきた． 今まで知らなかった点もいろいろとあり．
• ROC の receiver とは文字どおり受信者のこと． WWII の operations research の遺産． だからといってこんな意味不明な名前を長年つかわなくても．
• あたりまえのコトだが (といいつつ理解してなかったのだが)， (0, 0) から出発して (1, 1) にむかう ROC， スタート→ゴールにむけて linear predictor を小さくしていくかんぢ．
• 医者が偽陽性連発な検査をきらうのはコスト高になるから? しかし， 医者的には偽陰性のほうがマズいのではないかな? そして， ある種の個体群調査でソレ (偽陽性) がまずいのは …… 「好ましくない」方向にバイアスが入るから?
• けっきょく本日も岩波本作業には着手できずに， 1850 研究室発． 買いもの． 1920 帰宅． 晩飯．
• [今日の運動]
  • うんどう休養日．
• [今日の食卓]
  • 朝 (0820): 米麦 0.6 合． ゆでブロッコリーのひき肉・ショウガあんかけ．
  • 昼 (1300): 北部生協． 自宅からもってきた米麦 0.8 合． ゆでブロッコリーのひき肉・ショウガあんかけ． リンゴ．
  • 晩 (2100): 米麦 0.8 合． ニラ・シーフード炒飯． キャベツ・ニンジン・シイタケ・ホタテのスープ．

2010 年 02 月 02 日 (火)

• 0800 起床． ねむい． 朝飯． コーヒー． ぼーっと昨日の問題の検討のつづき． 0910 自宅発． 晴． 0930 研究室着．
• で， 13 時ごろまで， その昨日の問題の検討のつづきと， メイルかきでオワってしまった …… 一般性ありそうなところの抜粋．

  ROC や AUC は何だかよくわかっていなかったのですが，いろいろ勉強してみると，保全生物学的な文脈のもとでは，prediction が重要かつ偽陽性 (不在なのに在と主張) が偽陰性よりまずい，というのがひとつのまとめかたかも，という気がしてきました．

  つぎに ROC や AUC は何をしているのかと考えると，これはあてはまりの良さというより偽陽性のおこりやすさ (予測におけるまちがいの一部) を評価しているわけですが，保全生態学的な個体群サイズ予測では偽陽性は偽陰性より深刻な問題であると明確に指摘・主張するのもよいのではないかと思います． たとえば，個体群浮動 (population drift) が発生しうる状況では個体群サイズの過大予測は絶滅確率の予測をおおいにゆがませるといった指摘ですね (個体群浮動に関する論文には個体群サイズと絶滅確率の関係がいろいろとあると思います)．

  さて，ROC+AUC は偽陽性の存在を示すひとつの方法ですが，これだけでは個体群サイズ見積りに関してどのようなバイアスが生じているのかよくわかりません．そこで cross validation の結果を利用して，それぞれのモデルによる個体群サイズ (この場合だと合計「在」地点数の予測?) の予測のよしあしを明示的に示すのが良いかと思います．たとえば，このモデルだと個体群サイズを ○倍過大に推定し，その予測分布 (MCMC の結果を利用して生成) はこのようにずれていると示してみるのはどうでしょうか．

  これに対して，あてはまりの良さは必ずしも重要な点ではないように思いますので， deviance をつかった指標などは不要ではないでしょうか．あてはまりの良さにしつこくこだわる読者が多いような気もしますが，こういうヒトたちはだいたいにおいて (言うところの) 「説明力」と予測力の区別がついていないので，そのあたりを明確に記述すれば良いでしょう．

• …… などともっともらしく書いてはいるけれど， じつは一般性どころか， あまり正しくない ということが午後になってわかった．
• まず， ROC+AUC は偽陽性だけ着目しているわけではない， ということについて． こちらの 奥村さんの解説 がすごくわかりやすいのだが …… ようするに AUC なる確率の解釈は， 対象生物種の在・不在の問題をあつかっているときに， ある調査地点 i での 「在」 確率 p[i] が与えられているとすると，
  • ホントに在である地点の p[i] をランダムにひとつ選ぶ
  • ホントに不在である地点の p[j] をランダムにひとつ選ぶ
  としたときに， p[i] > p[j] となっている確率が AUC に等しい， ということになる． 完全識別できていれば p[i] > p[j] となる確率 (AUC) は 1 となる．
• ということで， p[i] > p[j] となる確率が低ければ低いほど， 偽陽性が多くなるとは言えるけれど， 同時にそれは偽陰性も増えるということだ． つまり， こういう識別の問題においては偽陽性が大きければ， 自動的に偽陰性も大きくなる， ということになる． これってあたり前のことなんだろうけど， ぜんぜん気づかなかった． マヌケだ．
• もうひとつは， 空間構造なしロジスティック回帰が 在確率の予測にバイアスを発生させるのでは， といった指摘なのだが …… これも正しくない といったことがわかってきた． 詳しくは明日以降にでも書くとして， ヒトことで書くならば， ナイーヴ平均場近似であっても場全体の平均「在」確率は それっぽく推定できる， とゆーことだ． 空間構造ありモデルはもっと局所的な構造を再現してみせるもの， ではないかな?
• …… といった 誤謬 には何も気づかぬまま， 平和に昼飯． 寒気が流入とかいってたけど， いい天気．
• …… うーむ， すでにアタマがばてたのか， メイルかき以外のぎょーむがススまん ……

[ぺたぺたブーム]

さいきん， 机のまわりで紙をぺたぺた貼るのがまいぶーむになっている． ディスプレイのうしろの紙は岩波本の目次で， その上に「これから書くこと」がぺたぺたされている． この紙片はぜんぜん減らず， 増えるばかり． 右はあれこれクリアホルダーなんだけど， 「耳」といいますか内容をあらわした紙片を 両面テイプとかではりつけている．

• 岩波本作業はあまりススまぬまま， 午前中に書いたあれこれの自分のまちがいに気づいて …… 1940 研究室発． 2000 帰宅． ちょっとうんどう． 晩飯．
• [今日の運動]
  • 腕立ふせ 5 × 3 回． ラジオ体操．
• [今日の食卓]
  • 朝 (0820): 米麦 0.6 合． キャベツ・ニンジン・シイタケ・ホタテのスープ．
  • 昼 (1320): 北部生協． 自宅からもってきた米麦 0.8 合． 生協豚汁とサラダバー． 394 円．
  • 晩 (2200): 米麦 0.8 合． ニラ・ホルモンの炒めもの． ヨーグルト．

2010 年 02 月 03 日 (水)

• 0700 起床． 朝飯． コーヒー． 0800 自宅発． 曇． 寒い． 本日の最低 & 最高気温 (予報) は -12°C & -8°C とのこと (あとで実測をみたら -12.6°C & -9.5°C だった!)． 0820 研究室着．
• 午前中は東大演習林の後藤さんとモデリング相談． トドマツ・アカエゾマツあれこれ．
• 後藤さんと昼飯にでて， この件は終了．
• 午後は岩波本作業に復帰を …… とのせつなる願いもムナしく， また別の問題にとらっぷされてしまった． モデル選択についての調査あれこれ． 以下， 断片的に．
• 「統計科学のフロンティア」 の「モデル選択」 の p.46 には情報量規準 TIC によるモデル選択と， cross validation によるそれは同等， といった説明がある．
• JAGS 開発者としても有名な Plummer さんの論文， Plummer (2008) では， 下のように expected deviance L^e(Y, Z) なる量が定義されているわけだが (Z は validation 用データ，Y はパラメーター推定用データ) …… cross validation によるモデル選択って， この L^e(Y, Z) が 最大 最小であるモデルを選べばよい? (cross validated な) 対数尤度関数を (事前分布は使わずに) 事後分布からのサンプルを使って積分 ……?

• …… てな検討でオワってしまった． 1830 研究室発． 1855 帰宅． ぐったり． 晩飯．
• [今日の運動]
  • ラジオ体操．
• [今日の食卓]
  • 朝 (0740): 米麦 0.6 合． 卵焼き．
  • 昼 (1210): 後藤さんと地環研前の生パスタ屋さんのレサン． 鶏キノコトマトソース．
  • 晩 (2120): 米麦 0.8 合． ハクサイ・ネギ・ミズナ・シイタケ・カジカの味噌鍋．

2010 年 02 月 04 日 (木)

• 0750 起床． ねむい． 朝飯． コーヒー． 0900 自宅発． 晴． 寒い日々がつづく． 0920 研究室着．
• ちょっとメイルかきあれこれ．
• 1030 より 研究室セミナー， 本日は宮田さんでアカエゾマツの炭素同位体なハナシ． 甲山さんの阿寒道楽から派生したハナシなんだが …… 気象データ→δ¹³C， 気象データ→直径成長， δ¹³C⇔直径成長， とゆーアヤしいとらいあんぐるを調べなさい， というハナシ． はてさて， どうなることやら ……
• ネット雑用あれこれ． よーやく 2/22-2/23 の A 棟帰還作業の情報がぽつぽつと． 下は会計係のゑくせるアートな新 A 棟地下図． 研究室雑用 wiki から参照できるように自分のディレクトリにほうりこんでみた (Pukiwiki のファイルアップロードはいまいち好きではないんで)．

• 仕事はススまぬまま， 昼飯． ははー， ろしあんなステルス戦闘機の試作 PAK-FA ですか …… しかもこの厳寒期にご近所のコムソモリスクナアムーレで? …… と思って 調べてみたら， スホーイの工場があるみたいね．
• また阿呆なことに時間をつかってしまった …… 本日， 生態学会東京大会のプログラムが届いたんだけど， それをぱらぱらとながめているうちに， いまさらながらに 3/20 (土) まで滞在 (6 泊!) するのではなく， 3/18 (木) の夜あたり (4 泊) で撤退したくなった． で， 旅行届け変更可能かどうかを確認してから， 航空券・宿泊の変更． 航空券は取り消しにまつわるあれこれ 3500 円ばかりとられたが (これは久保個人の損害)， 航空運賃そのものはなぜかしら (土曜から木曜に変えたから?) 1 万円ほど安くなった (これは公的な予算の節約) …… うーむ， 「変更するかも」といった予感があったのか， キャンセル料がまだしも安い航空券かっていたようだ． 自分ではそのあたり覚えていない． 宿泊のほうはキャンセル料とられなかったけど， 宿泊数減少で 1.8 万円ほどコスト圧縮．
• こんどから生態学会大会参加の旅の手配をする前に， オンライン版のプログラムをよくよく読んでおこう． といっても， 2011 年 3 月は札幌大会なのだが．
• 15 時から 1.5 時間ほど， 塩寺さんとキナバルデータ解析相談． いかにも植物生理学的な割算値地獄だったので， 問題点を整理するために， いったん階層ベイズモデルを設計して (実際に推定計算をやったわけではない)， データをみていくとだいたいどうなりそうかわかった． 現時点で私がしたうけぎょーむできる状況ではないので， とりあえず 「脱」わりざん値な offset 項 わざとかの説明 ……
• といった説明しているあいだにも， いろいろこんさるメイルが蓄積してしまう …… この季節っていったいどうなってんの?
• ほんのちょっとだけお手伝いしたカドフシアリ論文が accept された． 三浦・宮崎師弟ぱわー， おそるべし． アリといえば， 当方の石狩浜アリけんかモデル論文もさっさと再投稿しなくてはいかんのだが ……
• どうも最近は夕方以降は特にぐったりぎみだなー， と思いつつ撤退． 1850 研究室発． 1910 帰宅． ぐったり． 晩飯．
• [今日の運動]
  • うーむ …… なんで不調なんだろう?
• [今日の食卓]
  • 朝 (0820): 米麦 0.6 合． コマツナ・モヤシのあえもの．
  • 昼 (1320): 北部生協． 自宅からもってきた米麦 0.8 合． 生協サラダバー． 250 円．
  • 晩 (2200): 米麦 0.8 合． コンニャク・鶏レバの煮物．

2010 年 02 月 05 日 (金)

• 0800 起床． ふちょー． ぼー． 朝飯． コーヒー． 0910 自宅発． 曇． 0930 研究室着．
• 岩波本作業を再開したいところだが …… 2/9 (火) の公開講座投影資料の修正をやっといたほうがよさそうだ． 日曜日・月曜日はごたごたしそうなので．
• 伊庭さんの投影資料をよく読んでみる． けっこういろいろ改訂されているなぁ．
• 伊庭さんのと比べると， 私のハナシっていまいちおもしろくないような …… いまさらながら．
• ぼちぼちと修正してみる．
• 昼飯．
• 13 時から 2 時間ばかり， サロベツ M2 三人衆の発表練習会に …… みなさん， ぐったりおつかれですなぁ ……
• こちらもぐったり．
• 公開講座の修正作業のつづき． なかなか進捗しない ……

• 完了しなかった． なにごとにつけ作業がおそい …… 1920 研究室発． 1940 帰宅． 晩飯．
• [今日の運動]
  • うんどうしてない ……
• [今日の食卓]
  • 朝 (0800): シリアル． ヨーグルト．
  • 昼 (1230): 研究室． セヴンイレヴンのサンドイッチ．
  • 晩 (2130): 米麦 0.8 合． サトイモ・ニンジン・ネギ・豚肉の味噌汁．

2010 年 02 月 06 日 (土)

• 0830 起床 …… 朝飯． 洗濯． 怠業． うーむ．
• 昼飯． 1400 自宅発． 雪． すごい寒気． 1420 研究室着．
• メイルかきだけでオワってしまった． 1700 研究室発． 寒い． 1720 帰宅． 晩飯．
• 公開講座資料の修正作業． 能率わるい．
• どーにもアタマ痛いのでコーヒーのんでみる．
• 真夜中にとりあえず版アップロード． 明日から東京出張なのに何も旅じたくしてない．
• 岩波本作業はここでストップしたまま． 三週間以上なにもできてないかな? やれやれ．

• [今日の運動]
  • らぢをたいそう ……
• [今日の食卓]
  • 朝 (0900): シリアル． ヨーグルト．
  • 昼 (1300): シリアル． ヨーグルト．
  • 晩 (1900): 米麦 0.8 合． ダイコン・ネギ・煮干の味噌汁．

2010 年 02 月 07 日 (日)

• 0630 起床． 0725 自宅発． 気温はやや高い． -2 度くらい． 0735 JR 札幌着． Suica 一万円追加． これでしばらくはだいじょうぶだろう …… 今回もいれて今月・来月だけで東京往復 3 回だからなあ． いやはや． 0750 同発． ねむい． ぐったり．
• ぐったりしてるうちに快速エアポートが動かなくなった． 恵庭でポイント不転換事故．
• 30 分おくれで． 0855 新千歳空港． いそいで搭乗口をめざすが …… JAL の場合，保安ゲイト通過後に領収書発行システムなし． どうなることやら． 搭乗口での搭乗券発行もなし． いやはや．

[JAL 502 便 B777-200]

昨晩から千歳もけっこう積雪したよーで． 四割ぐらい空席．

• 飛行機は 15 分おくれで taxing 開始． さらに離陸まで 15 分かかる． 機内では勉強会予習．
• 15 分ぐらい遅れて， 1105 羽田着． 着陸時にゆれて， かなりきぶん悪くなった． 寝不足のときにそうなりがち．
• 1115 京急発． 1136 品川で JR に乗り換え． 1152 発のはずの中央線快速にのって， 1155 発． 1232 立川着． 車中ではひたすらぐったり． 胃のあたりがむかむかする．
• 多摩モノレイルには乗らず， 立川まで歩く． 北海道より 10 度ぐらい温かい， ちょうどよい天気． 統数研 にむけててけてけと早足で歩いているうちに， 回復してきた．
• 13 時すぎから統数研で， 深谷君・島谷さん主催の勉強会．
• 18 時すぎまで， 群集の中立理論の数理モデルの勉強． 低温研の平尾くんが，

  Rosindell, Wong & Etienne　(2008) A coalescent approach to spatial neutral ecology.Ecological Informatics 3:259-271

  そして深谷君が

  Etinenne & Alonso (2007)　Neutral Community Theory: How Stochasticity and Dispersal-Limitation Can Explain Species Coexistence. Journal of Statistical Physics 128 485-510

  の途中まで (つづきは明日)． どちらもたいへん勉強になった． このあたりまた後日にでも．
• 19 時ごろから島谷さんお気に入り (らしい?) 立川駅ちかくの 居酒屋 で懇親会． ビンビールケイスを逆さまにして合板をおいた楽しい店でした．
• 2120 ごろ立川駅北モスバーガーむかいの 宿 着． この宿の近所の 銭湯 の割安チケットもつくプランなので， チケットもらって風呂． なかなかよかった．
• [今日の運動]
  • うんどう不足ですなあ．
• [今日の食卓]
  • 朝 (0650): 米麦 0.6 合． ダイコン・ネギ・煮干の味噌汁．
  • 昼 (1300): 勉強会会場でこんびにおにぎり．
  • 晩 (1930): 立川駅ちかくの居酒屋． なぜか関西ふうな (なのだろうか? 白かったのだが) おでんがおいしかった． 皆さん酔っぱらってた．

2010 年 02 月 08 日 (月)

• 0740 起床． 東京滞在二日目． 朝飯． コーヒー．
• 岩波本相談のメイルかき ……

  1. データ解析に不可欠な統計モデル 
  ・以降の章の内容などが修正されないままになっている
  
  2. 確率分布と統計モデルの最尤推定 
  ・内容の入れかえなどの修正がまだ十分ではない
  
  3. 一般化線形モデル (GLM) 
  ・全体のまとめなどが新しい内容を反映していない
  
  4. GLMのモデル選択と検定
  ・まだまだ大幅な修正が必要
  ・モデル選択についてはあまりくわしく書けないかも
  ・検定については久保の意見みたいなものを大幅削除する
  
  5. GLMあれこれ 
  ・以前の GLM 章の内容の「雑」な部分を並べただけの状態
  
  6. 一般化線形混合モデル:観測できないランダム効果のモデル化 
  ・もうちょっと他の章とのつながりを説明する
  
  7. マルコフ連鎖モンテカルロ法と GLM のベイズ化 
  ・説明かきがまだまだ必要
  ・簡単な Gibbs sampler 例題を追加する必要あり
  ・WinBUGS の使い方はこちらに移動する
  
  8. GLMMの階層ベイズモデル化 
  ・WinBUGS 関係を第 7 章に移動
  ・個体差 + 場所差のような例題を追加
  
  9. 空間構造のある階層ベイズモデル 
  ・ここまでの章の変更を反映する
  

  …… ただし， ここに書いてたあれこれは， 夕方の相談では「無理なんじゃないの」 という方向になったわけだが．
• 0900 宿発． 統数研へ．

[多摩モノレイル]

…… の下をてけてけと歩くのは， なかなかよい． 雪もなく， 春みたいだ．

• 0930 から 16 時ごろまで， 群集勉強会のつづき．
• 16 時から伊庭さんと岩波本相談． 第 7 章 MCMC 章をどうするか， 問題． そして第 3 章の再りすとら．
• 立川駅ちかくで晩飯． その後， 伊庭さんお気に入りの オリオン書房 (立川北口店) に連れていってもらった …… たしかにこの本屋はすごい． 立川市， なかなかあなどれん ……
• 2040 宿に帰着．

[今日も銭湯]

なかなか良いところです． 写真は 2F ののんびりするコーナー (1F にもあり)． 漫画喫茶ふう， 店のヒトの趣味が反映されたジグソーパズルやフィギュアも多数展示． 不思議な空間です． この銭湯と連携している今回の 宿， なかなかよいかも (再来週は別の宿を予約ずみなのだが)．

• 明日は 公開講座 なのでさっさと寝なくては …… と思いつつ， なかなか眠れん．
[今日の運動]
• うんどう不足な日々 ……
• [今日の食卓]
  • 朝 (0820): 宿でこんびにサンドイッチ．
  • 昼 (1310): セミナー室でこんびに握り飯．
  • 晩 (1900): 立川駅北の高島屋の 中華料理屋 で海老炒飯・小籠包など．

2010 年 02 月 09 日 (火)

• 0730 起床． ねむい． 朝飯． コーヒー． 0905 宿発． 曇． 0920 研究室着．
• 10 時から 公開講座 「マルコフ連鎖モンテカルロ法の基礎と実践 (2010年版)」． 昨年と同じく，
  1. (10:00-11:00) 伊庭: イジング模型と２変量正規分布を使ったアルゴリズムの導入
  2. (11:00-12:30) 久保: ベイジアンモデリングとMCMC
  3. (13:30-14:30) 久保: R と WinBUGS の使い方
  4. (14:30-16:00) 伊庭: 数学的基礎・関連する話題など
  といった内容で．

[公開講座の会場]

進行も昨年 (ぎょーむ日誌 2009-02-23) みたいなかんじで． むろん， 伊庭さんも私も説明内容は手直ししたけどね． それから前回は準備ばてでへろへろしてたけど (上記リンク参照)， 今回はそのあたり気をつけていたので合計 2 時間少しを， なんとかのりきることできました． 参加者の皆さん (100 名ぐらい) は …… 年齢層たかし． まあ， 学割あるとはいえ， 有料だしね． 性比はさらにすごい偏りだけど．

• で， 今回も 16 時すぎぐらいに終了． 質疑応答などでいろいろ気づいたのだけど， またそれは後日にこのへんに書いてみる．
• 講座終了後， なぜかしら JR 立川ちかくでお茶 …… 今から考えると， この状況はなかなかすごくて， 伊庭さん・島谷さん・ほか野郎 3 名が エミリーフローゲ という名前の， 公開講座とはまったく逆の方向に性比が偏った店に突入して， コーヒーとケイキをいただきつつ， 何だかがやがやしている， といった光景が現出してしまって， ですね ……
• 往路に空港途中で立ち往生したので， 私は皆さんよりひと足さきにその店から逃げだして， 1743 立川発. 1818 新宿． 1822 同発． 1841 品川． 1842 同発. 東京では鉄道のポイントが低温・豪雪で動作しなくなる， といった事故はあまり発生しないようで， なんの問題もなく， 1907 羽田着． 時間あるので晩飯．
• 2020 JAL 549 便 B767-300 に搭乗． 9 割ぐらいの搭乗率? 10 分おくれの 2040 taxing 開始． 2205 新千歳着陸． 雪はだいぶ少なくなっている． 2231 同発 …… この時間は快速エアポートが運行されていない． 各駅なんで， 15 分ぐらいよけーに時間かかるんだよね．
• 2322 札幌着． Suica の残額から憶測すると， 今回の新千歳往復 (2080 円) + 東京うろうろで合計 5000 円ぐらい使ったことになるなぁ …… 札幌は当然ながら寒い． -3°C ぐらい． 2340 帰宅． さすがにぐったり．
• [今日の運動]
  • うんどう不足な日々．
• [今日の食卓]
  • 朝 (0800): 宿でこんびにサンドイッチ．
  • 昼 (1240): 講習会会場の講堂でこんびに握り飯．
  • 晩 (1940): 羽田空港第一ターミナルのイタリアンふうな店で， ミネストローネとカルボナーラ． 第二ターミナルに入っている店とは別なんだけど， 値段その他はまああまり変わらないといいますか …… 火曜日の夜ということで， 店内がすいててよかったけれど．

2010 年 02 月 10 日 (水)

• 0820 起床． まだねむい． 朝飯． コーヒー． 0920 自宅発． 曇． 0940 研究室着．
• 今日は雑用とメイルかきぐらいしかススまなかった．
• 昨日の MCMC 公開講座 で気づいたこと． やはり 100 名も集まるので， いろいろな人がいる (昨年の参加者アンケートによると職業もさまざま)．
• で， 今回もいろいろ質問いただいたのだが， やはりこれだけの人たちにいろいろ伝えるのは難しい， ということ．
• まあ， 私の説明がまずくて， 「ちょっとよくわからなかった」 という質問は予想されたものであったのだが． たとえば
  • ロジスティック回帰の意味がよくわからない
  • 階層ベイズモデルとそうでないベイズモデルの区別がよくわからない
  といったあたりはよくある質問なので， 質疑応答で簡単に対処できるし， まあ想定の範囲内なのだが ……
• やはり中にはいろいろなヒトがいるもので， 「Gibbs sampler とか自分で実装しているんだけど， 問題が難しくなるとしんどい …… WinBUGS ならそのあたりうまくいくのか， と考えて受講してみたけどよくわからない」 といった受講者も 1-2 名 (今回は少なくとも 2 名) いるんだよね． WinBUGS でなぜうまくいくのかよくわからない， ってのは当然で私がちゃんと説明してないからだ．
• で， 去年のぎょーむ日誌 などよくよく読み直してみると， 前回にも同じ質問うけてるんだよね (2009-02-25) …… ただし， これは講座終了後にメイルでの質問だったわけだが．
• で， 答えとしては， WinBUGS user manual をよく読んでください， ということになる． HTML 版マニュアル の MCMC method section みてもらうとわかるのだが …… 皆さん， おそらく私などよりよほど勉強熱心なのだけど， なぜかしら adaptive rejection sampler とか slice sampler とか意外と知られてないんだよね．
• まあ， 来年もそういう「Gibbs sampler 自作派，まじめに共役事前分布つかう流」 のヒトたちが来てくださるような気もするので， それにそなえてスライドを追加しておきますか．

• ということで， 公開講座のペイジ からダウンロードできる講座資料が日々更新されていく ……
• 公開講座に関してはまだいろいろあるけど， 今日はこのあたりで．
• 1940 研究室発． 買いもの． 2010 帰宅． 晩飯．
• [今日の運動]
  • 朝とかにうんどうしたほうがいいのかな ……
• [今日の食卓]
  • 朝 (0900): 米麦 0.6 合． 卵焼き．
  • 昼 (1310): 北部生協． 自宅からもってきた米麦 0.8 合． 生協サラダバー． 240 円． 今日からここ北部生協食堂は工事のため閉鎖．
  • 晩 (2050): 米麦 0.8 合． ネギ卵炒飯． 乾燥ダイコン葉・豆腐の味噌汁． ヨーグルト．