Proc GEMMOD

proc genmod data= ;

class ;

model y = x / type3 ;

repeated subject= /  ;

lsmeans x / diff=control ('ref') cl adjust=Dunnett(or tukey);

run;


Parameters estimated with maximum likelihood methods

Repeated statement: specifies the covariance structure of multivariate responses and iterative fitting algorithm for GEE model fitting.