《直接输出科研论文统计表格的SAS宏》中sas宏的疑问-经管之家官网！

人大经济论坛-经管之家 收藏本站
搜索
经济学 管理学 金融学 统计学
您当前的位置> 软件培训>>

SAS软件培训

>>

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

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

发布：dahufa123 | 分类：SAS软件培训

搜索

关于本站

人大经济论坛-经管之家：分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

获取电子版《CDA一级教材》

完整电子版已上线CDA网校，累计已有10万+在读~ 教材严格按考试大纲编写，适合CDA考生备考，也适合业务及数据分析岗位的从业者提升自我。

完整电子版已上线CDA网校，累计已有10万+在读~ 教材严格按考试大纲编写，适合CDA考生备考，也适合业务及数据分析岗位的从业者提升自我。

市场营销

人力资源

国际贸易

TOP热门关键词

专题页面精选

北大刘建蒙教授2010年在《中国卫生统计》上发表的文章《直接输出科研论文统计表格的SAS宏》文后附录了一个较完整的sas宏（省略了数值变量的统计），我将其抄下来在sas中运行发现有错误，能告诉我我哪错了还是程序有错 ...
免费学术公开课,扫码加入
北大刘建蒙教授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_copycode($('code_Okp'));">复制代码
    「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
    免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
    涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
    来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
    经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
    扫描下方二维码下载并注册APP
本文关键词：

本文论坛网址：https://bbs.pinggu.org/thread-1125632-1-1.html

上一篇 | 【问题】那位大侠能解释下边际收益递减规 ...
下一篇 | 提高创意产业潜力认知——国家创意产业政 ...
SAS软件培训  精彩帖子推荐 更多

您可能感兴趣的文章

本站推荐的文章

人气文章

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

本文链接网址：https://bbs.pinggu.org/jg/ruanjianpeixun_sasruanjianpeixun_1125632_1.html

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责；
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
学术培训 |认证培训 |注册会计师 |市场营销 |人力资源 |国际贸易 |学报
数据分析师 人大经济论坛 大学 专业 手机版

为做大做强人大经济论坛，本站接受风险投资商咨询，请联系（010-62719935）
京ICP备11001960号 京ICP证090565号 京公网安备1101084107号 论坛法律顾问：王进律师 知识产权保护声明 免责及隐私声明
邮箱：service@pinggu.org
合作咨询电话：(010)62719935 广告合作电话：010-68456523 13661292478 刘老师
投诉电话：(010)68466864 不良信息处理电话：13661292478

联系客服
值班时间：工作日（9:00--18:00）