# NOT RUN {
library(stargazer)
lm1 <- lm(mpg ~ hp, data = mtcars)
lm2 <- lm(mpg ~ hp + wt, data = mtcars)
# Use default "HC2" standard errors
stargazer(lm1, lm2,
se = starprep(lm1, lm2),
p = starprep(lm1, lm2, stat = "p.value"),
omit.stat = "f")
# NB: We remove the F-stat because stargazer only can use original F-stat
# which uses classical SEs
# Use default "CR2" standard errors with clusters
stargazer(lm1, lm2,
se = starprep(lm1, lm2, clusters = mtcars$carb),
p = starprep(lm1, lm2, clusters = mtcars$carb, stat = "p.value"),
omit.stat = "f")
# Can also specify significance levels and different standard errors
stargazer(lm1, lm2,
ci.custom = starprep(lm1, lm2, se_type = "HC3", alpha = 0.1, stat = "ci"),
omit.stat = "f")
# }
Run the code above in your browser using DataLab