ec_estimates_SIGMA function

Obtain posterior mean estimates of upper level covariance