ぎょーむ日誌 2006-03-08

2006 年 03 月 08 日 (水)

• 0900 起床． 今日から東京出張なんだが． といっても千歳離陸は 12 時間後 (ふむ， この時点では出発時刻をちゃんと記憶している)． コーヒー． 朝飯． 洗濯． すでに出張状態なんでだらだらと． しかし出張の準備 (といっても何もやることないんだけど) めんどう． たしか札幌-東京間は 2-3 時間ぐらい時差があったはずなんで (そう考えると日本列島も長大だね)， 出張すると生活周期がマシになるかも． 西ということはちょうど遅れる方向で都合がよろしい．
• 1110 自宅発． 雪． しかし歩道を被覆してた雪はかなりなくなってる． 1125 研究室着．
• 出張雑事のつづき． もってくものとか， 訪問先の電話番号確認とか．
• 生態学会新潟大会ポスター賞審査員ひきうけてしまった． ことわる理由をでっちあげるのがめんどうになったので． 3/25 (土) か． どわすれしそうでコワい．
• 窓の外， すごい吹雪 …… ダメだ， もう吹雪を見ても 「ああ，まだまだ冬だ」 という気分になれない．
• 返信とどこおりメイルはかなりへった． うっかりひきうけてしまった査読あるな． これは東京からの帰路とかで読んで， 翌日に報告を提出， とするしかあるまい．
• 昼飯． 苫小牧駒大付属高の甲子園出場辞退 (の原因となったうかつというか「北海道的」な卒業生たち) にご立腹の高校野球ファン大学院生弁ずる スポーツマン倫理に関するキビしいご意見を拝聴する．
• そして研究室のメイリングリストに甲山さんが雑用メイルを． 出張好きの甲山さんは現在 Pune にいる．
• 母子里 Gibbs 林冠モデル， 林冠生産力推定計算の布石として各個葉 …… と書いて計算まちがいに気づいた． ええい， ともかく MCMC 計算にはまったく影響しないんだけど (なぜかといえば MCMC 計算で重要なのは観測と対応のつく 光センサーで測定された光量であり， 個葉はそれをホンモノっぽくするための光遮蔽物にすぎないから) …… しかしながら， 個葉ごとの光合成速度とか計算するときに重要になる 「この位置にいてこの角度でついてる葉はどれぐらい光子を受け取ってるか」 の計算ですね．
• で， しくじったのは葉それぞれがてんでばらばらな方向むいてるのに， それを失念していた． つまり 「どの葉も水平においたとすればどれぐらいの受光量あるか」 なんてのを計算してんだよね． うかつだ．
• コンピューター内のダケカンバ林分は 30cm × 60cm × 60cm の大きさの計算区画「ハコ」に分割されてて， その中のそれぞれの葉はハコ内で 10 分割されてる計算セクターのうち 5 つに所属していることになっており， そのセクターごとに決まる光透過が光回廊を形成して …… と面倒な世界の中で正しい計算を実現するには， それなりの集中力が必要とされる．
• …… などとうだうだしてると， 私が「ウラ金づくりのためのカラ出張」といった違法行為を画策してるんでは， なる深刻な疑念にとらわれてしまったらしい女性大学院生が，
  「久保さんまだいるんですか， 東京いかないんですか?」
  と「せかす」ふりをして「悪事からは手をひけ」と牽制してくるんで，
  「もうすぐ出ます， ちょっとバグとりで追いつめられてるだけです」
  と答えてみると， ホントに追いつめられたような心境になってしまった．
• しかし人間は追いつめられると一時的に短期記憶容量が拡大することあるらしく …… これはバグとりにはひどく重宝される資源であり， そのおかげでバグ特定できた …… 光が直交する平面の法線ヴェクトルを規格化するの忘れてた． 阿呆阿呆だ．
• 5000 MCMC step だけ試験運転． 11 秒で終了． ざっとみたところ悪くなさそう． よーし …… make save & scp ... で一式を手もとの計算機つまり A803 (かとーオフィス) ThinkPad に転送． A801 の Dell デスクトップ機にはひきつづき 250 万 MCMC step の試験運転とサンプリングを命じておく． これは約 87 分後に終了するはず．
• 計算結果は東京からみることにして出発準備． 1810 研究室発． 雪はやんで星が見える． 1816 JR 札幌駅着． ANA チェックイン機にクレジットカードいれる． 航空券と領収書がでてくる． 1825 快速エアポートにのって同発． 帰宅時間なんでけっこう乗客おおい．
• いやはや， どたばたと出発してしまったけど， なンか忘れものとかしてないよな， とごそごそ調べてるうちに「大失敗」 いや「まぬけ失敗」に気づいた …… 出発が一時間ほど早すぎた． このまま汽車で運ばれていくと空港で二時間弱まつことになる． なんとゆーべきか， 追いつめられると人間の長期記憶ってのは錯乱したりするものなんですね ……
• 1910 なぜかしら 5 分ほど遅れて新千歳空港着． ここまで 1040 円． もーどんどん遅れてください． 期せずして空港で晩飯くうことに． 空港 3F で蕎麦をくってると， 私の搭乗する便が機材やりくりがどーのこーので 出発が 15 分ほど遅れるとの館内放送． もーどんどん遅れてください． とりあえず検問ぬける．
• 時間あるので試験運転 5000 MCMC step 計算の葉データしらべる． 個葉分布は興味ぶかい．

[光センサー vs 個葉]

縦軸は個葉の受け取ってる光量， 横軸はその個葉から見て最寄り (下方) にある光センサーでの受光量． 最初の 5000 MCMC step なので葉セット数はおよそ 4000 枚． 具体的には上述「ハコ」内の葉とそこに一個だけ設置されてる光センサーの対応である． ここからわかることは， ハコ内の葉たちと光センサーの値には相関あるけれど， ぢつはけっこうずれてるよというもの．

• なぜ同じハコ内でもこれほどばらつくのか， と言えばそれは 「葉が向いてる方角がばらばらだから」 ということになる．

[林冠上部の図，再掲]

これは昨日示した図なんだけど …… この Gibbs 林冠では葉の向きはばらばらである． ただし， これは本物ダケカンバ葉 3 枚ぶんぐらいの面積をもつ円盤葉である． この円盤葉の面積は自由に設定できる． そして着色は「同じハコの光センサーでの受光量」 (上の図の横軸) を表現している．

• 光センサーはどれも水平に設置されており， 直上 & 東西南北ななめ上からくる光を効率よく受けとる． しかしながら， 個葉は極端な場合「直立」してたりして， 同じ場所にある光センサーにくらべて 半分未満の光量しか受け取ってないんだよね．

[個葉受光量で着色]

これは東京についてから転送した 250 万 MCMC step 時点での林冠上部． 上の図との最大の違いは， 個葉受光量 (上の上の図の縦軸) で着色している． 林冠上部でも「暗い葉 (青)」が混じる．

• 光センサーより受光量の多い葉は? これは光センサーより上にあり， しかも水平に近い葉である．
• この「光センサーと近くの葉の受光量は別モノ」 は MCMC 計算の平衡状態ふきんでもかわらない．

[250 万 MCMC step 世界]

縦軸は個葉の受け取ってる光量， 横軸はその個葉から見て最寄り (下方) にある光センサーでの受光量． 葉セット数 14631 枚 (葉数でいうと 4.3 万枚ぐらい)． くらしっくな群落光合成モデルだと 赤線の上に個葉の受光量がならぶはず． この「個葉受光量のちがい」 が生産力のちがいにどう影響することになるのだろうか? …… バイアスありそうだ． 方向は計算してみないとわからん．

[観測 vs モデル]

念のためのせとくと …… MCMC 計算によって到達した Gibbs 林冠においては， 各光センサーの受光量 (縦軸はモデルで計算されたもの) は観測 (横軸) にかなり近いものになってる． これは「そうなるように」 MCMC 計算しているので， いわばあたりまえのことなんだけど． 個葉に関しては受光量の観測データないんで， 比較できない．

• 2105 ADO 26 便の B767-300 搭乗． こんな時間なのに搭乗率は九割以上． 明日おあいする伊庭さんの 「ベイズ統計と統計物理」 を再読してみると， 何だか入りこんでしまっていつ離陸したのかわからなかった． 2240 羽田空港着． 気温すごく高い．
• 京急の特急に乗って 2300 羽田空港発． 特急といっても蒲田まで各駅なんだが． 2310 京急蒲田着． ここまで 330 円． うわー …… 下車するとプラットホームにはスギ花粉仮面の男たちが． 駅の西口から出て， ややいかがわしい (←偏見) 蒲田の商店街を西に進む． 大田区内の京急の駅では線路に直交する方向に商店街が発達している．
• 12 月の出張 のときと同じ宿， 2322 東横イン蒲田東口着 (なぜ「東口」かといえば JR 蒲田からみれば 東 だから)． 保健所指導によって厨房区画設置工事中で朝の無料給食は中止だの， 区役所の点検で図面と実態のちがいがばれてしまって修正中だの， いろいろごたごたしてるよーで …… 二泊で 14280 円．
• 上の図を作るのに必要な計算結果を北大から転送する． 平衡状態での葉数 14600 枚? 今までより 1000 枚ぐらい増えたな …… つまり出発直前に修正した 「法線ヴェクトル規格化わすれ」 ばぐが他にも影響あたえてたわけだな． 三次元光計算は悪魔の迷宮だ．
• とかやってるうちにどんどん夜がふけていく． いかん， また時差が ……
• てきとーにきりあげて， ぬるい風呂にはいることにする． 時刻はすでに 2540．
• [今日の運動]
  • うんどう休養日ひたすらつづく
• [今日の食卓]
  • 朝 (0930): 米麦 0.6 合． ハクサイ・ネギ・ニラ・豆腐の煮物．
  • 昼 (1340): 研究室お茶部屋． 米麦 0.5 合． ハクサイ・ネギ・ニラ・豆腐の煮物．
  • 晩 (1930): 新千歳空港 3F の蕎麦屋 「一灯庵」 田舎ざるそば + 焼きおにぎり 840 円． 茶碗蒸しとそば湯つき．