Marcelo,
para especificar el tipo de heteroscedasticidad, deberias utilizar el argumento
weights de la función lme (en vez de "correlation").
Podria ser algo del tipo:
medidasrep<-lme(Cria_ab~Cria_ab0+ Trat+Dia+Trat*Dia,random=list(ncol=
pdIdent(~1)),weights= varIdent(form=~Dia), ...)
que permite que la varianza intra varie con el dia.
Un saludo. Olivier
De: "marcelo rodriguez" <rodriguezmarcelogaston en gmail.com>
Para: "Olivier Nuñez" <onunez en unex.es>
Enviados: Lunes, 2 de Octubre 2017 14:50:05
Asunto: Re: [R-es] Medidas Repetidas
Hola Olivier!!
Muchas gracias por tu respuesta. Estuve probando con el siguiente Scrpit:
medidasrep<-lme(Cria_ab~Cria_ab0+ Trat+Dia+Trat*Dia,random=list(ncol=
pdIdent(~1)),correlation=corAR1(form=~1|ncol),method="REML",na.action=na.omit,data=medidas,keep.data=FALSE)
Pero necesitaria poder modelar la heterogeneidad de varianzas ya que con la
autoregresiva de orden 1 no puedo. Si podes darme una mano con como indicarle
una no estructurada o alguna otra forma de modelar la falta de homocedasticidad,
te lo agradeceria.
Desde ya muchas gracias!!
Saludos!
El 29 de septiembre de 2017, 5:30, Olivier Nuñez < onunez en unex.es >
escribió:
Marcelo,
si "y" es la respuesta, "t" el tratamiento y "id"
el indice del insecto en tu base de datos "Insectos".
Un modelo que toma en cuenta la estructura de tus datos es puede escribir:
modelo=lme(y~t, data = Insectos, random = ~ id)
Un saludo. Olivier
----- Mensaje original -----
De: "marcelo rodriguez" < rodriguezmarcelogaston en gmail.com >
Para: r-help-es en r-project.org
Enviados: Lunes, 25 de Septiembre 2017 19:25:05
Asunto: [R-es] Medidas Repetidas
Buenos días a todos!
Quisiera saber si me pueden ayudar con un Script. Necesito hacer un ANOVA
para medidas repetidas. Sé que el modelo es el *lme,* pero no se indicarle
las medidas repetidas.
El conjunto de datos tiene como variable respuesta volumen ingresado y como
factores tiempo y tratamientos. El tiempo se midió a las 2,4,6 y 8 min para
los distintos tratamientos, son medidas repetidas porque se observó en los
distintos tiempos el mismo insecto.
Espero que se haya entendido lo que quiero realizar y que me puedan ayudar
con el Script para poder trabajarlo en R.
Muchas gracias!!
[[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
[[alternative HTML version deleted]]