R code for a regression table

Here is R code to make a table like this:

I wanted a regression table where:

• Outcome means are reported in the table
• N changes across columns
• Columns have different standard error adjustments, including clusters

Notes:

There are several pitfalls with standard packages.

• This simpler option works in stargazer but only if N does not change across columns:

  se=starprep(mod1, mod2, mod3, clusters = auto$foreign, se_type = "stata")

• lm_robust() does not work with stargazer (StackOverflow discussion)