Log_mix for missing categorical data

Check out the solution in this thread: Imputation of a 3 category covariate to model a binary outcome