Calculate Adjusted Event Rates

[mekono]

I'm trying to calculate the adjusted event rates for 2 patient cohorts, after adjusting the model (using proc glimmix) by propensity weights as well as some patient characteristics (age, sex, diabetes, smoking etc). 

After running my model (below) I'm not getting any values for the specific outcome.

Any advice on how to adjust my model to get the adjusted event rates, or is there another procedure I should use?

proc glimmix data= final_population;
class patient_group survival  gender (ref="0") daibetes (ref="0") smoking (ref="0") obesity  (ref="0") ;
model survival  (descending)=age gender daibetes smoking obesity propensity_score / dist=binary link=logit solution;
random intercept / subject=patient_group
run;

[SteveDenham · Accepted solution]

You will need to add an LSMEANS statement to get the estimates.  They will be reported as logits, so you will need to include the ILINK option.  If you want to calculate differences in risks, you will need to use the %nlmeans macro.

SteveDenham