Я запустил надежную модель смешанных эффектов с помощью команды rlmer, реализованной в пакете robustlmm. В пакет входит функция построения графиков plot.lmerMod для остаточного анализа, в результате получаются графики ggplot объектов. Как вы можете видеть из приведенного ниже примера, мне удалось добавить опорную линию на график зависимости QQ от остатков (см. Ниже). Однако я не уверен, как это сделать для графика QQ со случайными эффектами (RE). Через функцию getMe доступно несколько «компонентов». Я извлек условный режим из «сферической» переменной случайных эффектов , которая представлена `` u```. Линии на результирующем графике выглядят хорошо, но я действительно не уверен, что это правильно?
library(lme4)
library(robustlmm)
library(ggplot2)
fit <- rlmer(Reaction ~ Days + (Days|Subject), sleepstudy)
# Q-Q versus residuals
QQ1 <- plot.rlmerMod(fit, which=2)
QQ1[[1]]+
stat_qq_line(aes(sample=resid(fit)))
# Q-Q versus RE
RE <- getME(fit,name = c("u"))
QQ2 <- plot.rlmerMod(fit, which=3)
QQ2[[1]]+
stat_qq_line(aes(sample=RE2))
