KuboLog 2006-08-17

2006 年 08 月 17 日 (木)

0750 起床．朝飯．コーヒー．ノド問題は解決しつつあるように思える． 0920 自宅発．本日からひさしぶりに雨ふるはず …… という数日前の天気予報ハズれて晴．で，そのあてにならぬ天気予報によれば，また真夏日になるよーで．たしかに暑い． 0935 研究室着．

スリランカの争乱はぜんぜんおさまらないね．北端のジャフナ半島と北東部のトリンコマリーが主戦場，とのこと (地図)．島の南西部のコロンボ市内の爆弾テロはいまは一段落ついてるのかな．

また林冠モデル作文のつづき．昨日追加した結果，葉群の角度分布うんぬん，といったあたり．

先に revision をだしたアカマツ原稿の editing 請求． 41000 円．思えば不思議なコストだ …… とひっかかていてもしょうがないので，手紙を書いて九大事務に転送する．そちらで管理してもらってる科研費から支払ってもらうので．

ところでそのアカマツ原稿の現状は，といえば …… すでに二週間ちかくが経過したのに status: with editor ですか．あいかわらず「時間をかければかけるほど良い」なる反主流派ぽりしーを貫徹してやがるな．

午前中にすでに気温は 30°C 突破．

林冠モデル，角度うんぬんについて作文してるうちに気になったので，下のような図を描いてみる．やはり「高さ」によって葉群の角度分布がびみょーに異なる．しかしこのあたりは biological な意味があると考えるより，むしろ推定における偏りみたいなもの，という気がする …… そうか，このあたりを補正するべいづ推定を組みこめば，むしろ収束が速くなるかも!

ふーむ …… 上の図みてるうちに，高度別角度分布の図とか作らなくてもいいとは思うんだけど (あってもあまり意味がないから)，葉群垂直分布 (左上) は追加するべきだね．これまた「推定があまりうまくいってない」側面を示すものとして．

ということで論文用の図を描くべく， R の作図プログラムを改訂．

[タテ histgram]

これは hist() だけでは描けないので， hz <- hist(データ, ..., plot = FALSE) で結果を格納し， hz$breaks だの hz$counts を sapply() などで取りだしつつ rect() でハコを描く．
にしても地面ふきんに葉が多いのは，観測データと推定モデルの両方が「いまいち」であるから，というもの．このあたり解決したければ，かなり Bayesian 化をおしススめる必要がある．

他にも図を作り直したり，文章をそれにあわせて書き直したり． R で postscript() で出力するときに，どうもフォント問題か何かで問題が生じる場合がある …… といっても dvips 出力とかふつーのヒトには関係ない状況におけるものだが． R --vanilla < 作図プログラム で出力させれば何も問題ナシ．

お茶部屋でちょっと一休み．知床調査隊が本日もどってきて，その放出行動食であるリッツクラッカーなどかじり，お茶をのむ．高校野球ふぁんの大学院生への追従 (お茶部屋における私の社会的地位は最底辺に近いので)．
「本日も応援の甲斐あって，駒笘は逆転勝ちしちゃいましたね」
「です」 (当然です，の意)

夕刻．だんだん曇ってきた．気温低下中．よしよし．

results 節も discussion 節もぢりぢりと肥大化中．私はみょーなクセがあり，いちばん最後のパラグラフは何か「まとめ」じみたことを書かないと落ちつかない (体裁としてはそんなことはする必要ないんだが)．まとめみたいなもの，ということで抽象的になるんだけど，その作文がなかなかできない …… しかしながら，他の部分，文献とか読まずに書ける discussion はほぼでそろってきた．

1855 研究室発． 1910 帰宅． 2035 自宅発北大構内走．走ってる途中からぱらぱらと雨ふりはじめる． 2115 帰宅．体重 74.2kg．帰宅してしばらくすると本格的に降ってきた．慈雨．晩飯．

語学堪能なる院生から，昨日の quasi 夏風邪なる表現はなかなかよろしい，とのご講評をいただく．なンと学部生のときにラテン語の授業まで ……

[今日の運動]

北大構内走 40 分間．

[今日の食卓]

朝 (0820): クロワッサン．
昼 : 食ってない
晩 (2200): 米麦 0.8 合．昨日の麻婆豆腐ののこり．ネギ・卵の炒めもの．

ぎょーむ日誌 2006-08-17

2006 年 08 月 17 日 (木)