Ripleys K 绘图不正确?

Ripleys K not plotting correctly?

我正在使用包 spatstat 和 运行 Ripley's K (spatstat::Kest) 对我的点进行点模式分析,以查看是否存在任何聚类。然而,似乎并不是所有应该出现在图表 (kFem) 中的线都已绘制出来。例如,红线 (Kt运行s) 在 x=12 附近停止,而绿线 (Kbord) 根本没有出现。对于如何解释这一点以及是否存在错误,我将不胜感激。

这是我的书房window。它是一个不规则的形状,因为我正在分析沿着 t运行 扇形线的点模式。

这是我的点模式的密度图:

Kest 中不太可能(但并非不可能)存在导致此问题的简单错误,因为此特定功能已经过许多用户的深入测试。更有可能的是,您的观察结果 window 是不规则的,并且存在无法在所有距离上计算各种估计值的数学原因。请添加您的点模式的 plot/summary,以便我们了解观察结果 window(或者更好地访问观察结果 window)。

此外,要手动检查 K 函数的估计值,您可以将函数值 (fv) 对象转换为 data.frame 并打印它:

dat <- as.data.frame(kFem)
head(dat, n = 10)

更新:

您的 window 确实非常不规则,并解释了为什么它在远距离上没有产生一些修正。我猜你的横断面只有几米宽,你正在考虑长达 50m 的距离。边界校正只能计算距离不超过横断面的一半宽度。

使用 Kest 意味着您相信您的横断面是大 同质 点过程的一个子集(到处都具有相同的强度并且在整个 space).如果这是真的,那么 Kest 提供了对未知的真实齐次 K 函数的合理估计。但是,您提供了一个图,您将区域划分为高、中和低强度的部分,这与同质性假设不一致。与理论泊松线的偏差可能只是由于强度不均匀,而不是点之间的实际相关性。您可能应该只考虑远小于 50 的距离(您可以在调用 Kest 时设置 rmax)。