Thank you both. This is so helpful thank you. ... View more

Sorry,   This is my SAS code, but what I want is lsmeans for married by sex at age=0 but without putting the interaction in the model. And so far, that does not seem possible. PROC GLM DATA=AGE; CLASS MARRIED SEX; MODEL JOB_PRESTIGE = MARRIED SEX AGE/SS3 INT; LSMEANS MARRIED/ AT (age) = (0); RUN;   Thank you ... View more

I have pre post survey data on subjects. The outcome is 1(not at all) to 4 (very). I used the code with Proc Glimmix but I'm not sure how to translate these subject specific probabilities to say the probability of group 1 at time 1 being in Category 1.   Thank you so much for your help. Audrey   proc GLIMMIX data=gsdatav2 METHOD=LAPLACE; class RE_M_GCARDEV_3(REF='1') RE_M_GCARDEV_Event(REF='1') GGENDER(REF='1') Survey (REF='1') NIH_GROUP (REF='2') ANALYSIS_ID; model M_R_GPATHCONS_11(desc) = RE_M_GCARDEV_3 RE_M_GCARDEV_Event GGENDER Survey NIH_GROUP Survey*NIH_GROUP /dist=multi oddsratio (DIFF=last label ) solution; random int/SUBJECT=ANALYSIS_ID; output out=ds2 pred(noblup)=predpa stderr(noblup)=stderrmupa pred(ilink noblup)=predmupa stderr(ilink noblup)=stderrmupa; run; ... View more

Thank you so much. I just glanced through and it might do the trick, I will have to look at it closely. I’m not proficient at estimate statements so if that is what I need, wish me luck. Once again, thank you! ... View more