A robust heteroskedasticity consistent covariance matrix estimator