proc format; value yearfmt 1='Seventies' 2='Eighties' 3='Nineties'; value sexfmt 1='Not wrong' 2='Sometimes wrong' 3='Always wrong'; data gss; do i=1 to 3; do j=1 to 3; input ct @; premarital=i; decade=j; lxl=i*j; output; end; end; format premarital sexfmt. decade yearfmt.; datalines; 2445 3162 2650 1696 1676 1315 3226 2858 2228 ; proc freq; weight ct; table premarital*decade; run; proc genmod data=gss; class premarital decade; model ct=premarital decade lxl/link=log dist=poisson type3; run;