Addresses #646 (closed). There is some duplication of code for computing the log likelihood from the various SNR arrays. It might be good to fold this into a likelihood_from_snr_arrays function.
likelihood_from_snr_arrays