我有两个点模式(ppp)对象p1和p2。在p1和p2中分别有X和Y点。我在p1中拟合了一个ppm模型(以位置坐标为自变量),然后用它来预测p2中每个Y点的“强度”。
现在,我想在p2中获得在该点/区域发生事件的概率。如何将预测的强度用于此目的?
我可以使用Spatstat来做这件事吗?还有没有别的选择。
发布于 2018-12-13 08:54:40
强度是每单位面积的预期点数。在小区域(例如像素)中,您只需将强度乘以像素区域,即可获得像素中出现点的概率。
fit <- ppm(p1, .......)
inten <- predict(fit)
pixarea <- with(inten, xstep * ystep)
prob <- inten * pixarea
如果prob
值小于大约0.4,则此规则是准确的。
在较大的区域W
中,期望的点数是该区域上强度函数的积分:
EW <- integrate(inten, domain=W)
结果EW
是一个数字值,即W
中的预期总点数。为了得到至少一个点的概率,
P <- 1- exp(-EW)
您还可以使用带有参数interval="prediction"
的predict.ppm
计算点数的预测间隔。
https://stackoverflow.com/questions/53740217
复制相似问题