data ideology; input party $ ideology count @@; cards; R 1 30 R 2 46 R 3 148 R 4 84 R 5 99 D 1 80 D 2 81 D 3 171 D 4 41 D 5 55 ; proc logistic data = ideology order=data; class party /param = ref; freq count; model ideology = party /link=clogit scale=none aggregate; run; *Checking with nominal X2 and G2 methods (less powerful); proc freq; weight count; tables party*ideology/chisq; run;