如何在 spatstat 中以 1 km 分辨率进行点过程模型

Question

我想进行点过程分析，以对 1 平方公里分辨率的点模式进行建模。我有三层协变量作为像素图像命名； rd、pd 和 ras。 PPP 对象位于包含 [-3.1523926, -2.5752286] x [53.31128, 53.70412] 单位的多边形边界上

我已经尝试过 quadscheme 函数 Q <- quadscheme(data, method="grid", eps=1) 将 quad 类 Q 传递到 ppm 公式 model1 <- ppm(Q~ras/rd+pd)

它不能返回任何结果。任何帮助都非常感谢

Answer 1

你说ppm 命令没有返回任何结果。但是你使用了一个作业，model1 <- ppm(...)。这不会产生任何打印输出。您可以打印结果：尝试输入 model1 或 print(model1) 或 coef(model1)。

您说您想“以 1 平方公里的分辨率对点模式进行建模”。这是否意味着点模式数据是 (a) 四舍五入到最近 1 公里的坐标，(b) 像素网格中每个 1 公里像素内是否存在点的指标，(c) 以合理精度记录的空间坐标，但您想通过离散化到 1km 像素进行分析？

如果 (c) 是正确的，那么您可以使用 spatstat 函数 slrm 而不是 ppm。只需 slrm(data ~ ras + pd + rd) 即可。

quadscheme(data, method="grid", eps=1) 命令不会离散空间坐标，也不能帮助您实现目标 (c)。情况很复杂;我建议您根本不要使用quadscheme，因为这是一个相当高级的功能。相反，如果data 已经是"ppp" 类的点模式，则可以在调用模型拟合函数ppm 时使用data 而不是Q。

示例数据会有所帮助。