
Prezados, Estou utilizando da função lme (pacote nlme) para ajustar um modelo misto. Gostaria de saber porém se é possível flexibilizar a estrutura de correlação dos efeitos aleatórios para modelos fora da função "pdMat" e suas respectivas classes ( "pdClass"). Posso (vide código a seguir) propor uma estrutura bloco diagonal, com blocos na classe "corClasses" ou outras formas e/ou funções que eu mesmo descreveria? lme(Y~X1+X2-1, data=dados.frame, random=list(subject=pdBlocked(list(corARMA (~Z1+Z2-1), corARMA(~ Z3+Z4-1))))) Grato David

Até onde sei, a nlme::lme() não permite especificações de matrizes de covariância para o efeito aleatório, apenas para os erros (último termo do modelo). Para efeitos aleatórios acredito que só permita efeitos independentes. Porém, a lme4::lmer() está sendo desenvolvida para permitir mais especificação. Dê uma olhada na documentação correspondente. No caso de matrizes de covariância geradas por parentesco (pedigree, comuns em dados de programas genéticos) tem-se já algumas alternativas que não recordo o nome. À disposição. Walmes.

Obrigado pela atenção, Gostaria de ressaltar que, se não estou enganado, o pacote nlme::lme() permite algumas estruturas para a matriz de covariância do efeito aleatório, atribuídas pela função "pdMat" e suas respectivas classes ( "pdClass"). Já para o erro aleatório temos maiores possibilidades tanto para a estrutura de correlação como para a heterogeneidade, como por exemplo as classes provindas de "corClasses". O meu interesse portanto era, se possível, como faço para transpor algumas dessas estruturas, "corClasses", para a matriz de covariância do efeito aleatório. Vou dar um olhada no pacote que me sugeriu. Att. David Em 27 de maio de 2015 19:41, walmes . <walmeszeviani@gmail.com> escreveu:
Até onde sei, a nlme::lme() não permite especificações de matrizes de covariância para o efeito aleatório, apenas para os erros (último termo do modelo). Para efeitos aleatórios acredito que só permita efeitos independentes. Porém, a lme4::lmer() está sendo desenvolvida para permitir mais especificação. Dê uma olhada na documentação correspondente. No caso de matrizes de covariância geradas por parentesco (pedigree, comuns em dados de programas genéticos) tem-se já algumas alternativas que não recordo o nome.
À disposição. Walmes.
_______________________________________________ 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)
-
David Matta
-
walmes .