利用SAS软件制作报表的常用手法 |
您所在的位置：网站首页 › 统计描述分析怎么做表格模板 › 利用SAS软件制作报表的常用手法 |
利用SAS软件制作报表的常用手法
|
如何利用SAS 软件制作报表常用的技巧
1、先设计好你的报表格式 2、用SAS过程先初步得到这些结果 1. 比如可以用means过程、univariate过程得到一些描述统计量； 2. 可以用means过程后面加选项t来做t成组检验，可以用ttest过程来做成组或配对t检验。 而t检验的前途条件要求数据服从正态分布，如果不服从正态分布可以用非参数检验。 3、利用SAS的ODS中ods select选择你所需要的结果 可以先用ods trace来看一下输出的结果都有哪些，然后用ods select选择结果到相应的数据集中。4、为了直接输出得到上面的表格，可以有两种方法来做。 一种是ODS；一种是利用put在结果输出窗口得到它，再粘贴到word中。 参考代码如下： /* 该sas 文件旨在说明如何编写一次完成以后不用任何人手更改（包括更改日期、excel的格式等）的常规报表 参考了链接如下：http://bbs.pinggu.org/thread-1444265-1-1.html*/ options noxwait noxsync; options validvarname=any; /*该文件为固定报表模板，可以事先调整好单元格格式、字体颜色等*/ x '"C:\Users\Administrator\Desktop\报表模板.xlsx"'; /*设置各类时间,如test2012-05-16，就是用当天的时间、 sas程序运行的时间去得到这个"2012-05-16";另外一种是你要生成的excel文件是包含时间的，也在这里处理得出 */ data _null_; x = put(date()-1,yymmdd10.); /* 比如每日运行这个程序，处理前一日的文件，就是用date()-1 */ y =substr(x,6,2) || substr(x,9,2)||"b"; z =compress(input(substr(x,6,2),best8.)||"."||input(substr(x,9,2),best8.)); call symput('y',y); call symput('z',z); run; %let log=alla_&path; /* 我要处理的文件名就是 alla_0516b 这种形式 */ libname result "D:\test"; data temp1(compress=yes); set result.&log; /* 类似这样应用 */ run; /* 中间是你数据处理的过程，省略 */ filename r1 dde 'excel|[报表模板.xlsx]自定义表名1!r5c1:r60c6' ; /* 对某张表某些单元格进行写入 */ data _null_; set result1; file r1 notab linesize=2000; /* DDE默认空格为分隔符，如果一个变量中间有空格将会分开到两个单元格，用notab即可避免， linesize赋予一个足够大的值，则过长的变量不会错行 */ put date '09'x time '09'x source '09'x duration '09'x sql '09'x type; run; filename r1 dde 'excel|[报表模板.xlsx]自定义表名2!r5c1:r60c6' ; /* 继续写入下一张表 */ data _null_; set result2; file r1 notab linesize=2000; put date '09'x time '09'x source '09'x duration'09'x sql '09'x type; run; filename r1 dde "excel|system"; data _null_; file r1; put '[workbook.activate("自定义表名2")]'; /* 激活其中一张表 */ put '[row.height(0,"A1:A1",false,3)]'; /* 调整行高；类似这样的跟vba比较像 */ put '[workbook.activate("自定义表名1")]'; put '[row.height(0,"A1:A1",false,3)]'; x= compress('[save.as("D:\sql('||&path2||').xlsx")]'); /* 存储一个备份到某个路径，文件名为 sql(5.16).xlsx */ put x; y= compress('[save.as("E:\MailFile\sql('||&path2||').xlsx")]'); /* 存储到邮件文件夹，这样由邮件自动发送excel出去 */ put y; put '[quit]'; run; /*整个程序如上，然后txt写如下内容另存为bat文件，在windows-附件-系统工具-任务计划程序里面设置每日凌晨运行这个bat即可： D: cd: "D:\sas_program\sas\sasfoundation\9.2\" sas.exe -sysin "D:\thisprogramname.sas" -altlog "D:\test\log.log" */ /* 上述bat第二行是你sas程序的路径，第三行表示执行的sas程序的路径和名字，然后将sas运行的日志写入到log.log中， 以便事后查看日志 */ /*========================================================SAS 报表设计案例=====================================*/ /*四格表指标统计分析报表sas宏程序*/ %macro chisq (database= ,var= ,varfmt=,index=) ; proc freq data = &database. noprint ; tables &var. * group /out = TabFreq ; run ; data ma ( where = (group ="A") ) mb ( where= (group ="B") ) ; set TabFreq (keep = &var. group COUNT) ; where &var. = . ; run ; ods listing close ; ods output CrossTabFreqs = CrossTabFreqs chisq = chisqFishersExact=FishersExact ; proc freq data = &database. ; tables &var. * group /chisq exact expectednopercent norow ; run ; ods listing ; proc sort out = EXPECTED(keep = EXPECTED) data = CrossTabFreqs ( where = ( EXPECTED not=. ) ) ; by Expected ; run ; data T(keep = T) ; set EXPECTED ; if n = 1 then T = EXPECTED ; if T = . then delete ; run ; proc sort out = Frequency ( keep = Frequency) data = CrossTabFreqs ; by descending Frequency ; run ; data n (keep = n) ; set Frequency ; if n = 1 then n = Frequency ; if n = . then delete ;run ; run ; data chisq1 (keep =value1 Prob1 ) ; set chisq ( rename = (value = value1 Prob =Prob1) ) ; where (Statistic = ('Chi-Square') | Statistic = ('卡方') ) ; run ; data chisq2 (keep =value2 Prob2 ) ; set chisq ( rename = (value = value2 Prob =Prob2) ) ; where ( Statistic = ( 'Continuity Adj. Chi-Square') | Statistic =('连续校正卡方') ) ; run ; data Fishersexact (keep = cvalue1) ; set Fishersexact ; where (Name1 = ('XP2_FISH') ) ; run ; data st (keep = ma mb n T value1 Prob1 value2 Prob2 cvalue1&var. c) ; merge ma ( rename = (COUNT = ma) ) mb ( rename = (COUNT = mb) ) n Tchisq1 chisq2 Fishersexact ; &var. = - 99 ; if ma = . then ma = 0 ; if mb = . then mb = 0 ; format value1 6.2 value2 6.2 ; run ; data outa ( where = ( group ="A") ) outb (where = (group ="B") ); set CrossTabFreqs ; format ColPercent 6.2 ; run ; data outt (drop = table group _TYPE_ _TABLE_ Expected Missing); merge outa (rename = ( Frequency = A ColPercent = percentA) ) outb ( rename = ( Frequency = B ColPercent =percentB) ) ; if &var. = . then &var. = - 99; run ; proc sort ;by &var. ;run ; data _null_ ; file print notitle ; merge outt st ; if &var. = - 99 then do ; if n = 40 and 1 = 40 and T > = 5 then do ; put @3"&index. " @70 'chisq = 'value1 @83 'P='Prob1 ; end ; put @5 "例数(缺失) " @22 A '(' MA ') ' @47 B '(' MB ') '; end ; else do ; put @5 &var. '(%) ' @22 A '(' PercentA ') '@47 B'('PercentB ') '; format &var. &varfmt..; end ; run ; proc datasets; delete chisq chisq1 chisq2 Crosstabfreqs ExpectedFishersexact Frequency Ma Mb N Outa Outb Outt st stt t Tabfreq ; quit ; %mend chisq ; data jixian; input id group$ SEX; datalines; 1 A 1 2 A 2 3 A 2 4 A 2 5 B 1 6 B 2 7 B 1 8 B 2 ; run; proc format; value sexfmt 1='男性' 2='女性'; run; %chisq (database=jixian ,var=sex ,varfmt=sexfmt,index=性别)其出来的结果如下：
|
【本文地址】
公司简介
联系我们
今日新闻 |
点击排行 |
|
推荐新闻 |
图片新闻 |
|
专题文章 |