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 <[email protected]> Para: [email protected] 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 [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
