vcov_fixcoef function

Gets the variance-covariance matrix of fixed effects of a fitted model