Stochastic Variational Inference for Vine Copula Distributional Regression
Abstract
Structured additive distributional regression flexibly relates all parameters of a conditional response distribution to covariates, but multivariate extensions remain challenging when dependence is complex. We propose a multivariate structured additive distributional regression model based on regular vine copulas. Different vine edges may use different pair-copula families, while every pair-copula parameter may vary with covariates through a structured additive predictor. The model therefore accommodates heterogeneous marginal distributions together with pair-specific asymmetric, tail-dependent, and covariate-dependent dependence structures. For scalable inference, we develop a tree-wise stochastic variational inference procedure based on component-specific Gaussian variational approximations. Marginal models are estimated first, followed by pair-copula regressions sequentially along the vine trees. We also adapt sequential vine selection by fitting candidate covariate-dependent pair-copula regressions and using information criteria based on effective degrees of freedom for both family and tree selection. In simulations, the tree-wise estimator remains close to an oracle using the true recursive conditional inputs, whereas global refinement yields more concentrated approximations and lower frequentist coverage for upstream components. An application to six-dimensional meteorological data from the Netherlands yields a vine combining Gaussian and non-Gaussian pair copulas, with pronounced nonlinear spatial and temporal variation in dependence.
Related
- Joint Distributional Copula Regression for Insurance Claim Frequency and Severity via Boosting
- Variational Bayesian inference for logistic regression models for group testing data: Comparing a new hybrid variational inference approach to traditional methods.
- Adaptive mixture variational inference for spike-and-slab regression
- Amortized Structured Stochastic Variational Inference for Gaussian Process Latent Variable Models
- Quantile regression and copula-based analysis of risk and dependence in the South African stock market
BibTeX
@misc{callegher2026stochastic,
title = {{Stochastic Variational Inference for Vine Copula Distributional Regression}},
author = {Gianmarco Callegher and Thomas Kneib},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.37996},
eprint = {2609.37996},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.37996},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Stochastic Variational Inference for Vine Copula Distributional Regression
AU - Gianmarco Callegher
AU - Thomas Kneib
PY - 2026
DA - 2026/09/29
PB - arXiv
DO - 10.48550/arxiv.2609.37996
UR - https://arxiv.org/abs/2609.37996
AB - Structured additive distributional regression flexibly relates all parameters of a conditional response distribution to covariates, but multivariate extensions remain challenging when dependence is complex. We propose a multivariate structured additive distributional regression model based on regular vine copulas. Different vine edges may use different pair-copula families, while every pair-copula parameter may vary with covariates through a structured additive predictor. The model therefore accommodates heterogeneous marginal distributions together with pair-specific asymmetric, tail-dependent, and covariate-dependent dependence structures. For scalable inference, we develop a tree-wise stochastic variational inference procedure based on component-specific Gaussian variational approximations. Marginal models are estimated first, followed by pair-copula regressions sequentially along the vine trees. We also adapt sequential vine selection by fitting candidate covariate-dependent pair-copula regressions and using information criteria based on effective degrees of freedom for both family and tree selection. In simulations, the tree-wise estimator remains close to an oracle using the true recursive conditional inputs, whereas global refinement yields more concentrated approximations and lower frequentist coverage for upstream components. An application to six-dimensional meteorological data from the Netherlands yields a vine combining Gaussian and non-Gaussian pair copulas, with pronounced nonlinear spatial and temporal variation in dependence.
ER -