Computes credible intervals for Bayesian wavelet shrinkage