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! ")