I'd always just use Stata for this when possible since it's so simple. If you do have to use R though, you'd need to specify that the variance-covariance matrix be HC1 (which is one used by Stata). If mod is your model, use coeftest(mod, vcov = vcovHC(mod, type="HC1")). There are also...