[R-br] Matrizes de variâncias e covariâncias do pacote nlme

Tiago Souza Marçal tiagosouzamarcal em hotmail.com
Quarta Novembro 12 01:40:23 BRST 2014


Agradeço mais uma vez pelas valiosas recomendações Walmes. Vou verificar os pontos que você mencionou.
Mas fiquei com uma ultima dúvida. Eu consegui reproduzir as estimativas de componentes de variância da função lme4::lmer() com a função nlme::lme() usando o comando list. Mas não consegui inserir os argumento do termo random (lme) no termo corCompSymm() como mostrado a baixo.
int <- interaction(fa,fb)mod1 <- lme(y ~ fa, random=list(~1|fb,~1|int), correlation=corCompSymm(form=list(~1|fb,~1|int))) # Um erro é mostrado acho que é por causa do comando list(). Não encontrei outra forma de declarar os dois fatores simultaneamente para a função corCompSymm().
Desde já agradeço.
Att.
Tiago.
################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ################################################################# 

Date: Tue, 11 Nov 2014 21:49:30 -0200
From: walmeszeviani em gmail.com
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br]	Matrizes de variâncias e covariâncias do pacote nlme

Ao comparar 2 aplicativos certifique-se:
1) a documentação de ambos se refere ao mesmo modelo/estrutura/parametrização/suposição?
2) ambos estão utilizando o mesmo algorítimo de estimação, inclusive os mesmos critérios de parada/convergência?
3) simule dados de estrutura/modelo/parâmetros conhecida e estime pelos dois aplicativos.

Das 3 estruturas que você mencionou, eu nunca vi nada relacionado com a terceira em implementações em R.

À disposição.
Walmes.
​

_______________________________________________
R-br mailing list
R-br em 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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141112/0c7d0ec5/attachment.html>


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