ぎょーむ日誌 2011-10-20
2011 年 10 月 20 日 (木)
-
0620 起床．
ラジオ体操．
朝飯．
コーヒー．
弁当の準備．
0800 自宅発．
0820 研究室着．
-
またデータ解析こんさるめいるが
……
-
岩波統計モデル本
第 9 章・図解ギブスサンプリング章のかきなおしのつづき．
-
うー
-
うーな昼飯．
-
ぱらぱらと見なおしてるうちに
……
なんでこんなに書きなおしに苦闘してるかというと，
ギブスサンプリング用ソフトウェア (つまり WinBUGS とか)
を導入するから，
ギブスサンプリングのご利益みたいなものの章ということで，
わざわざ新しく分離独立させたわけだが
(そのために複数パラメーターの GLM のベイズモデル化なる
例題もでっちあげた)
……
-
今次の修正方針ではこういう汎用なソフトウェアに関しては，
ギブスサンプリングではなく MCMC サンプリングに修正ってことになっているので，
まあ，
章中のあちこちで「ぎぶすさんぷりんぐ」狩りをやってるんだけど
……
どうにも書きなおしがしんどい．
-
で，
いろいろ考えてるうちに，
「第 9 章 GLM のベイズモデル化とギブスサンプリング」
なんて章の題名はそもそもヤメてしまえばいいぢゃん，
と気づいた．
なーんで，
章名にギブスサンプリングとかいれちゃってるのかね，
ばかばかしい．
「第 9 章 GLM のベイズモデル化と事後分布の推定」
でいいでしょ．
汎用な MCMC サンプリングなソフトウェアにとって，
ギブスサンプリングなんてその一部にすぎないという位置づけになったわけだから．
-
そもそも，
今年 5 月末
の第 9 章のかきなおしの時点で，「図解ギブスサンプリング」
な部分は「9.5 複数パラメーターの MCMC サンプリング」
という section の「9.5.1 ギブスサンプリング: この章の例題の場合」
という subsection レヴェルにまでおしこまれているわけだし．
-
ということで，
ギブスサンプリングなるものがでてくるのはここだけ
(+ ゐんばぐすなる名前の解説)
に局限して，
章の他の部分はベイズ化ポアソン回帰の解説と
ゐんばぐすの使いかた説明に終始しておればよいということですな．
いやはや．
-
……
といった「ぎぶすさんぷりんぐしわけ」
で最初のめんどうかきなおしを突破したので，
ひとやすみ．
-
さて，
作業のつづき．
次のめんどうは WinBUGS 結果説明の用語の問題なのだが
……
-
Gelman の Bayesian Data Analysis でも読んで，
たまには勉強するか．
-
研究室ねっと雑用がまたわりこんでくる．
(2011-10-20 16:06:40)
-
convergence ってのがねえ．
-
とりあえず，
また図をひとつ追加することに．
R-hat 図．
-
しかし時間ぎれ．
1810 研究室発．
買いもの．
1840 帰宅．
晩飯の準備．
晩飯．
-
新しい図はこんなかんぢでどうかな?
図をかくと，
元気がでるね．
-
[今日の運動]
-
朝のラジオ体操．
30 秒ぶらさがり × 3 回
……うんどうぶそくぎみだな．
-
[今日の食卓]
- 朝 (0640):
シリアル．
ドライフルーツ．
ヨーグルト．
- 昼 (1200):
研究室．
米 0.7 合．
ニシン丸干し (カズノコつき) のグリル．
カボチャの煮物．
さつまあげ．
レタス．
- 晩 (1920):
ハクサイ・ニラ・ネギ・エノキダケ・ブナシメジ・ハンペン・豚肉のキムチ鍋．
卵雑炊．