deep copula modeling

This work addresses the lack of identifiability and interpretability in deep generative models when handling multivariate data with arbitrary marginal distributions. The authors propose an identifiable Bayesian deep generative Copula model based on a hierarchical binary latent variable structure. By employing a rank likelihood to decouple modeling of marginal distributions from dependence structure, the method achieves, for the first time, parameter identifiability in deep generative Copula models. A Bayesian rank selection prior with adaptive layer width is introduced to enhance flexibility. Theoretical analysis establishes posterior consistency under both continuous and mixed marginals, and the inference algorithm combines stochastic EM with maximum a posteriori estimation. Empirical results demonstrate superior performance in small-sample settings and reveal an interpretable hierarchical latent structure in personality survey data.

Existing spatial data models struggle to simultaneously capture complex extremal features such as tail dependence, asymptotic independence, and tail asymmetry. This work proposes a novel approach by introducing multivariate Pareto mixture distributions into a spatial copula framework, yielding a flexible model capable of jointly modeling both bulk and tail behaviors. The resulting formulation provides a unified treatment of the three aforementioned extremal dependence structures while preserving permutation asymmetry. Based on copula theory and maximum likelihood estimation, the method is validated through finite-sample simulations that demonstrate its computational feasibility and favorable parameter estimation performance. Empirical analysis of temperature data successfully reveals intricate tail structures, underscoring the model’s theoretical rigor and practical utility.

This work addresses probabilistic density estimation for high-dimensional complex data. We propose a neural Copula modeling framework that explicitly decouples marginal distributions from dependency structures—a departure from conventional joint modeling approaches. Integrating Copula theory with deep neural networks, our method employs a normalizing flow architecture to separately parameterize marginals and the high-dimensional copula function, enabling differentiable density evaluation, generative sampling, and tractable inference. Leveraging marginal-joint decoupled training and a differentiable mutual information estimator, the model achieves superior performance over kernel density estimation and state-of-the-art neural density estimators across diverse complex distributions. Empirical evaluation demonstrates its effectiveness in high-accuracy mutual information estimation and photorealistic data generation, validating both theoretical design and practical utility.

In clinical trials—particularly ophthalmology—homogeneity testing for bilateral proportion data is commonly constrained by pre-specified, inflexible dependence structures (e.g., independence, perfect positive/negative dependence), limiting interpretability and adaptability. This paper introduces the Clayton copula—a flexible, interpretable tool for modeling asymmetric lower-tail dependence—into bilateral proportion homogeneity testing for the first time, thereby eliminating reliance on a priori dependence assumptions. We propose three Clayton copula–based test statistics and rigorously evaluate them via Monte Carlo simulation, demonstrating well-controlled Type I error rates and superior statistical power. Furthermore, we validate the robustness and practical utility of our approach on two real-world ophthalmologic datasets. This work establishes a theoretically rigorous, computationally feasible, and clinically meaningful testing paradigm for bilateral proportion data, advancing both methodological foundations and applied biostatistical practice.

This paper addresses the challenge of identifying heterogeneous individual-level choice behaviors from macro-level aggregate selection data. To this end, it establishes, for the first time, a systematic theoretical linkage between ordered probit choice models and copula theory, mapping individual heterogeneity onto the structural form of copula functions. The authors propose an analytically tractable representation based on extreme-value theory, enabling unique and unbiased identification of both heterogeneity types and their mixing weights. Methodologically, the approach integrates copula modeling, extreme-value function analysis, and structural identification theory to derive a general closed-form extreme-value representation. This framework overcomes key limitations of conventional aggregate modeling—such as loss of behavioral granularity and identifiability constraints—thereby substantially improving the accuracy, interpretability, and structural fidelity of micro-behavioral inference. It introduces a novel paradigm for discrete choice analysis, behavioral econometrics, and multivariate dependence modeling.

This work addresses the inaccuracy of posterior approximations in variational inference caused by limited expressiveness of conventional variational families. To overcome this limitation, the authors propose a novel variational method that integrates wavelet basis representations with copula-based dependency structures. The approach leverages discrete wavelet transforms to flexibly capture complex marginal posterior densities and explicitly models high-dimensional dependencies among parameters using copula functions. Efficient inference is achieved through Monte Carlo estimation of the evidence lower bound (ELBO), automatic differentiation, and gradient-based optimization. As the first method to jointly employ wavelets and copulas in variational inference, it demonstrates competitive posterior mean estimates—on par with Markov chain Monte Carlo (MCMC)—and substantially improved uncertainty quantification over standard variational approaches across logistic regression, sparse linear models, and hierarchical models.

This work addresses the challenges of joint distribution modeling in probabilistic forecasting for irregular multivariate time series, where coupling between marginal distributions and dependency structures often introduces bias. To resolve this, we propose CopFITi, a novel model that decouples marginals from dependencies: it employs normalizing flows to flexibly model univariate marginal distributions and leverages a Gaussian mixture copula to capture complex multivariate dependencies. CopFITi is the first copula-based model for irregular multivariate time series that is explicitly constructed to satisfy marginal consistency by design. Experimental results demonstrate that CopFITi achieves state-of-the-art performance in joint density estimation, significantly improving both marginal calibration and overall predictive accuracy.

This work addresses the limitations of traditional D-vine copula fitting, which relies on greedy strategies prone to local optima and lacks support for interpretable local anomaly detection. The authors propose the first fully differentiable D-vine copula fitting framework, integrating beam search with gradient-based optimization to enhance global fit quality through multi-path exploration. Leveraging the hierarchical dependency structure inherent in D-vines, the method enables edge-level anomaly localization alongside global anomaly scoring. Furthermore, it incorporates Mondrian conformal prediction to provide statistically valid uncertainty quantification for local anomalies. Experimental results demonstrate that the proposed approach significantly outperforms existing methods across multiple benchmark and real-world datasets, achieving superior performance in both interpretability and anomaly detection accuracy.

Traditional geostatistical methods rely on second-order moments and Gaussian assumptions, which are inadequate for capturing non-Gaussian spatial dependence. This work addresses this limitation by leveraging Sklar’s theorem to introduce a copula-based framework that decouples marginal distributions from the spatial dependence structure. The authors systematically develop spatial copula models applicable at both fixed point sets and process levels, emphasizing Kolmogorov consistency to clarify distinctions between these two modeling paradigms. The framework is further extended to spatio-temporal settings and supports flexible marginal specifications. By integrating spatial statistics, copula theory, and stochastic processes, this study establishes a unified approach for modeling non-Gaussian spatial dependence, elucidating the relationships, strengths, and limitations of existing methodologies, thereby advancing both theoretical understanding and practical applications in the field.

This study addresses the challenge of modeling edge effects and dependence structures that evolve with covariates—such as age—in multivariate responses of mixed types. Existing approaches are often hindered by strong assumptions or insufficient flexibility. To overcome these limitations, this work proposes a Bayesian nonparametric framework that integrates adaptive spline-based marginal regression with a covariate-dependent Gaussian copula infinite mixture model. A probit stick-breaking process is introduced to flexibly capture the covariate-driven evolution of dependence patterns, avoiding restrictive global correlation matrix constraints. The method unifies heterogeneous response types and dynamic dependencies through varying-coefficient copula regression and employs Markov chain Monte Carlo algorithms for posterior inference. Simulation studies demonstrate its accuracy and robustness, while empirical analysis of the 2023 Behavioral Risk Factor Surveillance System (BRFSS) data reveals complex age-varying marginal and dependence structures in health outcomes.