
      # adjuntar ambos paquetes y cargar el conjunto de datos

      # obtener una descripción general del conjunto de datos

      # estimar el modelo de regresión simple

      # imprimir el resumen en la consola


      # adjuntar ambos paquetes y cargar el conjunto de datos
      library(AER)
      library(MASS)
      data("Boston")

      # obtener una descripción general del conjunto de datos
      summary(Boston)

      # o
      str(Boston)
      
      # o
      head(Boston)
      
      # estimar el modelo de regresión simple
      bh_mod <- lm(medv ~ lstat, data = Boston)
      
      # imprimir el resumen en la consola
      coeftest(bh_mod, vcov. = vcovHC)


      test_function("library", index = 1)
      test_function("library", index = 2)
      test_function("data")

      test_or({
      test_function("head")
      },{
      test_function("summary")
      },{
      test_function("str")
      })

      test_function("coeftest", args=c("x", "vcov."), eval = F)

      test_or({
        test_object("bh_mod")
      },{
        f <- ex() %>% override_solution("lm(Boston$medv ~ Boston$lstat)") %>% check_function("lm")
        f %>% check_arg("formula") %>% check_equal()
      },{
        f <- ex() %>% override_solution("attach(Boston);lm(medv ~ lstat)") %>% check_function("lm")
        f %>% check_arg("formula") %>% check_equal()
      })

      success_msg("Correcto. Observe que ambos coeficientes son muy significativos.")