Question

我想请求你的帮助，找出为什么当我使用函数信封时，我的参数不被接受，但定义了＆＃34;未使用的参数＆＃34;。我使用的数据是没有标记的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包引起的。请帮助我找到解决方案，因为我无法进行分析。

非常感谢，

马丁娜

Answer 1

论据collab.routing.route_provider: class: CollaborationBundle\Routing\RouteProvider arguments: - "@annotation_reader" - "@doctrine.orm.entity_manager" - "@module.provider" collab.routing.url_generator: class: CollaborationBundle\Routing\UrlGenerator应为rank。

在示例中，显着性级别和参数nrank之间的关系也不正确。对于双面测试，显着性级别为nrank，因此alpha = 2 * nrank/(nsim+1)。

您选择了 0.95 的显着性水平，但我认为您的意思是 0.05 。因此，对于nrank= alpha * (nsim+1)/2，您希望nsim=9999获得显着性水平为0.05的测试。

在这种应用中不需要如此大量的模拟（9999）。蒙特卡罗测试对nrank=0.05 * 10000/2 = 250的小值有效。在您的示例中，我通常会使用nsim和nsim=39。

见the spatstat book的第10章。

包络函数（spatstat） - 错误＆＃34;未使用的参数＆＃34;

1 个答案: