Version 1.1.0 Edited on October 9, 2022. Doctoral candidate Reinier van Linschoten kindly pointed out a mistake in my R code for \(V_B\), the between imputation variance. The blog post now includes the corrected workflow. What? If you’re in the know, you know there are three major ways to handle missing data: full-information maximum likelihood, multiple imputation, and one-step full-luxury1 Bayesian imputation. If you’re a frequentist, you only have the first two options.