获取与“iterMax”相关的错误`spatstat :: density`函数

Question

帮助错误“没有时间迭代超过iterMax;增加dtt或增加iterMax”for density函数在R

如何增加iterMax？

我使用的是spatstat库。

这是我的剧本：

d660＆lt; - density（unmark（control23_network.ppp），660）

线性网络上的My Point模式有：

1969分 具有20126顶点和21363行的线性网络 封闭窗口：矩形= [-87.63141，-87.55547] x [41.75817,41.83858]单位

它工作的一个例子是spatstat库中的样本数据集：

数据（芝加哥）

芝加哥

d60＆lt; - density（unmark（chicago），60）

积（D60）

D60

此示例数据集包含：

线性网络上的点图案 116点 多种类型，可能的类型：攻击，入室盗窃，卡特彼勒，损坏，抢劫，盗窃，非法侵入 具有338个顶点和503个线的线性网络 封闭窗口：矩形= [0.3894,1281.9863] x [153.1035,1276.5602]英尺

Answer 1

与网络规模相比，平滑带宽sigma似乎很大（660个单位）（封闭矩形的边长约为0.1个单位）。你尝试过像sigma = 0.001或类似的东西吗？ sigma的重要值注定会产生问题。帮助说：

  

计算时间很短，但随着西格玛的增加呈现平方增长。

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

最后，你的R和spatstat的版本是什么？ density.lpp的代码最近发生了很大的变化。

Answer 2

这是关于spatstat::density.lpp的问题，density用于线性网络上的点模式的方法（类lpp）。

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

此示例中合理的带宽sigma值介于0.001和0.05之间。