Computational methods
The joint posterior probability distributions of the model parameters described in the preceding sections are analytically intractable, and therefore computational Bayesian methods will be used for the data analyses. We will use Markov-chain Monte Carlo (MCMC) methods implemented in stan, a probabilistic programming language [