spatstat中线性网络上的复制点模式

时间:2018-10-03 03:56:41

标签: r spatstat

我有一个重复实验的点模式数据,其中每个重复实验中的点都被限制在相同的线性网络中(数据来自蛇的自行车道的每日调查:每天分别给出动物所在位置的单独点模式找到)。

我知道在ElementUI中可以同时将点过程拟合到多个点模式(使用spatstat),并且可以在线性网络上拟合点过程模型(使用mppm) ;可以同时做两个吗?据我所知,lppm将不接受mppm对象:还有另一种适合这种类型的模型的方法吗?

1 个答案:

答案 0 :(得分:1)

spatstat尚未完全支持此功能。

但是,您可以通过使用lpplinequad对象转换为正交方案,然后在超帧中使用这些正交方案而不是ppp对象来完成大部分操作。示例:

X1 <- spiders
X2 <- runiflpp(25, domain(spiders))
A <- linequad(X1)
B <- linequad(X2)
f <- function(x,y)x
H <- hyperframe(X=solist(A,B), Z=list(f,f))
fit <- mppm(X ~ Z, data=H)

mppm的大多数方法都可以正常工作;除了您无法模拟或预测拟合模型外,因为它不知道它应该在网络上。

如果有很长的ppp对象列表,则可以使用solapply( , linequad)来代替将点模式一一转换为正交方案。