Support fitting models using Bayes

Uncertainty quantification can be a challenge in error correction models using MLE.
We currently use MLE with the quadratic approximation. #4 partially addresses this, but maybe Stan can do an even better job and also support a wider range of models more easily.

One thing to consider is if this package should focus on MLE methods and supporting Bayesian models should be developed in parallel, but in a separate package.