Erro na análise de regressão
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.
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@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@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.
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: [cid:a1cd90f0-63e3-48e9-b9db-0aabefec4e46] 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@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@listas.c3sl.ufpr.br<mailto:r-br@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@listas.c3sl.ufpr.br<mailto:R-br@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.
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@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@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@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@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.
participantes (2)
-
Angélica Ricarte -
Cesar Rabak