library(AER)
data(CPS1985)
dummy_mod <- lm(wage ~ gender, data = CPS1985)
# compute robust standard errors and save them in rob_SEs
# compute robust standard errors and save them in rob_SEs
rob_SEs <- sqrt(diag(vcovHC(dummy_mod, type = "HC1")))
test_predefined_objects("dummy_mod")
test_object("rob_SEs")
success_msg("Nice! ")