<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Obrigado pela dica Walmes,<div><br></div><div>no entanto, estou trabalhando com um arquivo de dados extraído Searle et al. (1992) e citado por Marcelino & Lemma (2000). E quando comparo as estimativas utilizando a estrutura padrão do R (Componentes de variância - Variâncias heterogêneas) os resultados batem. Mas quando tento utilizar a opção de simetria composta os resultados divergem dos encontrados pelos autores. Portanto, acho que devo estar declarando algo errado. Os altores <span style="font-size: 12pt;">Marcelino & Lemma (2000) utilizaram o aplicativo SAS.</span></div><div><span style="font-size: 12pt;"><br></span></div><div><span style="font-size: 12pt;">Segue abaixo o CMR:</span></div><div><span style="font-size: 12pt;"><br></span></div><div><span style="font-size: 12pt;">data = </span>structure(list(fa = c(1, 1, 1, 1, 2, 2, 2, 2), fb = c(1, 1, 2, </div><div>3, 1, 2, 2, 3), rep = c(1, 2, 1, 1, 1, 1, 2, 1), y = c(7.9, 8.1, </div><div>6, 2, 8, 4.8, 7.2, 12)), .Names = c("fa", "fb", "rep", "y"), row.names = c(NA, </div><div>-8L), class = "data.frame")</div><div><span style="font-size: 12pt;"><br></span></div><div>require(nlme)</div><div><span style="font-size: 12pt;"><br></span></div><div>int = interaction(fa,fb) </div><div>mod <- lme(y ~ fa, random=list(~1|fb,~1|int))</div><div>summary(mod) # resultado bate<span style="font-size: 12pt;"> </span></div><div><span style="font-size: 12pt;"><br></span></div><div><span style="font-size: 12pt;">int <- interaction(fa,fb)</span></div><div>mod1 <- lme(y ~ fa, random=list(~1|fb,~1|int), correlation=corCompSymm())</div><div><span style="font-size: 12pt;">summary(mod1) # resultado não bate</span></div><div><span style="font-size: 12pt;"><br></span></div><div><span style="font-size: 12pt;">Gostaria de saber se eu posso obter todas as seguintes estruturas no pacote nlme(): não estruturada (UN), </span>Toeplitz<span style="font-size: 12pt;"> (T) e </span>Huynh-Feldt.</div><div><span style="font-size: 12pt;"><br></span></div><div><font size="3">Agradeço desde </font>já<font size="3"> por toda ajuda!!!</font></div><div><font size="3"><br></font></div><div><font size="3">Att.</font></div><div><font size="3"><br></font></div><div><font size="3">Tiago.</font></div><div><br><div>#################################################################</div><div> </div><div>Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES</div><div> </div><div>Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas</div><div> </div><div>################################################################# </div><br><br><div><hr id="stopSpelling">Date: Tue, 11 Nov 2014 06:26:37 -0200<br>From: walmeszeviani@gmail.com<br>To: r-br@listas.c3sl.ufpr.br<br>Subject: Re: [R-br] Matrizes de variâncias e covariâncias do pacote nlme<br><br><p dir="ltr">São várias. Veja as respectivas documentações. </p>
<p dir="ltr">help(corClasses)<br>
help(varClasses)</p>
<p dir="ltr">À disposição.<br>
Walmes. </p>
<div class="ecxgmail_quote">Em 10/11/2014 21:19, "Tiago Souza Marçal" <<a href="mailto:tiagosouzamarcal@hotmail.com">tiagosouzamarcal@hotmail.com</a>> escreveu:<br><blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex;">
<div><div dir="ltr">Boa noite pessoal,<div><br></div><div>gostaria de saber quais as estruturas de matrizes eu posso obter combinando os argumentos </div><div><br></div><div>correlation e weights nas funções lme e gls?? </div><div><br></div><div>Att.</div><div><br></div><div>Tiago. <br><br><div>#################################################################</div><div> </div><div>Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES</div><div> </div><div>Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas</div><div> </div><div>################################################################# </div></div> </div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div>
<br>_______________________________________________
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.</div></div> </div></body>
</html>