[R-br] Como fazer análise de relação peso-comprimento usando o R

Manassés P Nóbrega nobrega.manasses em gmail.com
Qui Fev 4 10:26:36 -02 2021


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 em 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 em usp.br <dtiezzi em usp.br>
>
> On 3 Feb 2021, at 23:35, Chiara Lubich por (R-br) <
> r-br em 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 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.
>
>
> _______________________________________________
> 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/20210204/5423314e/attachment.html>


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