
      library(AER)
      data(CPS1985)
      dummy_mod <- lm(wage ~ gender, data = CPS1985)


      # probar si la brecha de género es significativamente diferente de cero
      
      # utilizar errores estándar robustos


      coeftest(dummy_mod, vcov. = vcovHC, type = "HC0")
      
      # o
      linearHypothesis(dummy_mod, "genderfemale=0", vcov. = vcovHC, type = "HC0")


      test_predefined_objects("dummy_mod")
      test_or({
        test_function("coeftest", args = c("x", "vcov.")) 
      },{
        test_function("linearHypothesis", args = c("model", "vcov.", "hypothesis.matrix"))
      },{
        test_function("linearHypothesis", args = c("model", "hypothesis.matrix", "type"))
      })
      success_msg("¡Correcto! La hipótesis de que no existe una brecha salarial se rechaza a cualquier nivel.")