Simulation code used for the paper: Reference-based imputation methods based on conditional mean imputation (section 4)