[R-br] correlação robusta em dados longitudinais

Pedro Emmanuel Alvarenga Americano do Brasil emmanuel.brasil em gmail.com
Quinta Maio 29 14:17:50 BRT 2014


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 Brasil
http://blog.ipec.fiocruz.br/lapclin-chagas/
Curriculum Lattes:  http://lattes.cnpq.br/6597654894290806
ResearchGate.net: https://www.researchgate.net/profile/Pedro_Brasil2/
Instituto Nacional de Infectologia/Instituto de Pesquisa Clínica Evandro
Chagas
Fundação Oswaldo Cruz
Rio de Janeiro - Brasil
Av. Brasil 4365,
CEP 21040-360,
Tel 55 21 3865-9648
e-mail: pedro.brasil em ipec.fiocruz.br
e-mail: emmanuel.brasil em gmail.com

---Apoio aos softwares livres
www.zotero.org - gerenciamento de referências bibliográficas.
www.broffice.org ou www.libreoffice.org - textos, planilhas ou
apresentações.
www.epidata.dk - entrada de dados.
www.r-project.org - análise de dados.
www.ubuntu.com - sistema operacional
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140529/21b5880e/attachment.html>


Mais detalhes sobre a lista de discussão R-br