Bayesian Beta Regression, adapted for bounded discrete responses, commonly seen in survey responses.
Estimation is done via Markov Chain Monte Carlo sampling, using a Gibbs wrapper around univariate slice sampler
(Neal (2003) ), as implemented in the R package MfUSampler
(Mahani and Sharabiani (2017) ).