glmmPQL() 空间自相关输出的解释

机器算法验证 r 混合模式 解释 自相关 随机效应模型
2022-04-01 18:02:35

我正在使用随机效应和空间自相关对二项式数据进行建模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)

在此处输入图像描述

让我很困惑。

  1. 我如何解释summary()结果?

  2. 与残渣相比,如何处理nlme::Variogram()显示完全不同形状的输出。svgm/summary() 输出?

  3. 我做错了什么?

0个回答
没有发现任何回复~