<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Prof. Walmes,<br>
Me passe um comando com lme() por favor.<br>
Um abraço.<br>
André<br>
<br>
Em 3/26/aaaa 10:16, Walmes Zeviani escreveu:
<blockquote
cite="mid:CAFU=EkZZgkKj+28eLJaMCiR582dpTxa3k-M_UyoBzGJCHUZnxA@mail.gmail.com"
type="cite">
<div><font face="trebuchet ms,sans-serif">Para ter esses efeitos
você tem que usar o contraste tipo polinomial para o seu
fator. Tem diversas formas de fazer isso. Como você não deu um
CMR eu tô mandando um com lm() mesmo.</font></div>
<div><font face="trebuchet ms,sans-serif"><br>
</font></div>
<div>
<div><font face="'courier new', monospace">x <-
factor(rep(1:4,3))</font></div>
<div><font face="'courier new', monospace">y <- rnorm(x)</font></div>
<div><font face="'courier new', monospace"><br>
</font></div>
<div><font face="'courier new', monospace"># convertendo dentro</font></div>
<div><font face="'courier new', monospace">m0 <- lm(y~x,
contrast=list(x=contr.poly))</font></div>
<div><font face="'courier new', monospace">coef(m0)</font></div>
<div><font face="'courier new', monospace"><br>
</font></div>
<div><font face="'courier new', monospace">m0 <- lm(y~x)</font></div>
<div>
<font face="'courier new', monospace">coef(m0)</font></div>
<div><font face="'courier new', monospace"><br>
</font></div>
<div><font face="'courier new', monospace"># convertendo fora</font></div>
<div>
<font face="'courier new', monospace">contrasts(x) <-
contr.poly(nlevels(x))</font></div>
<div><font face="'courier new', monospace">m0 <- lm(y~x)</font></div>
<div><font face="'courier new', monospace">coef(m0)</font></div>
<div style="font-family:'trebuchet ms',sans-serif"><br>
</div>
<div style="font-family:'trebuchet ms',sans-serif">Lembrando que
usar polinômios ortogonais só tem sentido para **níveis
igualmente espaçados** do fator. Porque eles são/eram úteis?
Porque a matriz X'X fica ortogonal e dessa forma você faz
inferência para um efeito marginal aos demais (embora não faça
o mínimo de sentido julgar o termo linear depois do quadrático
etc). Na minha opinião os testes de hipótese sequênciais são
mais apropriados e não dependem dessa condição restritiva de
níveis igualmente espaçados.</div>
<div style="font-family:'trebuchet ms',sans-serif"><br>
</div>
<div style="font-family:'trebuchet ms',sans-serif">À disposição.</div>
<div style="font-family:'trebuchet ms',sans-serif">Walmes.</div>
<div style="font-family:'trebuchet ms',sans-serif">
<br>
</div>
</div>
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques
Zeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório
de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de
Estatística - Universidade Federal do Paraná</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41
3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361
3600) 1053 1173</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">e-mail: <a
moz-do-not-send="true" href="mailto:walmes@ufpr.br"
target="_blank">walmes@ufpr.br</a></span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">twitter:
@walmeszeviani</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">homepage: <a
moz-do-not-send="true" href="http://www.leg.ufpr.br/%7Ewalmes"
target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user
number: 531218</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
===========================
Aviso de confidencialidade:
===========================
- Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal n. 5.851, de 7 de dezembro de 1972, eh enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada eh ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco.
=====================
Confidentiality note:
=====================
- "This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure".
<br>
</body>
</html>