完全ランダム性(CSR: Complete Spatial Randomness, Diggle, 1983)とは、以下のとおり。
強度(単位面積当りの点の個数)が有界領域
A 内で変化しない
点どうしに相互作用はない
樹木データを樹木の種類別にプロットしてみる(下図)
→ 互いの種の生息を抑制する傾向が見られる
(図: Lansing データの点パターンを樹木の種類別にプロット)
di を i番目の点から一番近い点までの距離とする。点-点最近接距離の経験分布関数
( n
は A 内の点の個数)
を計算するのが、関数 Ghat。自動的に関数のプロットを行うので、作図したくなければ plot=F とする。
クラスタ型は y が小さいところですぐが大きくなるし、規則型では逆に
y が大きくなるまで
の値は大きくならない。理論分布(CSR)に対するプロットも可能だが、この場合、境界効果を無視するので、境界効果を取り入れたければ以下の手順を踏む。
1次元上データ(時系列など)の境界は点、すなわち0次元であったが、空間データの場合、境界が線、すなわち1次元になるため、もはやこの効果を無視することが出来ない。点パターンデータの場合、観測領域の外側にも同様に点が分布していることを考慮に入れなければならなくなる。
点パターンの作成には、関数 make.pattern を用いる。
ei を i番目の原点から一番近い点までの距離とする。原点-点最近接距離の経験分布関数を
( m
は原点の個数)
と定義する。この値は関数 Fhat を用いて計算することができる。原点どうしの距離、つまり
grid の幅は指定しなければ
となる。原点を座標指定して与えることもできる。
Fhat の解釈は Ghat と全く逆になり、x が大きいところでが大きくなればクラスタ型となる。
もし境界効果が無視できて、点の数が比較的多い場合、点パターンがCSRであれば、
G(y) = F(y) = 1 - exp( -πλy2 )
が成り立つ。ここで λ は強度。λ
を ( =
n/|A| ) で置き換えて、Fhat(y) を 1 - exp( -πλy2
) に対してプロットしてみて、直線になればCSRである。
ご意見、ご感想はNTTデータ数理システムS-PLUSグループ <splus-info@ml.msi.co.jp>