data a;
input id sex age race;
n=1;
cards;
1 1 2 2
2 1 2 1
3 2 3 1
4 2 3 1
;
run;
data b;
input id sex age race;
n=2;
cards;
1 1 4 1
2 2 4 2
3 1 5 3
;
run;
proc append base=a data=b;
run;
%macro freq(var);
proc freq data=a noprint;
tables n*&var./chisq;
output OUT=&var. chisq  ;
run;
quit;
data &var.;
set &var.;
var="&var.";
run;
%mend freq;
%freq(sex);
%freq(age);
%freq(race);
data r(keep=p var result);
set age sex race;
rename p_pchi=p;
if p<0.05 then result="构成比无差异";
else  result="构成比有差异";
run;

写的很简单
如要变量多 数据量大
变量不一样的时就rename吧
一样的话写个宏吧
sushe1527来完成

如果P值接近0.05而且有些格子数目为0 考虑用FISHER精确检验的P值
不过数据量大变量多的时候FISHER检验会很慢

这点小事你就一起写了吧

感谢crackman！我研究下