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 two-step penalization and shrinkage approach for binary response data that is jointly separated and correlated: The effects of social networks on diarrheal disease

Modeling of Correlated Cognitive Function and Functional Disability Outcomes with Bounded and Missing Data in a Longitudinal Aging Study

Bayesian Modeling Longitudinal Dyadic Data with Nonignorable Dropout, with Application to a Breast Cancer Study

Bayesian Latent-Class Mixed-Effect Hybrid Models for Dyadic Longitudinal Data with Non-Ignorable Dropouts.

Distributed Learning from Multi-Site Observational Health Data for Zero-Inflated Count Outcomes

A Fully Nonparametric Modelling Approach to Binary Regression

Modelling collinear and spatially correlated data

Statistical methods for longitudinal and clustered designs with binary responses

Binary regression analysis with network structure of respondent-driven sampling data

Regression Analysis of Dependent Binary Data for Estimating Disease Etiology from Case-Control Studies

Combining Cox Regressions Across a Heterogeneous Distributed Research Network Facing Small and Zero Counts

The impact of clustering binary data on relative risk towards a study of inferential methods

A Correlated Network Scale-up Model: Finding the Connection Between Subpopulations

Flexible Bayesian Modeling for Longitudinal Binary and Ordinal Responses

Variable selection for sparse Dirichlet-multinomial regression with an application to microbiome data analysis

A note on marginalization of regression parameters from mixed models of binary outcomes

Covariate Microaggregation for Logistic Regression: An Application for Analysis of Confidential Data

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

Advancing Interpretable Regression Analysis for Binary Data: A Novel Distributed Algorithm Approach

Joint mixed-effects models for causal inference in clustered network-based observational studies