Hello Steve, As discussed in this post, I finally got the right code for Proc Glimmix and I have been using it thus far. Only today when I started using it on much larger dataset (over 250000 records) . It is the same code as I have posted above. The Proc Glimmix worked fine on the 100000 cases. But on the new dataset witha around 250000 cases, the output I am getting is not complete. Here is the code.I am not sure why it is not yielding a complete. Any help is greatly appreciated. I am also attaching the output for your reference. proc glimmix data=library.nismi091011_gender5 ; *using oddsratio*/; class hospid gender racenew ZIPINC_QRTL_1 dm_all cm_htn_c cm_ANEMDEF cm_arth cm_bldloss cm_chf cm_chrnlung cm_coag cm_depress cm_drug cm_hypothy cm_liver cm_lymph cm_lytes cm_mets cm_neuro cm_obese cm_para cm_perivasc cm_pulmcirc cm_renlfail cm_wghtloss cararrhythmia; model died (event='1') = age gender racenew ZIPINC_QRTL_1 dm_all cm_htn_c cm_ANEMDEF cm_arth cm_bldloss cm_chf cm_chrnlung cm_coag cm_depress cm_drug cm_hypothy cm_liver cm_lymph cm_lytes cm_mets cm_neuro cm_obese cm_para cm_perivasc cm_pulmcirc cm_renlfail cm_wghtloss cararrhythmia /oddsratio dist=binary link=logit ddfm=bwsolution; random intercept / subject=hospid; where pcionly=1 and stemi=0; weight discwt; title 'Logi Reg in-hosp mortality vs gender in POST-PCI MI patients using "where" option with RACE +income for nonstemi with pcionly'; run; quit;
... View more