Restrictions of PCBNs for integration-free computations
Alexis Derumigny, Niels Horsman, Dorota Kurowicka
TL;DR
This paper addresses the computational bottleneck in PCBNs where density evaluation and simulation may require high-dimensional integration. It provides a complete graph-level characterization: integration-free computations are possible if and only if the DAG contains no active cycles or interfering v-structures, and offers an algorithm to construct parental orders that guarantee integration-free costs. The authors also establish the asymptotic normality of stepwise estimators for PCBNs under these restricted conditions and demonstrate favorable finite-sample behavior in simulations. The proposed approach enables scalable, flexible modeling with copula-based dependencies while controlling computational complexity, with estimation implemented via a stepwise rank-based estimating-equations framework. Overall, the work yields practical guidelines and theoretical guarantees for integrating PCBNs in high-dimensional settings.
Abstract
The pair-copula Bayesian Networks (PCBN) are graphical models composed of a directed acyclic graph (DAG) that represents (conditional) independence in a joint distribution. The nodes of the DAG are associated with marginal densities, and arcs are assigned with bivariate (conditional) copulas following a prescribed collection of parental orders. The choice of marginal densities and copulas is unconstrained. However, the simulation and inference of a PCBN model may necessitate possibly high-dimensional integration. We present the full characterization of DAGs that do not require any integration for density evaluation or simulations. Furthermore, we propose an algorithm that can find all possible parental orders that do not lead to (expensive) integration. Finally, we show the asymptotic normality of estimators of PCBN models using stepwise estimating equations. Such estimators can be computed effectively if the PCBN does not require integration. A simulation study shows the good finite-sample properties of our estimators.
