[R-br] Modelo de efeito misto com correlacao AR1

[Fabio Mathias Corrêa]

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 em centroin.com.br>
Para: r-br em 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 em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br