data a;
input subjid 字符4.   n02aa n02ab n02ae;
cards;
1154 1   0   1
1155 1   1   1
1156 1   0   0
1157 0   1   0
1158 1   1   0
1159 1   1   1
1160 1   1   1
;
run;

data bb;
        length all $32766;
        set a;
/*        all="";*/
        array num _numeric_;
        do over num;
/*                if num=1 then all=cats(vlabel(num),all);*/
                all=cats(vvalue(num),all);
        end;
run;

proc freq data=bb;
        table all/out=cc;
run;

data a;
input subjid   n02aa n02ab n02ae;
cards;
1154 1   0   1
1155 1   1   1
1156 1   0   0
1157 0   1   0
1158 1   1   0
1159 1   1   1
1160 1   1   1
;
run;


proc sql;
create table b as select count(unique subjid) as n, n02aa, n02ab, n02ae from a where  n02aa=1 & n02ab=1 & n02ae=1;
quit;

你要找的是关联规则分析吗，请看此例：
[https://blog.csdn.net/Tiaaaaa/article/details/58615565]

这个数据集中variable有很多 （300个左右），不知道该怎么写一段macro实现两两variable的组合呢？ 而不仅仅是给出的例子，我自己写了一段，有3w多个组合。。。运行到一半系统提示没有空间显示了

基本思路：transpose之后（注意排序），为1的数据catx起来，一个人做成一条，然后proc freq算频数