Hi, I am modelling hospital mortality data (dichotomous rsponse variable) using GLIMMIX at three levels: patient, doctor, hospital. I have over 3million patients, about 80,000 doctors and about 70 hospitals. My code looks like: proc glimmix data=mortality method=quad; class doctor hospital; model death(event=1)= patient_covariates hospital_covariates doctor_covariates doctor hospital/dist=b s; random intercept/ subject=hospital; random intecept /subject=doctor(hospital); run; I can not make this model to converge. I wonder whether I have too many levels on the 'doctor' variable which has so many lavels? Any help would be appreciated!
... View more