library(AER) data(CPS1985) dummy_mod <- lm(wage ~ gender, data = CPS1985) # calcular errores estándar robustos y guardarlos en rob_SEs # calcular errores estándar robustos y guardarlos en rob_SEs rob_SEs <- sqrt(diag(vcovHC(dummy_mod, type = "HC1"))) test_predefined_objects("dummy_mod") test_object("rob_SEs") success_msg("¡Excelente!")