These R codes implement the Bayesian methodology of Castelletti & Consonni (2021, Bayesian Analysis) for structure learning and causal inference in probit graphical models. Specifically: ...