[R-es] Medidas Repetidas

[Olivier Nuñez]

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]]