Bayesian analysis using Monte Carlo integration is a powerful method for univariate inference. This approach makes possible multiparameter flexibility within families of univariate distributions.

When two binary responses are measured for each study subject across time, it may be of interest to model how the bivariate associations and marginal univariate risks involving the two responses ...