我正在使用随机效应和空间自相关对二项式数据进行建模MASS::glmmPQL()。在不考虑空间自相关的情况下绘制模型拟合的残差半变异函数会导致我得到以下半变异函数:
因此,我设置了一个 corSpatial 对象来解释我的数据中的空间自相关。
correl = corSpatial(value = c(1911, 0.03), form = ~ry + rx,
nugget = TRUE,
fixed = FALSE, type = "spherical")
在使用上面介绍的对象检查了考虑空间自相关的模型的拟合后correl,我在调用时得到以下输出summary()(注意“日期”是随机效应):
Correlation Structure: Spherical spatial correlation
Formula: ~ry + rx | date
Parameter estimate(s):
range nugget
329.0934685 0.2516632
我使用变异函数再次可视化残差(?)nlme::Variogram()
plot(nlme::Variogram(fit5, form = ~rx + ry, data = d,
maxDist = 20000, breaks = c(50, 200, 400, 600, 1000,
1500, 2000, seq(2700,29700,by = 1000)),
robust = FALSE), smooth = TRUE, span = 0.65,
showModel = FALSE)
让我很困惑。
我如何解释
summary()结果?与残渣相比,如何处理
nlme::Variogram()显示完全不同形状的输出。svgm/summary() 输出?我做错了什么?

