O argumento correlation já é especificado para os efeitos aleatórios. Portanto, basta vc indicar o valor limite. modelo <- lme(valor~Gr+Idd+Escol,random=~Gr|Bloco,correlation=corAR1(0.8) ,data=base) Valeu!!! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 ----- Mensagem original ----- De: Bernardo Rangel Tura <tura@centroin.com.br> Para: r-br@listas.c3sl.ufpr.br Cc: Enviadas: Quarta-feira, 4 de Maio de 2011 19:50 Assunto: Re: [R-br] Modelo de efeito misto com correlacao AR1 On Wed, 2011-05-04 at 17:29 -0300, Walmes Zeviani wrote:
Tura,
O efeito aleatório você separa na sintaxe do argumento random=, assim random=~fator|bloco, penso (embora nunca fiz) que a mesma sintaxe se aplique no argumento correlation=.
À disposição. Walmes.
Walmes, Infelizmente não funciona modelo <- lme(valor~Gr+Idd+Escol,random=~Gr|Bloco,correlation=corAR1(, form=~Gr|Bloco) ,data=base) Error in Initialize.corAR1(X[[2L]], ...) : Covariate must have unique values within groups for corAR1 objects Você sabe o que pode estar acontecendo? -- []s Tura _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br