sand_vcov function

Compute sandwich estimate of variance-covariance matrix