抛砖引玉

proc sql;
create table a1 as
select
       distinct A.x as cat,
       count(A.x) as n,
           'x' as id
from a as A
group by a.x

union all

select
       distinct a.y as cat,
           count(a.y) as n,
           'y' as id
from a as A
group by a.y
;
quit;

%macro cnt;

proc sql noprint;
select name into: varlist separated by " " from dictionary.columns
where libname='WORK' and memname='A';
quit;

%let i=1;
%do %until(%scan(&varlist,&i)=%str( ));
proc sql;
create table a&i as
select "%scan(&varlist,&i)" as var, %scan(&varlist,&i) as cat,count(%scan(&varlist,&i)) as n
from a
group by 1,2;
quit;
%let i=%eval(&i+1);
%end;

%let cnt=%sysfunc(count(&varlist,%str( )));
data combine;
set
%do i=1 %to %eval(&cnt+1);
a&i
%end;
;
run;
%mend;
%cnt

data wanted;
set combine;
by notsorted var;
if not first.var then var='';
run;

我就想到简单的分开用proc univariate来做，再把结果set到一起

proc univariate data=a;
class x;
var x ;
output out=b
n=n;
run;

proc univariate data=a;
class y;
var y ;
output out=c
n=n;
run;

data b;
length var $1.;
set b;
var='x';
rename x=cat;
run;

data c;
length var $1.;
set c;
var='y';
rename y=cat;
run;

data d;
set b c;
run;

data b(keep=group cat);
set a;
array newvar _numeric_;
do over newvar;
  group=vname(newvar);
  cat=newvar;
  output;
end;
run;

proc tabulate data=b;
class group cat;
table group='var'*cat,n*f=2.;
run;