[R-es] intercepto en mixed-model

[Marcuzzi, Javier Rubén]

Si me permite, agrego a su pregunta

¿y si lo quito o veo el valor a 0?

library(nlme)
data(Orthodont)
plot(Orthodont)

m1<-lme(distance~age,data=Orthodont,random=~1|Subject)
m2<-lme(distance~I(age-11),data=Orthodont,random=~1|Subject)
m3<-lme(distance~age-1,data=Orthodont,random=~1|Subject)
m4<-lme(distance~I(age-0),data=Orthodont,random=~1|Subject)
fixef(m1)
fixef(m2)
fixef(m3)
fixef(m4)
anova(m1,m2,m3,m4)

Se que el modelo sin el intercepto es más "feo", pero a mi me sirve quitar 
el intercepto para visualizar más fácil los resultados, mirando el siguiente 
ejemplo (en otros modelos al quietar el intercepto el modelo es "malo")


# Pinheiro and Bates, pp. 222-247
# p. 222
options(contrasts = c("contr.treatment", "contr.poly"))
fm1BW.lme <- lme(weight ~ Time * Diet, BodyWeight, random = ~ Time)
#vs
fmjBW.lme <- lme(weight ~ Time * Diet-1, BodyWeight, random = ~ Time)
fixef(fm1BW.lme)
fixef(fmjBW.lme)
anova(fm1BW.lme, fmjBW.lme)

Javier

-----Mensaje original----- 
From: Jaime Otero Villar
Sent: Thursday, July 12, 2012 11:31 AM
To: r-help-es en r-project.org
Subject: [R-es] intercepto en mixed-model

Hola,

Esto es más una cuestión estadística que propiamente sobre R así que
pido disculpas de antemano. La pregunta es la siguiente:

En un modelo de efectos mixtos cuando a veces ocurre correlación entre
coeficientes una solución posible es centrar la variable predictora.
En el ejemplo siguiente tomado de Pinheiro & Bates (2000) se modela en
múltiples sujetos una distancia anatómica en función de la edad:

library(nlme)
data(Orthodont)
plot(Orthodont)

m1<-lme(distance~age,data=Orthodont,random=~1|Subject)

m2<-lme(distance~I(age-11),data=Orthodont,random=~1|Subject)

El modelo 2 centra la edad tomando el valor 11 como referencia, con lo
que el intercepto del modelo poblacional resultante sería la distancia
a la edad de 11 años. En este ejemplo las mediciones de distancia
fueron tomadas en cada sujeto a la misma edad. Pero que ocurriría si
esas mediciones hubieran sido hechas a edades muy diferentes en cada
sujeto y centráramos en vez de por 11 (que en este caso podría ser la
media de todas las edades) lo hiciéramos por la media de edad
particular de cada sujeto? Como se interpretaría ahora el intercepto?
Y tendría sentido esta segunda forma de centrar si la variable
dependiente fuera muy diferente en cada sujeto?

Muchas gracias

Saludos

Jaime

[[alternative HTML version deleted]]







_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es