get_robust_vcov function

Compute Robust Covariance Matrix (HC0-style)