ぎょーむ日誌 2005-03-14

2005 年 03 月 14 日 (月)

• 0740 起床． 朝飯． コーヒー． 0900 自宅発． 晴． 気温高い． 0915 研究室着．
• とりあえず， カーネルのかたちを少しかえて， 苫小牧 parasitoid 計算のつづきやらせてみる． たぶん， ほとんど変わらんだろうけど．
• 母子里したうけ準備のためのしたうけ作業， POV-Ray 出力をすこしなおして …… 結局， Tgif 上で文字など追加して HTML で説明文など追加して， アップロード． 発注者たる低温研の加藤さんに連絡メイル． 時刻は 1045． とりあえず， 一件落着ということにする．

• parasitoid の randomization 計算のほうは， まだ途中までしか結果が得られていないけれど …… それで図を描いてみると， やはり予想どーりというか今までとほとんど変わっていない． つまり こういう (ムシの写真が嫌いなヒトは見ないほうがよいかも) 小さいやつらは 樹木のあいだをぶんぶん飛びまわったりはしない， ということらしい． 苫小牧にいる leafminer parasitoid はせいぜい数ミリまで， とのことで．
• 問題は， だ …… いまやってる randomization は全樹木に関するものなんだけど， これをサイズによる条件つきの randomization にしたところで， よりよい結果が得られるのかどうか， というあたりで．
• それ以前にさー， 「条件つき randomization」 ってけっこう難しいよ …… プログラミング的に． データ構造から修正せんといかんのかしらん?
• などと， R の作図プログラムをひねりながら， うだうだ考えてるうちに， 計算がぜんぶ終ってしまった． 下の図は， season1 のカーネル幅は 10m． まあ， 妥当っぽくは見えるけれど． あいかわらず樹木間の影響はなさそう， なのか?

• さてさて． 次の一手が思いつかん． 条件つき randomization か， 樹高差 (非対称的?) を考慮した重みづけか， はたまた非線形な密度依存性か．
• 次のテを考えつかないので昼飯． なぜか東北大の鳴子の木村さんがみえたので， 楽しく雑談．
• 低温研の加藤さんから図の追加発注がきたので， 作図スクリプトをてきとーに作ってみる． まあ， 「展示用」 ならこれで OK だろう， とアップロード．

• 苫小牧 parasitoid 問題は解決方向が見えないので， とりあえず関係ない雑用書類かき．
• それが片づいた， と思ったとたんに， 作図やりなおし指令 …… といってもたんに「視点」をかえるだけなんだけど．

• 苫小牧からの統計解析こんさるなメイルへの回答かき， とか． 時刻はすでに 1545 かい．
• 進展ないまま一時間経過． えーい．
• その間に， 低温研の加藤さん (名前からわかるように名古屋人) のほうでは名古屋大での報告会の準備が終ったようで …… なになに明日は「せんとれあ」 (この名前もどうにかしてほしい) で「赤みそ」そふとくりーむ?? (正確には 八丁味噌ソフトクリーム のよーで) …… 誰か， 万博うかれ名古屋人たちの無軌道なる暴走を止めてくれ ……
• 苫小牧 parasitoid モデル， 策がつきてしまったんで， 試みに現象論的な非線形モデルを導入してみる． これは他樹木から影響される側の host 密度が高ければ， その影響が「うすめられる」という指数関数をくみこんでいる …… こういうモデリングはやや不本意なんだが． 利点としては， 「こういうのやらんといかんかなー」 と考えていた修正案，
  • 条件つき randomization
  • 樹高差の考慮
  などなどの効果を部分的ながら， 「代用」できている， ということだ． つまり， この計算をやることで， 他の修正案で得られるであろう結果が だいたい推定できるようになる， というあたり．
• それから計算プログラムの変更が最小限ですむ， というのもこの追いつめられた状況ではなかなかよい． 具体的には一行だけ変更すればすむんで．
• 計算開始は 1820． さて， これで「大木効果」をウスめた場合の影響とやらがはっきりするはずだ．
• 30 分ほどで season1 の計算結果はでた …… これだけをまず図にしてみると …… こりゃ， おもしろいな． もしこの推定結果が正しいとすると (ただし randomization 検定の結果はいづれも「ゆーい」ではないんだが)， 多くの樹種では 「他樹種の大木」がそばにあると寄生率は低下する， という傾向をあらわしているかのように見える． そしてこの調査地でいう 「他樹種の大木」とは ……

  	> tapply(d$DBH.cm, d$Sp, function(x) sum(x * x * pi * 0.25))
  	      Am       Ap      Cac       Mk       Mo       Qc       Tj 
  	 3426.14  1220.24   942.21   601.10   837.45 38863.24  3959.88 
  

  ほかならぬミズナラのことだ． こいつだけは「同樹種の大木」に parasitoid を「吸いとられる」 ような推定結果になっている．

• つまり apprent competition とは「逆」の効果が発生してる可能性あるわけだ．
• もうちょいこの方向で追及することに． 以降の season の計算を中断して， さらに「小さい樹木での寄生率はどうきまるか」 に焦点をあてた推定計算 & randomization やらせてみる． 具体的には距離加重密度の非線形係数を exp(-0.01 * density) から exp(-0.02 * density) と変更してみる．
• これを再計算させてから撤退． 1920 研究室発． 1900 帰宅． 体重 74.4kg． 重い． 晩飯．
• 「さらに小さい樹木に重心をおいた」 とでもいうべき parasitoid の推定計算， 新しい結果でた． まあ， ますます上のよーな傾向が明らか …… と言えるかな． 依然としてびみょーではあるけれど． ふーむ．

• これでだいたい 「大きな樹木 → 小さな樹木に及ぼす影響， として推定されるであろうパラメーター」 の動向 (あるいは randomization 法変更などの影響) はつかめてきた． しかし …… 次の一手がうかばん． そもそも， この傾向を調べるために使ったモデルってすごく現象論だしなぁ．
• おお， 粕谷さんみずから 異端審問官ぶりを公開 しておられる …… まわりに冒涜的な統計学つかいがいたら， こっそり密告しましょう．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0810): バゲット． ブロッコリー・タマネギ・ジャガイモ・ニンニク茎・ブナシメジ・豆腐のコーンシチュー．
  • 昼 (1250): 研究室お茶部屋． 米麦 0.7 合． ブロッコリー・タマネギ・ジャガイモ・ニンニク茎・ブナシメジ・豆腐のコーンシチュー．
  • 晩 (2040): 米麦 0.7 合． 麻婆豆腐． ハクサイ・マイタケ・煮干の味噌汁．