vcov-method function

Extract parameter variance covariance matrix