ぎょーむ日誌 2007-05-03

2007 年 05 月 03 日 (木)

• 0930 起床． コーヒー． 朝飯．
• 富良野トドマツややめんどうな確認計算， 午後すぎまでかかってしまった．
• いわゆる fat-tail な散布カーネルで観測されたパターンが説明されるか， 問題． ここまでのさまざまな R 作図から 「全個体共通の散布カーネルのわけないぢゃん」 と思っていたんだけど， 数量的にそれを示さねばならぬことに．
• ということで， 新しく「散布カーネル & 倒木よしあし重きつき resampling」 法を開発した (富良野調査地のトドマツは倒木の上で更新する; 倒木が更新適地)． R の sample() 関数をつかうと簡単にできる． 手順はこう:
  1. ある母樹 i を選ぶ
  2. その母樹 i からの全倒木への距離を計算する -> v.distance
  3. 各倒木 j の sampling weight を計算する -> v.weight
     • sampling weight ∝ K_i(Rij) * (その倒木上で観測された全子供数)
     • K_i(r) : Fig.6 破線の周辺散布カーネル
     • Rij : 母樹-倒木間距離
  4. 母樹 i の子供数 ni を全倒木からカーネル & 倒木良さ重みつきで 「散布」してその最大値を得る
     max(sample(v.distance, ni, prob = v.weight))
  5. 2.-4. を 1000 回くりかえして母樹 i における「最大散布距離」分布 をもとめる
• 結果そのイチ: 子供数最大のトドマツ A089 (発見された子供数 37) の場合． fat-tail な散布カーネルを仮定した場合， 「最遠方に飛ばされた子供までの距離」 は 200m 超になる確率が高い (しかし現実は 20m ぐらいだ) ということがわかる． つまり， このぎょーかいでよく使われてる 「いつでも fat-tail 散布カーネルで OK!」 てなふうに使うのはダメ， ということ (観測されたパターンを説明できない)．

[最大母子間距離の分布]

子供数最多 (n = 37) の A089 母樹の場合． 横軸は 1000 回の各試行における最大散布距離で，灰色の棒はそ の回数 (縦軸に対応)．縦の青線が A089 の子供の観察された最大散布距 離で，fat-tail な散布カーネルから期待されるよりぜんぜん子供がとんで いないことがわかる． 破線曲線はその周辺カーネルの概形．

• 結果そのニ: 上のトドマツ A089 も含む全 79 母樹の library(lattice) 作図 (→ 全体)． 先日の 母樹からみた倒木分布 という図も参考になるでしょう．

• 皆さんにメイルだしてひとやすみ．
• 10 日前にとつぜん交通事故でなくなった David Halberstam の War in a time of peace ペイパーバック版とどく (Amazon から)． Chapter は番号だけでそっけないんだけど， 詳細な索引がついてるので助かる．
• ネット雑用． セミナー情報 更新． 興味あるヒトは誰でも参加できます．
• なんとなく， いつのまにか晩飯．
• [今日の運動]
  • なんとなくうんどう休養日
• [今日の食卓]
  • 朝 (1020): イソップの巨大クロワッサン．
  • 昼 (1350): 博多ラーメン．
  • 晩 (2030): スパゲッティー． ミートソース．