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

1. /*建立数据集cancer,用于后续画图*/
   
2. data cancer;
   
3. input cause : &$20.  /*死亡原因，各种癌症等*/
   
4.       cases       /*各种原因的死亡数*/
   
5.       sex : &$20.;
   
6. cards4;
   
7. Lung Cancer  160000  female
   
8. Colorectal Cancer  50000  male
   
9. Breast Cancer  45000  female
   
10. Pancreatic Cancer  40000  male
    
11. Prostate Cancer  35000  female
    
12. Leukemia  30000  male
    
13. Lymphoma  28000  female
    
14. Liver Cancer  26000  male
    
15. Ovarian Cancer  24000  female
    
16. Esophageal Cancer  22000  male
    
17. Bladder Cancer  20000  female
    
18. Kidney Cancer  18000  male
    
19. ;;;;     /*注意“cause : &$20.”“cards4”“;;;;”用法！*/
    
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. /*Horizontal Bar Chart in a ROWLATTICE Layout*/
    
30. proc sgpanel data=cancer;
    
31. panelby sex / layout=rowlattice onepanel
    
32. novarname uniscale=column;
    
33. hbar cause / response=cases stat=sum;
    
34. rowaxis discreteorder=data
    
35. display=(nolabel);
    
36. colaxis display=(nolabel);
    
37. run;
    
38. 
    
39. ods _all_ close;
    
40. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. data heart;
   
2. input weight cholesterol systolic diastolic ageatdeath smoking;
   
3. cards;
   
4. 50 100 120 80 70 1
   
5. 60 120 130 70 60 0
   
6. 70 140 135 90 55 1
   
7. 80 180 125 100 50 1
   
8. ;
   
9. run;
   
10. 
    
11. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
12. ods listing gpath='C:\Users\Administrator\Desktop\test20180124';
    
13. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
14. /*ods graphics on / imagename="picture" outputfmt=pdf noborder;*/
    
15. ods graphics off;
    
16. /*图片名称picture，输出类型pdf，图片不需要边框*/
    
17. 
    
18. /*Paneled Graph of Independent Scatter Plots*/
    
19. proc sgscatter data=heart;
    
20. plot (weight cholesterol systolic
    
21. diastolic) * (ageatdeath) /
    
22. group=smoking columns=2
    
23. markerattrs=(size=10 symbol=circle)
    
24. legend=(notitle);
    
25. run;
    
26. 
    
27. ods _all_ close;
    
28. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
   
2. ods listing gpath='C:\Users\Administrator\Desktop\test20180124';
   
3. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
   
4. /*ods graphics on / imagename="picture" outputfmt=pdf noborder;*/
   
5. ods graphics off;
   
6. /*图片名称picture，输出类型pdf，图片不需要边框*/
   
7. 
   
8. /*Scatter Plot Matrix*/
   
9. proc sgscatter data= heart;
   
10. matrix weight systolic diastolic cholesterol/
    
11. ellipse=(alpha=0.05 type=predicted)
    
12. diagonal=(histogram normal);
    
13. run;
    
14. 
    
15. ods _all_ close;
    
16. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

支持分享

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Laboratory Panel of WBC Results (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180131
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH® SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20180131;
   
9. %mend;
   
10. 
    
11. /*建立数据集labs,用于后续画图*/
    
12. data labs;
    
13. /*label response = 'Mean Clinical Response'    label会体现在图的x、y轴label
    
14.       visit    = 'Visit';*/
    
15. input line :$20.    /*临床上关注的某个定量指标,即testname*/
    
16.       numlow
    
17.       numhi
    
18.       result
    
19.       visitnum :$20.       /*访视次数，第0天、第1天……*/
    
20.       patid
    
21.       panel :$20.;
    
22. cards;
    
23. Neutrophils(%) 40.5 75.0 48.60 PreRx 102 Hematology
    
24. Neutrophils(%) 40.5 75.0 48.60 week1 102 Hematology
    
25. Neutrophils(%) 40.5 75.0 48.60 week2 102 Hematology
    
26. Neutrophils(%) 40.5 75.0 48.60 week3 102 Hematology
    
27. Neutrophils(%) 40.5 75.0 48.60 week4 102 Hematology
    
28. Neutrophils(%) 40.5 75.0 48.60 week5 102 Hematology
    
29. Neutrophils(%) 40.5 75.0 48.60 week6 102 Hematology
    
30. Lymphocytes(%) 15.4 48.5 41.60 PreRx 102 Hematology
    
31. Basophils(%) 40.5 75.0 48.60 PreRx 102 Hematology
    
32. Eosinophils(%) 15.4 48.5 41.60 PreRx 102 Hematology
    
33. Monocytes(%) 15.4 48.5 41.60 PreRx 102 Hematology
    
34. WBC 15.4 48.5 41.60 PreRx 102 Hematology
    
35. ;
    
36. run;
    
37. 
    
38. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
39. ods listing gpath='C:\Users\Administrator\Desktop\test20180131';
    
40. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
41. /*ods graphics on / imagename="picture" outputfmt=pdf;     图片名称picture，输出类型pdf*/
    
42. ods graphics off;
    
43. 
    
44. /*Laboratory Panel of WBC Results*/
    
45. ods listing style=statistical;
    
46. proc sgpanel data=labs;
    
47.    panelby line /
    
48.    /*如果要指定顺序，就赋line变量为1、2、3再label*/
    
49.    onepanel uniscale=column
    
50.    novarname layout=rowlattice;
    
51.    refline numlow / label noclip;/*参考线label注释*/
    
52.    /*如果参考线在数据范围之外，也要显示*/
    
53.    refline numhi / label noclip;
    
54.    scatter x=visitnum y=result /
    
55.    markerattrs=(symbol=asterisk size=10);
    
56.    rowaxis display=(nolabel novalues noticks);
    
57.    /*rowaxis不要坐标注释、不要坐标值、不要坐标点*/
    
58.    colaxis offsetmax=0.1 display=(nolabel);
    
59. run;
    
60. 
    
61. ods _all_ close;
    
62. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. /*建立数据集labs,用于后续画图*/
   
2. data labs;
   
3. /*label response = 'Mean Clinical Response'    label会体现在图的x、y轴label
   
4.       visit    = 'Visit';*/
   
5. input visitnum     /*访视次数，第0天、第1天……*/
   
6.       labtest :$20.     /*临床上关注的某个定量指标*/
   
7.       drug :$20.
   
8.       pre
   
9.       result;
   
10. cards;
    
11. 1 alat drug_a 0.27216 0.64893
    
12. 1 bilirubin_total drug_a 1.28536 0.54600
    
13. 1 alk drug_a 0.27216 0.64893
    
14. 1 asat drug_a 1.28536 0.54600
    
15. 2 alat drug_a 0.27216 0.64893
    
16. 2 bilirubin_total drug_a 1.28536 0.54600
    
17. 2 alk drug_a 0.27216 0.64893
    
18. 2 asat drug_a 1.28536 0.54600
    
19. 3 alat drug_a 0.27216 0.64893
    
20. 3 bilirubin_total drug_a 1.28536 0.54600
    
21. 3 alk drug_a 0.27216 0.64893
    
22. 3 asat drug_a 1.28536 0.54600
    
23. ;
    
24. run;
    
25. 
    
26. /*自定义格式，用于刻度值的自定义显示*/
    
27. proc format;
    
28.    value visitnumfmt
    
29.       1="1 week"
    
30.       2="3 months"
    
31.       3="6 months";
    
32. run;
    
33. 
    
34. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
35. ods listing gpath='C:\Users\Administrator\Desktop\test20180131';
    
36. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
37. /*ods graphics on / imagename="picture" outputfmt=pdf;     图片名称picture，输出类型pdf*/
    
38. ods graphics off;
    
39. 
    
40. /*Panel of LFT Results*/
    
41. ods listing style=statistical;
    
42. proc sgpanel data=labs;
    
43.    format visitnum visitnumfmt.;
    
44.    panelby labtest visitnum /
    
45.    onepanel novarname layout=lattice;
    
46.    scatter x=pre y=result/ group=drug;
    
47.    refline 1 1.5 2 / axis=Y
    
48.    lineattrs=(pattern=dash);
    
49.    refline 1 1.5 2 / axis=X
    
50.    lineattrs=(pattern=dash);
    
51.    rowaxis integer min=0 max=4
    
52.    label='Study (/ULN)';
    
53.    colaxis integer min=0 max=4
    
54.    label='Baseline (/ULN) *';
    
55.    keylegend / title=" " noborder;/*不显示图例的题目，图例无边框*/
    
56. run;
    
57. 
    
58. ods _all_ close;
    
59. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. /*建立数据集labs,用于后续画图*/
   
2. data labs;
   
3. /*label response = 'Mean Clinical Response'    label会体现在图的x、y轴label
   
4.       visit    = 'Visit';*/
   
5. input Patient
   
6.       alat
   
7.       biltot
   
8.       alkph
   
9.       asat
   
10.       relday  /*days*/
    
11.       miny    /*dval*/
    
12.       sday;   /*sdays*/
    
13. cards;
    
14. 1 0.50000 0.40000 0.40000 0.50000 -25 . .
    
15. 1 1.36533 0.64286 0.68972 1.08498 0 -0.5 0
    
16. 1 2.18967 0.87118 0.97332 1.64750 25 -0.5 25
    
17. 2 0.50000 0.40000 0.40000 0.50000 -25 . .
    
18. 2 1.36533 0.64286 0.68972 1.08498 0 -0.5 0
    
19. 2 2.18967 0.87118 0.97332 1.64750 25 -0.5 25
    
20. 3 0.50000 0.40000 0.40000 0.50000 -25 . .
    
21. 3 1.36533 0.64286 0.68972 1.08498 0 -0.5 0
    
22. 3 2.18967 0.87118 0.97332 1.64750 25 -0.5 25
    
23. 4 0.50000 0.40000 0.40000 0.50000 -25 . .
    
24. 4 1.36533 0.64286 0.68972 1.08498 0 -0.5 0
    
25. 4 2.18967 0.87118 0.97332 1.64750 25 -0.5 25
    
26. ;
    
27. run;
    
28. 
    
29. /*自定义格式，用于刻度值的自定义显示*/
    
30. proc format;
    
31.    value Patientfmt
    
32.       1="Patient 5152: White Male Age 48; Drug: A"
    
33.       2="Patient 6416: White Male Age 64; Drug: A"
    
34.       3="Patient 6950: White Male Age 51; Drug: A"
    
35.       4="Patient 6969: White Male Age 48; Drug: B";
    
36. run;
    
37. 
    
38. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
39. ods listing gpath='C:\Users\Administrator\Desktop\test20180131';
    
40. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
41. /*ods graphics on / imagename="picture" outputfmt=pdf;     图片名称picture，输出类型pdf*/
    
42. ods graphics off;
    
43. 
    
44. /*Panel of LFT Results for At-Risk Patients*/
    
45. ods listing style=statistical;
    
46. proc sgpanel data=labs cycleattrs;
    
47.    format Patient Patientfmt.;
    
48.    panelby patient / novarname;
    
49.    series x=relday y=alat / markers
    
50.    lineattrs=(thickness=2px
    
51.    pattern=solid);
    
52.    series x=relday y=asat / markers
    
53.    lineattrs=(thickness=2px
    
54.    pattern=solid);
    
55.    series x=relday y=alkph / markers
    
56.    lineattrs=(thickness=2px pattern=solid);
    
57.    series x=relday y=biltot / markers
    
58.    lineattrs=(thickness=2px pattern=solid);
    
59.    band x=sday lower=miny upper=4.5 /
    
60.    transparency=0.8
    
61.    legendlabel='Trial Duration';
    
62.    refline 1 1.5 2 / axis=Y
    
63.    lineattrs=(pattern=dash);
    
64.    colaxis min=-50 max= 200;
    
65.    rowaxis label="Upper Limit Normal";
    
66. run;
    
67. 
    
68. ods _all_ close;
    
69. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. proc format;
   
2. value smokfmt
   
3. 1-999 = "Smoker"
   
4. other = "Non-Smoker";
   
5. value diasfmt
   
6. 0-90 = "Diastolic < 90 mmHg"
   
7. 91-95 = "Diastolic 90 - 95 mmHg"
   
8. 96-999 = "Diastolic > 95 mmHg";
   
9. run;
   
10. 
    
11. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
12. ods listing gpath='C:\Users\Administrator\Desktop\test20180131';
    
13. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
14. /*ods graphics on / imagename="picture" outputfmt=pdf;     图片名称picture，输出类型pdf*/
    
15. ods graphics off;
    
16. 
    
17. /*Two Aggregating Formats Used in a LATTICE Layout*/
    
18. ods listing style=statistical;
    
19. proc sgpanel data=sashelp.heart;
    
20.    format smoking smokfmt.
    
21.    diastolic diasfmt.;
    
22.    panelby smoking diastolic /
    
23.    novarname missing onepanel
    
24.    layout=lattice;
    
25.    histogram weight;
    
26.    density weight;
    
27.    rowaxis integer grid;
    
28.    colaxis grid label='Weight (in pounds)';
    
29. run;
    
30. 
    
31. ods _all_ close;
    
32. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: GROUPED BAR CHARTS (Clinical Trial Reporting)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180201
   
6.    From: Clinical Trial Reporting Using SAS/GRAPH? SG Procedures
   
7.          http://support.sas.com/resources/papers/proceedings09/174-2009.pdf
   
8.    History: 20180201;
   
9. %mend;
   
10. 
    
11. /*建立数据集freq,用于后续画图*/
    
12. data freq;
    
13. /*label response = 'Mean Clinical Response'    label会体现在图的x、y轴label
    
14.       visit    = 'Visit';*/
    
15. input trt
    
16.       pain :$10.
    
17.       percent percent.;
    
18. cards;
    
19. 1 0 12%
    
20. 1 1-2 14%
    
21. 1 3-4 16%
    
22. 2 0 12%
    
23. 2 1-2 14%
    
24. 2 3-4 16%
    
25. ;
    
26. run;
    
27. 
    
28. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
29. ods listing gpath='C:\Users\Administrator\Desktop\test20180201';
    
30. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
31. /*ods graphics on / imagename="picture" outputfmt=pdf;     图片名称picture，输出类型pdf*/
    
32. ods graphics off;
    
33. 
    
34. /*Stacked Group Bar Chart*/
    
35. ods listing style=analysis;
    
36. proc sgplot data=freq;
    
37.    vbar trt / response=percent group=pain;
    
38.    xaxis display=(nolabel);
    
39.    yaxis
    
40.    label='Cumulative Percent of
    
41.    Patients';
    
42.    keylegend /position=bottom
    
43.    title="Pain Score";
    
44. run;
    
45. 
    
46. ods _all_ close;
    
47. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
   
2. ods listing gpath='C:\Users\Administrator\Desktop\test20180201';
   
3. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
   
4. /*ods graphics on / imagename="picture" outputfmt=pdf;     图片名称picture，输出类型pdf*/
   
5. ods graphics off;
   
6. 
   
7. /*Adjacent Group Bar Chart*/
   
8. ods listing style=analysis;
   
9. footnote bold 'Pain Score';
   
10. proc sgpanel data=freq;
    
11.    panelby pain / layout=columnlattice /*layout根据不同的列值*/
    
12.    onepanel noborder novarname  /*cell不加边框*/
    
13.    colheaderpos=bottom; /*列头位置在最下面*/
    
14.    vbar trt / response=percent group=trt; /*画垂直条图*/
    
15.    rowaxis grid offsetmax=0.05;
    
16.    colaxis display=none;
    
17.    keylegend /position=right across=1;
    
18. run;
    
19. 
    
20. ods _all_ close;
    
21. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码

1. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
   
2. ods listing gpath='C:\Users\Administrator\Desktop\test20180201';
   
3. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
   
4. /*ods graphics on / imagename="picture" outputfmt=pdf;     图片名称picture，输出类型pdf*/
   
5. ods graphics off;
   
6. 
   
7. /*LATTICED GROUP BAR CHART*/
   
8. ods listing style=analysis;
   
9. proc sgpanel data=freq noautolegend;
   
10. panelby trt / layout=columnlattice
    
11. novarname;
    
12. vbar pain / response=percent group=trt;
    
13. rowaxis grid;
    
14. run;
    
15. 
    
16. ods _all_ close;
    
17. ods html;      /*最好将ods destination恢复至默认的html*/

复制代码