Amigos de R,
Estou conduzindo uma análise de dados longitudinais com a função nlme::lme.
POr enquanto o melhor modelo que tenho é expresso como:
f7 <- lme(Velocidade ~ lado + as.factor(dia) ,data=l,random = ~1|Id/Arteria,correlation=corAR1(.8,form=~dia))
Nessa função eu gostaria que a correlçao fosse considerada para os diferentes grupos dado para a variável lado, porque que graficamente me parece ser diferente. Mas eu não consegui colocar esse argumento dentro de corAR1. Mas esse nem é o problema que me preocupa mais.
A leitura fonte que tenho indica que sempre que possivel, deve-se utilizar uma correlaçao robusta, pois assim provavelmente evita-se o erro de assumir uma estrutura incorreta da correlação, apesar de eu ter experimentado várias.
Mas eu não achei no pacote nlme qualquer função/argumento que me levasse acreditar que uma correlçao robusta esteja disponível. Alguem saberia saber no lugar de corAR1 é possivel utilizar alguma coisa que funcionasse como uma correlação robusta?
Abraço,
Dr. Pedro Emmanuel A. A. do BrasilInstituto Nacional de Infectologia/Instituto de Pesquisa Clínica Evandro Chagas