ぎょーむ日誌 2007-12-26

2007 年 12 月 26 日 (水)

• 0820 起床． さいきん「電灯つけっぱなし寝」が増えてしまっている …… 朝飯． コーヒー． 0940 自宅発． 晴． いつのまにか歩道が雪で被覆されている． 0955 研究室着．
• 苫小牧 + 屋久島統合モデリングのつづき． 意味のない中央化をやってたことによーやく気づいたので削除 …… しかしあとになって， まあこのへんはどう「ずらし」ても結果は同じだし， 中央化してたほうがゐんばぐすの収束は良いのでそのままに． 昨日わずらわされていた「耐陰性」なる確率論的なグループわけは 初期値のおきかたにけっこう依存することがわかってきた． このあたり修正しつつ， 再計算をくりかえす．
• いろいろとややこしい推定結果がでてくるので， いまさらながら観測データ作図を工夫してみる． うーむ． なンというか， べいづな方法は強力なのでかかるノイズだらけの中から 何か傾向みたいなものをあらわす事後分布をすくいだしてくるわけだが ……

• …… 実際のところ， 統計モデリングにおいて重要なのは (とくにモデルが link 関数ナシの最単純な線形モデルでない場合は)， 個々のパラメーターのゆーい差ではなく， むしろ説明変数 (input) のとりうる範囲における 応答変数 (output) の変化のしかた， つまり統計モデルの全体的な挙動そのもの に注意せよ …… ってことなんだよね．
• どうしたものかと思いつつ昼飯． 昼飯．
• ともかく推定計算全体の中でうまくいってない部分 (原理的には識別不可能な個体差 + 測定誤差のきりわけ) ははっきりしてるので， これを何とかごまかす方法を考えないといかんな． 試行錯誤． 計算まち時間には， 他の M2 こんさるとかメイル書きとか．
• 家族の一人が数ヶ月前にくも膜下出血になった， という観測データが得られてしまったので， 私のアタマの中にも動脈瘤が存在する 条件つき確率 の odds はいっきに数十倍に高まった …… ということで文科省共済組合の脳ドック (アタマの中の動脈瘤の有無も調べてくれる) の 組合員割引 みたいなものが使えないかしらん， と先日にもうしこんでみたんだが …… どうやら脳ドックは人間ドックとセットじゃないとだめみたいだね．

  札幌厚生病院にてご予約を受けておりますが、
  ご希望頂きました脳ドッグは、人間ドッグを受診した上で、オプション受診になりま
  すので、補助金の対象外となります。
  また、札幌厚生病院は人間ドッグが本年度は予約不可となっていますので、申し訳ご
  ざいませんが、別の医療機関で人間ドッグのお申込を頂きたいと存じます。
  

  ということで， めんどうになったので次の定期健康診断 (6 月ごろか?) のときに (おそらく次回も必要ないんだけど) 人間ドックを受診してそのついでに脳の断面診断もやってもらうことに …… ということで， それまでは血圧が急上昇するような状況に私を追いこまないでください．
• なお上のメイルの 「脳ドッグ」「人間ドッグ」はそれぞれ脳ドック・人間ドックの 北海道方言 (なのだろう)， イヌとは特に関係ないことを注記しておく．
• 苫小牧 + 屋久島統合モデリング， 階層ベイズの MCMC 計算を終わらせるための試行錯誤あれこれ． 「耐陰性」なるラベリングを「皆さんのおおまかな意見」 にもとづいて「耐陰性あり」の assign を潜在変数をもちいて表現する， という魅力的なアイデアは今回は却下することに …… というのも， ホントに「耐陰性」あり・なしどちらとも決めがたい樹種というのが いくつか存在するとわかったんで．
• そのかわり先験的な (というほどたいしたものなのか?) 耐陰性 3 水準 (-1, 0, +1) なる数量を平均とする事前分布から 「耐陰性値」なるウソにウソを重ねた 事後分布で表現する安易なモデリングに変える …… これで収束は少しマシになった．
• だいたいこのモデルでうまくいくだろう， というのがわかったので 4000 MCMC step 計算を命じて撤退． 1940 研究室発． 2000 帰宅． 晩飯の準備． ちょっとうんどう． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0900): 米麦 0.6 合． ワカメ・豆腐の味噌汁．
  • 昼 (1330): 研究室お茶部屋． 米麦 0.6 合． もくずスープ．
  • 晩 (2250): 米麦 0.8 合． ワカメ・豆腐の味噌汁． カスベ 電磁波酒蒸し． キャベツ・ニンジン・キュウリ・タマネギのサラダ．