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

classic Classic list List threaded Threaded
4 messages Options
Reply | Threaded
Open this post in threaded view
|

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

R-br mailing list

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
[hidden email]
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.
Reply | Threaded
Open this post in threaded view
|

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

R-br mailing list
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 <[hidden email]>:

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
[hidden email]
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.


_______________________________________________
R-br mailing list
[hidden email]
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.
Reply | Threaded
Open this post in threaded view
|

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

R-br mailing list

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)

Currículo Lattes




De: Cesar Rabak <[hidden email]>
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 <[hidden email]>:

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
[hidden email]
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.


_______________________________________________
R-br mailing list
[hidden email]
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.
Reply | Threaded
Open this post in threaded view
|

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

R-br mailing list
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 <[hidden email]>:

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)

Currículo Lattes




De: Cesar Rabak <[hidden email]>
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 <[hidden email]>:

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
[hidden email]
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.



_______________________________________________
R-br mailing list
[hidden email]
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.