Hi everyone,

I hope you'll be able to help me with this question. I have multiple years of cross-sectional data for examining how a categorical independent variable x relates to a binary outcome y. I had some missing data, so I imputed the missing observations following appropriate procedures. I'd like to plot means and confidence intervals for the outcome variable over the different categories of the independent variable and over time while accounting for the imputed dataset. In other words, my code looks something like this:

HTML Code:
mi svyset id [pweight = svy_weight] 

mi estimate, post: mean i.y if x==0, over(svy_year)
estimates store x0

mi estimate, post: mean i.y if x==1, over(svy_year)
estimates store x1

mi estimate, post: mean i.y if x==2, over(svy_year)
estimates store x2

coefplot  x1 x2 x3, vert
Unfortunately, the resulting coefplot includes the mean values for when y==0 and when y==1. I'd like to plot only the means + cis for when y==1. Is there a way to do that? I've tried other commands (ciplot and cibar) but I can't seem to make them work while accounting for the multiply imputed data and survey weights.

Any guidance would be most appreciated!