🤖 AI Summary
Density evaluation and random sampling in Pair-Copula Bayesian Networks (PCBNs) typically require high-dimensional numerical integration, leading to prohibitive computational costs.
Method: We first derive necessary and sufficient conditions on the DAG structure under which exact density evaluation and simulation can be performed *without any numerical integration*. We then propose a polynomial-time algorithm to automatically identify “integration-free parent orderings” and develop a parameter estimation framework based on stepwise estimating equations.
Contribution/Results: We rigorously prove the asymptotic normality of the resulting estimators. By leveraging Pair-Copula decomposition, DAG topological constraints, and conditional independence analysis, our approach preserves statistical accuracy while drastically improving computational efficiency. Monte Carlo simulations demonstrate favorable finite-sample properties, significantly enhancing the practicality and scalability of PCBNs.
📝 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.