fast_spike_slab_beta function

Compute marginal posterior estimates for beta-spike-and-slab prior