R (spatstat) : Built-In Functions Not Evaluating to Desired Limit

标签 r function plot ppp spatstat

所以我遇到的问题很基本,但我似乎无法解决它。我一直在使用 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/

相关文章:

r - 通过订购数据框来产生奇怪的行为

C函数没有返回值

python - 如何将传奇元素彼此相邻放置?

r - 如何根据 R 中的位置向量替换字符串向量中的字符?

r - 如果循环中出现错误,则创建 "n"的向量并继续

r - 如果R中数据帧为空,如何编写if/else语句

python - 是否可以将 Python Seaborn 与 plotly 集成

ios - 添加每个 UITableViewCell 单元格中每次出现的键

python - 非常基本的带参数的函数

javascript - 无法使用 d3 创建曼哈顿图