如何恢复glmer模型的残差

Question

我正在研究一个带有 glmer 包 lme4 和 year 作为随机效应（截距）的逻辑混合模型，以获取空间自相关（SAC） .我想使用 Moran's I 指数测试模型残差上的 SAC，但我无法恢复模型残差的系列。如何从 glmer 类中恢复任何类型的残差以执行空间自相关测试？

Answer 1

所有 merMod 对象都有一个 residuals() 方法。

假设您拟合的模型对象名为 mymod，那么将残差打印到控制台是 residuals(mymod)。

要将结果保存到您的数据框中

mydata$mod1resid <- residuals(mymod)