Formula method for regression with HC4 Heteroscedasticity-consistent standard errors

Source rallfun first.