<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Cesar Rabak e demais colegas, boa tarde.</p>
<p><br>
</p>
<p>Agradeço o retorno e peço desculpas pela demora em responder esta postagem.</p>
<p><br>
</p>
<p>Respondendo os questionamentos:</p>
<p><br>
</p>
<p><b>- </b><span style="font-size:12.8px"><b>Por que então que você constrói (especifica) o seu termo Error como PARC:REP?</b></span></p>
<p>R: o modelo da anova que eu utilizei [mod<-aov(var~PARC*SubP+Error(PARC:REP))] retirei de outras postagens feitas no fórum e também o vi em outras consultas, portanto, não sei o porquê do "Error como PARC:REP<span style="font-size:12.8px"></span>" (e gostaria
de entender), mas entendo que seria o erro referente às parcelas ('Erro a' do resultado da anova). Em todas as minhas consultas o modelo para parcela subdividida é sempre este, mas comigo não funcionou (ou melhor, funcionou, mas apareceu uma mensagem de erro,
como mostrado anteriormente). Entendo que se eu não inserir o termo do 'Error', a análise será gerada considerando um esquema fatorial.</p>
<p><br>
</p>
<p>- <b></b><span style="font-size:12.8px"><b>Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer (uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar em conta na sua ANOVA?</b></span></p>
<p><span style="font-size:12.8px"><span style="font-size: 12pt;">R</span><span style="font-size: 12pt;">: como dito anteriormente, o erro referente à Parcela (Erro a), visto que no resultado da anova devo ter dois Erros (a e b).</span></span></p>
<p><span style="font-size:12.8px"><span style="font-size: 12pt;"><br>
</span></span></p>
<p><span style="font-size:12.8px"><span style="font-size: 12pt;">- <span style="font-size:12.8px">
<b>Outra perguntinha: você <i>entende</i> que ao especificar o termo erro você reduziu o número de pontos para ajuste para dez e esse é o número de fatores PARC vezes o número de parcelas SubP?</b></span></span><span style="font-size: 12pt;"></span></span><br>
</p>
<p>R: não consigo compreender isto.</p>
<p><br>
</p>
<p>Se você, ou alguém aqui do fórum, puder me ajudar, agradeço muito, estou ansiosa para resolver este problema. Se eu estiver errada em algum ponto, por favor, corrijam-me. Se o problema está no modelo, qual seria o que eu poderia utilizar para parcela subdividida
(PS) (tenho um experimento em DIC, com parcelas subdividas (2 x 5) e 4 repetições)?
<br>
</p>
<p>Em suma, meu objetivo é gerar o resíduo (que não sei se utilizei o comando correto, visto que é uma PS) para poder realizar o teste de Shapiro-Wilk (SW), pois o quadro da anova que utilizo é aquele do pacote 'ExpDes.pt' (e, infelizmente, para o caso de PS
este pacote não apresenta o resultado do teste de SW). Assim, se houver uma solução ainda mais simples para que possa atingir este objetivo, peço, por gentileza, que me apresentem.</p>
<p><br>
</p>
<p>Muito obrigada.</p>
<p>Aguardo retorno.</p>
<p><br>
</p>
<p>Att.,</p>
<p>Angélica.<br>
</p>
<div style="color: rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>De:</b> Cesar Rabak <cesar.rabak@gmail.com><br>
<b>Enviado:</b> sexta-feira, 9 de setembro de 2016 16:14<br>
<b>Para:</b> Angélica Ricarte<br>
<b>Assunto:</b> Re: [R-br] Parcela Subdividida</font>
<div> </div>
</div>
<div>
<div dir="ltr">Olá Angélica,
<div><br>
</div>
<div>CMR (que está descrito no FAQ deste grupo) é "código mínimo reproduzível".</div>
<div><br>
</div>
<div>A ideia é permitir que a gente refaça seus passos e acelere a ajuda.</div>
<div><br>
</div>
<div>Seu código acima não é exatamente reproduzível porque esta linha:</div>
<div><br>
</div>
<div>> <span style="font-size:12.8px">PS<-read.xlsx(file.choose(),</span><span style="font-size:12.8px"><wbr>sheetName="ParcSub_DIC",h=T) </span></div>
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">diz respeito a um arquivo que está na sua máquina. Contudo, como você imprimiu o total do dataframe importado co o comando seguinte, está quase!</span></div>
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">OK, então você tem uma resposta "var" para o fator PARC e aninhado em PARC você tem quatro subparcelas, com quatro repetições por medida.</span></div>
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">Por que então que você constrói (especifica) o seu termo Error como PARC:REP?</span></div>
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer (uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar em conta na sua ANOVA?</span></div>
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">Outra perguntinha: você <i>entende</i> que ao especificar o termo erro você reduziu o número de pontos para ajuste para dez e esse é o número de fatores PARC vezes o número de parcelas SubP?</span></div>
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">HTH</span></div>
<div><span style="font-size:12.8px">--</span></div>
<div><span style="font-size:12.8px">Cesar Rabak</span></div>
<div><span style="font-size:12.8px"><br>
</span></div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2016-09-08 14:27 GMT-03:00 Angélica Ricarte <span dir="ltr">
<<a href="mailto:angelica.ricarte@hotmail.com" target="_blank">angelica.ricarte@hotmail.com</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<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>
______________________________<wbr>______________________________<wbr>_<br>
> PS<-read.xlsx(file.choose(),<wbr>sheetName="ParcSub_DIC",h=T);<wbr>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.<wbr>factor(SubP)<br>
[1] TRUE<br>
<br>
> mod<-aov(var~PARC*SubP+Error(<wbr>PARC:REP));summary(mod)<br>
<br>
Error: PARC:REP<span class=""><br>
Df Sum Sq Mean Sq F value Pr(>F) <br>
</span>PARC 1 6.453 6.453 80.67 0.0706 .<br>
Residuals 1 0.080 0.080 <br>
<span class="">---<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>
</span>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>
<span class="">---<br>
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<br>
<br>
</span>> resíduo<-residuals(mod);<wbr>resíduo<br>
NULL<br>
______________________________<wbr>______________________________<wbr>_______<br>
Aguardo retorno.<br>
<br>
Obrigada.<br>
<br>
Att.,<br>
Angélica.<br>
<br>
<div>
<hr>
From: <a href="mailto:cesar.rabak@gmail.com" target="_blank">cesar.rabak@gmail.com</a><br>
Date: Tue, 6 Sep 2016 21:49:25 -0300<span class=""><br>
Subject: Re: [R-br] Parcela Subdividida<br>
</span>To: <a href="mailto:angelica.ricarte@hotmail.com" target="_blank">angelica.ricarte@hotmail.com</a>;
<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>
<div>
<div class="h5"><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><br>
<div>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 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@n4<wbr>.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>2016-09-04 1:22 GMT-03:00 Angélica Ricarte via R-br <span dir="ltr"><<a rel="nofollow">[hidden email]</a>></span>:<br>
</div>
</div>
<blockquote style="border-left:2px solid #cccccc; padding:0 1em">
<div>
<div>
<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-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-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">[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">[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/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>
<br>
<hr size="1" noshade="" color="#cccccc">
<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-t<wbr>p4666613p4666617.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" target="_blank" style="font:9px serif">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" target="_blank">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/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="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-g<wbr>uia</a>) e forneça código mínimo reproduzível.<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</body>
</html>