get_var_EB function

Using simulation to obtain the asymptotic variance-covariance matrix of gamma_EB, package corpcor and MASS are required