Copula Graphical Models for Heterogeneous Mixed Data

TL;DR

This paper introduces a novel graphical modeling approach for heterogeneous, multi-group data with mixed discrete and continuous variables, leveraging Gaussian copulas to accurately recover underlying network structures.

Contribution

It develops a multi-group graphical model that accounts for data heterogeneity and mixed types using Gaussian copulas and fused penalties, improving structure recovery over existing methods.

Findings

Better graph structure recovery in simulations

Effective handling of mixed data types

Application to ecological data yields new insights

Abstract

This article proposes a graphical model that handles mixed-type, multi-group data. The motivation for such a model originates from real-world observational data, which often contain groups of samples obtained under heterogeneous conditions in space and time, potentially resulting in differences in network structure among groups. Therefore, the i.i.d. assumption is unrealistic, and fitting a single graphical model on all data results in a network that does not accurately represent the between group differences. In addition, real-world observational data is typically of mixed discrete-and-continuous type, violating the Gaussian assumption that is typical of graphical models, which leads to the model being unable to adequately recover the underlying graph structure. The proposed model takes into account these properties of data, by treating observed data as transformed latent Gaussian data, by means of the Gaussian copula, and thereby allowing for the attractive properties of the Gaussian distribution such as estimating the optimal number of model parameter using the inverse covariance matrix. The multi-group setting is addressed by jointly fitting a graphical model for each group, and applying the fused group penalty to fuse similar graphs together. In an extensive simulation study, the proposed model is evaluated against alternative models, where the proposed model is better able to recover the true underlying graph structure for different groups. Finally, the proposed model is applied on real production-ecological data pertaining to on-farm maize yield in order to showcase the added value of the proposed method in generating new hypotheses for production ecologists.