Programming the statistical procedures from SAS

estimating hospital average using glimmix

Reply
Occasional Learner
Posts: 1

estimating hospital average using glimmix

Dear all,

I have a question regarding estimating hospital effect using random intercept logistic regression.

For example, I have crude mortality rate and would like to know hospital variation on mortality.

 

proc glimmix;

class hospital;

model death=x1-x5/dist=binary link=logit ddfm=bw s;

random intercept/subject=hospital solution type=vc;

output out=glimmixout pred(blup ilink)=predprob
                         pred(noblup ilink)=predprob_pa;

run;

 

How can I get predicted hospital mortality or adjusted hospital mortality?

Can I use the random intercept for each hospital (using the random effect plus the estimate of the intercept) to compare the adjusted mortaliy across hospitals? This would assume that all covariates x1-x5 are at centered values, right?

 

Or should I use the predprob for each patient and sum them for each hospital and get an average probablity for each hospital?

 

Any thoughts about this?

 

Thanks a lot.

vivian

 

 

 

 

Ask a Question
Discussion stats
  • 0 replies
  • 184 views
  • 0 likes
  • 1 in conversation