Necesitaría extraer los valores ajustados de un modelo usando lmer pero solo para la porción fija del modelo y no para porción aleatoria. Busqué en muchos post pero no logro encontrar una solución. Cualquier ayuda es muy bienvenida.
Por ejemplo:
library("lme4")
fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
fitted(fm1) # Asi obtenemos los valores ajustados
La función fitted(fm1)
me da los valores ajustados pero incluye la porción fija y aleatoria del modelo.
model.matrix(<modelo>, type = "fixed") %*% fixef(<modelo>)
, donde obviamente<modelo>
es tu variable del modelo.