ぎょーむ日誌 2009-01-04
2009 年 01 月 04 日 (日)
-
0820 起床．
うーむ．
朝飯．
コーヒー．
怠業．
うんどう．
体重 71.8 kg．
うんどう不足で体重かわらず，
といったところか?
-
1150 自宅発．
曇．
地下鉄北 12 条駅ちかくの 100 円ショップでいろいろ買う．
いろいろな工具 (けっこう品質高そう) とか
金属メジャー (ストッパーつき) とか計算機ソフトバッグとか
何もかも 100 円で買えてしまう不思議さよ
……
1220 研究室着．
昼飯．
お茶部屋の流しのそうじ．
-
昨年末から続いているシラカンバ・ホオノキ一次枝解析，
統計モデリングは昨日だいたいのところ終了したので，
本日はその細かい手なおし．
長さの単位を cm から m に直したり，
決め打ちと申しますか「主観的」な測定誤差
(これは文字どおりの意味で，測定したときの誤差)
のケタ数の調節とか
……
に時間をとられてしまう．
-
計算まち時間に雑用いろいろ．
研究室 server 機の backup 状況のチェック，
ついでに
sudo /sbin/fsck.ext3 -C -y /dev/hdb[13]．
-
シラカンバ・ホオノキ計算，
細かい改良をほどこして，
まあこんなところかなといった結果が得られた．
BUGS code
と 事後分布表．
結果を格納した R オブジェクトが肥大化するのをイヤがって，
下の図のごとくケチってるわけだが
(枝・樹木ごとの RGR だけ)，
これがアダになってまた再計算，
とゆー可能性も残されている．
6000 MCMC step で 3000 秒ぐらいだったかな?
-
まあ，
やっぱりこの一次枝の成長・死亡・発生モデルがすっきり作れたのは，
(implicit ながら)
資源の流れってのを想定しながら統計モデリングできたからだよなぁ．
で，
そのベイスになっているのは「2007 年から 2008 年にかけての変化」
という時間変化データであった，
と．
見かけとはうらはらに，
動態に関しては意外と個体差・枝差が少なかったこともさいわいしました，
ってかんぢでしょうか．
-
で，
この M2 下うけぢょぶを終了させるためには，
計算説明・結果説明の報告書をかかねばならぬのだが
……
いつもながら「めんどうだなぁ」気分がつのって，
なかなか進捗せず．
とはいえ，
BUGS code と結果を格納している R の
bugs
オブジェクト (bugs class オブジェクトの
概要)
をぽいっとわたしてオワり，
ではすまされないからなぁ
……
1915 研究室発．
雪．
寒い．
氷点下．
1930 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0850):
シリアル．
ヨーグルト．
- 昼 (1250):
研究室お茶部屋．
米麦 0.8 合．
タマネギ・ニンジン・シイタケ・豆腐の炒めもの．
リンゴ．
- 晩 (2100):
コンニャクの煮物．
キャベツ・ニンジン・ブロッコリー・ソーセイジの酒蒸し．
ヨーグルト．