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

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Distribution Plot (using template and sgrender procedure)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180229
   
6.    From: Getting Started with the Graph Template Language in SAS:
   
7.          Examples, Tips, and Techniques for Creating Custom Graphs
   
8.    History: 20180229;
   
9. %mend;
   
10. 
    
11. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
12. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180229' image_dpi=100;*/
    
13. ods listing gpath='C:\Users\Administrator\Desktop\test20180229';
    
14. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
15. /*ods graphics on / imagename="picture" outputfmt=pdf width=4in height=2.5in;*/
    
16. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
17. /*图片名称picture，输出类型pdf*/
    
18. ods graphics off;
    
19. 
    
20. /*Distribution Plot*/
    
21. proc template;
    
22. define statgraph Fig_4_7_2;
    
23. begingraph;
    
24. entrytitle 'Distribution of Mileage';
    
25. layout overlay;
    
26.    histogram mpg_city / binaxis=false;
    
27.    densityplot mpg_city / name='n' lineattrs=graphfit
    
28.       legendlabel='Normal';
    
29.    densityplot mpg_city / kernel() name='k' lineattrs=graphfit2
    
30.       legendlabel='Kernel';
    
31.    discretelegend 'n' 'k' / location=inside across=1
    
32.       halign=right valign=top;
    
33. endlayout;
    
34. endgraph;
    
35. end;
    
36. run;
    
37. proc sgrender data=sashelp.cars(where=(type ne 'Hybrid'))
    
38. template=Fig_4_7_2;
    
39. run;
    
40. 
    
41. ods _all_ close;
    
42. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Survival Plot (using template and sgrender procedure)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180229
   
6.    From: Getting Started with the Graph Template Language in SAS:
   
7.          Examples, Tips, and Techniques for Creating Custom Graphs
   
8.    History: 20180229;
   
9. %mend;
   
10. 
    
11. ods trace on;
    
12. ods output Survivalplot=GTL_GS_SurvivalPlotData;
    
13. ods graphics / reset imagename='Survival-LifeTest';
    
14. proc lifetest data=sashelp.BMT plots=survival(atrisk=0 to 2500 by 500);
    
15.    time T * Status(0);
    
16.    strata Group / test=logrank adjust=sidak;
    
17.    run;
    
18. ods trace off;
    
19. proc format;
    
20.   value aml
    
21.   3 = 'AML-Low'
    
22.   2 = 'AML-High'
    
23.   1 = 'All';
    
24. run;
    
25. 
    
26. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
27. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180229' image_dpi=100;*/
    
28. ods listing gpath='C:\Users\Administrator\Desktop\test20180229';
    
29. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
30. /*ods graphics on / imagename="picture" outputfmt=pdf width=4in height=2.5in;*/
    
31. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
32. /*图片名称picture，输出类型pdf*/
    
33. ods graphics off;
    
34. 
    
35. /*Survival Plot*/
    
36. proc template;
    
37. define statgraph Fig_4_7_3;
    
38. begingraph;
    
39. entrytitle 'Product-Limit Survival Estimates';
    
40. layout overlay;
    
41.    stepplot x=time y=survival / group=stratum
    
42.       lineattrs=(pattern=solid) name='s';
    
43.    scatterplot x=time y=censored /
    
44.       markerattrs=graphdatadefault(symbol=plus) name='c';
    
45.    scatterplot x=time y=censored / group=stratum
    
46.       markerattrs=(symbol=plus);
    
47.    innermargin;
    
48.    blockplot x=tatrisk block=atrisk / class=stratumnum
    
49.       display=(values label) valueattrs=(size=8)
    
50.       labelattrs=(size=8);
    
51.    endinnermargin;
    
52.    discretelegend 'c'/ location=inside halign=right valign=top;
    
53.    discretelegend 's' / valueattrs=(size=8);
    
54. endlayout;
    
55. endgraph;
    
56. end;
    
57. run;
    
58. proc sgrender data=GTL_GS_SurvivalPlotData template=Fig_4_7_3;
    
59. format stratumnum aml.;
    
60. run;
    
61. 
    
62. ods _all_ close;
    
63. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

楼主好厉害，一直在分享自己的知识，谢谢谢谢！

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Row and Column Data Lattice
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180305
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180305;
   
10. %mend;
    
11. 
    
12. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
13. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
14. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
15. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
16. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
17. /*图片名称picture，输出类型pdf*/
    
18. ods graphics off;
    
19. 
    
20. /*Row and Column Data Lattice*/
    
21. proc template;
    
22. define statgraph Fig_5_2_1;
    
23. begingraph;
    
24. entrytitle 'Mileage by Horsepower by Type and Origin';
    
25. layout datalattice rowvar=origin columnvar=type /
    
26.       headerlabeldisplay=value/*head显示*/
    
27.       rowaxisopts=(griddisplay=on)
    
28.       columnaxisopts=(griddisplay=on);
    
29.    layout prototype;
    
30.       scatterplot x=horsepower y=mpg_city / datatransparency=0.7
    
31.          markerattrs=(symbol=circlefilled size=11);
    
32.    endlayout;
    
33. endlayout;
    
34. endgraph;
    
35. end;
    
36. run;
    
37. proc sgrender template=Fig_5_2_1
    
38. data=sashelp.cars(where=(type not in ('Hybrid', 'Truck')));
    
39. run;
    
40. 
    
41. ods _all_ close;
    
42. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Column Data Lattice with Union-All Data Range
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180305
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180305;
   
10. %mend;
    
11. 
    
12. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
13. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
14. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
15. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
16. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
17. /*图片名称picture，输出类型pdf*/
    
18. ods graphics off;
    
19. 
    
20. /*Column Data Lattice with Union-All Data Range*/
    
21. proc template;
    
22. define statgraph Fig_5_2_2;
    
23. begingraph;
    
24. entrytitle 'Mileage by Type and Origin';
    
25. layout datalattice columnvar=origin / headerlabeldisplay=value
    
26.       rowaxisopts=(griddisplay=on offsetmin=0)
    
27.       columnaxisopts=(display=(ticks tickvalues)
    
28.       tickvalueattrs=(size=7));
    
29.    layout prototype / cycleattrs=true;
    
30.       barchart x=type y=mpg_city / stat=mean dataskin=gloss
    
31.          datatransparency=0.3 barlabel=true;
    
32.    endlayout;
    
33. endlayout;
    
34. endgraph;
    
35. end;
    
36. run;
    
37. proc sgrender template=Fig_5_2_2
    
38. data=sashelp.cars(where=(type not in ('Hybrid', 'Wagon')));
    
39. format mpg_city 4.1; run;
    
40. 
    
41. ods _all_ close;
    
42. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Column Data Lattice with Union Data Range
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180305
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180305;
   
10. %mend;
    
11. 
    
12. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
13. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
14. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
15. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
16. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
17. /*图片名称picture，输出类型pdf*/
    
18. ods graphics off;
    
19. 
    
20. /*Column Data Lattice with Union Data Range*/
    
21. proc template;
    
22. define statgraph Fig_5_2_3;
    
23. begingraph;
    
24. entrytitle 'Mileage by Type and Origin';
    
25. layout datalattice columnvar=origin / headerlabeldisplay=value
    
26.       columndatarange=union/*关键！*/
    
27.       rowaxisopts=(griddisplay=on offsetmin=0)
    
28.       columnaxisopts=(display=(ticks tickvalues)
    
29.       tickvalueattrs=(size=6));
    
30.    layout prototype;
    
31.       barchart x=type y=mpg_city / stat=mean dataskin=gloss
    
32.          datatransparency=0.3 barlabel=true;
    
33.    endlayout;
    
34. endlayout;
    
35. endgraph;
    
36. end;
    
37. run;
    
38. proc sgrender template=Fig_5_2_3
    
39. data=sashelp.cars(where=(type not in ('Hybrid', 'Wagon')));
    
40. format mpg_city 4.1;
    
41. run;
    
42. 
    
43. ods _all_ close;
    
44. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Proportional Data Lattice
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180305
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180305;
   
10. %mend;
    
11. 
    
12. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
13. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
14. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
15. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
16. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
17. /*图片名称picture，输出类型pdf*/
    
18. ods graphics off;
    
19. 
    
20. /*Proportional Data Lattice*/
    
21. proc template; /*--SAS 9.4--*/
    
22. define statgraph Fig_5_2_4;
    
23. begingraph;
    
24. entrytitle 'Mileage by Type and Origin';
    
25. layout datalattice columnvar=origin / headerlabeldisplay=value
    
26.       columnweight=proportional/*关键！*/
    
27.       rowaxisopts=(griddisplay=on offsetmin=0)
    
28.       columnaxisopts=(display=(ticks tickvalues)
    
29.       tickvalueattrs=(size=7));
    
30.    layout prototype;
    
31.       barchart x=type y=mpg_city / stat=mean dataskin=gloss
    
32.          barlabel=true;
    
33.    endlayout;
    
34. endlayout;
    
35. endgraph;
    
36. end;
    
37. run;
    
38. proc sgrender template=Fig_5_2_4
    
39. data=sashelp.cars(where=(type not in ('Hybrid', 'Wagon')));
    
40. format mpg_city 4.1; run;
    
41. 
    
42. ods _all_ close;
    
43. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Data Lattice with Side Bar
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180305
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180305;
   
10. %mend;
    
11. 
    
12. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
13. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
14. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
15. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
16. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
17. /*图片名称picture，输出类型pdf*/
    
18. ods graphics off;
    
19. 
    
20. /*Data Lattice with Side Bar*/
    
21. proc template;
    
22. define statgraph Fig_5_2_5;
    
23. begingraph;
    
24. entrytitle 'Mileage by Type and Origin';
    
25. layout datalattice columnvar=origin / headerlabeldisplay=value
    
26.       rowaxisopts=(label='Mileage' griddisplay=on offsetmin=0)
    
27.       columnaxisopts=(display=(ticks tickvalues)
    
28.       tickvalueattrs=(size=7));
    
29.    layout prototype / cycleattrs=true;
    
30.       barchart x=type y=mpg_city / stat=mean dataskin=gloss
    
31.          datatransparency=0.3 name='c' legendlabel='City';
    
32.       barchart x=type y=mpg_highway / stat=mean dataskin=gloss
    
33.          barwidth=0.5 datatransparency=0.3 name='h'
    
34.          legendlabel='Highway';
    
35.    endlayout;
    
36.    sidebar / spacefill=false;/* if true，框框会变得很长！*/
    
37.       discretelegend 'c' 'h';
    
38.    endsidebar;
    
39. endlayout;
    
40. endgraph;
    
41. end;
    
42. run;
    
43. proc sgrender template=Fig_5_2_5
    
44. data=sashelp.cars(where=(type not in ('Hybrid', 'Truck', 'Wagon')));
    
45. run;
    
46. 
    
47. ods _all_ close;
    
48. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Row Data Lattice (using template and sgrender procedure)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180305
   
6.    From: Getting Started with the Graph Template Language in SAS:
   
7.          Examples, Tips, and Techniques for Creating Custom Graphs
   
8.    History: 20180305;
   
9. %mend;
   
10. 
    
11. /*建立数据集GTL_GS_Labs,用于后续画图*/
    
12. /*--Fig 5.2.6 Labs--*/
    
13. data GTL_GS_Labs;
    
14.   pi=constant('PI');
    
15.   do Lab='Lab-1', 'Lab-2', 'Lab-3';
    
16.     amp=100*ranuni(3);
    
17.         phase=pi*ranuni(2);
    
18.     do day=1 to 300 by 5;
    
19.           Value=amp*(1+0.1*(sin(day*pi/180+phase)+0.2*ranuni(2)));
    
20.           output;
    
21.         end;
    
22.   end;
    
23. run;
    
24. /*proc print;run;*/
    
25. 
    
26. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
27. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
28. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
29. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
30. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
31. /*图片名称picture，输出类型pdf*/
    
32. ods graphics off;
    
33. 
    
34. /*Row Data Lattice*/
    
35. proc template;
    
36. define statgraph Fig_5_2_6;
    
37. begingraph;
    
38. entrytitle 'Lab Value by Day';
    
39. layout datalattice rowvar=lab / headerlabeldisplay=value/*默认为lab=...*/
    
40.       rowaxisopts=(label='Value' griddisplay=on
    
41.       tickvalueattrs=(size=7))
    
42.       columnaxisopts=(display=(ticks tickvalues)
    
43.       tickvalueattrs=(size=7))
    
44.       rowdatarange=union rowgutter=5;/*根据单个cell，rowgutter单元间隔，5应该是px*/
    
45.    layout prototype / cycleattrs=true;/*cycleattrs指定精美样式*/
    
46.       seriesplot x=day y=value / lineattrs=(thickness=2);
    
47.    endlayout;
    
48. endlayout;
    
49. endgraph;
    
50. end;
    
51. run;
    
52. proc sgrender template=Fig_5_2_6 data=GTL_GS_Labs;
    
53. run;
    
54. 
    
55. ods _all_ close;
    
56. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Data Lattice with Multi-Page Output (using template and sgrender procedure)
   
3.    Software： SAS 9.4 TS Level 1M2
   
4.    Name: justsoso
   
5.    Date: 20180305
   
6.    From: Getting Started with the Graph Template Language in SAS:
   
7.          Examples, Tips, and Techniques for Creating Custom Graphs
   
8.    History: 20180305;
   
9. %mend;
   
10. 
    
11. /*建立数据集GTL_GS_Labs,用于后续画图*/
    
12. /*--Fig 5.2.7 Labs2--*/
    
13. data GTL_GS_Labs2;
    
14.   pi=constant('PI');
    
15.   do Lab='Lab-1', 'Lab-2', 'Lab-3', 'Lab-4', 'Lab-5', 'Lab-6';
    
16.     amp=100*ranuni(4);
    
17.         phase=2*pi*ranuni(3);
    
18.     do day=1 to 300 by 5;
    
19.           Value=amp*(1+0.1*(sin(day*pi/180+phase)+0.2*ranuni(5)));
    
20.           output;
    
21.         end;
    
22.   end;
    
23. run;
    
24. 
    
25. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
26. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
27. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
28. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
29. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
30. /*图片名称picture，输出类型pdf*/
    
31. ods graphics off;
    
32. 
    
33. /*Data Lattice with Multi-Page Output*/
    
34. proc template;
    
35. define statgraph Fig_5_2_7;
    
36. begingraph;
    
37. entrytitle 'Lab Value by Day';
    
38. layout datalattice rowvar=lab / headerlabeldisplay=value
    
39.       rowaxisopts=(label='Value' griddisplay=on
    
40.       tickvalueattrs=(size=7))
    
41.       columnaxisopts=(display=(ticks tickvalues)
    
42.       tickvalueattrs=(size=7))
    
43.       rowdatarange=union rows=3 rowgutter=5;/*rows=3限制单元高度*/
    
44.    layout prototype / cycleattrs=true;
    
45.       seriesplot x=day y=value /
    
46.       lineattrs=graphdata2(thickness=2 pattern=solid);
    
47.    endlayout;
    
48. endlayout;
    
49. endgraph;
    
50. end;
    
51. run;
    
52. proc sgrender template=Fig_5_2_7 data=GTL_GS_Labs2;
    
53. run;
    
54. 
    
55. ods _all_ close;
    
56. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码