Como fazer análise de relação peso-comprimento usando o R
Boa noite, Gostaria de saber como fazer a relação peso-comprimento para peixes utilizando o R. Normalmente, faço essa análise utilizando o Statistica e estou tentando atualmente fazer no R. Logaritmizei mesmo dados de peso e comprimento e fiz análise utilizando o seguinte comando (conforme um colega me passou): lwr<-read.table("clipboard",sep="\t", header=T, dec=".") attach(lwr) lwr lpeso<-log(peso) lcomprimento<-log(comprimento) lmlwr<-lm(lpeso~lcomprimento) summary(lmlwr) e obtive o seguinte resultado: Call: lm(formula = lpeso ~ lcomp) Residuals: Min 1Q Median 3Q Max -0.5495 -0.1027 -0.0046 0.1154 0.3636 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -3.73920 0.06033 -61.98 <2e-16 *** lcomp 3.18886 0.03571 89.31 <2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.1874 on 46 degrees of freedom Multiple R-squared: 0.9943, Adjusted R-squared: 0.9941 F-statistic: 7975 on 1 and 46 DF, p-value: < 2.2e-16 No entanto, onde vejo a variação (mínimo - máximo) dos valores de a e b? Acredito que esteja usando o comando errado. Se puderem ajudar, agradeço Muito obrigada desde já
Você quer estimar o índice de correlação entre o peso e comprimento e quer o IC95%, isso? A função lm() que você usou ajusta um modelo linear entre as duas variáveis. Daniel ---------------------------------------------------------------- Daniel Tiezzi, MD, PhD Oncologia / Mastologia Professor Associado - Livre Docente Departamento de Ginecologia e Obstetrícia Setor de Mastologia e Oncologia Ginecológica Faculdade de Medicina de Ribeirão Preto - USP Tel.: 16 3602-2488 https://github.com/dtiezzi <https://github.com/dtiezzi> http://danieltiezzi.pro.br <http://danieltiezzi.pro.br/>e-mail: dtiezzi@usp.br
On 3 Feb 2021, at 23:35, Chiara Lubich por (R-br) <r-br@listas.c3sl.ufpr.br> wrote:
Boa noite,
Gostaria de saber como fazer a relação peso-comprimento para peixes utilizando o R. Normalmente, faço essa análise utilizando o Statistica e estou tentando atualmente fazer no R. Logaritmizei mesmo dados de peso e comprimento e fiz análise utilizando o seguinte comando (conforme um colega me passou): lwr<-read.table("clipboard",sep="\t", header=T, dec=".") attach(lwr) lwr lpeso<-log(peso) lcomprimento<-log(comprimento) lmlwr<-lm(lpeso~lcomprimento) summary(lmlwr)
e obtive o seguinte resultado: Call: lm(formula = lpeso ~ lcomp)
Residuals: Min 1Q Median 3Q Max -0.5495 -0.1027 -0.0046 0.1154 0.3636
Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -3.73920 0.06033 -61.98 <2e-16 *** lcomp 3.18886 0.03571 89.31 <2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.1874 on 46 degrees of freedom Multiple R-squared: 0.9943, Adjusted R-squared: 0.9941 F-statistic: 7975 on 1 and 46 DF, p-value: < 2.2e-16
No entanto, onde vejo a variação (mínimo - máximo) dos valores de a e b? Acredito que esteja usando o comando errado.
Se puderem ajudar, agradeço
Muito obrigada desde já _______________________________________________ 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.
Bom dia, Chiara. Como o Daniel comentou, você ajustou um modelo de regressão linear aos seus dados. Se sua pergunta: "onde vejo a variação (mínimo - máximo) dos valores de a e b?", significa que você quer construir um intervalo de confiança para os parâmetros do seu modelo de regressão, basta usar o comando confint(lmlwr ). Um abraço, Manassés Em qua., 3 de fev. de 2021 às 23:51, Daniel Guimarães Tiezzi por (R-br) < r-br@listas.c3sl.ufpr.br> escreveu:
Você quer estimar o índice de correlação entre o peso e comprimento e quer o IC95%, isso? A função lm() que você usou ajusta um modelo linear entre as duas variáveis.
Daniel
---------------------------------------------------------------- Daniel Tiezzi, MD, PhD Oncologia / Mastologia Professor Associado - Livre Docente Departamento de Ginecologia e Obstetrícia Setor de Mastologia e Oncologia Ginecológica Faculdade de Medicina de Ribeirão Preto - USP Tel.: 16 3602-2488 https://github.com/dtiezzi http://danieltiezzi.pro.br e-mail: dtiezzi@usp.br <dtiezzi@usp.br>
On 3 Feb 2021, at 23:35, Chiara Lubich por (R-br) < r-br@listas.c3sl.ufpr.br> wrote:
Boa noite,
Gostaria de saber como fazer a relação peso-comprimento para peixes utilizando o R. Normalmente, faço essa análise utilizando o Statistica e estou tentando atualmente fazer no R. Logaritmizei mesmo dados de peso e comprimento e fiz análise utilizando o seguinte comando (conforme um colega me passou): lwr<-read.table("clipboard",sep="\t", header=T, dec=".") attach(lwr) lwr lpeso<-log(peso) lcomprimento<-log(comprimento) lmlwr<-lm(lpeso~lcomprimento) summary(lmlwr)
e obtive o seguinte resultado: Call: lm(formula = lpeso ~ lcomp)
Residuals: Min 1Q Median 3Q Max -0.5495 -0.1027 -0.0046 0.1154 0.3636
Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -3.73920 0.06033 -61.98 <2e-16 *** lcomp 3.18886 0.03571 89.31 <2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.1874 on 46 degrees of freedom Multiple R-squared: 0.9943, Adjusted R-squared: 0.9941 F-statistic: 7975 on 1 and 46 DF, p-value: < 2.2e-16
No entanto, onde vejo a variação (mínimo - máximo) dos valores de a e b? Acredito que esteja usando o comando errado.
Se puderem ajudar, agradeço
Muito obrigada desde já _______________________________________________ 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.
_______________________________________________ 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.
Cesar, muito obrigada pelos esclarecimentos. Farei o cálculo do IC (pro modelo não linear) por *bootstraping, *como você indicou. Obrigada!
participantes (3)
-
Chiara Lubich -
Daniel Guimarães Tiezzi -
Manassés P Nóbrega