如何用SAS画统计图，这是我的学习分享！

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

准备花时间深入学习SAS的画图功能，主要方式是看到好的统计图反过来研究具体的SAS程序，我会不定期更新学习笔记。

另外找了一些资料，供大家参考学习。
http://robslink.com/SAS/Home.htm
https://blogs.sas.com/content/topic/data-visualization/

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Dose Response Plot (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20171129
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH® SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20171129;
   
9. %mend;
   
10. 
    
11. /*建立数据集response,用于后续画图*/
    
12. data response;
    
13. label response = 'Mean Clinical Response'    /*label会体现在图的x、y轴label*/
    
14.       visit    = 'Visit';
    
15. input response    /*临床上关注的某个定量指标*/
    
16.       visit       /*访视次数，第0天、第1天……*/
    
17.       trt @@;     /*处理因素（药物种类）*/
    
18. cards;
    
19. 9.4 0 1 9.35 1 1 8.1 2 1 6.3 3 1 4 4 1 2.2 5 1 1.6 6 1 1.3 7 1 0.7 8 1
    
20. 0.8 9 1 0.65 10 1 9.35 0 2 9.37 1 2 8.7 2 2 8.2 3 2 7.8 4 2 4.4 5 2
    
21. 2 6 2 1.9 7 2 1.5 8 2 1.4 9 2 1.2 10 2
    
22. ;
    
23. run;
    
24. 
    
25. /*自定义格式，用于刻度值的自定义显示*/
    
26. proc format;
    
27.    value visitfmt
    
28.       0 = "Baseline"   1 = "Day 1"   2 = "Day 2"    3 = "Day 3"
    
29.       4 = "Day 4"      5 = "Day 5"   6 = "Day 6"    7 = "Day 7"
    
30.       8 = "Day 8"      9 = "Day 9"   10 = "Day 10";
    
31.    value trtfmt
    
32.       1 = "Super Drug" 2 = "Old Drug";
    
33. run;
    
34. 
    
35. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
36. ods listing gpath='C:\Users\Administrator\Desktop\test20171129';
    
37. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
38. ods graphics on / imagename="picture" outputfmt=pdf;     /*图片名称picture，输出类型pdf*/
    
39. proc sgplot data=response;
    
40.    format visit visitfmt. trt trtfmt.;    /*这里用到前面的自定义格式*/
    
41.    series x=visit y=response / group=trt markers;
    
42.    /*指定x、y轴变量，group指定分组变量，节点处画markers*/
    
43.    refline 1 / axis=x;     /*x=1处画参考线*/
    
44.    keylegend / position=topright across=1 location=inside;
    
45.    /*画图例，位置右上，一列，在frame里面*/
    
46.    xaxis grid;    /*画x轴网格线*/
    
47.    yaxis grid;    /*画y轴网格线*/
    
48. run;
    
49. ods _all_ close;
    
50. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

扫码加我 拉你入群

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

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

分享0 收藏54 回帖

关键词：如何用 统计图 destination Proceedings Procedures

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Cancer Mortality Plot (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20171130
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH? SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20171130;
   
9. %mend;
   
10. 
    
11. /*建立数据集cancer,用于后续画图*/
    
12. data cancer;
    
13. input cause : &$20.  /*死亡原因，各种癌症等*/
    
14.       deaths;        /*各种原因的死亡数*/
    
15. cards4;
    
16. Lung Cancer  160000
    
17. Colorectal Cancer  50000
    
18. Breast Cancer  45000
    
19. Pancreatic Cancer  40000
    
20. Prostate Cancer  35000
    
21. Leukemia  30000
    
22. Lymphoma  28000
    
23. Liver Cancer  26000
    
24. Ovarian Cancer  24000
    
25. Esophageal Cancer  22000
    
26. Bladder Cancer  20000
    
27. Kidney Cancer  18000
    
28. ;;;;     /*注意“cause : &$20.”“cards4”“;;;;”用法！*/
    
29. run;
    
30. 
    
31. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
32. ods listing gpath='C:\Users\Administrator\Desktop\test20171130';
    
33. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
34. ods graphics on / imagename="picture" outputfmt=pdf;     /*图片名称picture，输出类型pdf*/
    
35. /*按变量deaths降序*/
    
36. proc sort data=cancer;
    
37.    by descending deaths;
    
38. run;
    
39. proc sgplot data=cancer;
    
40.    dot cause / response=deaths;     /*y轴cause，x轴deaths*/
    
41.    yaxis discreteorder=data display=(nolabel);     /*y轴顺序按照排序后的data？没有label*/
    
42.    xaxis label='Number of Deaths';     /*指定x轴label*/
    
43. run;
    
44. ods _all_ close;
    
45. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Plasma Concentration Profile (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20171201
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH? SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20171201;
   
9. %mend;
   
10. 
    
11. /*建立数据集pk,用于后续画图*/
    
12. data pk;
    
13. input time conc vx vy vlabel :$20.
    
14. vx2 vy2 vlabel2 :$20.
    
15. vx3 vy3 vlabel3 :$20.;
    
16. cards;
    
17. 0.25 0.00 3.5 3 AUC=12.03 5 6.76 Cmax=6.76 0.75 8 Tmax=0.75
    
18. 0.25 1.00 . . . . . . . . .
    
19. 0.50 3.50 . . . . . . . . .
    
20. 0.75 6.7 . . . . . . . . .
    
21. 1 5.5 . . . . . . . . .
    
22. 1.5 4.5 . . . . . . . . .
    
23. 1.75 3.8 . . . . . . . . .
    
24. 2 3 . . . . . . . . .
    
25. 2.5 2.5 . . . . . . . . .
    
26. 3 2 . . . . . . . . .
    
27. 4 1 . . . . . . . . .
    
28. 6 0.6 . . . . . . . . .
    
29. 8 0.1 . . . . . . . . .
    
30. ;
    
31. run;
    
32. 
    
33. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
34. ods listing gpath='C:\Users\Administrator\Desktop\test20171201';
    
35. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
36. ods graphics on / imagename="picture" outputfmt=pdf noborder;
    
37. /*图片名称picture，输出类型pdf，图片不需要边框*/
    
38. proc sgplot data=pk noautolegend;   /*不需要图例*/
    
39.    band x=time lower=0 upper=conc;  /*色带，沿着x轴刷，从y=0向上到y=conc*/
    
40.    series x=time y=conc / markers   /*指定time为x轴变量，conc为y轴变量*/
    
41.       markerattrs=graphdata1 (symbol=circlefilled);  /*用markers画出小圆点*/
    
42.    xaxis offsetmin=0  /*x轴上最小值与y轴的距离比例*/
    
43.       label='Time relative to dose (hours)';
    
44.    yaxis max=8.0 offsetmin=0  /*y轴最大值为8*/
    
45.       label='Plasma concentration [ng/mL]';
    
46.    vector x=vx y=vy / xorigin=2 yorigin=1.5  /*画箭头，起点是(vx,vy)，终点(2,1.5)*/
    
47.       lineattrs=(color=black)
    
48.       datalabel=vlabel arrowshape=filled;  /*箭头注释，箭头形状*/
    
49.    vector x=vx2 y=vy2 / xorigin=0  /*画参考线，起点是(vx2,vy2)，终点(0,6.76)*/
    
50.       yorigin=6.76 lineattrs=
    
51.       (pattern=mediumdash color=black)  /*线条属性*/
    
52.       datalabel=vlabel2 noarrowheads;
    
53.    vector x=vx3 y=vy3 / xorigin=.75
    
54.       yorigin=0 lineattrs=
    
55.       (pattern=shortdash color=black) datalabel=vlabel3 noarrowheads;
    
56. run;
    
57. ods _all_ close;
    
58. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Mean Percent Change from Baseline (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20171202
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH? SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20171202;
   
9. %mend;
   
10. 
    
11. /*建立数据集response,用于后续画图*/
    
12. data response;
    
13. infile cards truncover;
    
14. input trt week perc lo hi end;
    
15. cards;
    
16. 1 0 0
    
17. 1 3 8 7.5 8.5
    
18. 1 6 7 6.5 7.5
    
19. 1 9 8 7.5 8.5
    
20. 1 12 9 8.5 9.5 9
    
21. 1 15 6 5.5 6.5 4.5
    
22. 2 0 0
    
23. 2 3 2 1.5 2.5
    
24. 2 6 -1 -1.5 -0.5
    
25. 2 9 2 1.5 2.5
    
26. 2 12 3 2.5 3.5 3
    
27. 2 15 4 3.5 4.5 4.5
    
28. ;
    
29. run;
    
30. 
    
31. /*自定义格式，用于图例的自定义显示*/
    
32. proc format;
    
33.    value trtfmt
    
34.       1 = "Placebo"   2 = "Treetment 10mg";
    
35. run;
    
36. 
    
37. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
38. ods listing gpath='C:\Users\Administrator\Desktop\test20171202';
    
39. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
40. ods graphics on / imagename="picture" outputfmt=pdf noborder;
    
41. /*图片名称picture，输出类型pdf，图片不需要边框*/
    
42. proc sgplot data=response;
    
43.    format trt trtfmt.;  /*使用自定义格式*/
    
44.    band y=week lower=12.1 upper=15 /  /*画色带，x=12.1到x=15*/
    
45.    transparency=0.8 fillattrs=graphdata1;  /*透明度0.8，指定填充属性*/
    
46.    scatter x=week y=perc / group=trt   /*画散点图，下T为lo，上T为hi*/
    
47.       yerrorlower=lo yerrorupper=hi
    
48.       markerattrs=(symbol=circlefilled)
    
49.    name="scat";   /*这里是指定图例？*/
    
50.    series x=week y=perc / group=trt  /*连线*/
    
51.       lineattrs=(pattern=solid);
    
52.    series x=week y=end / group=trt  /*连线*/
    
53.       lineattrs=(pattern=shortdash)
    
54.       markers markerattrs=(symbol=circle);   /**/
    
55.    xaxis integer values=(0 to 15 by 3) /*设置x轴*/
    
56.       label="Weeks in Treatment";
    
57.    yaxis max=10 label="Percent Change";   /*设置y轴*/
    
58.    refline 0;
    
59.    refline 13.5 / axis=x
    
60.       label="|-- Washout --|"
    
61.       labelloc=outside labelpos=min    /*纵向min*/
    
62.       lineattrs=(thickness=0px);    /*不要显示这条线*/
    
63.    keylegend "scat" / title="" noborder;  /*设置图例*/
    
64. run;
    
65. ods _all_ close;
    
66. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Clinical Study Design Plot (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180123
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH? SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20180123;
   
9. %mend;
   
10. 
    
11. /*建立数据集response,用于后续画图*/
    
12. data txsize;
    
13. input (c500 t500 c100 t100 c50 t50 c25 t25) (: percent3.);/*此方式有一定简化作用*/
    
14. cards;
    
15. 10% 5% 12% 10% 18% 18% . .
    
16. 20% 7% 20% 12% 20% 19% . .
    
17. 30% 9% 30% 17% 30% 22% 30% 29%
    
18. 40% 10% 40% 18% 40% 27% 40% 35%
    
19. 50% 10% 50% 20% 50% 29% 50% 39%
    
20. 60% 9% 60% 20% 60% 30% 60% 41%
    
21. 70% 9% 70% 20% 70% 29% 70% 41%
    
22. 80% 8% 80% 19% 80% 28% 80% 40%
    
23. 90% 7% 90% 18% 90% 27% 90% 39%
    
24. 95% 6% 95% 17% 95% 26% 95% 38%
    
25. ;
    
26. run;
    
27. 
    
28. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
29. ods listing gpath='C:\Users\Administrator\Desktop\test20180123';
    
30. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
31. ods graphics on / imagename="picture" outputfmt=pdf noborder;
    
32. /*图片名称picture，输出类型pdf，图片不需要边框*/
    
33. ods escapechar='~';  /*The ODS escape character is set to ~*/
    
34. proc sgplot data=txsize noautolegend;    /*不要自动画图例*/
    
35.    reg x=c500 y=t500 / curvelabel="n=500"  /*回归线标注"n=500"*/
    
36.       curvelabelpos=max  /*max表示回归线最右边*/
    
37.       degree=2    /*degree取值1-10越小越光滑*/
    
38.       nomarkers   /*原始点的位置不要画出来*/
    
39.       lineattrs=(thickness=2px  /*回归线粗细，2个像素*/
    
40.       pattern=solid);  /*实线*/
    
41. 
    
42.    reg x=c100 y=t100 / curvelabel="n=100"
    
43.       curvelabelpos=max degree=2 nomarkers
    
44.       lineattrs=(thickness=2px
    
45.       pattern=solid);
    
46. 
    
47.    reg x=c50 y=t50 / curvelabel="n=50"
    
48.       curvelabelpos=max degree=2 nomarkers
    
49.       lineattrs=(thickness=2px
    
50.       pattern=solid);
    
51. 
    
52.    reg x=c25 y=t25 / curvelabel="n=25"
    
53.       curvelabelpos=max degree=2 nomarkers
    
54.       lineattrs=(thickness=2px
    
55.       pattern=solid);
    
56. 
    
57.    refline 0.5 /axis=x transparency=0.5;  /*参考线x=0.5，透明度0.5*/
    
58.    xaxis label='Response rate in the control group' /*x轴label*/
    
59.       values=(0.0 to 1.0 by 0.1) grid; /*x轴刻度*/
    
60.    yaxis label='Anticipated smallest
    
61.       treatment effect (~{unicode delta_u} vs control)'
    
62.       /*注意~{unicode delta_u}*/
    
63.        min=0.0 max=0.5 grid;  
    
64.    inset ("~{unicode alpha} ="="5%"
    
65.       "~{unicode beta} ="="20%") /
    
66.       position=topright noborder
    
67.       textattrs=graphtitletext;
    
68. run;
    
69. ods _all_ close;
    
70. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Demographic Profile in a LATTICE Layout (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180124
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH? SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20180124;
   
9. %mend;
   
10. 
    
11. /*建立数据集demog,用于后续画图*/
    
12. data demog;
    
13.    infile cards missover;  /*这句不可少*/
    
14.    input (gender race trt) (:$20.) age @;  /*这种用法很巧妙！★★★★★*/
    
15.    do while(not missing(age));
    
16.       output;
    
17.       input age @;
    
18.    end;
    
19. cards;
    
20. male white control 10 20 20 20 30 25 10 4 6 2
    
21. male white experimental 12 5 2 54 3 43 12
    
22. male black control 12 3 5 10 5 32 4 3 30 3 20
    
23. male black experimental 5 45 6 13 22 54 12
    
24. female white control 12 45 45 1 2 3 54 2 1 2 35
    
25. female white experimental 5 23 5 2
    
26. female black control 12 50 25 84 30 15 8 15
    
27. female black experimental 2 35 4 3 45 40
    
28. ;
    
29. run;
    
30. 
    
31. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
32. ods listing gpath='C:\Users\Administrator\Desktop\test20180124';
    
33. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
34. /*ods graphics on / imagename="picture" outputfmt=pdf noborder;*/
    
35. ods graphics off / imagename="picture" outputfmt=pdf noborder;
    
36. /*图片名称picture，输出类型pdf，图片不需要边框*/
    
37. 
    
38. proc sgpanel data=demog;
    
39.    panelby gender race / layout=lattice;
    
40.    hbox age / category=trt datalabel;  /*datalabel,异常点标注*/
    
41.    rowaxis display=(nolabel);  /*行坐标不标注*/
    
42.    /*colaxis display=(nolabel); */
    
43. run;
    
44. 
    
45. ods _all_ close;
    
46. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. /*Demographic Profile in a PANEL Layout*/
   
2. proc sgpanel data=demog;
   
3.    panelby gender race / layout=panel columns=4;
   
4.    hbox age / category=trt;
   
5.    rowaxis display=(nolabel);
   
6. run;

复制代码

1. /*建立数据集demog,用于后续画图*/
   
2. data hyp;
   
3. label rawp='rawp p-value' adjp='weight adjusted p-value';
   
4. do trt='TreatmentA', 'TreatmentB';
   
5.    do hyp='H11', 'H12', 'H21', 'H22';
   
6.       input rawp adjp;
   
7.       output;
   
8.    end;
   
9. end;
   
10. cards;
    
11. 0.025 0.028
    
12. 0.004 0.022
    
13. 0.01 0.029
    
14. 0.006 0.029
    
15. 0.025 0.028
    
16. 0.007 0.06
    
17. 0.01 0.025
    
18. 0.006 0.025
    
19. ;
    
20. run;
    
21. 
    
22. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
23. ods listing gpath='C:\Users\Administrator\Desktop\test20180124';
    
24. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
25. /*ods graphics on / imagename="picture" outputfmt=pdf noborder;*/
    
26. ods graphics off;
    
27. /*图片名称picture，输出类型pdf，图片不需要边框*/
    
28. 
    
29. /*Comparison of p-values in a COLUMNLATTICE Layout*/
    
30. proc sgpanel data=hyp;
    
31.    panelby trt / novarname
    
32.       layout=columnlattice;
    
33.    series x=hyp y=rawp / markers
    
34.       markerattrs=(size=6pt);
    
35.    series x=hyp y=adjp / markers
    
36.       markerattrs=(size=6pt symbol=square);
    
37.    rowaxis label='P-value';
    
38.    colaxis label='Hypothesis';
    
39.    refline 0.05 /
    
40.    lineattrs=(pattern=shortdash);
    
41. run;
    
42. 
    
43. ods _all_ close;
    
44. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

好贴，必须要顶，楼主的经验总结很宝贵。

建议补充不同图形的适用范围，比如时间反应关系，时间剂量效应……

384390412 发表于 2018-1-24 15:53
建议补充不同图形的适用范围，比如时间反应关系，时间剂量效应……

好建议 ，适当时候统一补充