<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'>Boa tarde, Cezar Rabak.<br><br>Entendi melhor, porém ainda não por completo.<br><br>Mas, o meu problema agora é que não consigo mais obter o mesmo resultado. Estou importando exatamente os mesmos dados e utilizando exatamente os mesmos comandos, porém não consigo mais o mesmo resultado que enviei no fórum dia 04/09/16. Quando realizo a ANOVA o resultado gerado está sendo diferente (e está errado) daquele que apresentei anteriormente (este está correto). Gostaria de saber como posso obter o "resíduo" para inserir no comando para o teste de Shapiro Willk: "shapiro.test(resíduo)".<br><br> Seguem abaixo os comandos utilizados (acho que é o 'CMR' a que você se referiu) e o dados, alguém poderia me ajudar a encontrar o meu erro e mostrar uma solução? Já procurei muito, mas não encontro.<br>_____________________________________________________________<br>> PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T);PS<br> PARC SubP REP var<br>1 C 30 1 12<br>2 C 30 2 13<br>3 C 30 3 14<br>4 C 30 4 12<br>5 C 45 1 18<br>6 C 45 2 17<br>7 C 45 3 16<br>8 C 45 4 17<br>9 C 60 1 22<br>10 C 60 2 21<br>11 C 60 3 24<br>12 C 60 4 23<br>13 C 75 1 24<br>14 C 75 2 23<br>15 C 75 3 22<br>16 C 75 4 24<br>17 C 90 1 21<br>18 C 90 2 20<br>19 C 90 3 21<br>20 C 90 4 19<br>21 S 30 1 22<br>22 S 30 2 21<br>23 S 30 3 24<br>24 S 30 4 23<br>25 S 45 1 21<br>26 S 45 2 22<br>27 S 45 3 23<br>28 S 45 4 21<br>29 S 60 1 20<br>30 S 60 2 18<br>31 S 60 3 19<br>32 S 60 4 20<br>33 S 75 1 17<br>34 S 75 2 15<br>35 S 75 3 16<br>36 S 75 4 15<br>37 S 90 1 12<br>38 S 90 2 13<br>39 S 90 3 12<br>40 S 90 4 11<br>> attach(PS)<br><br>> SubP<-as.factor(SubP);is.factor(SubP)<br>[1] TRUE<br><br>> mod<-aov(var~PARC*SubP+Error(PARC:REP));summary(mod)<br><br>Error: PARC:REP<br> Df Sum Sq Mean Sq F value Pr(>F) <br>PARC 1 6.453 6.453 80.67 0.0706 .<br>Residuals 1 0.080 0.080 <br>---<br>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<br><br>Error: Within<br> Df Sum Sq Mean Sq F value Pr(>F) <br>PARC 1 1.7 1.67 1.535 0.226 <br>SubP 4 108.9 27.23 25.076 6.64e-09 ***<br>PARC:SubP 4 496.9 124.22 114.418 < 2e-16 ***<br>Residuals 28 30.4 1.09 <br>---<br>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<br><br>> resíduo<-residuals(mod);resíduo<br>NULL<br>___________________________________________________________________<br>Aguardo retorno.<br><br>Obrigada.<br><br>Att.,<br>Angélica.<br><br><div><hr id="stopSpelling">From: cesar.rabak@gmail.com<br>Date: Tue, 6 Sep 2016 21:49:25 -0300<br>Subject: Re: [R-br] Parcela Subdividida<br>To: angelica.ricarte@hotmail.com; r-br@listas.c3sl.ufpr.br<br><br><div dir="ltr">Angélica,<div><br></div><div>Como você deu nomes "crípticos" para seu dataframe, e não postou um CMR, dá só para fazer uma 'interpretação' das mensagens.</div><div><br></div><div>Aparentemente você quer a resposta PS[,4] para a interação entre as variáveis (grupos[?]) PS[,1] e PS[,2] e considera a combinação dos fatores PS[,1] e PS[,3] como medidas repetidas e tem 28 casos ao todo.</div><div><br></div><div>Além do mais você considera o fator PS[,1] como resposta e Error() ou seja parte da especificação das medidas repetidas.</div><div><br></div><div>A mensagem sobre o modelo Error() ser singular ocorre quando você tenta ajustar mais parâmetros que seus dados permitem...</div><div><br></div><div>Se você continuar com dificuldades, você precisa explicar melhor o que deseja fazer ou consultar uma obra sobre ANOVA com medidas repetidas.</div><div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div><div><br></div><div><br></div></div><div class="ecxgmail_extra"><br><div class="ecxgmail_quote">2016-09-05 21:05 GMT-03:00 Angélica Ricarte via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span>:<br><blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex;">
<div><div dir="ltr">Prezado Cesar,<br><div><div dir="ltr"><br>Não consegui compreender a mensagem de erro e nem o porquê que minha matriz é nula. Você pode ajudar?<br><br>Agradeço e aguardo retorno.<br><br>Obrigada.<br>Att.,<br>Angélica<br><br><div><hr>Date: Mon, 5 Sep 2016 12:52:41 -0700<br>From: <a href="mailto:ml-node%2Bs2285057n4666617h13@n4.nabble.com" target="_blank">ml-node+s2285057n4666617h13@<wbr>n4.nabble.com</a><br>To: <a href="mailto:angelica.ricarte@hotmail.com" target="_blank">angelica.ricarte@hotmail.com</a><br>Subject: Re: [R-br] Parcela Subdividida<span><br><br>
<div dir="ltr">Angélica,<div><br></div><div>Você entendeu a mensagem de erro?</div><div><br></div><div>Não é possível analisar uma matriz nula....</div><div><br></div><div><br></div></div></span><div><br><div><div><div class="h5">2016-09-04 1:22 GMT-03:00 Angélica Ricarte via R-br <span dir="ltr"><<a rel="nofollow" target="_blank">[hidden email]</a>></span>:<br></div></div><blockquote style="border-left:2px solid #cccccc;padding:0 1em;"><div><div class="h5">
<div><div dir="ltr">Prezados,<br><br>Estou tentando realizar uma análise de variância para um experimentos em parcelas subdivididas em DIC, mas, apesar de ser gerado um resultado, aparece uma mensagem de erro. Quando tento realizar o cálculo dos resíduos também dá erro e, consequentemente, não consigo fazer o teste de normalidade. Alguém pode me ajudar? Vejam como realizei as análises e os seus resultados:<br><pre style="font-family:'Lucida Console';font-size:16pt !important;border:none;word-break:break-all;white-space:pre-wrap !important;line-height:25px;color:rgb(0,0,0);font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;text-align:-webkit-left;text-indent:0px;text-transform:none;word-spacing:0px;background-color:rgb(255,255,255);"><font style="font-size:10pt;" size="2"><span style="color:blue;"><font color="#000000">#ANOVA:</font><br>> mod<-aov(PS[,4]~PS[,1]*PS[,2]+<wbr>Error(PS[,1]:PS[,3]));summary(<wbr>mod)
</span></font><font style="font-size:10pt;" size="2"><span style="color:rgb(197,6,11);">Warning message:
</span></font><font style="font-size:10pt;" size="2"><span style="color:rgb(197,6,11);">In aov(PS[, 4] ~ PS[, 1] * PS[, 2] + Error(PS[, 1]:PS[, 3])) :</span></font><font style="font-size:10pt;" size="2"><span style="color:rgb(197,6,11);">
</span></font><font style="font-size:10pt;" size="2"><span style="color:rgb(197,6,11);"> modelo Error() é singular
</span></font><font style="font-size:10pt;" size="2">
Error: PS[, 1]:PS[, 3]
Df Sum Sq Mean Sq F value Pr(>F)
PS[, 1] 1 8.1 8.100 11.3 0.0152 *
Residuals 6 4.3 0.717
---
Signif. codes:
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
PS[, 2] 4 108.9 27.23 24.94 3.02e-08 ***
PS[, 1]:PS[, 2] 4 496.9 124.23 113.79 3.09e-15 ***
Residuals 24 26.2 1.09
---
Signif. codes:
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1</font></pre>#Resíduos:<br><pre style="font-family:'Lucida Console';font-size:16pt !important;border:none;word-break:break-all;white-space:pre-wrap !important;line-height:25px;color:rgb(0,0,0);font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;text-align:-webkit-left;text-indent:0px;text-transform:none;word-spacing:0px;background-color:rgb(255,255,255);"><font style="font-size:10pt;" size="2"><span style="color:blue;white-space:pre-wrap;">> </span></font><font style="font-size:10pt;" size="2"><span style="color:blue;">residuals(mod)
</span></font><font style="font-size:10pt;" size="2">NULL<br>#Teste de Shapiro-Wilk:<br></font><font style="font-size:10pt;" size="2"><span style="color:blue;white-space:pre-wrap;">> </span></font><font style="font-size:10pt;" size="2"><span style="color:blue;">res<-residuals(mod)
</span></font><font style="font-size:10pt;" size="2"><span style="color:blue;white-space:pre-wrap;">> </span></font><font style="font-size:10pt;" size="2"><span style="color:blue;">shapiro.test(res)
</span></font><font style="font-size:10pt;" size="2"><span style="color:rgb(197,6,11);">Error: is.numeric(x) is not TRUE</span></font><br></pre><br>Obrigada.<br>Att.,<br>Angélica.<br><br> </div></div>
<br>______________________________<wbr>_________________<br>
R-br mailing list<br>
</div></div><a rel="nofollow" target="_blank">[hidden email]</a><span><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" target="_blank">https://listas.inf.ufpr.br/cgi<wbr>-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" target="_blank">http://www.leg.ufpr.br/r-br-g<wbr>uia</a>) e forneça código mínimo reproduzível.<br></span></blockquote></div><br></div>
<br>______________________________<wbr>_________________
<br>R-br mailing list
<br><a rel="nofollow" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forne�a c�digo m�nimo reproduz�vel.
<br>
<br>
<hr color="#cccccc" noshade="" size="1">
<div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif;">
<div style="font-weight:bold;">If you reply to this email, your message will be added to the discussion below:</div>
<a href="http://r-br.2285057.n4.nabble.com/R-br-Parcela-Subdividida-tp4666613p4666617.html" target="_blank">http://r-br.2285057.n4.nabble.<wbr>com/R-br-Parcela-Subdividida-<wbr>tp4666613p4666617.html</a>
</div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;line-height:1.5em;">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=YW5nZWxpY2EucmljYXJ0ZUBob3RtYWlsLmNvbXwzMzU3OTgyfDEwODg4MzAyMzk=" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble:email.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble:email.naml-instant_emails%21nabble:email.naml-send_instant_email%21nabble:email.naml" rel="nofollow" style="font:9px serif;" target="_blank">NAML</a>
</div></div> </div></div> </div></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div> </div></body>
</html>