Calculate Adjusted Event Rates

mekono — Wed, 20 Jan 2021 14:53:56 GMT

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;

Re: Calculate Adjusted Event Rates

SteveDenham — Thu, 21 Jan 2021 13:13:41 GMT

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

topic Calculate Adjusted Event Rates in Statistical Procedures

Calculate Adjusted Event Rates

Re: Calculate Adjusted Event Rates