Hi, I'm conducting a mediation analysis with a binary outcome using proc causalmed. Can this procedure account for repeated observations (i.e., clustering) similar to a "repeated subject" statement in proc genmod? If so, where/how do I place this code? If not, does anyone know how to do this type of analysis?   Thanks!    ... View more

I have a competing risks model using phreg and am trying to find the group means through lsmeans. I'm not sure I understand the results; the estimates are incredibly small (or negative if I include control variables in the model) compared to the actual group means when I run the same model through proc glimmix as a linear regression model. The outcome is hospital length of stay. Here is some of the code:   proc phreg data=cif_icp_crani; class group / param=glm order=internal ref=last; model los_adj*outcome(0) = group / rl eventcode=1; format group groupf.; hazardratio group / diff=pairwise;; lsmeans group / means exp cl; run; proc glimmix data=cif_icp_crani; class group; model los_adj = group; lsmeans group; run;   Can someone please explain to me what the lsmeans output in a phreg model actually means and why the estimates are so different than what the actually group means are? Thank you. ... View more