[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