ぎょーむ日誌 2010-11-01

2010 年 11 月 01 日 (月)

• どーにもよくねむれなくて， 0600 起床． ねむれないけどねむい． 朝飯の準備． 朝飯． コーヒー． 0700 自宅発． 曇． 0720 研究室着．
• じゅぎょーの準備．
• 9 時から本日の じゅぎょー， 本日は MCMC とベイズを強引に導入する 岩波統計モデル本 第 8 章 …… やっぱりここは GLM ベイズ化を分離するだけでなく， 図の修正その他も必要だわ ……
• こういう第 8 章的な内容を北大 M1 じゅぎょー (参加者のたいはんは M1 ではないが) でやるのは初めてなので …… ペイスがつかめない， 実演もうまくいかない， 口先でごまかそうとしてコトバ数が多くなる …… といったばてばて状態で 1015 終了．
• 予定は変更． もともとは次回は階層ベイズ & 最終回は空間ベイズの予定だったんだけど， 次回は第 9 章だけをねたに WinBUGS をもうちょっと丁寧にやる， 最終回は階層ベイズモデルをゆっくりやって全体の復習， 時間があれば空間統計についてもちょろっと．
• ぐったり．
• スンモヌさんの Mac のネットワーク設定してたら， 午前終了．
• 昼飯． 雨． ぐったり．
• うう ……
• ということで， 次回 (次の水曜日は祝日なので来週月曜日) は BUGS ちゅーとりある的な内容， その準備をせんといかん．
• これは一見よけーな仕事が増えたかのように見えるが， じつは一ヶ月後に九州大学理学部でそういう講習会じみたことをやるので， 今回の準備は使いまわせる …… はず．
• いつまでもゐんばぐすばかり， とゆーのもなンなので， ちかごろの JAGS のしあがりぐあいでも …… と調べようとしたんだけど， install.packages(rjags) がエラーにはばまれてインストールできん!
• …… じつはこれはばかばかしいミスで， Ubuntu 10.04 用の jags (1.04) がインストールされていたためだった． sudo apt-get remove jags してから， JAGS の source forge から JAGS-2.1.0.tar.gz をダウンロードして， 手動でインストールすれば， rjags の最新版はインストールできる．
• ここまでで， またぐったりしてしまったので， JAGS 調査はススまず．
• ねむいよーねむいよー
• 院生からの質問． 95% 区間を示す箱ヒゲ図はどう書けばよいか …… 先日 も同じよーな質問あったんだけど， R の boxplot() とかの「ヒゲ」の両端は近似的な 95% 区間を示していて (Tukey 式) ホントの 95% 区間ではない．
• いつもお世話になってる R Graphical Manual (RGM) 検索して， てきとーに library(fields) の bplot() とかでいいぢゃない …… とかいいかげんな返信したんだけど， これの style="quantile" 指定は 95% ではなくて 90% 区間だった!
• よりよい解の例としては， 質問してきた院生自身がみつけてきた library(scapeMCMC) の plotQuant() 関数． mcmc オブジェクトが主な対象だが (ふつーの data.frame でも使えるらしい)， まあ WinBUGS な出力の作図をしたいよーで …… はい， そうですね． いやー， ぼくはこんなのは自分で ぱっぱっと作図関数を組み合わせて目的の図を つくっちゃう からねえ， と負けおしみを述べてみる．

• まけおしみといえば， 私が先日あらしコメントを書きこんだ 牧さん blog で， ながのさんというかたがより良い解を示しておられた． そーだ， which.max() ってあったよな …… こんどからあらしコメントするとき前に， 間瀬さんの R 間数本を確認することにしよう ……
• あまりぎょーむは進捗せづ …… 1800 研究室発． 小雨． 買いもの． 1820 帰宅． 晩飯．
• うー
• [今日の運動]
  • 小雨 & ねむくてうんどう休養日．
• [今日の食卓]
  • 朝 (0740): 米 0.5 合． ゴボウ・タラのあらの煮物．
  • 昼 (1200): 研究室． 食パン． ハム． バナナ． リンゴ．
  • 晩 (1900): 米 0.5 合． ダイコン・ニンジン・春雨のカレースープ． ゴボウ・タラのあらの煮物． カボチャの煮物．