data smoking; input city $ smoking $ yes no @@; cases = yes + no; smoking1 = 0; if smoking = "Yes" then smoking1 = cards; Beijing Yes 126 100 Harbin Yes Beijing No 35 61 Harbin No Shanghai Yes 908 688 Zhengzhou Yes Shanghai No 497 807 Zhengzhou No Shenyang Yes 913 747 Taiyuan Yes Shenyang No 336 598 Taiyuan No Nanjing Yes 235 172 Nanchang Yes Nanjing No 58 121 Nanchang No ; run; 1; 402 121 182 72 60 11 104 21 308 215 156 98 99 43 89 36 proc genmod data = smoking; class city ; model yes/cases = city smoking1 / dist=bin link=logit ; run