《直接输出科研论文统计表格的SAS宏》中sas宏的疑问

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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

1. 
   
2. /* 调用宏前的准备*/
   
3. /* 定义分类变量的类别*/
   
4. /* 定义数值变量的缺失值*/
   
5. /* 定义基线特征变量标签*/
   
6. /* 定义宏变量*/
   
7. /* 说明:
   
8. Dataset指拟分析的数据集;
   
9. Group指分组变量;
   
10. Cat_var指分类变量名;
    
11. Cat_index 指基线表中该分类变量的顺次;
    
12. Cat_n分类变量的总数;
    
13. Num_var指数值变量名;
    
14. Num_index指基线表中该数值变量的顺次;
    
15. Num_n数值变量的总数;
    
16. */
    
17. %LET Dataset=ss;
    
18. %LET Group=s;
    
19. %LET Cat_var=%str(AA);
    
20. %LET Cat_index=%str(1);
    
21. %LET Cat_n=1;
    
22. %LET Num_var=%str(Age height weight fee);
    
23. %LET Num_index=%str(2 3 4 5);
    
24. %LET Num_n=4;
    
25. /* 宏主体及解释*/
    
26. /* 宏名称为Baseline, 由4步生成*/
    
27. %MACRO Baseline;
    
28.   
    
29.   /* 步骤1. 生成分类变量的输出要素*/
    
30.   /* 生成频数数据集*/  
    
31.   
    
32.   %MACRO Stepa;
    
33.     DATA Freq_a;
    
34.      SET &Dataset;
    
35.       Treat=(&Group=1)*1+(&Group=2)*0;
    
36.        Control=(&Group=1)*0+(&Group=2)*1;
    
37.     RUN;
    
38.    
    
39.     PROC MEANS DATA=Freq_a NOPRINT;
    
40.      CLASS &Cat_var;
    
41.       VAR Treat Control; OUTPUT OUT=Freq_b SUM=/AUTONAME;
    
42.        WAYS 1;
    
43.     RUN;
    
44.    
    
45.     %GLOBAL Total_treat Total_control;
    
46.    
    
47.     PROC SQL NOPRINT;
    
48.      SELECT COUNT(*) INTO:Total_treat FROM &Dataset WHERE &Group=1;
    
49.       SELECT COUNT(*) INTO:Total_control FROM &Dataset WHERE &Group=2;
    
50.     QUIT;
    
51.    
    
52.     DATA Freq_c;
    
53.      SET Freq_b;
    
54.       LENGTH Variable Level $200.;
    
55.        Treat=CATS(Treat_sum,"(", PUT((Treat_sum/&Total_treat)*100, 8.1),")");
    
56.         Control=CATS(Control_sum,"(",PUT((Control_sum/&Total_control)*100, 8.1),")");
    
57.          %DO i=1 %TO &CAT_n;
    
58.            IF %SCAN(&Cat_var,&i) THEN DO;
    
59.              Variable="%SCAN(&Cat_var,&i)";
    
60.               Index=%SCAN(&Cat_var,&i);
    
61.                Level=PUT(%SCAN(&Cat_var,&i),%SCAN(&Cat_var,&i).);
    
62. /*
    
63. 此处sas会提示错误，
    
64. “
    
65. NOTE: 由宏函数“SCAN”生成行。
    
66. 1     AA
    
67.       --
    
68.       48
    
69. ERROR 48-59: 输出格式 AA 没有找到或无法加载。
    
70. ”
    
71. */
    
72.            END;
    
73.          %END;
    
74.     RUN;
    
75.    
    
76.     /* 输出卡方值*/
    
77.     ODS LISTING CLOSE;
    
78.      ODS OUTPUT ChiSq=Freq_d;
    
79.      
    
80.     PROC FREQ DATA=&Dataset;
    
81.      TABLE &Group*(&Cat_var)/Chisq;
    
82.     RUN;
    
83.    
    
84.     ODS LISTING;
    
85.    
    
86.     DATA Freq_e(KEEP=Variable Prob);
    
87.      SET Freq_d;
    
88.       Variable=SCAN(Table,-1,"");
    
89.        WHERE Statistic="Chi-Square" or Statistic="卡方";
    
90.     RUN;
    
91.    
    
92.     /* 合并卡方与频数数据集*/
    
93.     PROC SORT DATA=Freq_c;
    
94.      BY Variable;
    
95.      
    
96.     PROC SORT DATA=Freq_e;
    
97.      BY Variable;
    
98.     RUN;
    
99.    
    
100.     DATA Freq_f;
     
101.      MERGE Freq_c Freq_e;
     
102.       BY Variable;
     
103.        KEEP Variable Level Treat Control Index Prob;
     
104.     RUN;
     
105.    
     
106.   %mend;
     
107.   
     
108.   %Stepa;
     
109.   
     
110.   /*步骤2. 生成数值变量的输出要素*/
     
111.   
     
112.   %MACRO Stepb;
     
113.    DATA Mean_f;RUN;/*此处建一空数据集代替*/
     
114.     /* 参考步骤1, 首先将T 检验的结果以均数(标准差) 、方差齐性检验和t检验3个数据集存储备用, 进而生成数值变量的输出要素*/
     
115.   %mend;
     
116.   
     
117.   %Stepb;
     
118.   
     
119.   /* 步骤3. 合并分类变量与数值变量输出要素, 并按顺序变量排序*/
     
120.   
     
121.   %MACRO Stepc;
     
122.     DATA Final;
     
123.      SET Freq_f Mean_f;
     
124.     RUN;
     
125.    
     
126.     PROC SORT DATA=Final;
     
127.      BY Index;
     
128.     RUN;
     
129.   %mend;
     
130.   
     
131.   %Stepc;
     
132.   
     
133.   /* 步骤4: 将输出元素以RTF类型输出并存至C盘根目录下的Example.rtf文件*/
     
134.   %MACRO Stepd;
     
135.     Options Nodate Nonumber Orientation=Portrait;
     
136.      ODS LISTING CLOSE;
     
137.       ODS RTF FILE="C:\Example.rtf" bodytitle;
     
138.        TITLE FONT="宋体" HEIGHT=10pt "表1 干预和对照组基线特征分布";
     
139.    
     
140.     PROC REPORT DATA=Final NOWD STYLE(REPORT)={BACKGROUND=WHITE RULES=GROUP FRAME=HSIDES CELLWIDTH=60%}
     
141.       STYLE(HEADER)={FONT_WEIGHT=LIGHT BACKGROUND=WHITE FONT_SIZE=10pt};
     
142.      COLUMN Variable Level Treat Control Prob;
     
143.       DEFINE Variable/GROUP ORDER=DATA "基线特征" STYLE(COLUMN)={JUST=LEFT} STYLE(HEADER)={JUST=LEFT};
     
144.        DEFINE Level/DISPLAY " " STYLE(COLUMN)={JUST=LEFT};
     
145.         DEFINE Treat/DISPLAY "治疗(n=%left(&total_treat))" STYLE(COLUMN)={JUST=RIGHT} STYLE(HEADER)={JUST=RIGHT};
     
146.        DEFINE Control/DISPLAY "对照(n=%left(&total_control))" STYLE(COLUMN)={JUST=RIGHT} STYLE(HEADER)={JUST=RIGHT};
     
147.           DEFINE Prob/GROUP "P" STYLE(COLUMN)={JUST=CENTER} STYLE(HEADER)={FONT_STYLE=ITALIC};
     
148.            COMPUTE BEFORE Variable/STYLE(LINES)={FONT_SIZE=2pt};
     
149.              LINE " ";
     
150.            ENDCOMP;
     
151.             FORMAT Variable $Variable.;
     
152. /*
     
153. 此处sas会提示错误，
     
154. “
     
155. ERROR: 输出格式 $VARIABLE 没有找到或无法加载。
     
156. 
     
157. ”
     
158. */
     
159. 
     
160.     RUN;
     
161.    
     
162.      ODS RTF CLOSE;
     
163.     ODS LISTING;
     
164.   %mend;
     
165.   
     
166.   %Stepd;
     
167. %mend Baseline;
     
168. %Baseline;

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏4 回帖

关键词：SAS宏 论文统计 科研论文 统计表 orientation 论文 SAS SAS 科研 表格

also want to know!!

看起来是你少了 $variable 格式的定义
你需要一个proc format来定义 $variable

ntsean 发表于 2011-6-27 09:47
看起来是你少了 $variable 格式的定义
你需要一个proc format来定义 $variable

是的，应该定义aa的格式。在数据步后，proc format;value aa 1='one' 2='two';run;

userzht 发表于 2014-7-24 09:58
是的，应该定义aa的格式。在数据步后，proc format;value aa 1='one' 2='two';run;

请问是在哪里添加？第几行呢？

没看懂

ntsean 发表于 2011-6-27 09:47
看起来是你少了 $variable 格式的定义
你需要一个proc format来定义 $variable

您好，请问在哪加？

yezi1107 发表于 2015-2-27 11:02
您好，请问在哪加？

proc format；
value aa 1='one' 2='two';
run;/*level变量所在的data步之上就行*/

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).);

天穹下 发表于 2015-3-3 13:29
proc format；
value aa 1='one' 2='two';
run;/*level变量所在的data步之上就行*/

非常感谢您，我去试试~

您好！最近也在研究这篇文章的SAS代码，请问代码中的步骤2（/*步骤2. 生成数值变量的输出要素*/），您有扩展代码吗？