如何从混合模型（lmer）中提取学生化残差？

Question

这是我的模特，

MY = lmer(MY ~ Energy*Protein + (1|Sequence/Cow) + (1|Period), data = CH4CO2)

是否有可以提取学生化残差的函数？

Answer 1

以下是我根据定义所做的事情，

StudentResid = function(fit){
  res = residuals(fit)
  H = hatvalues(fit)
  sigma = summary(fit)$sigm
  sres = sapply(1:length(res), function(i) res[[i]]/(sigma*sqrt(1-H[[i]])))
  test = abs(sres) > 3
  ifelse(any(test), which(test), NA)
}

StudentResid(MY)

在这里，我返回具有相对较大的学生化残差的观察结果，将其视为异常值。看起来不错？