4

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.

4
  • 2
    No estoy familiarizado con esto, pero por lo que leí rápidamente, entiendo que podría ser algo así: model.matrix(<modelo>, type = "fixed") %*% fixef(<modelo>), donde obviamente <modelo> es tu variable del modelo. Commented el 3 ene. 2020 a las 20:52
  • 1
    Chequea esto: r.789695.n4.nabble.com/… Commented el 3 ene. 2020 a las 20:53
  • Muchisimas gracias! Funcionó perfecto!
    – Calfiao
    Commented el 3 ene. 2020 a las 21:00
  • 1
    Me alegro, edito un poco tu pregunta y dejo la respuesta. Saludos. Commented el 3 ene. 2020 a las 21:02

1 respuesta 1

4

Según uno de los autores del paquete lme4, Douglas Bates, en este hilo, la forma de calcular esto sería haciendo:

model.matrix(fm1, type = "fixed") %*% fixef(fm1)

o en su defecto, también:

getME(fm1, "X")  %*% fixef(fm1)
1
  • Una nueva consulta. Hay forma de usar los coeficientes de los factores fijos sobre un nuevo dataset para predecir nuevos valores?
    – Calfiao
    Commented el 7 may. 2020 a las 18:19

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.