En mucho manuales se recomienda trabajar con variables centradas, ya que 
los Betas resultantes no están correlacionados. Comprueba

summary(m1)$varFix ; summary(m2)$varFix

Si el diseños no esta perfectametne balanzeado (estimaciones a 
diferentes edades), se introduce ruido en la estimación de los errores, 
pero no afecta significativamente  a los efecto fijos. Puedes 
comprobarlo, con esta modificación de tu código

Orthodont$age2 <-  jitter( Orthodont$age)
m3<-lme(distance~I(age2-mean(age2)),data=Orthodont,random=~1|Subject)
summary(m3))$varFix

Por los demás los tres modelos son prácticamente equivalentes.

En diseños longitudinales con tiempos de medición irregulares de un 
proceso de crecimiento (pendientes distintas de cero), plantear central 
en la media de los individuos creo que no tiene ningún sentido.

Saludos

El 12/07/2012 16:31, Jaime Otero Villar escribió:
> 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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es


-- 
+--------------------------------------------------------------
| Francisco J. Viciana Fernández
| Coordinador del Registro de Población
| Servicio de Estadísticas Demográficas y Sociales
| Instituto de Estadística y Cartografía de Andalucía
| Leonardo Da Vinci, nº 21. Isla de La Cartuja.
| 41071 SEVILLA.
| franciscoj.viciana@juntadeandalucia.es
| Tlf.: +(34) 95 503 38 21
+--------------------------------------------------------------


	[[alternative HTML version deleted]]

