ぎょーむ日誌 2007-02-12

2007 年 02 月 12 日 (月)

• 0720 起床． 朝飯． コーヒー． 洗濯しつつ計算試行錯誤をススめる．
• 木金土日月と雪山ラッセル連続アタック & 連続撤退的な 「豊平川アミノ酸濃度の階層ベイズモデリング」 試行錯誤， よーやくオワリが見えてきた ……
• 昼飯の準備しつつしあげの計算 12000 MCMC step とりぷるチェイン 620 秒 …… まあ， 収束とかまだまだがたがたしてるところあるけれど， 各種ややこしいばらつき錯綜している この観測データセットをシメあげて 人間が引っぱり出せる情報はこんなもんでしょう． 計算終了． 昼飯．

• 今回の解析あれこれ．
  • 水槽効果とでもいうべき random effects があちこちに入っている (図を描くとよくわかる)
  • さらに， 測定機器にかける前の処理によっては 測定値がすごく下ることがあり， その変化はかなり大きい (これも統計モデルに組みこむ試行錯誤が必要であった)
  • 以上のような状況なので 水槽ごとのアミノ酸組成の時間変化， みたいな細かい変化は抽出できないはず …… なンだけど， この組成変化をあらわすような パラメーターが必要である (なぜそうなのか理由はよくわからないんでけど MCMC 計算における混交ぐあいがかなり改善される …… 「動かないパラメーター」をゆさぶるバネとして 機能している?)
  • 今回も多変量確率分布の事前分布は使わなかった …… アミノ酸濃度間の分散共分散 (17 + 17 * (17 - 1) / 2 = 153 個) ぜんぶを推定できるデータではなさそうだったんで …… というか， 多変量事前分布の使いどころをまだまだ理解できていない
• …… てなコトを考えつつ 1340 自宅発． 吹雪だけど晴． 気温は低い． 1355 研究室着．
• 脱力したせいか仕事ススまん． 窓の外は猛烈な吹雪． と思ったらいつのまにか晴れていたり．
• アミノ酸問題に没頭というかかまけていてこの数日間を費してしまったので …… やるべきことは多々あるんだけど， とりあえず 2/15 (木) 〆切の査読原稿にまずはとりくんでみる． こんな状況でどうして査読なんか引き受けるんだ， と思われるかもしれませんが， 世の中なかなか断わりにくい依頼というものがあって， ですね ……
• 先日 30 分ばかりざっと見て (このやや奇妙な) 計算の核心部分は理解できたんで， 本日はなぜかくのごとき状況になってしまったのか， あたりを解明する読解．
• …… するつもりが， 1.5 時間ほど熱帯林解析ハナシ． このハナシ， 次に何をやるにしても計算が格段にたいへんになりそう．
• 査読原稿よみ …… ああ， (私のところにまわってくるだけあって) コレもけっこうハズレだ． つっこみどころ満載． アタマくらくらする． さてさて， どうしたもんだか．
• 1920 研究室発． 1940 帰宅． 体重 68.6kg． 晩飯の準備． 晩飯．
• ちょっとデータ解析こんさるメイルかき， とか．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0800): 米麦 0.6 合． 卵焼き．
  • 昼 (1310): 明太子スパゲッティー．
  • 晩 (2130): 米麦 0.7 合． 納豆． ネギ・シイタケ・豆腐・ワカメの煮物．