r - 在 R 的同一图中用两个样本模拟雷区

Question

我试图通过在同一个地块中绘制两个泊松分布的样本来模拟雷区，一个比另一个具有更高的强度和更小的面积。这是雷区，另一个只是被视为点的噪音（石头、洞、金属）。我无法让 R 在轴上绘制具有相同单位的点。无论我做什么，这些点都会覆盖整个情节，即使我只希望 X 点覆盖四分之一的情节。我的 R 代码如下：

library(spatstat)
Y = rpoispp(c(5),win=owin(c(0,10),c(0,10)))
X = rpoispp(c(10),win=owin(c(0,5),c(0,5)))

如果你能帮助我，请告诉我。

Answer 1

我的猜测是你正在做类似的事情：

> plot(Y)
> plot(X)

绘制点。

这样做的问题是，该类的绘图函数的默认行为ppp（这是rpoispp函数返回的）是创建一个仅包含其点的新绘图。所以第二次plot调用基本上会删除第一个图，并在不同比例的窗口中绘制自己的点。add=TRUE您可以通过设置第二个绘图的选项来覆盖此行为。所以代码

> plot(Y)
> plot(X, add=TRUE, cols="red")

应该给你类似的东西：

查看文档 ( help(plot.ppp)) 以获得更多解释和其他美化情节的选项。