KuboLog 2006-11-10

2006 年 11 月 10 日 (金)

0745 起床． 0755 自宅発．曇． 0810 研究室着．朝飯．コーヒー．

とりあえず昨晩生成したマレイシア調査用 PDF ファイルのうち，必要最小限のものだけ選抜して圧縮．それでも 1.4MB あるか．で，ダウンロードできるところにおく．そうか時差が一時間あったな．

ぎょーむ日誌とか昨日から放置してるメイル返信かきあれこれとか．また査読しやがれメイルとかもあるし．

R のこれまたよくある質問．グラフの軸で単位とかかきたいときのベキ乗 (上つき) とかはどう書けばよいか．ぎょーむ日誌でも何度か紹介してるけど，たとえば expression("density (g cm" ^ -3 * ")") とすれば，

となる． expression() の () 内に「数字だけでなく文字列をも項とする疑似数式みたいなものを書く，二項演算子でつなぎつつ」と考えるところがポイントである．すると * は「かける」をあらわし expression() はかけ算記号が省略した表記 ("x" * "y" はかけ算記号を省略して "xy" と表記する，といったルールになっている) で graphic device 上に描くという対応関係が理解できるだろう．

また院生のマシンに R のインストール & Tinn-R のインストール．

アップロードしたあれこれに関して，マレイシアからはやりなおし命令もなく「問題なし」との連絡が． IR (いんたーなしょなるれすきゅー) ぎょーむ，無事に終了，か?

12/7 (木) - 12/10 (日) 東京出張の宿の手配．昨年 12 月の R ゆーざー会のときは，羽田から至近の蒲田に泊まった，というのも南麻布には土曜日の朝に移動すればよかったので．しかし今回は金曜日朝の通勤通学時間帯にも移動せんといかんわけで …… 蒲田から都心への経路はスシづめ状態だろう．ということで軟弱にも統数研から至近距離に宿をとる …… 一泊 9450 円! うーむ …… 蒲田あたりだと 6900 円ぐらいなンだが．航空券も予約．

竹中さんから 5th International Workshop on Functional-Structural Plant Models (FSPM, 2007 年 11 月 4-9 開催) の案内いただいたわけだが …… いやはや，要旨提出は来年の 3/19 か．ぱいぷ樹木の改良改善にかまけているヒマあるんだったら，アリ科研費しがらみ研究と未提出宿題あれこれやりなさいよ，という「借金」状況だよなあ …… ああ，にゅーじーらんど行ってみたかった．

だったらアリ関係のメイルでも書けよ，ということになるが …… とりあえずすでに午後だし，気力もつきたので昼飯．

MacOS X のヒトたちも Darwine とかインストールすれば WinBUGS が動かせる (??) だろうから R2WinBUGS とか使えるようになるのではなかろうか ……? 特に Intel Mac 上ではかなり高速に動作するような気がするんだけど．

といったことを考えつつメイルなど書いたり．しかしアリ関連メイルのほうは書く気合いが足りない．

なンかいろいろ雑用とか．

そして小林さんの新アイデアでアカマツ (野外データのほう) 解析やりなおし …… うーむ，やはり「個体差」とか考慮せんといかんのか?

[光合成速度 vs 葉重量]

まあ，最尤推定値は「まがる」わけですが． AIC で比較すると，単なる直線回帰が勝つんだよね．そしてわれながらバカげていることやってみたんだが …… 「一個体から一データしか取れてないけど『個体差』推定しろ」とかいう命令だしても， R の lmer() や nlme() はなかなかカシこくて (というか当然ながらというか) 「そんなのできるか」と却下されてしまった．これまた順当なことにというか，処理ごとにセンをあてはめる，というのもパラメーター使いすぎで AIC よくない．

1900 研究室発． 1920 帰宅．体重 72.4kg．晩飯．すぐ寝てしまった．

[今日の運動]

うんどうふそく

[今日の食卓]

朝 (0820): 研究室お茶部屋．クロワッサン．
昼 (1320): 研究室お茶部屋．食パン．
晩 (2040): 米麦 0.7 合．ダイコン・豆腐・ワカメ・煮干の味噌汁．

ぎょーむ日誌 2006-11-10

2006 年 11 月 10 日 (金)