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" <rodriguezmarcelogas...@gmail.com> Para: "Olivier Nuñez" <onu...@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 < onu...@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" < rodriguezmarcelogas...@gmail.com > Para: r-help-es@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@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es