data poli;
input gender $ party $ ideology$ count ;
datalines;
fe do vl 44
fe do sl 47
fe do mo 118
fe do sc 23
fe do vc 32
fe re vl 18
fe re sl 28
fe re mo 86
fe re sc 39
fe re vc 48
ma do vl 36
ma do sl 34
ma do mo 53
ma do sc 18
ma do vc 23
ma re vl 12
ma re sl 18
ma re mo 62
ma re sc 45
ma re vc 51
;
run;
proc print data=poli;
run;
proc logistic data=poli order=data;
class party gender/param=ref;
freq count;
model ideology = party gender /link=clogit aggregate scale=none;
output out = porb predprobs=I;
run;
proc logistic data=poli order=data;
class party gender/param=ref;
freq count;
model ideology = party /link=clogit aggregate scale=none;
output out = prob predprobs=I;
run;
proc freq data=prob noprint;
WEIGHT count;
tables party*ip_vl*ip_sl*ip_mo*ip_sc*ip_sl*ip_vc/
list nocum nopercent out=test;
run;
No comments:
Post a Comment