data Speech; input Religion Racist Communist Militarist Speech Suicide Score; /* Religion here refers to an anti-religious speech. 1 indicates */ /* support of free speech in a given context, 0 indicates opposition */ /* We won't analyze the Suicide data here */ datalines; 1 1 1 1 571 105 4 0 1 1 1 11 0 3 1 0 1 1 64 1 3 0 0 1 1 19 0 2 1 1 0 1 67 4 3 0 1 0 1 18 0 2 1 0 0 1 18 4 2 0 0 0 1 30 3 1 1 1 1 0 66 10 3 0 1 1 0 14 1 2 1 0 1 0 36 3 2 0 0 1 0 39 3 1 1 1 0 0 63 62 2 0 1 0 0 53 16 1 1 0 0 0 49 444 1 0 0 0 0 270 724 0 ; run; /* Use order=data so that 1's are coded as the lower level of the data */ /* and 0's are used as the reference level */ proc genmod data=speech order=data; class score religion racist communist militarist; model speech=score religion racist communist militarist/dist=poi link=log; run;