owin - 等高線マップをプロットし、各等高線レベルのマーク (x、y に関連付けられている) の平均を表示する方法

Question

ポイントパターンX（「木の位置」など）の強度関数から生成されたspatstatの等高線マップがあります。この点パターンの各 x、y 座標は、対応する 3 番目のベクトル (「木の直径」など) でマークされます。-->cf image (もちろん木を表す縦線は省略可能)

等高線の各レベルのマーク (直径) の平均を異なる色で表示したいと考えています。提案？ありがとう！

Answer 1

一種のノンパラメトリック回帰を効果的に求めています。

これは、関数を使用した簡単な計算であり、データセットrhohatで実証されていlongleafます。

最初に強度推定値を計算します:Z <- density(longleaf)画像を生成しますZ。次に、コマンドZの呼び出しで共変量として扱います。rhohat

f <- rhohat(unmark(longleaf), Z)

と

g <- rhohat(unmark(longleaf), Z, weights=marks(longleaf)).

ここで比率を取り、h <- eval.fv(g/f)それをプロットしplot(h)ます。これは、推定平均樹木直径を森林密度の関数として示しています。この関数 h を Z の元の等高線に適用するには、最初にhtrue 関数に変換し、H <- as.function(h)次に を評価hZ <- eval.im(H(Z))し、最後に を行いplot(hZ)ます。