Regresión Lineal Simple

(Ejemplo extraido del Libro: Bioestadistica, Autor: Daniels Cap9-ejercicio=22)

Los datos muestran los valores del consumo de metilmercurio y la concentración total de mercurio en la sangre de 12 idividuos expuestos al metilmercurio por consumir peces contaminados.

metilmercurio=c(180,200,230,410,600,550,275,580,105,250,460,650)
totalmercurio=c(90,120,125,290,310,290,170,375,70,105,205,480)
tabla=data.frame(metilmercurio,totalmercurio)
tabla
##    metilmercurio totalmercurio
## 1            180            90
## 2            200           120
## 3            230           125
## 4            410           290
## 5            600           310
## 6            550           290
## 7            275           170
## 8            580           375
## 9            105            70
## 10           250           105
## 11           460           205
## 12           650           480

Estimación de la recta de regresión, coeficientes, R^2 y ANOVA

modelo2=lm(formula=totalmercurio~metilmercurio)
modelo2
## 
## Call:
## lm(formula = totalmercurio ~ metilmercurio)
## 
## Coefficients:
##   (Intercept)  metilmercurio  
##      -20.5791         0.6407
summary(modelo2)
## 
## Call:
## lm(formula = totalmercurio ~ metilmercurio)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -69.164 -36.413   5.319  23.462  84.094 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -20.57906   30.66189  -0.671    0.517    
## metilmercurio   0.64075    0.07373   8.691 5.66e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 46.37 on 10 degrees of freedom
## Multiple R-squared:  0.8831, Adjusted R-squared:  0.8714 
## F-statistic: 75.53 on 1 and 10 DF,  p-value: 5.66e-06
anova(modelo2)
## Analysis of Variance Table
## 
## Response: totalmercurio
##               Df Sum Sq Mean Sq F value   Pr(>F)    
## metilmercurio  1 162392  162392  75.531 5.66e-06 ***
## Residuals     10  21500    2150                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Gráfico de dispersión

Relación entre el consumo de mercurio de metil(ug Hg/día) y Mercurio en toda la sangre(ng/g

plot(metilmercurio,totalmercurio)
abline(modelo2,col="pink")