Modelos multinível com variância complexa
Boa tarde, pessoal! Tive um curso introdutório aos modelos multinível, que utilizou quase exclusivamente o MLwiN. Procurei treinar em paralelo o ajuste dos modelos com lme4, o que não foi difícil para modelos com variância simples. Mas, quando o negócio chegou na variância complexa, não consegui encontrar correspondência no lme4 do que me foi ensinado no MLwiN. Nem sei se não seria necessário utilizar nlme. Um exemplo de eventual necessidade seria em algum estudo longitudinal, onde a variabilidade das medidas (por exemplo, altura; primeiro nível) de cada pessoa (segundo nível) poderia mudar ao longo do tempo. Outro exemplo seria verificar se a variabilidade da variável de reposta é diferente conforme alguma variável explanatória dicotômica ou não. Onde encontro material sobre como fazer isso no R? Obrigado, [1]Leonardo Ferreira Fontenelle References 1. http://lattes.cnpq.br/9234772336296638
Leonardo, Há alguns anos eu tive que trabalhar num modelo longitudinal multinível em que usei o pacote nlme (dê uma olhada em pdMatrix). Acho que há outras opções como MCMCglmm, mas eu não o conheço. De cabeça, eu me lembro de ter consultado o livro do Pinheiro e Bates. Daniel 2014-09-27 14:25 GMT-03:00 Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br>:
Boa tarde, pessoal!
Tive um curso introdutório aos modelos multinível, que utilizou quase exclusivamente o MLwiN. Procurei treinar em paralelo o ajuste dos modelos com lme4, o que não foi difícil para modelos com variância simples. Mas, quando o negócio chegou na variância complexa, não consegui encontrar correspondência no lme4 do que me foi ensinado no MLwiN. Nem sei se não seria necessário utilizar nlme.
Um exemplo de eventual necessidade seria em algum estudo longitudinal, onde a variabilidade das medidas (por exemplo, altura; primeiro nível) de cada pessoa (segundo nível) poderia mudar ao longo do tempo. Outro exemplo seria verificar se a variabilidade da variável de reposta é diferente conforme alguma variável explanatória dicotômica ou não.
Onde encontro material sobre como fazer isso no R?
Obrigado,
Leonardo Ferreira Fontenelle
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Obrigado pela resposta, Daniel. Pelo que entendi passando o olho, o que eu preciso é especificar uma varFunc para o argumento nlme. É isso mesmo? -- Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638 Em Sáb 27 set. 2014, às 14:57, Daniel Marcelino escreveu:
Leonardo,
Há alguns anos eu tive que trabalhar num modelo longitudinal multinível em que usei o pacote nlme (dê uma olhada em pdMatrix). Acho que há outras opções como MCMCglmm, mas eu não o conheço. De cabeça, eu me lembro de ter consultado o livro do Pinheiro e Bates.
Daniel
2014-09-27 14:25 GMT-03:00 Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br>:
Boa tarde, pessoal!
Tive um curso introdutório aos modelos multinível, que utilizou quase exclusivamente o MLwiN. Procurei treinar em paralelo o ajuste dos modelos com lme4, o que não foi difícil para modelos com variância simples. Mas, quando o negócio chegou na variância complexa, não consegui encontrar correspondência no lme4 do que me foi ensinado no MLwiN. Nem sei se não seria necessário utilizar nlme.
Um exemplo de eventual necessidade seria em algum estudo longitudinal, onde a variabilidade das medidas (por exemplo, altura; primeiro nível) de cada pessoa (segundo nível) poderia mudar ao longo do tempo. Outro exemplo seria verificar se a variabilidade da variável de reposta é diferente conforme alguma variável explanatória dicotômica ou não.
Onde encontro material sobre como fazer isso no R?
Obrigado,
Leonardo Ferreira Fontenelle
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Eu não tenho o livro aqui além de fazer uns 3 anos que não ajusto um modelo multinível. Mas eu acho que é por aí mesmo. Por exemplo, considerando que a variância dentro de um grupo aumente com a altura mais ou menos de forma linear, é possível criar um componente disso: vf <- varFixed(~ altura) vf <- initialize( altura, data=data) Daniel 2014-09-27 15:13 GMT-03:00 Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br>:
Obrigado pela resposta, Daniel.
Pelo que entendi passando o olho, o que eu preciso é especificar uma varFunc para o argumento nlme. É isso mesmo?
-- Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638
Em Sáb 27 set. 2014, às 14:57, Daniel Marcelino escreveu:
Leonardo,
Há alguns anos eu tive que trabalhar num modelo longitudinal multinível em que usei o pacote nlme (dê uma olhada em pdMatrix). Acho que há outras opções como MCMCglmm, mas eu não o conheço. De cabeça, eu me lembro de ter consultado o livro do Pinheiro e Bates.
Daniel
2014-09-27 14:25 GMT-03:00 Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br>:
Boa tarde, pessoal!
Tive um curso introdutório aos modelos multinível, que utilizou quase exclusivamente o MLwiN. Procurei treinar em paralelo o ajuste dos modelos com lme4, o que não foi difícil para modelos com variância simples. Mas, quando o negócio chegou na variância complexa, não consegui encontrar correspondência no lme4 do que me foi ensinado no MLwiN. Nem sei se não seria necessário utilizar nlme.
Um exemplo de eventual necessidade seria em algum estudo longitudinal, onde a variabilidade das medidas (por exemplo, altura; primeiro nível) de cada pessoa (segundo nível) poderia mudar ao longo do tempo. Outro exemplo seria verificar se a variabilidade da variável de reposta é diferente conforme alguma variável explanatória dicotômica ou não.
Onde encontro material sobre como fazer isso no R?
Obrigado,
Leonardo Ferreira Fontenelle
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Obrigado! -- Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638 Em Sáb 27 set. 2014, às 21:50, Daniel Marcelino escreveu:
Eu não tenho o livro aqui além de fazer uns 3 anos que não ajusto um modelo multinível. Mas eu acho que é por aí mesmo. Por exemplo, considerando que a variância dentro de um grupo aumente com a altura mais ou menos de forma linear, é possível criar um componente disso: vf <- varFixed(~ altura) vf <- initialize( altura, data=data)
Daniel
2014-09-27 15:13 GMT-03:00 Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br>:
Obrigado pela resposta, Daniel.
Pelo que entendi passando o olho, o que eu preciso é especificar uma varFunc para o argumento nlme. É isso mesmo?
-- Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638
Em Sáb 27 set. 2014, às 14:57, Daniel Marcelino escreveu:
Leonardo,
Há alguns anos eu tive que trabalhar num modelo longitudinal multinível em que usei o pacote nlme (dê uma olhada em pdMatrix). Acho que há outras opções como MCMCglmm, mas eu não o conheço. De cabeça, eu me lembro de ter consultado o livro do Pinheiro e Bates.
Daniel
2014-09-27 14:25 GMT-03:00 Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br>:
Boa tarde, pessoal!
Tive um curso introdutório aos modelos multinível, que utilizou quase exclusivamente o MLwiN. Procurei treinar em paralelo o ajuste dos modelos com lme4, o que não foi difícil para modelos com variância simples. Mas, quando o negócio chegou na variância complexa, não consegui encontrar correspondência no lme4 do que me foi ensinado no MLwiN. Nem sei se não seria necessário utilizar nlme.
Um exemplo de eventual necessidade seria em algum estudo longitudinal, onde a variabilidade das medidas (por exemplo, altura; primeiro nível) de cada pessoa (segundo nível) poderia mudar ao longo do tempo. Outro exemplo seria verificar se a variabilidade da variável de reposta é diferente conforme alguma variável explanatória dicotômica ou não.
Onde encontro material sobre como fazer isso no R?
Obrigado,
Leonardo Ferreira Fontenelle
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (2)
-
Daniel Marcelino -
Leonardo Ferreira Fontenelle