Thank you so much for your help! I am using SAS 9.3 so I will give that a try and see how it goes. I basically have a variable age that I have grouped into ranges using a format statement, and I would like to have a plot of age group on the x-axis vs. mortality rate in that group on the y-axis. So my code for producing a frequency table (without output) would look something like: proc freq data=work.temp;      tables age*dead_yn;      format age age_groups.; run; I will give this a try and let you know how it works out! ... View more

The hard part is interpreting anything like an HL stat in light of the clustering.  I would suggest doing a within-hospital lack of fit test for each hospital, as well as one overall that essentially ignored the clustering.  If the latter shows a lack of fit, it might then be quickly identified as being due to a specific hospital.  I think all of these tests would have to be obtained from PROC LOGISTIC, first with a by hospital statement, and then without. You might be able to take the within hospital p values as data for a generalized linear model with a beta distribution, and use the sample sizes as weights.  This might provide a better pooled representation than the "ignore the clustering" approach. Good luck. Steve Denham ... View more