Calculate Covariance Matrix from a linear model fitted with lm() using different estimators
Useful links