<div dir="ltr">OK, Angélica.<div><br></div><div>Veja que seu plot indica uma varição não linear com relação à sobrevivência, ela é máxima para um certo valor de concentração do óleo, enquanto sua curva e uma regressão linear faz entender que o aumento do óleo diminui monotonicamente a taxa de sobrevivência.</div><div><br></div><div>SE for feita uma regressão linear, obtem-se com seus dados o seguinte:</div><div><br></div><div><div><div><font face="monospace, monospace">> summary(lm(S ~ Oleo))</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Call:</font></div><div><font face="monospace, monospace">lm(formula = S ~ Oleo)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Residuals:</font></div><div><font face="monospace, monospace"> Min 1Q Median 3Q Max </font></div><div><font face="monospace, monospace">-10,003 -5,716 -2,858 2,858 15,719 </font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Coefficients:</font></div><div><font face="monospace, monospace"> Estimate Std. Error t value Pr(>|t|) </font></div><div><font face="monospace, monospace">(Intercept) 89,997 3,781 23,804 4,68e-15 ***</font></div><div><font face="monospace, monospace">Oleo -5,716 4,042 -1,414 0,174 </font></div><div><font face="monospace, monospace">---</font></div><div><font face="monospace, monospace">Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Residual standard error: 10,1 on 18 degrees of freedom</font></div><div><font face="monospace, monospace">Multiple R-squared: 0,1,<span style="white-space:pre"> </span>Adjusted R-squared: 0,05 </font></div><div><font face="monospace, monospace">F-statistic: 2 on 1 and 18 DF, p-value: 0,1744</font></div></div><div><br></div><div><font face="monospace, monospace">> anova(lm(S ~ Oleo))</font></div><div><font face="monospace, monospace">Analysis of Variance Table</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Response: S</font></div><div><font face="monospace, monospace"> Df Sum Sq Mean Sq F value Pr(>F)</font></div><div><font face="monospace, monospace">Oleo 1 204,2 204,2 2 0,1744</font></div><div><font face="monospace, monospace">Residuals 18 1837,8 102,1 </font></div></div><div><br></div><div><br></div><div>Quando você chama dic, ele força a variável tratamento (mesmo neste caso sendo uma variável numérica contínua) a ser vista pelo SW que chama aov() em dic, como uma variável categórica não ordenada, e como dentro de cada categoria não há variabilidade nenhuma no seus dados você acaba com isto:</div><div><br></div><div><div><font face="monospace, monospace">> summary(aov(S ~ factor(Oleo)))</font></div><div><font face="monospace, monospace"> Df Sum Sq Mean Sq F value Pr(>F) </font></div><div><font face="monospace, monospace">factor(Oleo) 3 2042 680,7 4,407e+30 <2e-16 ***</font></div><div><font face="monospace, monospace">Residuals 16 0 0,0 </font></div><div><font face="monospace, monospace">---</font></div><div><font face="monospace, monospace">Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1</font></div></div><div><br></div><div>O código acima redunda o óbvio que se saberia lendo os seus dados, para cada valor de conc. de óleo há um valor de sobrevivência!</div><div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-06-04 21:16 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 dir="ltr">
<div id="m_4323487977333733041divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p style="margin-top:0;margin-bottom:0">Olá, Cesar Rabak. </p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">O experimento foi realizado em DIC com adição de doses de extrato de óleo essencial de capim santo em ração de peixes. Assim, o tratamento é "Oleo" e a variável é a taxa de sobrevivência (S).</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Os dados plotados:</p>
<p style="margin-top:0;margin-bottom:0"><img size="12523" id="m_4323487977333733041img112010" style="max-width:99.9%" src="cid:a1cd90f0-63e3-48e9-b9db-0aabefec4e46"><br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Eu queria entender porque não dá certo a análise de regressão e o teste de médias utilizando o pacote ExpDes.pt, sendo que para outras variáveis deu tudo certo. O que me levou a acreditar que o problema está no fato de
não haver variância e o CV=0%. Mas não sei explicar o porquê.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<div id="m_4323487977333733041Signature">
<div id="m_4323487977333733041divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p style="text-align:center"><b></b></p>
<span id="m_4323487977333733041ms-rterangepaste-start"></span>
<p style="text-align:left;font-family:Calibri,Arial,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">
<b>Angélica Ricarte da Silva Batista</b></p>
<p style="font-family:Calibri,Arial,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px;text-align:center">
</p>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">
<div style="text-align:left"><span style="font-size:11pt">Mestre em Agronomia/ Graduada em Engenharia Agronômica (Universidade Federal do Vale do São Francisco - UNIVASF)</span></div>
<span style="font-size:11pt"></span></div>
<span style="font-size:11pt"></span>
<p style="font-family:Calibri,Arial,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px;text-align:center">
<span style="font-size:11pt"></span></p>
<span style="font-size:11pt"></span>
<p style="text-align:left;font-family:Calibri,Arial,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">
<span style="font-size:11pt"><a href="http://lattes.cnpq.br/0601801978475226" class="m_4323487977333733041OWAAutoLink" id="m_4323487977333733041LPNoLP" target="_blank"></a><a href="http://lattes.cnpq.br/0601801978475226" class="m_4323487977333733041OWAAutoLink" id="m_4323487977333733041LPNoLP" target="_blank">Currículo Lattes</a></span><a href="http://lattes.cnpq.br/0601801978475226" id="m_4323487977333733041LPNoLP" target="_blank"><span style="font-size:11pt"></span></a></p>
<span id="m_4323487977333733041ms-rterangepaste-end"></span>
<p style="text-align:center"><a href="http://lattes.cnpq.br/0601801978475226" id="m_4323487977333733041LPNoLP" target="_blank"></a></p>
</div>
</div>
<br>
<br>
<div style="color:rgb(0,0,0)">
<hr style="display:inline-block;width:98%">
<div id="m_4323487977333733041divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>De:</b> Cesar Rabak <<a href="mailto:cesar.rabak@gmail.com" target="_blank">cesar.rabak@gmail.com</a>><br>
<b>Enviado:</b> segunda-feira, 4 de junho de 2018 22:09<br>
<b>Para:</b> Angélica Ricarte; a lista Brasileira oficial de discussão do programa R.<br>
<b>Assunto:</b> Re: [R-br] Erro na análise de regressão</font>
<div> </div>
</div><div><div class="h5">
<div>
<div dir="ltr">Angélica,
<div><br>
</div>
<div>O tratamento seria a variável <i>Oleo</i> e a resposta, a var. <i>S</i> ?</div>
<div><br>
</div>
<div>Você já experimentou plotar esses dados?</div>
<div><br>
</div>
<div>Qual é o seu pressuposto para chamar dic do pacote ExpDes.pt, dessa maneira?</div>
<div><br>
</div>
<div>--</div>
<div>Cesar Rabak</div>
<div><br>
</div>
</div>
<div class="m_4323487977333733041x_gmail_extra"><br>
<div class="m_4323487977333733041x_gmail_quote">2018-06-04 20:15 GMT-03:00 Angélica Ricarte via R-br <span dir="ltr">
<<a href="mailto:r-br@listas.c3sl.ufpr.br" id="m_4323487977333733041LPlnk303370" class="m_4323487977333733041OWAAutoLink" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span>:<br>
<blockquote class="m_4323487977333733041x_gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_4323487977333733041x_m_2229806041982709220divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top:0;margin-bottom:0"><span style="font-size:11pt">Prezados, boa noite.</span><br>
</p>
<div style="color:rgb(0,0,0)">
<div class="m_4323487977333733041x_m_2229806041982709220BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="m_4323487977333733041x_m_2229806041982709220PlainText"><br>
Estou com a seguinte dúvida, se puderem me ajudar, agradeço demais:<br>
<br>
Utilizando o pacote "ExpDes.pt", fiz a seguinte análise:<br>
>dados<br>
Oleo S<br>
1 0.0 85.71<br>
2 0.0 85.71<br>
3 0.0 85.71<br>
4 0.0 85.71<br>
5 0.0 85.71<br>
6 0.5 85.71<br>
7 0.5 85.71<br>
8 0.5 85.71<br>
9 0.5 85.71<br>
10 0.5 85.71<br>
11 1.0 100.00<br>
12 1.0 100.00<br>
13 1.0 100.00<br>
14 1.0 100.00<br>
15 1.0 100.00<br>
16 1.5 71.42<br>
17 1.5 71.42<br>
18 1.5 71.42<br>
19 1.5 71.42<br>
20 1.5 71.42<br>
> dic(dados[,1],dados[,2],quali=<wbr>F,sigF=0.05)<br>
------------------------------<wbr>------------------------------<wbr>------------<br>
Quadro da analise de variancia<br>
------------------------------<wbr>------------------------------<wbr>------------<br>
GL SQ QM Fc Pr>Fc<br>
Tratamento 3 2042 680.68 4.4071e+30 1.5357e-239<br>
Residuo 16 0 0.00 <br>
Total 19 2042 <wbr> <br>
------------------------------<wbr>------------------------------<wbr>------------<br>
CV = 0 %<br>
<br>
------------------------------<wbr>------------------------------<wbr>------------<br>
Teste de normalidade dos residuos <br>
Valor-p: 5.640644e-07 <br>
ATENCAO: a 5% de significancia, os residuos nao podem ser considerados normais!<br>
------------------------------<wbr>------------------------------<wbr>------------<br>
<br>
------------------------------<wbr>------------------------------<wbr>------------<br>
Teste de homogeneidade de variancia <br>
valor-p: NaN <br>
Error in if (pvalor.hvar < 0.05) { : <br>
valor ausente onde TRUE/FALSE necessário<br>
> dic(dados[,1],dados[,2],quali=<wbr>F,sigF=0.05)<br>
<br>
Porque não consigo obter um resultado para a análise de regressão? Seria porque o CV é igual a zero e porque não há variância? Se sim, poderiam me dar uma explicação sobre a relação entre a variância e a regressão? Também tentei fazer o mesmo considerando os
tratamentos qualitativos (quali=T) para obter um teste de médias, porém aparece a mesma mensagem de erro. OBS: para outras variáveis (que possuem variância e CV diferente de zero) deu certo a análise.<br>
<br>
Grata, aguardo retorno.<br>
Att.,<br>
Angélica Batista.<br>
<br>
</div>
</span></font></div>
</div>
</div>
</div>
<br>
______________________________<wbr>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" id="m_4323487977333733041LPlnk952576" class="m_4323487977333733041OWAAutoLink" 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" id="m_4323487977333733041LPlnk939596" class="m_4323487977333733041OWAAutoLink" 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" id="m_4323487977333733041LPlnk955088" class="m_4323487977333733041OWAAutoLink" 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>
</div>
</blockquote></div><br></div>