2da-clase-practica.R

# Segunda Clase Práctica: Estimacion de MRLS

# Librerías
library(readxl)
library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(e1071)

## Warning: package 'e1071' was built under R version 4.4.1

# Importamos la base de datos
EAWE21 <- read_excel("EAWE21.xlsx")

# Hacemos MRLS entre ingresos y nivel de escolaridad
modelo = lm(EAWE21$EARNINGS ~ EAWE21$S)

# Graficamos la relacion y la estimación
plot(EAWE21$S,EAWE21$EARNINGS,
     main = "Relación entre Ingresos y Nivel de escolaridad",
     xlab = "Nivel de escolaridad",
     ylab = "Ingresos")
abline(coef = coef(modelo), col = "blue", lwd = 2)

# Test de hipótesis para la pendiente
resumen_modelo = summary(modelo)
p_value = resumen_modelo$coefficients["EAWE21$S", "Pr(>|t|)"]

cat("El p-valor para la estimación de la pendiente muestra significancia al", p_value ,"% con lo que se puede rechazar la hipótesis nula de que beta es = 0. \n")

## El p-valor para la estimación de la pendiente muestra significancia al 2.579002e-11 % con lo que se puede rechazar la hipótesis nula de que beta es = 0.

# Interpretación de los resultados del modelo
b_coef = coef(modelo)
b1 = b_coef[2]

Y_20 = b_coef[1] + 25*b_coef[2]

cat("Según la estimación realizada, por cada año adicional de estudios es esperable que el ingreso aumente en", b1, "dólares por hora. De esta manera, se esperaría que alguien con 25 años de estudios cobre alrededor de", Y_20, "dólares por hora\n")

## Según la estimación realizada, por cada año adicional de estudios es esperable que el ingreso aumente en 1.265712 dólares por hora. De esta manera, se esperaría que alguien con 25 años de estudios cobre alrededor de 32.40749 dólares por hora