import matplotlib.pyplot as plt import pymc3 as pm import aesara.tensor as at import arviz as az import scipy.stats as st import numpy as np az.style.use(‘arviz-whitegrid’) Multivariate bernoulli distribution Sources
https://doi.org/10.1016/0047-259X(90)90084-U
https://arxiv.org/pdf/1206.1874
Motivation I have recently been introduced to Bayesian framework for statistical analysis with PyMC3 package for Probabilistic Programming in Python 3.
I came across to a data given to me by my friend, which had multiple columns of binary data (responses to various questions with answers in binary format, yes or no).