implement estimation of correlated error matrix

[stephens999]

this requires

• ability to add a covariance matrix during set_mash_data
• changing of calc_lik_matrix in likelihood.R
• changing of compute_posterior_arrays in posterior.R

[stephens999]

also need to consider changing the ed wrapper to allow for covariances. (maybe a separate issue)

[stephens999]

this is done