Causal discovery algorithms aim at untangling complex causal relationships from data. Here, we study causal discovery and inference methods based on staged tree models, which can represent complex and asymmetric causal relationships between categorical variables. We provide a first graphical representation of the equivalence class of a staged tree, by looking only at a specific subset of its underlying independences. We further define a new pre-metric, inspired by the widely used structural intervention distance, to quantify the closeness between two staged trees in terms of their corresponding causal inference statements. A simulation study highlights the efficacy of staged trees in uncovering complexes, asymmetric causal relationships from data, and real-world data applications illustrate their use in practical causal analysis.

Bayesian networks are a class of models that are widely used for the diagnosis, prediction, and risk assessment of complex operational systems. Multiple approaches, as well as implemented software, now guide their construction via learning from data or expert elicitation. However, current software only includes minimal functionalities to explore the assumptions, quality of fit, and sensitivity to learned parameters of a constructed Bayesian network. Here, we illustrate the usage of the bnmonitor R package: the first comprehensive software for model-checking of a Bayesian network. An applied data analysis using bnmonitor is carried out over a medical dataset to illustrate the use of its wide array of functions.

Sensitivity analysis in probabilistic discrete graphical models is usually conducted by varying one probability at a time and observing how this affects output probabilities of interest. When one probability is varied, then others are proportionally covaried to respect the sum-to-one condition of probabilities. The choice of proportional covariation is justified by multiple optimality conditions, under which the original and the varied distributions are as close as possible under different measures. For variations of more than one parameter at a time and for the large class of discrete statistical models entertaining a regular monomial parametrisation, we demonstrate the optimality of newly defined proportional multi-way schemes with respect to an optimality criterion based on the I-divergence. We demonstrate that there are varying parameters' choices for which proportional covariation is not optimal and identify the sub-family of distributions where the distance between the original distribution and the one where probabilities are covaried proportionally is minimum. This is shown by adopting a new geometric characterization of sensitivity analysis in monomial models, which include most probabilistic graphical models. We also demonstrate the optimality of proportional covariation for multi-way analyses in Naive Bayes classifiers.

The accuracy of probability distributions inferred using machine-learning algorithms heavily depends on data availability and quality. In practical applications it is therefore fundamental to investigate the robustness of a statistical model to misspecification of some of its underlying probabilities. In the context of graphical models, investigations of robustness fall under the notion of sensitivity analyses. These analyses consist in varying some of the model's probabilities or parameters and then assessing how far apart the original and the varied distributions are. However, for Gaussian graphical models, such variations usually make the original graph an incoherent representation of the model's conditional independence structure. Here we develop an approach to sensitivity analysis which guarantees the original graph remains valid after any probability variation and we quantify the effect of such variations using different measures. To achieve this we take advantage of algebraic techniques to both concisely represent conditional independence and to provide a straightforward way of checking the validity of such relationships. Our methods are demonstrated to be robust and comparable to standard ones, which can break the conditional independence structure of the model, using an artificial example and a medical real-world application.