r - 在没有重叠点的情况下绘制 spatstat 中的残差

标签 r spatial spatstat

这可能非常简单,但我花了很多时间试图弄清楚但没有任何运气,也许任何人都可以帮助我。

我通过 spatstat 中的 ppm() 函数拟合了点模式模型(下面的可重现代码),当我绘制残差时,它会自动在残差图像上绘制点,这使得很难看到任何东西。有谁知道如何避免它?

代码:

library(spatstat)
pattern <- rpoispp(300)
cov <- rnoise(rgen = rnorm, dimyx=32, mean=2, sd=1, w = pattern$window)
fit <- ppm(pattern ~ cov)
res <- residuals.ppm(fit, type = "raw")
plot(res, how = "imagecontour")

最佳答案

数据点没有被“过度绘制”在残差之上:残差测量包括每个数据点的质量“原子”,以及平滑的密度,因此情节是正确的。

如果问题是因为代表原子的符号太大而看不到细节,那么您可以使用参数之一markscalemaxsize 将传递给 plot.ppp

话又说回来,如果有很多数据点,最好只平滑残差度量。如果 res 是您计算的残差度量,则尝试 plot(Smooth(res))。有关详细信息,请参阅 Smooth.msr 的帮助。

如果您确实需要提取残差度量的平滑密度分量,您可以遵循 Ege 的建议,或者使用 with.msr。例如

with(res, Smooth(qlocations %mark% density)) 

给出表示残差测量的连续分量的图像。

这些注释仅适用于原始残差,其中所有原子的质量都为 1。对于其他类型的残差,原子的质量不相等,显示它们变得更加重要。

关于r - 在没有重叠点的情况下绘制 spatstat 中的残差,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36544600/

相关文章:

r - 在 Azure Databricks 上安装 rgdal 和 rgeos

math - ℝ³ -> ℕ 有限数量值的映射

r - {spatstat }'s `nearestsegment()` 是否会产生不精确的结果?

spatial - 用作 as.linnet 对象的数据框对象

optimization - R:FAST 多变量优化包?

r - 将 Matlab/R 中的模型转换为 C++/Java

database - ORA-20085 : Different SRID's

r - 在 R 中将样方计数绘制为强度图

r - 旋转 dplyr 的 tibble 摘要的简单方法

data-structures - 移动物体的空间数据结构?