Monte Carlo 空间隔离测试

Monte Carlo test of spatial segregation

我无法解释我从 spatstat 中的 segregation.test 方法得到的结果。但是,我有三个不同的点模式 A、B、C,我想证明 C 和 B 是相关的,而 A 和 B 不是。您可以在这张图片中看到内核对强度的估计:

但是在 R 中使用 spatstat 包计算这个我总是得到相同的 p 值,尽管测试统计量 T 不同......这怎么可能?在这种情况下,检验统计量 T 意味着什么?为什么我得到完全相同的 p 值?

我希望你能帮助我做这个 Monte Carlo 测试的错误。

检验统计量T的含义在帮助文件中解释的很清楚。你看过了吗?

?segregation.test

在 Monte Carlo 测试中没有隔离的原假设下,数据模式和模拟模式是可交换的。 p 值是根据观察到的模式的检验统计量在模式总数中的排名来计算的。在这两种情况下,您都提供了观察到的数据具有最极端的分离统计量 T,并且 p 值为 1/26 = 0.03846。

要了解详细信息,请查看提到的帮助文件以及 spatstat book 的第 10 章和第 14 章。 (不幸的是,其中 none 是免费样本章节。)

编辑:检验统计量 T 是隔离程度的度量。如果这些点是随机标记的,它往往接近于 0,如果标记分离得很好,它往往在数字上“很大”。由于没有“大”的概念,Monte Carlo p 值用于判断观察到的 T 是否大到我们应该拒绝随机标记的原假设。