R（spatstat）：未计算到所需极限的内置函数

所以我遇到的问题很基本，但我似乎无法解决它。我一直在使用spatstat的F、G、J、K和L函数，我想将一个点过程的F函数绘制在与另一个点过程的F函数相同的绘图上。这很容易，只是我的一个点过程将函数计算为r=20，另一个计算为r=15（“r”是x轴上的自变量）。我希望它们都被评估为相同的r值，最好是20。我该怎么做？本质上，这就是问题所在

首先，我创建包含函数值的数据框：

mydata.Fest <- Fest(mydata)

不计算为与此行代码相同的“r”（自变量）值：

plot(mydata.Fest)

plot(mydata.Fest, xlim=c(0,20), ylim=c(0,0.8))

这两个图都是从xlim=c（0,20）和ylim=c（0,0.8）开始的，但出于某种原因，第二个图在x=15时停止计算函数。考虑到它们是从同一个对象绘制的，我不知道它们为什么不生成相同的绘图

作为旁注，我遇到了另一个问题：当使用F、J和K函数时，我发现调整限制（特别是自变量的限制）会导致R意外崩溃。如果你有这个问题，或者知道解决方案，请告诉我

---

谢谢

这有点奇怪。抱歉，这不是一个答案，但我认为最好让其他人留下一个你描述的行为的可复制的例子来调查

我过去所做的是制作一个一致的容器向量，用来评估分离的点模式

library(spatstat)
data(lansing)

mydata1 <- lansing[lansing$marks == "blackoak",]
mydata2 <- lansing[lansing$marks == "hickory",]

my_r <- (0:40*.002)

mydata1.Fest2 <- Fest(mydata1, r = my_r)
mydata2.Fest2 <- Fest(mydata2, r = my_r)

plot(mydata1.Fest2)
plot(mydata2.Fest2, add = TRUE)

mydata1.Fest2$r
mydata2.Fest2$r

---

我假设这不是Fest的预期行为，但我不太熟悉这项技术，也不知道它与Ripley的K有什么不同，所以我肯定会弄错，这是预期行为，原因我不知道。我还遇到了一个意外的崩溃，与您在试图设置输出设备上的轴限制时描述的类似，但我无法重现它。所以，也许你知道了什么，或者我们都疯了

关于这个问题，我联系了Adrian Baddeley博士（“spatstat”的创建者），他告诉我，当add=TRUE时，这个示例演示了plot.fv中的一个bug。这将在spatstat版本1.28-1中修复

这是一个常见问题。spatstat.org常见问题页面对此进行了讨论

spatstat中的摘要函数（如Fest和Kest）具有自变量的“推荐范围”，该范围可能小于计算函数的值范围。建议的范围是函数估计在统计上可靠的范围。默认情况下，函数仅按建议的限制打印。（这是空间统计的标准做法。如果我们不这样做，那么大多数图都会看起来很奇怪，因为函数值会放大，可靠信息会缩小，只占图左下角的一小部分。）

要检查r的可用范围和推荐范围，请打印函数对象（只需键入其名称）并查看最后几行输出

要控制打印限制，请使用xlim和ylim

有关更多详细信息，请参阅研讨会说明www.csiro.au/resources/pf16h.html

---

Adrian Baddeley-软件包作者

安迪，谢谢你的意见。我已经对函数进行了更多的实验，问题出现在不考虑移位方向的情况下。我目前不知道如何解决这个问题。@MikeZ，这可能是一个bug，我会向r-sig-geo组发送一条消息，看看他们是否能提供更好的帮助。只要确保像我在这里做的那样，使示例可复制，并在解决问题时更新此线程。

mydata1.Kest <- Kest(mydata1)
my_k <- mydata1.Kest$r
mydata2.Kest <- Kest(mydata2, r = my_k)

plot(mydata1.Kest)
plot(mydata2.Kest, add = TRUE)