A continuación el código para cargar los datos y para mostrar las 10 primeras observaciones de la base de datos llamada Medifis, en total tenemos 27 observaciones.

datos8 <- read.table("C:/Users/juanv/OneDrive/Semestre 5/PYEM/Semana 4/DatosT2.txt", header=TRUE, sep='')
head(datos8, 10)
##    Estatura L.Pie
## 1       159    36
## 2       164    39
## 3       172    38
## 4       167    37
## 5       164    36
## 6       161    38
## 7       168    39
## 8       181    43
## 9       183    41
## 10      158    36

Se crea un diagrama de dispersión que muestre la relación entre las dos variables, se usan las siguientes instrucciones.

library(ggplot2)

ggplot(datos8, aes(x = Estatura, y = L.Pie)) +
  geom_point(color = "blue", size = 3) +
  labs(title = "Diagrama de Dispersión entre Estatura y Longitud", x = "Estatura", y = "Longitud Del Pie") +
  theme_minimal()

mod1 <- lm(L.Pie ~ Estatura, data=datos8)
mod1 # Para imprimir el objeto mod1
## 
## Call:
## lm(formula = L.Pie ~ Estatura, data = datos8)
## 
## Coefficients:
## (Intercept)     Estatura  
##     -4.9638       0.2605

Para graficar la recta de regresión:

plot(datos8$Estatura, datos8$L.Pie, 
     main = "Regresión Lineal: L.Pie vs Estatura", 
     xlab = "Estatura", 
     ylab = "Longitud del Pie",
     col = "blue", pch = 19)

abline(mod1, col = "darkgreen", lwd = 2)  # Agrega la línea de regresión