Dear SAS users, I would like to perfom a multilevel multinomial logistic regression analysis on a dataset with missing values. I have a problem with proc mianalyze. My outcome has 3 levels ("area_tert": 1, 2, 3). While for each imputed data set, I obtain estimates for each covariate at each outcome level, when I combine the estimates, i.e. using proc mianalyze, I only obtain one estimate for each covariate, and not at each outcome level (2 and 3) and don't understand these pooled estimates! Could you please help me? Many thanks! Below is my syntax: title "analysis per imputed dataset "; proc glimmix data=imputed method=RSPL noclprint ; class id; model area_tert (order=freq ref='1')=age gender Education / link=glogit dist=multinomial solution; nloptions tech=newrap; random Intercept / subject=id group=area_tert; by _imputation_; ods output ParameterEstimates=gparms FitStatistics=gfitstats CovParms=randomeffect; run; title "combined results - fixed effects"; proc mianalyze parms=gparms; modeleffects Intercept age gender Education; stderr Intercept age gender Education;*/ run; title "combined results - random effect"; proc mianalyze data=Randomeffect; modeleffects estimate ; stderr stderr; run; Best regards,
... View more