[R-br] Erro na análise de regressão

Cesar Rabak cesar.rabak em gmail.com
Segunda Junho 4 21:51:22 -03 2018


OK, Angélica.

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.

SE for feita uma regressão linear, obtem-se com seus dados o seguinte:

> summary(lm(S ~ Oleo))

Call:
lm(formula = S ~ Oleo)

Residuals:
    Min      1Q  Median      3Q     Max
-10,003  -5,716  -2,858   2,858  15,719

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept)   89,997      3,781  23,804 4,68e-15 ***
Oleo          -5,716      4,042  -1,414    0,174
---
Signif. codes:  0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1

Residual standard error: 10,1 on 18 degrees of freedom
Multiple R-squared:    0,1, Adjusted R-squared:   0,05
F-statistic:     2 on 1 and 18 DF,  p-value: 0,1744

> anova(lm(S ~ Oleo))
Analysis of Variance Table

Response: S
          Df Sum Sq Mean Sq F value Pr(>F)
Oleo       1  204,2   204,2       2 0,1744
Residuals 18 1837,8   102,1


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:

> summary(aov(S ~ factor(Oleo)))
             Df Sum Sq Mean Sq   F value Pr(>F)
factor(Oleo)  3   2042   680,7 4,407e+30 <2e-16 ***
Residuals    16      0     0,0
---
Signif. codes:  0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1

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!

HTH
--
Cesar Rabak


2018-06-04 21:16 GMT-03:00 Angélica Ricarte <angelica.ricarte em hotmail.com>:

> Olá, Cesar Rabak.
>
>
> 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).
>
>
> Os dados plotados:
>
>
>
> 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ê.
>
>
> *Angélica Ricarte da Silva Batista*
>
> Mestre em Agronomia/ Graduada em Engenharia Agronômica (Universidade
> Federal do Vale do São Francisco - UNIVASF)
>
> <http://lattes.cnpq.br/0601801978475226>Currículo Lattes
> <http://lattes.cnpq.br/0601801978475226>
> <http://lattes.cnpq.br/0601801978475226>
>
> <http://lattes.cnpq.br/0601801978475226>
>
>
> ------------------------------
> *De:* Cesar Rabak <cesar.rabak em gmail.com>
> *Enviado:* segunda-feira, 4 de junho de 2018 22:09
> *Para:* Angélica Ricarte; a lista Brasileira oficial de discussão do
> programa R.
> *Assunto:* Re: [R-br] Erro na análise de regressão
>
> Angélica,
>
> O tratamento seria a variável *Oleo* e a resposta, a var. *S* ?
>
> Você já experimentou plotar esses dados?
>
> Qual é o seu pressuposto para chamar dic do pacote ExpDes.pt, dessa
> maneira?
>
> --
> Cesar Rabak
>
>
> 2018-06-04 20:15 GMT-03:00 Angélica Ricarte via R-br <
> r-br em listas.c3sl.ufpr.br>:
>
> Prezados, boa noite.
>
> Estou com a seguinte dúvida, se puderem me ajudar, agradeço demais:
>
> Utilizando o pacote "ExpDes.pt", fiz a seguinte análise:
> >dados
>    Oleo      S
> 1   0.0  85.71
> 2   0.0  85.71
> 3   0.0  85.71
> 4   0.0  85.71
> 5   0.0  85.71
> 6   0.5  85.71
> 7   0.5  85.71
> 8   0.5  85.71
> 9   0.5  85.71
> 10  0.5  85.71
> 11  1.0 100.00
> 12  1.0 100.00
> 13  1.0 100.00
> 14  1.0 100.00
> 15  1.0 100.00
> 16  1.5  71.42
> 17  1.5  71.42
> 18  1.5  71.42
> 19  1.5  71.42
> 20  1.5  71.42
> > dic(dados[,1],dados[,2],quali=F,sigF=0.05)
> ------------------------------------------------------------------------
> Quadro da analise de variancia
> ------------------------------------------------------------------------
>            GL   SQ     QM         Fc       Pr>Fc
> Tratamento  3 2042 680.68 4.4071e+30 1.5357e-239
> Residuo    16    0   0.00
> Total      19 2042
> ------------------------------------------------------------------------
> CV = 0 %
>
> ------------------------------------------------------------------------
> Teste de normalidade dos residuos
> Valor-p:  5.640644e-07
> ATENCAO: a 5% de significancia, os residuos nao podem ser considerados
> normais!
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------
> Teste de homogeneidade de variancia
> valor-p:  NaN
> Error in if (pvalor.hvar < 0.05) { :
>   valor ausente onde TRUE/FALSE necessário
> > dic(dados[,1],dados[,2],quali=F,sigF=0.05)
>
> 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.
>
> Grata, aguardo retorno.
> Att.,
> Angélica Batista.
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180604/67a775ad/attachment-0001.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: pastedImage.png
Tipo: image/png
Tamanho: 12523 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180604/67a775ad/attachment-0001.png>


Mais detalhes sobre a lista de discussão R-br