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.

Structure estimation of binary graphical models on stratified data: application to the description of injury tables for victims of road accidents

On a general structure for hazard-based regression models: An application to population-based cancer research

A bayesian hierarchical spatial model for dental caries assessment using non-gaussian markov random fields.

Learning Undirected Graphical Models with Structure Penalty

Classifying victim degree of injury in road traffic accidents: a novel stacked DCL-X approach

Estimating heterogeneous graphical models for discrete data with an application to roll call voting

Applying Machine Learning Methods and Models to Explore the Structure of Traffic Accident Data

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

Bivariate Extreme Value Modeling for Road Safety Estimation.

Enhanced spatial modeling on linear networks using Gaussian Whittle-Matérn fields

Doubly structured sparsity for grouped multivariate responses with application to functional outcome score modeling

A probabilistic Bayesian inference model to investigate injury severity in automobile crashes

Applications of Generalized Additive and Bayesian Hierarchical Models for Areal Safety Analysis: Case Study of an Urban Multimodal Transportation System in Chicago, Illinois

Hierarchical approaches for flexible and interpretable binary regression models

Bayesian Lesion Estimation with a Structured Spike-and-Slab Prior

Bayesian Structure Learning in Multi-layered Genomic Networks

A Bayesian multivariate spatial approach for illness-death survival models

Using hierarchical Bayesian binary probit models to analyze crash injury severity on high speed facilities with real-time traffic data

Ensemble-based Model Selection for Imbalanced Data to Investigate the Contributing Factors to Multiple Fatality Road Crashes in Ghana.

Construction Safety Risk Modeling and Simulation

Insurance pricing with hierarchically structured data: An illustration with a workers' compensation insurance portfolio