
      library(foreign)
      library(dplyr)
      data_URL <- "https://github.com/Jeshua-Romero-Guadarrama/Econoalgoritmia/tree/main/data/fastfood.dta?raw=true"
      dat <- read.dta(data_URL)
      dat <- dat %>% 
        mutate(FTE = nmgrs + empft + (0.5 * emppt),
        FTE2 = nmgrs2 + empft2 + (0.5 * emppt2))
      dat_NJ <- subset(dat, state == 1)
      dat_PA <- subset(dat, state == 0)
      reg_dat <- data.frame(
            rbind(
              data.frame(id = dat$sheet, 
                         chain = dat$chain,
                         state = dat$state,
                         empl = dat$FTE,
                         D = 0),
              data.frame(id = dat$sheet,
                         chain = dat$chain,
                         state = dat$state,
                         empl = dat$FTE2,
                         D = 1)))


      # estimar el modelo de regresión
      
      # obtener un resumen sólido


      # estimar el modelo de regresión
      emp_mod <- lm(empl ~ D, data = reg_dat, subset = state == 1)

      # obtener un resumen sólido
      library(AER)
      coeftest(emp_mod, vcov. = vcovHC, type = "HC1")


      ex() %>% check_predefined_objects("reg_dat")
      ex() %>% check_library("AER")
      ex() %>% check_object("emp_mod") %>% check_equal(undefined_msg = "¡Asegúrate de no eliminar `emp_mod`!")
      check_or(
          ex() %>%
          check_function(., "lm") %>% {
            check_arg(., "formula") %>% check_equal(incorrect_msg = "Parece que el argumento `formula` en `lm()` es incorrecto.")
            check_arg(., "data") %>% check_equal(incorrect_msg = "Parece que el argumento `data` en `lm()` es incorrecto.")
          },
          ex() %>%
          override_solution("emp_mod <- lm(dat_NJ$empl ~ dat_NJ$D)") %>% 
          check_function("lm") %>% check_arg("formula") %>% check_equal(incorrect_msg = "Parece que el argumento `data` en `lm()` es incorrecto."))
          check_or(
            ex() %>% check_function("coeftest") %>% {
            check_arg(., "x") %>% check_equal()
            check_arg(., "vcov.") %>% check_equal()
            },
          ex() %>% 
          override_solution("emp_mod <- lm(dat_NJ$empl ~ dat_NJ$D); coeftest(emp_mod, vcov. = vcovHC, type = 'HC1')") %>%
          check_function("coeftest") %>% {
          check_arg(., "x") %>% check_equal()
          check_arg(., "vcov.") %>% check_equal()
            })
           success_msg(msg = "¡Bien hecho! La estimación del coeficiente en "D", así como el estadístico t y, por lo tanto, también el valor p, coinciden con los resultados informados por "t.test()". No se puede rechazar la hipótesis de que el coeficiente de "D" sea cero.")