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. el 3 ene. 2020 a las 20:52
  • 1
    Chequea esto: r.789695.n4.nabble.com/… el 3 ene. 2020 a las 20:53
  • Muchisimas gracias! Funcionó perfecto!
    – Calfiao
    el 3 ene. 2020 a las 21:00
  • 1
    Me alegro, edito un poco tu pregunta y dejo la respuesta. Saludos. 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
    el 7 may. 2020 a las 18:19

Tu Respuesta

Al pulsar en “Publica tu respuesta”, muestras tu consentimiento a nuestros términos de servicio, política de privacidad y política de cookies

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