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.

A Note on Ising Network Analysis with Missing Data

Simulation-Based Performance Evaluation of Missing Data Handling in Network Analysis

A Note on the Robustness of a Full Bayesian Method for Nonignorable Missing Data Analysis

Network Imputation for a Spatial Autoregression Model with Incomplete Data

Tutorial: Assessing the impact of nonignorable missingness on regression analysis using Index of Local Sensitivity to Nonignorability.

Model-based Inference for Two-Stage Cluster Samples Subject to Nonignorable Item Nonresponse

Bayesian Model Selection for High-Dimensional Ising Models, With Applications to Educational Data

An Ising Similarity Regression Model for Modeling Multivariate Binary Data

The Analysis of Social-Science Data with Missing Values

Full-semiparametric-likelihood-based inference for non-ignorable missing data

A Note on Estimating Network Dependence in a Discrete Choice Model

Nonstandard Conditionally Specified Models for Nonignorable Missing Data.

A new method for constructing networks from binary data

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

A unified framework of analyzing missing data and variable selection using regularized likelihood

On Missing Data Imputation for IRB Models

Simulation-based Sensitivity Analysis for Non-ignorable Missing Data

Missing at Random or Not: A Semiparametric Testing Approach

Semiparametric Inference for Estimating Equations with Nonignorably Missing Covariates

Maximum likelihood estimation of missing data probability for nonmonotone missing at random data