点パターンのモデリング。空間点過程が
ここで、2次の強度(Gatrell et al., 1995) は、点どうしの空間的依存関係を示す尺度である。
定常仮定では、λ は有界領域 A のどこでも一定。関数 intensity は強度を求める関数。方法は以下の4つから選択して引数 method に与える。
点どうしの関係が、空間の中でどのように変化しているか記述する量
K(d) = λ-1 E [ イベント数 ≦ 任意のイベントの距離 d ] (λ は強度、 E は期待値)
を K関数と呼ぶ。
K関数の強みは、いくつかのモデルに対して K(d) の理論値が分かっていることである。例えばCSRでは K(d) = πd2 である。クラスタ型であれば、小さな d に対して K(d) > πd2 であろうし、規則型であれば K(d) < πd2 になるだろう。
ここでは n は A 内の点の数、 |A| は A の面積、di,j は i番目の点と j番目の点との距離、 wi,j は点 j を中心にして j を通る円の、領域 A に含まれている部分の面積の割合。 Id(di,j) は指示関数で、 di,j ≦ d ならば1、その他は0。関数 Khat でK関数の推定値を求めることができる。
CSRの K関数は2次関数なので、1次関数にした
の方が分かりやすい。これを L関数と呼ぶ。 CSR の L関数は L(d) = d となる。
Ghat(d) のときと同じように、Khat(d)のCSRの検定(に類するもの)を、全く同じ方法で行える。
ご意見、ご感想はNTTデータ数理システムS-PLUSグループ <splus-info@ml.msi.co.jp>