ぎょーむ日誌 2007-01-17

2007 年 01 月 17 日 (水)

• 0650 起床． 朝飯． コーヒー． 0850 自宅発． 晴． 0905 研究室着．
• 屋久島データ教材化作業つづく．
• (甲山さんとの会議にきてた) 地環研起学専攻の Y せんせーと 7F トイレ (この建物は奇数階が男性用) における会話．
  「久保君，何才だっけ」
  「37 です」
  「もうそんなトシだったのか」
  「すみません」
  「いやいや，おたがいさまだ」
• そろそろ昼飯でも …… というタイミングで 雪野さん Mac とらぶるしゅーてぃんぐ依頼． GLP 事務書類 (例によって例のごとくわーどファイル) れたーへっどのロゴがばけます …… しょうがないので， NeoOffice をインストールしてその場をしのぐ． この NeoOffice， うわさどおりむちゃくちゃ挙動が重い． あとで考えてみたら， 雪野さん PowerBook に入ってる MS おひすが古かったんだろうな． えーと Mac 版の「最新」は office2004 とかいうやつか． といったことも連絡．
• そろそろ昼飯でも …… というタイミングで 松田さんから random effects きわめて大かつ線形ならざる 駒ヶ岳データ解析プログラム質問． やっぱモデル選択やっても同じのよーで， で一件落着?
• へろへろと昼飯調達の旅． 1530 昼飯． ふう． 食後の休憩と称して， 粕谷さん馬車馬ぶろぐ ［統計］なんでもノンパラメトリクス にコメントしてみる …… 「あらし」ではありません．
• へー， Ecology and Ecological Modeling in R - A special volume of the Journal of Statistical Software か …… しかし何も貢献できそうにない私はまだまだ実力不足 ……
• 1800 よーやくにして， 牛原さん屋久島データ教材化作業の一部， library(grid) わざ濫用の作図プログラム全面かきなおしが終了した ……

• えー? またまたいつもと同じ図ぢゃん， などとあなどるなかれ．
  • 自作 MCMC 計算 (Metropolis-Hastings 法) Perl コードではなく， R2WinBUGS 出力に対応
  • ベイズモデルの改良 (以前のはよけーな「個体差」が含まれてた)
  • nested 事後分布作図関数の汎用化 (以前はこのあたりつくづくヒドかった)
  …… という， まあ「見た目」ではわからん改良あれこれだったんだが ……
• つぎ． シュート伸長の休眠 (dormancy) と同時に二度伸び (lammas) を考慮する， つまり二項分布から三項分布に変えた統計モデル． WinBUGS だと何も考えずにあっさりと実現できる …… 初期パラメーターの調節とかは注意ぶかくやらんといかんけど． モデルを複雑にすると当然収束まで時間もかかる． これはステップ数で 400 → 1000， 私の非力な Thinkpad X31 で計算させると 800 秒ぐらい．

• へろへろと帰ろうとすると素早く院生に迎撃される． 矢澤さんと「樹冠頂部の破損を考慮した樹高モデル (もちろん「破損」の事後分布を推定する Bayesian な)」の相談 というか雑談．
• そしてお茶部屋に入るとすでに原稿をかなり書きススめてる 松田さんから統計モデルまわりの説明のかきかた質問． 例文をでっちあげてからふらふらと撤退 …… しかし建物をでた瞬間に random effects の説明がまるっきりぬけてたことに気づいた．
• 2110 研究室発． 2125 帰宅． 体重 69.6kg． 晩飯の準備． 晩飯．
• 梅干しとヨーグルトって意外と …… いいくみあわせ?
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0810): 米麦 0.6 合． 卵焼き．
  • 昼 (1530): 研究室お茶部屋． 北大生協のりまき．
  • 晩 (2240): 米麦 0.7 合． チンゲンサイ・ショウガ・マイタケのカレー．