所以我遇到的问题很基本,但我似乎无法解决它。我一直在使用 spatstat 的 F、G、J、K 和 L 函数,我想将一个点过程的 F 函数绘制在与另一个点过程的 F 函数相同的图上。这很容易,只是我的一个点过程将函数计算为 r = 20,另一个计算为 r = 15('r' 是 x 轴上的自变量)。我希望它们都被评估为相同的 r 值,最好是 20。我该怎么做?本质上,这就是问题所在。
首先我创建包含函数值的数据框:
mydata.Fest <- Fest(mydata)
然后,出于某种原因,这行代码:
plot(mydata.Fest)
不求出与这行代码相同的“r”(自变量)值:
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 意外崩溃。如果您遇到过这个问题,或者知道解决方案,请告诉我。
谢谢!
最佳答案
这是一个常见问题解答。它在 spatstat.org FAQ 页面上进行了讨论。
spatstat 中的汇总函数(例如 Fest 和 Kest)对自变量有一个“推荐范围”,该范围可能比计算函数的值范围更短。推荐范围是函数估计在统计上可靠的范围。默认情况下,函数仅绘制到推荐的限制。 (这是空间统计中的标准做法。如果我们不这样做,那么大多数图都会看起来很奇怪,因为函数值会爆炸,可靠的信息会按比例缩小,只占据左下角的一小部分情节。)
要检查 r 的可用和推荐范围,打印函数对象(只需键入其名称)并查看输出的最后几行。
要控制绘图限制,请使用 xlim 和 ylim。
有关详细信息,请参阅研讨会笔记 www.csiro.au/resources/pf16h.html
Adrian Baddeley - 包作者
关于R (spatstat) : Built-In Functions Not Evaluating to Desired Limit,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/10933637/