vcov_global_loading function

Get the covariance estimates for the global factor loadings