Computes "M" and "Sigma" matrices for the sandwich estimator of variance-covariance matrix.
Useful links