We don't actually use Gaussian models for any of our Bayesian tests =). For continuous data, we use log-normal models, and the next blog post will be about how to make all that work.
And for intervals: yup, we just use sampling. But it's super fast because all you need are independent samples. For more details, see this first post:
It does a huge amount better.