Bayesian multiple imputation for large-scale categorical data with structural zeros