我在(R 库)中创建了一个无错误的 ppmspatstat
以及相关的残差。它基于大约 25,000 个数据点。但是,当我尝试使用 plot.msr 绘制残差时,我收到以下消息:
Error: all(sigma > 0) is not TRUE
关于这意味着什么以及我可以做些什么来解决问题的任何建议
谢谢
到目前为止,原始海报 (OP) 的响应不是很灵敏,因此尚不完全清楚导致此错误的原因。对于遇到此问题的其他人,可能的解决方法是在调用之前调用具有特定平滑带宽augment.msr
的-object 。例如,如果包含-object(例如 -object 的残差):msr
plot.msr
x
msr
ppm
x <- augment.msr(x, sigma=s)
plot(x)
其中s
是一个合适的平滑带宽。
希望问题在下一版本的 spatstat 中得到解决(可能会在 2014 年 10 月发布),但如果没有 OP 的帮助,就不可能知道。