收到与 "iterMax" 相关的错误 `spatstat::density` 函数

getting error `spatstat::density` function related to "iterMax"

R

density 函数的错误 "No of time iterations exceeds iterMax; increase dtt or increase iterMax" 帮助

如何增加 iterMax?

我正在使用库 spatstat。

这是我的脚本:

d660 <- 密度(未标记(control23_network.ppp),660)

我在线性网络上的点模式有:

1969分 具有 20126 个顶点和 21363 条线的线性网络 封闭 window:矩形 = [-87.63141, -87.55547] x [41.75817, 41.83858] 个单位

它工作的一个例子是在带有示例数据集的 spatstat 库中:

数据(芝加哥)

芝加哥

d60 <- 密度(未标记(芝加哥),60)

剧情(d60)

d60

这个样本数据集有:

线性网络上的点模式 116分 多类型,可能的类型:攻击、入室盗窃、汽车盗窃、损坏、抢劫、盗窃、侵入 具有 338 个顶点和 503 条直线的线性网络 封闭 window:矩形 = [0.3894, 1281.9863] x [153.1035, 1276.5602] 英尺

与网络规模(封闭矩形的边长约为 0.1 个单位)相比,您的平滑带宽 sigma 似乎很大(660 个单位)。您是否尝试过 sigma = 0.001 或类似的东西? sigma的大值注定会出问题。帮助说:

Computational time is short, but increases quadratically with sigma.

这可能不是根本问题,但至少在做很多其他事情之前值得尝试。

最后,您的 R 和 spatstat 是什么版本? density.lpp 的代码最近变化很大。

这是一个关于spatstat::density.lpp的问题,density线性网络点模式的方法(class lpp)。

选择的带宽sigma太大了。封闭矩形的宽度约为 0.1 个单位,而带宽为 sigma=660 个单位。这将需要大量的算法迭代,因此它拒绝了。

在此示例中,合理的带宽值 sigma 应介于 0.001 和 0.05 个单位之间。