ぎょーむ日誌 2010-03-02

2010 年 03 月 02 日 (火)

• 0710 起床． 朝飯． コーヒー． 0850 自宅発． 晴． 今日も真冬日， とゆー予報． 0900 研究室着．
• 深谷君メイルに触発されて， 手もとのベイズ本あれこれを調べてみる．
• 小山さんとサロベツやちパワー解析の 階層ベイズモデルの詳細を再検討する作業のつづき． 全パラメーターのうちひとつだけ 「収束がおそい」 ものがあって， それは (私のなかばやけくそな formulation によって定義された) tau[1] なるばらつき (== 分散の逆数) パラメーターで …… これを全種共通にしてるんだよね．
• そのあたりに無理があるのか， こいつの収束がおそい …… 調べてみると， この tau[1] によってばらつきをコントロールされてるパラメーターたちが， 数千ステップに一度「全パラメーターほとんどゼロ」 といった状態におちこむんだよね． そうすると tau[1] は 1.0E+3 から 1.0E+4 といった範囲の値になる (分散の逆数なので)． で， いったんこうなると， またゼロ近くまでもどってくるのに時間がかかって …… そう， あたかも収束してないかのように見えるんだよね． ただし必ずゼロ近辺に戻るので，収束が遅いと解釈すべきなのだろう．
• とゆーことで抜本的な解説策としては， 今より計算時間を 10-100 倍ほど長くしてやればめでたしめでたしなのだが …… うーむ， 私自身もあまり非力な PC に長時間計算させるのは好きではないし， まあなんといいますか， 私の一存では冷汗の流出量を制御できない問題でもあり ……
• しょうがないので， tau[1] の事前分布 (つまり超事前分布) を dgamma(1.0E-4, 1.0E-4) から dgamma(1.0E-2, 1.0E-2) と変更してみたら， (みかけ上) 収束が速くなった． まあ， この変更によって影響をこうむるのはどのあたりだろう， と R で作図などしてココロの平和を追求してみたのだが …… tau[1] > 1.0E-2 では「見た目」はそんなにかわらないんじゃないの， ということで (tau[1] は 1.0E+0 から 1.0E+1 ぐらいのケタ)， この設定でいきますか．
• ここまでの顛末を要約すると， 「無情報」 (non-informative) ということになっている事前分布は， たしかにカタチを多少ひねくっても全体的な傾向はほとんど変えないけれど …… ただし， 収束のおそいばらつきパラメーターが 「あたかもさっさと収束している」 かのように見せかけることができたりしますよ， とゆーことで …… うーむ．
• 昨日の CANON LBP5400 (in A710 院生部屋) のふちょー問題のつづき． とりあえず， K トナーカートリッジを交換してみたら， と再生トナーカートリッジ屋さんが助言してくれたので， そうしてみたらうまく動作してるように見える． 問題の K トナーカートリッジは持ち帰って調査するとのこと．
• 生態学会から何やらとどく．

[東京大会 ESJ57 の名札など]

毎年問題になる「初日の受付ぎょーむ」 その他を回避するために事前に名札・領収書などが 郵送されてきた． これはたいへん便利なのだが …… 当日， 何割ぐらいのヒトがこれを忘れるだろうか? 私などはそういったど忘れ危険性が高いので， ESJ57 関係あれこれをまとめてファイリングしてみる ……

• 昼飯．
• またちょっと部屋のかたづけ． 机まわり．

[耐震つっかい棒の流用]

このように壁際に設置して ……

[机の前の本棚を増設]

さきほどの棒とファイリングキャビネット (左下) に三段やすもの本棚をのせてみた． 先日 Homac で買った別の耐震つっかい棒で上からも おさえる …… ただしここは天井板がぜいじゃくなので， 震度 6 とかにはぜんぜん耐えられないだろうけど．
ともあれ， ここ A 棟最上階は地震に敏感で， 札幌市だと 2 年に一度ぐらいある震度 3 の地震でも， ここでは震度 4 ぐらいのゆれになるので， そのときに被害を多少は軽減できるかな， ぐらいの対策．

• で， またベイズ本しらべたり， メイルを書いたり ……
• A801 院生部屋のプリンター， CANON LBP3410 もトナーぎれ．
• そして A801 で配線工事．
• 1910 研究室発． 買いもの． 1930 帰宅． 晩飯の準備． 晩飯．
• えーい， なかなか仕事に復帰できない日々だ．
• そしてまた新しい雑用が累積することになってしまった． なんでだろうね． 根気よくぢりぢりと攻略する以外に， 方法はない．
• [今日の運動]
  • うんどうもできていないし
• [今日の食卓]
  • 朝 (0730): 米麦 0.6 合． ダイコン・ニンジン・タマネギ・豚肉の味噌汁．
  • 昼 (1320): とつぜん蕎麦を食べたくなったので， 北 13 条門前の「さか田」のあられかきあげ蕎麦．
  • 晩 (2030): 米麦 0.8 合． ダイコン・ニンジン・タマネギ・豚肉の味噌汁． コマツナ・卵の炒めもの． ヨーグルト．