mirror_vcov function

Create a properly symmetric variance covariance matrix