Abstract:Background: Availability of linked biomedical and social science data has risen dramatically in past decades, facilitating holistic and systems-based analyses. Among these, Bayesian networks have great potential to tackle complex interdisciplinary problems, because they can easily model inter-relations between variables. They work by encoding conditional independence relationships discovered via advanced inference algorithms. One challenge is dealing with missing data, ubiquitous in survey or biomedical datasets. Missing data is rarely addressed in an advanced way in Bayesian networks; the most common approach is to discard all samples containing missing measurements. This can lead to biased estimates. Here, we examine how Bayesian network structure learning can incorporate missing data. Methods: We use a simulation approach to compare a commonly used method in frequentist statistics, multiple imputation by chained equations (MICE), with one specific for Bayesian network learning, structural expectation-maximization (SEM). We simulate multiple incomplete categorical (discrete) data sets with different missingness mechanisms, variable numbers, data amount, and missingness proportions. We evaluate performance of MICE and SEM in capturing network structure. We then apply SEM combined with community analysis to a real-world dataset of linked biomedical and social data to investigate associations between socio-demographic factors and multiple chronic conditions in the US elderly population. Results: We find that applying either method (MICE or SEM) provides better structure recovery than doing nothing, and SEM in general outperforms MICE. This finding is robust across missingness mechanisms, variable numbers, data amount and missingness proportions. We also find that imputed data from SEM is more accurate than from MICE. Our real-world application recovers known inter-relationships among socio-demographic factors and common multimorbidities. This network analysis also highlights potential areas of investigation, such as links between cancer and cognitive impairment and disconnect between self-assessed memory decline and standard cognitive impairment measurement. Conclusion: Our simulation results suggest taking advantage of the additional information provided by network structure during SEM improves the performance of Bayesian networks; this might be especially useful for social science and other interdisciplinary analyses. Our case study show that comorbidities of different diseases interact with each other and are closely associated with socio-demographic factors.

Bayesian Inference of Networks Across Multiple Sample Groups and Data Types

Bayesian learning of multiple directed networks from observational data

Bayesian dynamic network modelling: an application to metabolic associations in cardiovascular diseases

Bayesian Structure Learning in Multi-layered Genomic Networks

Advances to Bayesian network inference for generating causal networks from observational biological data

A Full Bayesian Approach to Sparse Network Inference Using Heterogeneous Datasets

Inferring Interaction Networks from Multi-Omics Data.

A modeling framework for detecting and leveraging node-level information in Bayesian network inference

Bayesian mixed model inference for genetic association under related samples with brain network phenotype

Bayesian model-based clustering for populations of network data

NExUS: Bayesian simultaneous network estimation across unequal sample sizes

A Bayesian Approach to Joint Estimation of Multiple Graphical Models

Bayesian Approaches for Large Biological Networks

Treatment of missing data in Bayesian network structure learning: an application to linked biomedical and social survey data

Bayesian inference of sample-specific coexpression networks

Accounting for network noise in graph-guided Bayesian modeling of structured high-dimensional data

Bidirectional Inference Networks: A Class of Deep Bayesian Networks for Health Profiling

Bayesian inference for multiple Gaussian graphical models with application to metabolic association networks

Bayesian Joint Modeling of Multiple Brain Functional Networks

Clustering of heterogeneous populations of networks

Use of Bayesian networks to dissect the complexity of genetic disease: application to the Genetic Analysis Workshop 17 simulated data