1

我想要一些有关 spatstat 上的 Kest 功能的帮助。我正在尝试为我的模拟数据计算 de Ripley's 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

谢谢你，我真的很感谢你的帮助！

4

2 回答 2

1

我无法重现您的问题。这是我生成的K图。Kfunc.jpg

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

于 2014-06-29T11:44:43.573 回答
1

这与Kestor无关ppp。

请注意，错误消息是在讨论段，而不是点。

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

首先，更新所有软件包是明智之举。

于 2016-06-27T02:32:02.007 回答