ぎょーむ日誌 2009-05-28

2009 年 05 月 28 日 (木)

• 0700 起床． 朝飯． コーヒー． 洗濯． 一次元空間上の確率場についてうだうだと考えてみる． 0930 自宅発． 晴． 川むかいの北 12 条郵便局に転送してもらった 定額給付金の紙束をうけとる (札幌市は全国の自治体の中でこの紙束を郵送するのが もっとも遅い)． まあ， 15 兆円の補正予算の一部はのちのちわれわれが支払わされるので (納税者ひとりあたり 20-30 万円ぐらい?)， 1.2 万円ぐらいは返してもらってもいいでしょう， ということで． 0955 研究室着．
• 昨日のつづきで， また一次元 GRF 問題で試行錯誤してみる． しかしあまり進捗なく，

  	x[i] ~ dnorm(q[i], tau)
  	q[i] <- (x[LeftRight[i, 1]] + x[LeftRight[i, 2]]) * w
  

  としたときに
  • w が 0.5 より小さい正の値であれば (そんなにいろいろ調べたわけではないが)， まあだいたい同じような結果に収束する
  • w = 0.49 など 0.5 に近づけていくと， 徐々に収束が遅くなる
  • w = 0.5 は「臨界点」みたいな値で， この値だと一気に収束が遅くなる (あるいは永久に収束しない?)
• 雪野さんとプリンター相談， とか．
• で， 夜まで藤沼君と相談しながら Pasoh ブタ-Clidemia 空間統計モデリング． 段階的に試行錯誤しながら， いろいろな要因をほうりこんでみる．
• とりあえず， ひととーりあれこれを考慮したモデリングができるようになった …… とはいえ， これにはまだまだ改善すべき点がある， とあとになって気づいた． (BUGS code， 事後分布表)． 計算時間は 1300 秒ぐらい．

• 1915 研究室発北大構内走． こんな統計モデリングでよいのだろうか， と検討しつつ． 1955 研究室もどる． ばてた． やはり他人にあれこれ説明しながら， 試行錯誤していたのでアタマがクタびれてた．
• 2030 研究室発． 買いもの． 2050 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • 北大構内走 40 分間． ストレッチ．
• [今日の食卓]
  • 朝 (0730): 米麦 0.6 合． ニラ卵炒飯．
  • 昼 (1240): 研究室お茶部屋． 食パン． ダイコン・ニンジン・ブナシメジ・鶏肉のコンソメ煮物．
  • 晩 (2200): 米麦 1.0 合． ダイコン・ニンジン・豚モツのカレー． ヨーグルト．