关于spatstat的Ripley K函数的错误

Question

我想在spatstat上使用Kest函数提供一些帮助。我试图为我的模拟数据计算de Ripley的K.当我尝试运行该功能时，R会向我显示一些错误，即：

  

错误：7段不完全位于窗口内。

但是当我绘制地图时，所有点都完美地显示在观察窗口内。 我使用以下代码：

library(spatstat)
library(sp)
library(maptools)

pemt.area <- readShapePoly("Fragment/Fragment_SIRGAS.shp")
pemt.area <- as(pemt.area, "owin")
plot(pemt.area)

tree.simu <- read.table("simulation.txt", h=T)
tree.simu
tree.simu.ppp<- ppp(x=tree.simu$X,y=tree.simu$Y,window=pemt.area)
plot(tree.simu.ppp)

ripley <- Kest(tree.simu.ppp, correction="isotropic", domain=pemt.area, var.approx=TRUE, ratio=TRUE)
plot(ripley, main = "K function")

我试着查看小插图，但我没有成功解决这个问题。

以下链接包含数据文件。如果您还需要其他信息，请告诉我：https://drive.google.com/folderview?id=0B_cWshGkkAzjbHBnemNuMzh1Q1k&usp=sharing

谢谢，非常感谢您的帮助！

Answer 1

我无法重现你的问题。这是我生成的K图。

就我而言：OSX，R 3.1.0。我建议的是确保您的软件包是最新的，并仔细检查您使用的文件是否与链接中的文件相同。

Answer 2

这与Kest或ppp无关。

请注意，错误消息是关于细分，而不是积分。

这是来自spatstat函数psp的错误消息，它会创建线段的空间模式。该错误表示多边形数据不符合要求。最有可能在致电readShapePoly或as.owin期间发生。要跟踪它，您可以在错误报告后立即键入traceback()，这将显示错误发生的位置。

首先，更新所有包是明智的。