Hello all   My current work involves profiling of hospital performance regarding a specific binary quality indicator (Patients are nested within hospitals, i.e. 2-level data)   Performance will be adjusted for various risk factors, e.g. age, diagnosis etc.   My intention is then to use a random hospital intercept binary logit model (a fixed effect model  - estimating intercepts for each hospital using dummies - could perhaps  be an alternative ) ,    As a risk adjusted measure of hospital performance I would use so called "predictive margins"; i.e.:   1. Estimate the random intercept model including the empirical bayesian estimates (BLUPs) of the hospital specific intercepts 2. Assume all patients belong to hospital 1 - values on risk factors remain unchanged 3. calculate the predicted probability using the estimated fixed parameters and BLUP for hospital 1 for each patient and average over patients   Do 2 and 3 for every hospital.   Questions: 1. This is a build in command in stata (the margins command), but I guess i could obtain the same results using "store" in glimmix and proc plm? 2. the margins command in Stata uses the delta method to obtain standard errors of the predictive margins. Is there some way I could get the standard errors in SAS without having to bootstrap (time consuming since > 100 hospitals)   Thanks in advance.   Anders   ... View more