ぎょーむ日誌 2009-01-14

2009 年 01 月 14 日 (水)

• 0820 起床． うーむ． 0845 自宅発． 雪． 0900 研究室着． 朝飯． コーヒー．
• いろいろメイルかきとか． 窓の外は吹雪になったり， いきなり晴れたり．
• とあるデータ解析依頼 (やや急ぎ系?) があったので「蟻」解説かきはいったんほうりだして， そちらにとりくんでみる． なんとならば， 植物の重量分配問題であり， これはいままでさんざん苦闘してきたので， その経験と作図プログラム (とくに昨年の修論したうけでいろいろ開発したもの) を使えばあっという間にカタづく， と楽観していたため．
• この憶測は部分的にはあたった． ベイズモデリングや MCMC 計算そのものは， それこそあっという間にかたづいた …… しかし作図でツマづいてしまった． これは昨年の修論したうけ解析がかなり奇怪な設定 (わかりやすく言えば甲山さん的な設定) になっていて， それにあわせた読解困難な作図プログラムになっていたためである． 作図プログラムから読解困難な部分を削除して， わかりやすいものにする， という過程は意外と技術的に困難なモノだった …… 作図プログラム改訂が終了したらとっくに 13 時すぎてた． うう． やれやれ． がっくし．
• 発注者むけの簡単な解説メイル書きにさらに 20 分ほど費やしてしまった． まあ， 結果はわかりやすいものでしたよ． BUGS code， 事後分布表． 脱「なんでもかんでもあろめとりー」うんどうを推進しますかね．

核心部は地上部の総重量 (exp(log.totalw[i])) のうち p[i] を葉
(y[i]) にまわし，残りを x[i] にまわす，という部分です．

	# weight allocation between root and shoot
	x[i] <- exp(log.totalw[i]) * (1.0 - p[i])
	y[i] <- exp(log.totalw[i]) * (p[i])
	logit(p[i]) <- logit.p[i]
	logit.p[i] ~ dnorm(mean.logit.p[i], tau.tree)
	mean.logit.p[i] <- (
		bb[1] + bs[1, Spc[i]]
		+ (bb[2] + bs[2, Spc[i]]) * (log.totalw[i] - Mean.log.totalw)
	)

このようなモデリングによって，いいかげんな因果関係を仮定する
アロメトリーを回避できます．

• 昼飯．
• お茶部屋で宮田さんと苫小牧樹木光合成議論．
• 来月末の「北海道若手生態学研究会」 (なかなかスゴい改名ですな) の 案内ペイジ が アナウンスされる． 誰でも参加できます．
• ぐったりしてないで， 「蟻」解説かきをすすめなくては ……
• スウェーデンに帰国された Annika さんから 「北大のメイル転送設定どうやるの?」 質問が …… いつまでたっても複雑怪奇なままの北大ネット世界をしばし探索すると …… よーやく， わざわざ https://www.hines.hokudai.ac.jp:20000 から login して設定変更すればよい， とわかった． きわめてわかりにくい経路で到達した． そして将来的にはまた別の場所に移転するであろう， と．
• 解説かき， あまりススまぬまま撤退． 1910 研究室発． かなりの降雪． 1925 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • うーむ ……
• [今日の食卓]
  • 朝 (0910): 研究室． 食パン．
  • 昼 (1400): 研究室お茶部屋． 米麦 0.8 合． フリーズドライ野菜スープ． リンゴ．
  • 晩 (2030): 米麦 0.8 合． ホウレンソウ・豚肉の炒めもの． ヨーグルト．