Abstract
Computer simulation has played a pivotal role in analyzing alternative organ allocation strategies in transplantation. The current approach to producing cohorts of organ donors and candidates for individual-level simulation requires directly re-sampling retrospective data from a transplant registry. This historical data may reflect outmoded policies and practices as well as systemic inequities in candidate listing, limiting contemporary applicability of simulation results. We describe the development of an alternative approach for generating synthetic donors and candidates using hierarchical Bayesian network probability models. We developed two Bayesian networks to model dependencies among 10 donor and 36 candidate characteristics relevant to waitlist survival, donor-candidate matching, and post-transplant survival. We estimated parameters for each model using Scientific Registry of Transplant Recipients (SRTR) data. For 100 donor and 100 candidate synthetic populations generated, proportions for each categorical donor or candidate attribute, respectively, fell within one percentage point of observed values; the interquartile ranges (IQRs) of each continuous variable contained the corresponding SRTR observed median. Comparisons of synthetic to observed stratified distributions demonstrated the ability of the method to capture complex joint variability among multiple characteristics. We also demonstrated how changing two upstream population parameters can exert cascading effects on multiple relevant clinical variables in a synthetic population. Generating synthetic donor and candidate populations in transplant simulation may help overcome critical limitations related to the re-sampling of historical data, allowing developers and decision makers to customize the parameters of these populations to reflect realistic or hypothetical future states.