SAS中文论坛

标题: 《直接输出科研论文统计表格的SAS宏》中sas宏的疑问 [打印本页]

作者: shiyiming 时间: 2011-6-27 06:39
标题: 《直接输出科研论文统计表格的SAS宏》中sas宏的疑问
《直接输出科研论文统计表格的SAS宏》中sas宏的疑问
北大刘建蒙教授2010年在《中国卫生统计》上发表的文章《直接输出科研论文统计表格的SAS宏》文后附录了一个较完整的sas宏（省略了数值变量的统计），我将其抄下来在sas中运行发现有错误，能告诉我我哪错了还是程序有错误？

[code:23jfjdst]

/* 调用宏前的准备*/
/* 定义分类变量的类别*/
/* 定义数值变量的缺失值*/
/* 定义基线特征变量标签*/
/* 定义宏变量*/
/* 说明:
Dataset指拟分析的数据集;
Group指分组变量;
Cat_var指分类变量名;
Cat_index 指基线表中该分类变量的顺次;
Cat_n分类变量的总数;
Num_var指数值变量名;
Num_index指基线表中该数值变量的顺次;
Num_n数值变量的总数;
*/
%LET Dataset=ss;
%LET Group=s;
%LET Cat_var=%str(AA);
%LET Cat_index=%str(1);
%LET Cat_n=1;
%LET Num_var=%str(Age height weight fee);
%LET Num_index=%str(2 3 4 5);
%LET Num_n=4;

/* 宏主体及解释*/
/* 宏名称为Baseline, 由4步生成*/

%MACRO Baseline;

  /* 步骤1. 生成分类变量的输出要素*/
  /* 生成频数数据集*/

  %MACRO Stepa;
DATA Freq_a;
   SET &Dataset;
   Treat=(&Group=1)*1+(&Group=2)*0;
   Control=(&Group=1)*0+(&Group=2)*1;
RUN;

PROC MEANS DATA=Freq_a NOPRINT;
   CLASS &Cat_var;
   VAR Treat Control; OUTPUT OUT=Freq_b SUM=/AUTONAME;
   WAYS 1;
RUN;

%GLOBAL Total_treat Total_control;

PROC SQL NOPRINT;
   SELECT COUNT(*) INTO:Total_treat FROM &Dataset WHERE &Group=1;
   SELECT COUNT(*) INTO:Total_control FROM &Dataset WHERE &Group=2;
QUIT;

DATA Freq_c;
   SET Freq_b;
   LENGTH Variable Level $200.;
   Treat=CATS(Treat_sum,"(", PUT((Treat_sum/&Total_treat)*100, 8.1),")");
      Control=CATS(Control_sum,"(",PUT((Control_sum/&Total_control)*100, 8.1),")");
      %DO i=1 %TO &CAT_n;
         IF %SCAN(&Cat_var,&i) THEN DO;
         Variable="%SCAN(&Cat_var,&i)";
            Index=%SCAN(&Cat_var,&i);
            Level=PUT(%SCAN(&Cat_var,&i),%SCAN(&Cat_var,&i).);

/*

此处sas会提示错误，

“

NOTE: 由宏函数“SCAN”生成行。
1    AA
   --
   48

ERROR 48-59: 输出格式 AA 没有找到或无法加载。
”

*/
         END;
      %END;
RUN;

/* 输出卡方值*/
ODS LISTING CLOSE;
   ODS OUTPUT ChiSq=Freq_d;

PROC FREQ DATA=&Dataset;
   TABLE &Group*(&Cat_var)/Chisq;
RUN;

ODS LISTING;

DATA Freq_e(KEEP=Variable Prob);
   SET Freq_d;
   Variable=SCAN(Table,-1,"");
   WHERE Statistic="Chi-Square" or Statistic="卡方";
RUN;

/* 合并卡方与频数数据集*/
PROC SORT DATA=Freq_c;
   BY Variable;

PROC SORT DATA=Freq_e;
   BY Variable;
RUN;

DATA Freq_f;
   MERGE Freq_c Freq_e;
   BY Variable;
   KEEP Variable Level Treat Control Index Prob;
RUN;

  %mend;

  %Stepa;

  /*步骤2. 生成数值变量的输出要素*/

  %MACRO Stepb;
DATA Mean_f;RUN;/*此处建一空数据集代替*/
/* 参考步骤1, 首先将T 检验的结果以均数(标准差) 、方差齐性检验和t检验3个数据集存储备用, 进而生成数值变量的输出要素*/
  %mend;

  %Stepb;

  /* 步骤3. 合并分类变量与数值变量输出要素, 并按顺序变量排序*/

  %MACRO Stepc;
DATA Final;
   SET Freq_f Mean_f;
RUN;

PROC SORT DATA=Final;
   BY Index;
RUN;
  %mend;

  %Stepc;

  /* 步骤4: 将输出元素以RTF类型输出并存至C盘根目录下的Example.rtf文件*/
  %MACRO Stepd;
Options Nodate Nonumber Orientation=Portrait;
   ODS LISTING CLOSE;
   ODS RTF FILE="C:\Example.rtf" bodytitle;
   TITLE FONT="宋体" HEIGHT=10pt "表1 干预和对照组基线特征分布";

PROC REPORT DATA=Final NOWD STYLE(REPORT)={BACKGROUND=WHITE RULES=GROUP FRAME=HSIDES CELLWIDTH=60%}
   STYLE(HEADER)={FONT_WEIGHT=LIGHT BACKGROUND=WHITE FONT_SIZE=10pt};
   COLUMN Variable Level Treat Control Prob;
   DEFINE Variable/GROUP ORDER=DATA "基线特征" STYLE(COLUMN)={JUST=LEFT} STYLE(HEADER)={JUST=LEFT};
   DEFINE Level/DISPLAY " " STYLE(COLUMN)={JUST=LEFT};
      DEFINE Treat/DISPLAY "治疗(n=%left(&total_treat))" STYLE(COLUMN)={JUST=RIGHT} STYLE(HEADER)={JUST=RIGHT};
   DEFINE Control/DISPLAY "对照(n=%left(&total_control))" STYLE(COLUMN)={JUST=RIGHT} STYLE(HEADER)={JUST=RIGHT};
      DEFINE Prob/GROUP "P" STYLE(COLUMN)={JUST=CENTER} STYLE(HEADER)={FONT_STYLE=ITALIC};
         COMPUTE BEFORE Variable/STYLE(LINES)={FONT_SIZE=2pt};
         LINE " ";
         ENDCOMP;
         FORMAT Variable $Variable.;

/*

此处sas会提示错误，

“

[color=#FF0000]ERROR: 输出格式 $VARIABLE 没有找到或无法加载。                                                                                                                                     [/color]


”

*/

RUN;

   ODS RTF CLOSE;
ODS LISTING;
  %mend;

  %Stepd;
%mend Baseline;

%Baseline;
[/code:23jfjdst]

作者: shiyiming 时间: 2011-6-28 22:40
标题: Re: 《直接输出科研论文统计表格的SAS宏》中sas宏的疑问
定义你的$VARIABLE format。

欢迎光临 SAS中文论坛 (http://mysas.net/forum/)