【问题标题】:Ripleys K not plotting correctly?Ripley's K 绘图不正确?
【发布时间】:2018-11-08 21:00:54
【问题描述】:

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

这是我的学习窗口。这是一个不规则的形状,因为我正在分析沿样带线的点模式。

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

【问题讨论】:

  • 我已经根据这些新信息更新了我的答案。希望对你有帮助。结论:没有错误,只是观察窗口非常拉长,无法计算更远距离的一些估计值。

标签: sp spatstat


【解决方案1】:

Kest 中不太可能(但并非不可能)存在导致此问题的简单错误,因为此特定功能已被许多用户密集测试。您更有可能有一个不规则的观察窗口,并且有一个数学原因导致无法在所有距离上计算各种估计值。请添加您的点模式的绘图/摘要,以便我们了解观察窗口(或者更好地访问观察窗口)。

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

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

更新:

您的窗口确实非常不规则,并解释了为什么它没有在远距离产生一些修正。我猜你的样带只有几米宽,你正在考虑的距离可达 50m。边界校正只能针对不超过样带一半宽度的距离进行计算。

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

【讨论】:

  • 好的,谢谢,这很有帮助。澄清一下,High、Medium 和 Low 指的是属于一个特定标记的点的密度。采样是在整个样带线上均匀进行的,我想看看属于该特定标记的点是否聚集在一起。你有什么不同的方法推荐吗?
猜你喜欢
  • 2015-07-01
  • 2020-07-21
  • 1970-01-01
  • 2016-07-20
  • 2018-07-09
  • 2020-06-01
  • 2012-01-16
  • 2014-01-19
  • 2012-07-01
相关资源
最近更新 更多