Dear members, I am analyzing a cross-sectional big dataset includes serval CVDs and demographic data and want to plot a plot like the attachment. Our dataset is post-imputated, Hence, I am using a code: mi estimate: “mixed sbpa age if medication!=1||country:, reml emonly” However after this model, It seems “predict plot1, fitted” doesn’t work as usual. It can use if I drop the mi estimate commend. I wanna plot this graph with 95%CI fitted in the line and without the scartter of true value. Could you please kindly give me some advice? Many Thanks for your help!