第一次写按照SOC-PT分级的ADAE表格，在简书上记录下自己的思路，避免遗忘。
首先看一下AE表格的格式，发现第一列中包含了不同的层级。我的思路是分别对不同层级进行计数，最后汇总起来。
1654159395(1).png
读取数据后，output一次，用来计算不分级的事件数和例数。
data adae1;
set ads.adae;
where SAFFL="是";
if AEBODSYS="" then AEBODSYS="未编码";
if AEDECOD="" then AEDECOD="未编码";
run;
data adae2;
set adae1;
output;
if saffl="是" then aebodsys="all";
output;
run;
之后通过排序去重，生成第一列的模板。并通过retain语句生成对模板排序用的序号。之后第一次合并dummy1_1 dummy2_1，将ord1添加到dummy2_1中，使得AEDECOD标记上了所对应的AEBODSYS。再将AEDECOD以及AEBODSYS添加到一列中，按照ORD1 ORD2排序后，就可以得到我们最终需要的模板了。
proc sort data=adae2 out=dummy1(keep=aebodsys) nodupkey sortseq=linguistic(locale=zh_CN collation=pinyin numeric_collation=on);
by aebodsys;
run;
proc sort data=adae1 out=dummy2(keep=aedecod aebodsys) nodupkey sortseq=linguistic(locale=zh_CN collation=pinyin numeric_collation=on);
by aebodsys aedecod;
run;
data dummy1_1;
set dummy1;
retain ord1 0;
by aebodsys;
ord1=ord1+1;
run;
data dummy2_1;
set dummy2;
retain ord2 0;
by aebodsys aedecod;
ord2=ord2+1;
if first.aebodsys then ord2=2;
run;
proc sql;
create table dummy_1 as select a.*,b.ord1 from dummy2_1 a left join dummy1_1 b on a.aebodsys=b.aebodsys;quit;
data dummy;
set dummy_1(keep=aedecod ord1 ord2 rename=(aedecod=name)) dummy1_1(keep=aebodsys ord1 rename=(aebodsys=name));
if ord2=. then ord2=1;
proc sort;
by ord1 ord2;
run;
之后我使用的是freq步来进行计数，得到对应AEBODSYS AEDECOD的频数。之后按照分组，整理成最终需要的样子。
proc freq data=adae2(where=(aebodsys="ALL")) noprint;
table aebodsys*trtan/out=count1(drop=percent) nopercent;
run;
proc freq data=adae2(where=(aebodsys^="ALL")) noprint;
table aebodsys*trtan/out=count2(drop=percent) nopercent;
run;
proc freq data=adae2(where=(aebodsys="ALL")) noprint;
table aedecod*trtan/out=count3(drop=percent) nopercent;
run;
%macro filter(indata=,outdata=,var=);
data temp1 temp2;
set &indata.;
if trtan=1 then output temp1;
if trtan=2 then output temp2;
run;
proc sort data=temp1(rename=(count=count1));
by &var.;
run;
proc sort data=temp2(rename=(count=count2));
by &var.;
run;
data &outdata.;
merge temp1 temp2;
by &var.;
if count1=. then count1=0;
if count2=. then count2=0;
rename &var.=col1;
drop trtan;
run;
%mend;
%filter(indata=count1,outdata=freq1,var=aebodsys);
%filter(indata=count2,outdata=freq2,var=aebodsys);
%filter(indata=count3,outdata=freq3,var=aedecod);
data freq;
set freq1-freq3;
proc sort;
by col1;
run;
对于例数，这个方案的定义是每个人每个SOC-PT发生的不良事件只计算一次。
所以排序得到每个受试者每个SOC-PT唯一不良事件记录后，按照上面的步骤得到例数。最后将例数除以各组受试者的数目，或者例数的百分比。最后将整理好的数据output成RTF文件即可。
proc sort data=adae2 out=adae3 nodupkey;
by subjid aebodsys aedecod ;
run;
proc freq data=adae3(where=(aebodsys="ALL")) noprint;
table aebodsys*trtan/out=count4(drop=percent) nopercent;
run;
proc freq data=adae3(where=(aebodsys^="ALL")) noprint;
table aebodsys*trtan/out=count5(drop=percent) nopercent;
run;
proc freq data=adae3(where=(aebodsys="ALL")) noprint;
table aedecod*trtan/out=count6(drop=percent) nopercent;
run;
%filter(indata=count4,outdata=freq4,var=aebodsys);
%filter(indata=count5,outdata=freq5,var=aebodsys);
%filter(indata=count6,outdata=freq6,var=aedecod);
data freq_;
set freq4-freq6;
proc sort;
by col1;
run;
data freq_all;
merge freq freq_(rename=(count1=count3 count2=count4));
by col1;
rename col1=col1_;
run;
data fin;
length col1-col5 $200;
set freq_all;
col1=col1_;
if count1=0 then col3="0";
else if count1>0 then col3=put(count1,3.);
if count2=0 then col2="0";
else if count2>0 then col2=put(count2,3.)||" ("||strip(put(100*count2/&bign1.,5.1))||"）";
if count3=0 then col5="0";
else if count3>0 then col5=put(count3,3.);
if count4=0 then col4="0";
else if count4>0 then col4=put(count4,3.)||" ("||strip(put(100*count4/&bign2.,5.1))||"）";
keep col1-col5;
run;
proc sql noprint;
create table final as select a.*,b.* from dummy a left join fin b on a.name=b.col1;
quit;
proc sort data=final;
by ord1 ord2;
run;
因为是第一次做，可能还存在一些问题，打算之后发现了再回来修正。
