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



我想做一个点过程分析,以1平方公里的分辨率对点模式进行建模。我有三层协变量作为像素图像命名;RD,PD和RAS。PPP 对象位于包围 [-3.1523926, -2.5752286] x [53.31128, 53.70412] 单位的多边形边界上

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

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

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

您说您想"在 1 公里正方形分辨率上对点图案进行建模"。这是否意味着点模式数据是 (a( 四舍五入到最接近的 1km 的坐标,(b( 像素网格中每个 1km 像素内是否存在点的指标,(c( 以合理的精度记录的空间坐标,但您希望通过离散化为 1km 像素来分析?

如果 (c( 是正确的,那么您可以使用spatstat函数slrm而不是ppm。只是slrm(data ~ ras + pd + rd)应该工作。

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

示例数据会有所帮助。

最新更新