自动生成标准的统计表格（配合ODS）的macro

job · 发表于 2008-12-26 12:37:14

注册后推荐绑定QQ，之后方才可以使用下方的“用QQ帐号登录”。

您需要登录才可以下载或查看，没有账号？立即注册

x

wmqy2004@acrp

程序在sas9.13中运行通过。

%macro comdts(j,k,dataname);
%do i=&j %to &k;
&dataname&i
%end;
%mend comdts;

%macro comparison(data=,vars=,vg1=,vg2=,vg3=,vg4=,class=,par=1,format=6.1,total=1,where=);
/*扫描宏参数中包含的应变量名称，并将其赋值给一系列宏变量*/
ods select none;
%let i=1;
%do %while(%scan(&vars,&i) ne );%let var&i=%scan(&vars,&i);%let i=%eval(&i+1);%end;
/*获取应变量个数*/
%let vnum=%eval(&i-1);
/*获取指定分组变量的水平数*/
proc freq data=&data;tables &class/out=freout;run;
data _null_;set freout end=last;if last then call symput('lvnum',compress(_n_));
run;
/*分组变量水平数小于2时给出警告信息，并跳出宏运行进程*/
%if &lvnum<2 %then %do;
%put Warning: Level number of class variable is less than 2.;
%goto exitmacro;
%end;

/*数据描述*/
%do i=1 %to &vnum;
      proc means data=&data n nmiss mean std min max lclm uclm median maxdec=2;
         output out=tmpmean&i n=n nmiss=nmiss mean=mean std=std min=min max=max lclm=lclm uclm=uclm median=median;class &class;%if &where ne %then %do;where &where;%end;
         var &&var&i;
      run;

      data tmpmean&i;set tmpmean&i;
         nnmiss=compress(put(n,6.0)||"("||left(put(nmiss,6.0))||")");
         ms=put(mean,&format)||"+"||left(put(std,&format));
         minax=put(min,&format)||"~"||left(put(max,&format));
         luclm=put(lclm,&format)||"~"||left(put(uclm,&format));
         med=put(median,&format);
      run;

%if &lvnum=2 %then %do;

/* t-test */;
      %if &par=1 %then %do;
         ods output ttests=ttests equality=equality;
         proc ttest data=&data;var &&var&i;class &class;
            %if &where ne %then %do;where &where;%end;
         data test2;merge ttests equality;by variable;
         data test21;set test2;where probf>0.05 and variances='Equal';length stat $ 5; stat="";
         data test22;set test2;where probf<0.05 and variances='Unequal';length stat $ 5;stat="*";
         data test(keep=variable pvalue);merge test21 test22;
         length variable $10.;variable="&&var&i.";
            pvalue=left(put(probt,6.3))||stat;
         run;
      %end;
/*秩和检验*/;
      %if &par^=1 %then %do;
         proc npar1way wilcoxon data=&data;
         %if &where ne %then %do;where &where;%end;
         ods output kruskalwallistest=tmp1;var &&var&i;class &class;run;
         data test(keep=variable pvalue);set tmp1(rename=nvalue1=p);
         where label1='Pr > Chi-Square';methods='**';
         pvalue=left(put(p,6.3))||methods;length variable $10.;variable="&&var&i.";run;
      %end;
/*描述数据*/;
/*给出两组之合计值*/;
      %if &total=1 %then %do;
         proc iml;reset noname linesize=100;
         use tmpmean&i;setin tmpmean&i;
         read all var {nnmiss ms minax luclm med} where(&class=&vg1) into nn1;
         read all var {nnmiss ms minax luclm med} where(&class=&vg2) into nn2;
         read all var {nnmiss ms minax luclm med} where(&class=.) into nn;
         nn11=nn1`;nn12=nn2`;nnt=nn`;
         aaa=nn11||nn12||nnt;bbb={"N(miss)","Mean+SD","Min~Max","95% CI(L~U)","Median"};
         ab=bbb||aaa;create mydata&i from ab;append from ab;quit;
      %end;
data mydata1&i;merge test mydata&i;run;
data mydata1&i;retain variable col1 col2 col3 col4 pvalue;set mydata1&i;
run;
/*不给出两组之合计值*/;

      %if &total^=1 %then %do;
         proc iml;reset noname linesize=100;
         use tmpmean&i;setin tmpmean&i;
         read all var {nnmiss ms minax luclm med} where(&class=&vg1) into nn1;
         read all var {nnmiss ms minax luclm med} where(&class=&vg2) into nn2;
         nn11=nn1`;nn12=nn2`;
         aaa=nn11||nn12;bbb={"N(miss)","Mean+SD","Min~Max","95% CI(L~U)","Median"};
         ab=bbb||aaa;create mydata&i from ab;append from ab;quit;
      %end;
data mydata1&i;merge test mydata&i;run;
data mydata1&i;retain variable col1 col2 col3 pvalue;set mydata1&i;
run;
%end;

%if &lvnum>2 %then %do;

/* ANOVA test */;
      %if &par=1 %then %do;
         proc glm data=&data;class &class;model &&var&i=&class;
            ods output overallanova=tmpanova diff=tmpf(rename=(_1=t12 _2=t13 _3=t23 p1=p12 p2=p13 p3=p23));
            lsmeans &class/tdiff pdiff;%if &where ne %then %do;where &where;%end;
         run;

         data t12(keep=pvalue);set tmpf;if _n_=1;rename t13=t1;
            p1=left(put(p13,6.3));pvalue=p1||'*';*Group1:group2;
         run;
         data t13(keep=pvalue);set tmpf;if _n_=1;rename t23=t2;
            p2=left(put(p23,6.3));pvalue=p2||'**';*Group1:Group3;
         run;
         data t23(keep=pvalue);set tmpf;if _n_=2;rename t23=t3;
            p3=left(put(p23,6.3));pvalue=p3||'***';*Group2:Group3;
         run;
         data tmpanova(keep=variable pvalue);set tmpanova;if _n_=1;
         length variable $10.;variable="&&var&i.";length pvalue $10.;
            pvalue=left(put(probf,6.3));run;
         data test;set tmpanova t12 t13 t23;run;
      %end;

/*秩和检验*/;
      %if &par^=1 %then %do;
         proc npar1way wilcoxon data=&data;%if &where ne %then %do;where &where;%end;
         ods output kruskalwallistest=tmp1;var &&var&i;class &class;run;

         data test(keep=variable pvalue);set tmp1(rename=nvalue1=p);
         where label1='Pr > Chi-Square';methods='**';
         pvalue=left(put(p,6.3))||methods;length variable $10.;variable="&&var&i.";run;
      %end;

/*描述数据*/;
/*给出三组的合计值*/;
      %if &total=1 %then %do;
         proc iml;reset noname linesize=100;
         use tmpmean&i;setin tmpmean&i;
         read all var {nnmiss ms minax luclm med} where(&class=&vg1) into nn1;
         read all var {nnmiss ms minax luclm med} where(&class=&vg2) into nn2;
         read all var {nnmiss ms minax luclm med} where(&class=&vg3) into nn3;
         read all var {nnmiss ms minax luclm med} where(&class=.) into nn;
         nn11=nn1`;nn12=nn2`;nn13=nn3`;nnt=nn`;
         aaa=nn11||nn12||nn13||nnt;bbb={"N(miss)","Mean+SD","Min~Max","95% CI(L~U)","Median"};
         ab=bbb||aaa;create mydata&i from ab;append from ab;quit;
      %end;
data mydata1&i;merge test mydata&i;run;
data mydata1&i;retain variable col1 col2 col3 col4 col5 pvalue;set mydata1&i;
run;

/*不给三组的合计值*/;
      %if &total^=1 %then %do;
         proc iml;reset noname linesize=100;
         use tmpmean&i;setin tmpmean&i;
         read all var {nnmiss ms minax luclm med} where(&class=&vg1) into nn1;
         read all var {nnmiss ms minax luclm med} where(&class=&vg2) into nn2;
         read all var {nnmiss ms minax luclm med} where(&class=&vg3) into nn3;
         nn11=nn1`;nn12=nn2`;nn13=nn3`;
         aaa=nn11||nn12||nn13;bbb={"N(miss)","Mean+SD","Min~Max","95% CI(L~U)","Median"};
         ab=bbb||aaa;create mydata&i from ab;append from ab;quit;
      %end;

data mydata1&i;merge test mydata&i;run;
data mydata1&i;retain variable col1 col2 col3 col4 pvalue;set mydata1&i;
run;
%end;%end;
data mydata;set %comdts(1,&vnum,mydata1);
ods select all;
proc print data=mydata noobs;run;
%exitmacro:
%mend comparison;
*应用举例：%comparison(data=aaa,vars=height weight age hcy1 sbp1 dbp1 sbp2 dbp2 sbp4 dbp4
,vg1=1,vg2=2,class=sex,par=1,format=6.1,total=1);
*%comparison(data=aaa,vars=age hcy1 hcy4 hcy8 sbp1 dbp1 sbp2 dbp2 weight height,vg1=1,vg2=2,
vg3=3(分组变量的值),class=group,par=1(参数检验),format=6.1,total=1(有合计值));

lhenghui · 发表于 2008-12-31 18:00:17

学习了，谢过lz！

shanshan1002 · 发表于 2009-2-13 21:53:05

先保存了，慢慢研究。谢谢了！

wmqy2004 · 发表于 2009-2-16 22:58:27

版主，这应该是我此前编写的宏程序啊，曾发表在其它的论坛上

xinwq · 发表于 2009-7-20 10:28:49

鲜见clinical trial 中使用proc iml的，楼主大概搞学术太久了吧？

账号		自动登录	找回密码
密码			立即注册

[分享] 自动生成标准的统计表格（配合ODS）的macro

注册后推荐绑定QQ，之后方才可以使用下方的“用QQ帐号登录”。