搜索
人大经济论坛 › 附件下载

附件下载

所在主题：
文件名:  直接输出科研论文统计表格的SAS宏.pdf
资料下载链接地址: https://bbs.pinggu.org/a-929434.html
附件大小:
238.63 KB   举报本内容
北大刘建蒙教授2010年在《中国卫生统计》上发表的文章《直接输出科研论文统计表格的SAS宏》文后附录了一个较完整的sas宏（省略了数值变量的统计），我将其抄下来在sas中运行发现有错误，能告诉我我哪错了还是程序有错误？

  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;

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

  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;

  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;

  45. %GLOBAL Total_treat Total_control;

  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;

  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;

  76. /* 输出卡方值*/
  77. ODS LISTING CLOSE;
  78. ODS OUTPUT ChiSq=Freq_d;

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

  84. ODS LISTING;

  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;

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

  96. PROC SORT DATA=Freq_e;
  97. BY Variable;
  98. RUN;

  100. DATA Freq_f;
  101. MERGE Freq_c Freq_e;
  102. BY Variable;
  103. KEEP Variable Level Treat Control Index Prob;
  104. RUN;

  106. %mend;

  108. %Stepa;

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

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

  117. %Stepb;

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

  121. %MACRO Stepc;
  122. DATA Final;
  123. SET Freq_f Mean_f;
  124. RUN;

  126. PROC SORT DATA=Final;
  127. BY Index;
  128. RUN;
  129. %mend;

  131. %Stepc;

  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 干预和对照组基线特征分布";

  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 没有找到或无法加载。

  157. ”
  158. */

  160. RUN;

  162. ODS RTF CLOSE;
  163. ODS LISTING;
  164. %mend;

  166. %Stepd;
  167. %mend Baseline;
  168. %Baseline;
复制代码


    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载，请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品，拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知，将积极的采取必要措施；同时，本站也将在技术手段和能力范围内，履行版权保护的注意义务。
(如有侵权，欢迎举报)
二维码

扫码加我 拉你入群

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

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

京ICP备16021002号-2 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2025-12-29 09:00