Multiple imputation for two-level hierarchical models with categorical variables and missing at random data