vcov-methods function

Calculate variance-covariance matrix for a model