
      library(AER)
      library(MASS)
      data("Boston")
      attach(Boston)
      bh_mod <- lm(medv ~ lstat, data = Boston)
      R2_res <- summary(bh_mod)$r.squared


      # realizar la regresión
      
      # obtener un resumen robusto de coeficientes
       
      # comparar ambos coeficientes de determinación


      # realizar la regresión
      mod <- lm(medv ~ lstat + crim + age, data = Boston)

      # obtener un resumen robusto de coeficientes
      coeftest(mod, vcov. = vcovHC)

      # comparar ambos coeficientes de determinación
      R2_unres <- summary(mod)$r.squared
      R2_unres < R2_res


      test_or({
        ex() %>% check_function("lm") %>% check_result()
      },{
        ex() %>% override_solution("lm(Boston$medv ~ Boston$lstat + Boston$crim + Boston$age)") %>% check_function("lm") %>% check_result()
      },{
        ex() %>% override_solution("lm(medv ~ lstat + crim + age)") %>% check_function("lm") %>% check_result()
      })

      test_or({
        test_function("coeftest")
      },{
        f <- ex() %>% override_solution("coeftest(lm(Boston$medv ~ Boston$lstat + Boston$crim + Boston$age), vcov. = vcovHC)") %>% check_function("coeftest") 
        f %>% check_arg("x") %>% check_equal()
      },{   
        f <- ex() %>% override_solution("coeftest(lm(medv ~ lstat + crim + age), vcov. = vcovHC)") %>% check_function("coeftest") 
        f %>% check_arg("x") %>% check_equal()
      })

      test_student_typed("vcov. = vcovHC")

      test_predefined_objects("R2_res")
      test_object("R2_unres")

      test_or({
        test_student_typed("R2_unres > R2_res")
      },{
        test_student_typed("R2_unres < R2_res")
      },{
        test_student_typed("R2_res < R2_unres")
      },{
        test_student_typed("R2_res > R2_unres")
      })

      success_msg("Correcto. Recuerde que R2 generalmente crece cuando se agrega un regresor al modelo, pero esto no significa necesariamente que el modelo mejore sobre un modelo con menos regresores en todos los casos.")