ぎょーむ日誌 2006-03-14

2006 年 03 月 14 日 (火)

• 0920 起床． 朝飯． コーヒー． 0950 自宅発． 曇． 1000 研究室着．
• 1020 から 1530 ごろまで 苫小牧 の平尾君と geoRglm と苫小牧 leafminer 寄生率問題に適用する相談 …… というか平尾君に geoRglm についていろいろと教えていただいた． この package については 先週 東京で牧さんにいろいろと教えていただいたんだけど， まだまだ私が理解できてない部分が多々あったと思い知らされたわけで ……
  • こういう geoR 系は variogram の理解というか， モデルと variogram の対応づけが必要である …… いままで「variogram なんて単なる『めやす』みたいなもんだし」 と重視してなかったんだけど， geoR 系の空間モデリングはこのあたりに (良くも悪くも) べったりと依存しているのである
  • たとえば空間相関の関数型は cov.spatial おぶぢぇくとの cov.model で指定するんだけど， ここで選べるのは variogram 世界でよくつかわれてるものである …… default は matern (Mátern だからまてるん?)
  • variogram 依存の欠陥として， 空間相関パラメーター phi の prior として平均ゼロの事前分布を選べない (はず …… あるいは強引に指定してもよい?)
  • さらに動物の分布拡大など「動いてる分布」 へのモデルあてはめを部分的に支援する trend.spatial おぶぢぇくとがある …… default は「平ら」な cte， それだけでなく座標値の多項式も使える
  • DSC 2003 という会議の proceedings に ``geoR'' Ribeiro, ``geoRglm'' Christensen たちの解説記事あり， 参考になる
  • (追記 …… うっかり書き忘れていた) geoRglm 駆使した論文も教えてもらった:

    J. van de Kassteele, A. L. M. Dekkers, A. Stein, G. J. M. Velders. 2005. Model-based geostatistical interpolation of the annual number of ozone exceedance days in the Netherlands. Stochastic Environmental Research and Risk Assessment 19:173-183. (URL)

• さらに先週は完全に誤解してたんだけど …… binom.krige.bayes や pois.krige.bayes でパラメーターのベイズ推定をやっているとき， Gaussian Random Field (GRF) の値 (s) は「観測地点」でのみ計算されており， 平面をうめつくす格子点でやっているわけではない (これは location など指定する conditional simulation mode)
• このあたり， まだまだ私が理解できてないんだが
  • GRF は「すかすか」でも 「空間びっしり」でも挙動は同じなのか?
  • 境界条件とか気にしなくてもいいのか?
  といったあたりが気になっている． 専門家にきいたほうがいいのかな．
• 平尾君発表はなんとかなりそうだし， 村上さん対策もいろいろと考えて本日は終了．
• なぜかこういう時期にかぎってメイルとか出さないといけなかったり．
• そして A 棟 8F の「真の権力者たち」にアタマ下げつつあれこれと根回し …… 来年度からの座席配分に関する件で．
• 母子里 Gibbs 林冠 MCMC 計算， 香川大の小林さんからのご指摘で気づいたんだけど， 林冠「下」の葉はもうちょい消せるな． というのも「葉っぱ以外による光減衰」があるので (つまり現状ではそれを葉っぱたちが「代行」している)．
• ついでに「見てくれ」も良さそうだからということで， 仮想葉面積も少し大きくしてみる (ダケカンバ葉 3 枚ぶん → 6 枚ぶん)． 下の図は途中経過， 150 万 MCMC step における林冠． まあ， 多少は「下」がすっきりしてる， かな? ……

• A 棟 8F hub が取り換えられたので 「ポート利用許可」とかも全部変更になった． ただしこれは管理者がわが全部やってくれた．
• 北大ネットのメイルサーヴァー， 3/8 (水) から「迷惑メール縮減」とやらのフィルタリングやってるんだが， これによって「届くはずのメイルが削除された」 疑惑が発生中 (しかも日によって削除されたり削除されなかったり?)． 「メイルだしたのに返信ない」 といったことにココロあたりある人はお手数ですが教えてください (とりあえずメイルで?)．
• 2025 研究室発． 2045 帰宅． 晩飯の準備． うんどう． 体重 72.6kg． 晩飯．
• [今日の運動]
  • エアロバイク 50 分間．
• [今日の食卓]
  • 朝 (0930): ヨーグルト．
  • 昼 (1230): 平尾君と北大百周年記念会館 1F のきゃら亭で昼定食． シーフードオムライス， サラダ・コーヒーつき． 680 円．
  • 晩 (2230): 米麦 0.8 合． タマネギ・ニンジン・タラすりみのカレー．