Caro walmes,<br>Mais uma vez meu muito obrigado realmente era o que eu desejava fazer.<br>até a próxima<br><br><div class="gmail_quote">Em 23 de setembro de 2011 11:03, Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><font face="trebuchet ms,sans-serif">Fernando,<br><br>Para isso você precisa saber qual a definição dele de intercepto-inclinação geral. Penso que seja aquela obtida usando a restrição soma zero nas estimativas, mas pode ser outra. Para conseguir essas estimativas você nem precisa de pacote extra, é sou definir corretamente o tipo de contraste (ou restrição) que você quer. Evite enviar dados em anexo ou colados na mensagem. Isso não deixa o seu código reproduzível com um copia e cola. Eu teria que fazer download, blá blá, isso enrola e eu prefiro gerar dados do que ter que todo esse trabalhão. Para facilitar a vida de quem tá disposto a te ajudar, passe o dput() do seus dados, que assim com um copia e cola reproduzímos o seu código. Mostro isso no final.<br>
<br><span style="font-family:courier new,monospace">da <- expand.grid(estudo=gl(5,1), x=1:10)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">da$y <- round(rnorm(nrow(da),0,1),2)</span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">m0 <- lm(y~estudo*x, data=da)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">coef(m0)</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">m1 <- lm(y~-1+estudo/x, data=da) # mesmo espaço coluna do anterior</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">coef(m1) # são os interceptos e inclinações para cada estudo</span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">m2 <- lm(y~estudo*x, data=da,</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace"> contrast=list(estudo=contr.sum))</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">coef(m2)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">coef(m2)[c(1,6)] # intercepto e inclinação geral</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">coef(m2)[2:5] # interceptos desvios do geral para estudos de 1 a 4</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">-sum(coef(m2)[2:5]) # estudo 5</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">coef(m2)[7:10] # inclinações desvios da geral para estudos de 1 a 4</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">-sum(coef(m2)[7:10]) # estudo 5</span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">summary(m2)</span><br style="font-family:courier new,monospace"><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace"># o resultado dessa função desse ver enviado para um</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace"># código reproduzível com copia e cola</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">dput(da)</span><br style="font-family:courier new,monospace">
<br>À disposição.<br><font color="#888888">Walmes.<br><br clear="all"></font></font><div><div></div><div class="h5"><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: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a></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 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 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>
</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>