R 的 spatstat 中出现错误：点因位于外部而被拒绝

Question

我尝试获取威斯康星州的地图，并尝试对空间数据点进行克拉克-埃文斯测试。但是，当我导入威斯康星州地图并尝试将其转换为窗口时，我遇到了错误：

警告消息：244 个点因超出范围而被拒绝 指定窗口

我使用了以下代码：

wisc.map <- map("state", "wisconsin", fill=T)
w <- owin(poly=data.frame(x=rev(wisc.map$x), y=rev(wisc.map$y)))

运行后出现错误

spat = as.ppp(data, w)

很明显，我的所有数据点都在威斯康星州，如下图所示：

我在这里放了一些示例数据：

  latitude longitude
1 43.30716 -89.43137
2 43.30388 -89.42823
3 43.31845 -89.44372
4 43.33766 -89.44634
5 43.50815 -89.52010
6 43.34179 -89.01988

知道发生了什么吗？

谢谢

---

P.S 我想做克拉克埃文斯测试，因为我想检查我的数据中是否存在聚类模式。我想确保从数据中获得的热点具有统计显着性。我查了很多文献，但没有找到直接的方法来测试它，所以我必须使用Clark Evans Test...有人知道有更好的测试来检查聚类模式或热点显着性吗？

Answer 1

这不是错误，而是警告，因为输出清楚地表明了这一点。它告诉您，您提供给

ppp

的点中有 244 个位于您声称是观察窗口的多边形区域之外。它们可能非常接近边缘，而逼近威斯康星州的多边形不够详细。

注意：

spatstat

包当前不处理球面几何，因此在纬度和经度数据上使用

spatstat

函数没有意义。过去使用

spTransform

包中的

sp

来完成投影，但现在使用更现代的简单功能包

sf

可能会更方便，它具有

st_transform

功能。

旁注： 为什么要对此数据进行 Clark Evans 检验？它检查该点模式的最近邻距离是否相似，以及在同一观察窗口上具有相同预期点数的同质泊松过程。因此，参考过程将在整个威斯康星州均匀分布与您的数据大致相同数量的点，并且您的数据仅在威斯康星州西南部有点。