I am trying to find the adjusted means for values following multiple imputation using the mice package for a ANCOVA model (zScorePALAdj~Grouping+age+sex+education_years).
I have had to impute some missing values for a model, and then used the mi.anova function to provide summary output, that looks like the below:
mi.anova(mi.res=imputed_Data, formula="zScorePALAdj~Grouping+age+sex+education_years", type = 3)
I would like to get the adjusted means for the zScorePALAdj for the different groups from the 'grouping' variable.
I understand that emmeans would allow this if it were not pooled imputed data, for instance: emmeans(model, "Grouping"). However it does not seem to be able to use the output from the mi.anova model.
I was wondering if anyone can point me in the right direction as I am drawing a blank.
