如何從混合模型（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)

在這里，我返回具有相對較大的學生化殘差的觀測值（考慮為離群值）。 看起來不錯？