0

我想请你帮忙找出为什么当我使用函数信封时，我的参数不被接受，而是定义了“未使用的参数”。我使用的数据是没有标记的 ppp，我想用模拟数据和我的数据创建一个 L 函数图。这是我的 ppp 数据的代码：

 map2008MLW = ppp(xy2008_BNGppp$x, xy2008_BNGppp$y, window = IoM_polygon_MLWowin)

接着：

L2008 = Lest(map2008MLW,correction="Ripley")
OP = par(mar=c(5,5,4,4))
plot(L2008, . -r ~ r, ylab=expression(hat("L")), xlab = "d (m)"); par(OP)
L2008$iso  = L$iso  - L$r
L2008$theo = L$theo - L$r

所需的模拟次数

n = 9999

要显示的所需 p 显着性水平

 p  = 0.05

而此时的包络函数似乎不是很开心：

    EL2008 = envelope(map2008MLW[W], Lest, nsim=n, rank=(p * (n + 1)))
    Error in envelope(map2008MLW[W], Lest, nsim = n, rank = (p * (n + 1))) : 
  unused arguments (nsim = n, rank = (p * (n + 1)))

这似乎是一个通用错误，我不确定它是由包 spatstat 引起的。请帮助我找到解决方案，因为我无法继续进行分析。

非常感谢，

玛蒂娜

4

1 回答 1

0

论据rank应该是nrank。

此外，示例中的显着性水平和论点之间的关系nrank也不正确。对于双边检验，显着性水平为alpha = 2 * nrank/(nsim+1)，所以nrank= alpha * (nsim+1)/2。

您选择了0.95的显着性水平，但我假设您的意思是0.05。因此，nsim=9999您希望nrank=0.05 * 10000/2 = 250获得显着性水平为 0.05 的测试。

在这种应用中，不需要如此大量的模拟（9999）。蒙特卡洛检验在 的值较小时有效nsim。在您的示例中，我通常会使用nsim=39and nrank=1。

请参阅spatstat 书的第 10 章。

于 2015-12-31T00:44:27.103 回答