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

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Time Axis Features
   
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. /*Time Axis Features*/
    
21. proc template;
    
22. define statgraph Fig_6_2_5;
    
23. begingraph;
    
24. entrytitle 'Monthly Stock Prices by Company';
    
25. layout overlay / xaxisopts=( display=(ticks tickvalues)
    
26.       timeopts=(interval=semimonth minorticks=true))
    
27.       y2axisopts=(label='Price' griddisplay=on);
    
28.    highlowplot x=date high=high low=low / close=close
    
29.       yaxis=y2 lineattrs=graphdata2(thickness=3) ;
    
30. endlayout;
    
31. endgraph;
    
32. end;
    
33. run;
    
34. ods listing style=htmlblue;
    
35. proc sgrender template=Fig_6_2_5
    
36. data=sashelp.stocks(where=(date > '01jan2004'd));
    
37. where stock eq 'IBM';
    
38. run;
    
39. 
    
40. ods _all_ close;
    
41. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Discrete Legend
   
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 trace on;
    
13. ods output Survivalplot=GTL_GS_SurvivalPlotData;
    
14. ods graphics / reset imagename='Survival-LifeTest';
    
15. proc lifetest data=sashelp.BMT plots=survival(atrisk=0 to 2500 by 500);
    
16.    time T * Status(0);
    
17.    strata Group / test=logrank adjust=sidak;
    
18.    run;
    
19. ods trace off;
    
20. 
    
21. proc format;
    
22.   value aml
    
23.   3 = 'AML-Low'
    
24.   2 = 'AML-High'
    
25.   1 = 'All';
    
26. run;
    
27. 
    
28. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
29. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
    
30. ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
    
31. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
32. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
33. /*图片名称picture，输出类型pdf*/
    
34. ods graphics off;
    
35. 
    
36. /*Discrete Legend*/
    
37. proc template;
    
38. define statgraph Fig_6_3_1;
    
39. begingraph;
    
40. entrytitle 'Product-Limit Survival Estimates';
    
41. legendItem type=marker name='c' / label='Censored'
    
42. markerattrs=(symbol=plus);
    
43. layout overlay;
    
44.    stepplot x=time y=survival / group=stratum
    
45.       lineattrs=(pattern=solid) name='s';
    
46.    scatterplot x=time y=censored / group=stratum
    
47.       markerattrs=(symbol=plus);
    
48.    innermargin;/*这个部分还需要进一步学习*/
    
49.       blockplot x=tatrisk block=atrisk / class=stratumnum
    
50.          display=(values label) valueattrs=(size=8)
    
51.          labelattrs=(size=8);
    
52.    endinnermargin;
    
53.    discretelegend 'c'/ location=inside halign=right valign=top;
    
54.    discretelegend 's' / valueattrs=(size=8);
    
55. endlayout;
    
56. endgraph;
    
57. end;
    
58. run;
    
59. proc sgrender data=GTL_GS_SurvivalPlotData template=Fig_6_3_1;
    
60. format stratumnum aml.;
    
61. run;
    
62. 
    
63. ods _all_ close;
    
64. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Continuous Legend
   
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. data GTL_GS_HeatmapParm;
    
13.   drop pi val;
    
14.   pi=constant('Pi');
    
15.   do x=1 to 20;
    
16.     do y=1 to 10;
    
17.           value=sin(x*pi/10) + cos(y*pi/10);
    
18.           val=ranuni(2);
    
19.       group=ifc(val < 0.5, 'M', 'F');
    
20.           output;
    
21.         end;
    
22.   end;
    
23. run;
    
24. /*proc print data=GTL_GS_HeatmapParm; 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. /*Continuous Legend*/
    
35. proc template;
    
36. define statgraph Fig_6_3_2;
    
37. begingraph;
    
38. entrytitle "Basic Heat Map";
    
39. layout overlay / xaxisopts=(display=(ticks tickvalues))
    
40.       yaxisopts=(display=(ticks tickvalues));
    
41.    heatmapparm x=x y=y colorresponse=value / name='a';
    
42.    continuouslegend 'a';
    
43. endlayout;
    
44. endgraph;
    
45. end;
    
46. run;
    
47. proc sgrender data=GTL_GS_HeatmapParm template=Fig_6_3_2;
    
48. run;
    
49. ods _all_ close;
    
50. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Titles and Footnotes
   
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. /*Titles and Footnotes*/
    
21. proc template;
    
22. define statgraph Fig_6_4_1;
    
23. begingraph;
    
24. entrytitle halign=left textattrs=(color=blue) 'Item 1:'
    
25.    halign=center textattrs=graphtitletext 'Mileage by Type';
    
26. entryfootnote halign=left 'Excluding Hybrids';
    
27. layout overlay / xaxisopts=(display=(ticks tickvalues))
    
28.       yaxisopts=(griddisplay=on);
    
29.    barchart x=type y=mpg_city / stat=mean fillattrs=graphdata3
    
30.       dataskin=gloss;
    
31. endlayout;
    
32. endgraph;
    
33. end;
    
34. run;
    
35. proc sgrender template=Fig_6_4_1
    
36. data=sashelp.cars(where=(type ne 'Hybrid'));
    
37. run;
    
38. 
    
39. ods _all_ close;
    
40. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Entries
   
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. /*Entries*/
    
21. proc template; /*--SAS 9.4--*/
    
22. define statgraph Fig_6_4_2;
    
23. begingraph / subpixel=on;
    
24. entrytitle 'Car Statistics';
    
25. layout overlay / yaxisopts=(griddisplay=on)
    
26.       xaxisopts=(griddisplay=on);
    
27.    scatterplot x=horsepower y=mpg_city / filledoutlinedmarkers=true
    
28.       dataskin=gloss markerattrs=(size=13 symbol=circlefilled)
    
29.       markerfillattrs=graphdata2;
    
30.    layout gridded / rows=2 order=columnmajor border=true
    
31.          autoalign=(topright);
    
32.       entry textattrs=(size=10) halign=right 'N = ';
    
33.       entry textattrs=(size=10) halign=right {unicode SIGMA}
    
34.          {sub '1'} ' =';
    
35.       entry textattrs=(size=10) halign=left
    
36.          eval(strip(put(n(mpg_city),4.0)));
    
37.       entry textattrs=(size=10) halign=left
    
38.          eval(strip(put(stddev(mpg_city),4.1)));
    
39.    endlayout;
    
40.    regressionplot x=horsepower y=mpg_city / degree=2
    
41.       lineattrs=(color=black thickness=4);
    
42. endlayout;
    
43. endgraph;
    
44. end;
    
45. run;
    
46. proc sgrender template=Fig_6_4_2 data=sashelp.cars; run;
    
47. 
    
48. ods _all_ close;
    
49. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

公卫的道友？学习了，很有用，一直想学sas绘图，回头好好看看分析资料

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Discrete Attributes Map – Line Attributes
   
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. /*--Fig_7_2_1 Group Series Discrete Attr Map--*/
    
13. data GTL_GS_Series_Dis_Attr_Map;
    
14. format Date Date9.;
    
15. do i=0 to 334 by 30;
    
16. date='01jan2009'd+i;
    
17. if mod (i, 30) =0 then freq=1; else freq=0;
    
18. Drug='Drug A'; Val = 16+ 3*sin(i/90+0.5) + 1*sin(3*i/90+0.7); output;
    
19. Drug='Drug C'; Val =  6+ 3*cos(i/90+0.5) + 1*sin(3*i/90+0.7); output;
    
20. end;
    
21. run;
    
22. 
    
23. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
24. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180308' image_dpi=300;*/
    
25. ods listing gpath='C:\Users\Administrator\Desktop\test20180308';
    
26. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
27. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
28. /*图片名称picture，输出类型pdf*/
    
29. ods graphics off;
    
30. 
    
31. /*Discrete Attributes Map – Line Attributes*/
    
32. proc template; /*--SAS 9.4--*/
    
33. define statgraph Fig_7_2_1;
    
34. dynamic _title;
    
35. begingraph / subpixel=on;/*subpixel=on使连接线更加光滑*/
    
36. discreteattrmap name='drug';
    
37.    value 'Drug A' / lineattrs=(color=red pattern=solid);
    
38.    value 'Drug B' / lineattrs=(color=green pattern=solid);
    
39.    value 'Drug C' / lineattrs=(color=blue pattern=solid);
    
40. enddiscreteattrmap;
    
41. discreteattrvar attrvar=drugname var=drug attrmap='drug';/*attrvar=drugname用于group*/
    
42. entrytitle 'Response by Time and Treatment';
    
43. layout overlay / yaxisopts=(griddisplay=on label='Response')
    
44.       xaxisopts=(griddisplay=on display=(ticks tickvalues));
    
45.    seriesplot x=date y=val / group=drugname curvelabel=drug name='a'
    
46.       smoothconnect=true lineattrs=(thickness=2);/*smoothconnect=true使连接线更加光滑*/
    
47.    discretelegend 'drug' / type=line;
    
48. endlayout;
    
49. endgraph;
    
50. end;
    
51. run;
    
52. 
    
53. proc sgrender data=GTL_GS_Series_Dis_Attr_Map template=Fig_7_2_1; run;
    
54. 
    
55. ods _all_ close;
    
56. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Discrete Attributes Map – Fill Attributes
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180308
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180308;
   
10. %mend;
    
11. 
    
12. /*--7_2_2 HIGHLOW Discrete Attr Map--*/
    
13. data GTL_GS_highlow_Dis_Attr_Map;
    
14. length cap $ 12;
    
15. input drug $ low high cap $;
    
16. datalines;
    
17. A      10 20 NONE
    
18. A      30 60 NONE
    
19. B      20 35 NONE
    
20. B      50 75 NONE
    
21. C      30 90 FILLEDARROW
    
22. ;
    
23. run;
    
24. 
    
25. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
26. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180308' image_dpi=300;*/
    
27. ods listing gpath='C:\Users\Administrator\Desktop\test20180308';
    
28. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
29. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
30. /*图片名称picture，输出类型pdf*/
    
31. ods graphics off;
    
32. 
    
33. /*Discrete Attributes Map – Fill Attributes*/
    
34. proc template;
    
35. define statgraph Fig_7_2_2;
    
36. dynamic _title;
    
37. begingraph;
    
38.    discreteattrmap name='drug';
    
39.       value 'A' / fillattrs=(color=lilac);
    
40.       value 'B' / fillattrs=(color=powderblue);
    
41.       value 'C' / fillattrs=(color=gold);
    
42.    enddiscreteattrmap;
    
43.    discreteattrvar attrvar=drugname var=drug attrmap='drug';
    
44.    entrytitle 'Response by Time and Treatment';
    
45.    layout overlay / yaxisopts=(griddisplay=on label='Response')
    
46.          xaxisopts=(griddisplay=on display=(ticks tickvalues));
    
47.       highlowplot y=drug high=high low=low / type=bar highcap=cap
    
48.          group=drugname outlineattrs=(color=black) name='a';
    
49.       discretelegend 'a';
    
50.    endlayout;
    
51. endgraph;
    
52. end;
    
53. run;
    
54. proc sgrender data=GTL_GS_highlow_Dis_Attr_Map template=Fig_7_2_2;
    
55. run;
    
56. 
    
57. ods _all_ close;
    
58. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Range Attributes Map – Discrete Colors
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180308
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180308;
   
10. %mend;
    
11. 
    
12. data GTL_GS_HeatmapParm;
    
13.   drop pi val;
    
14.   pi=constant('Pi');
    
15.   do x=1 to 20;
    
16.     do y=1 to 10;
    
17.           value=sin(x*pi/10) + cos(y*pi/10);
    
18.           val=ranuni(2);
    
19.       group=ifc(val < 0.5, 'M', 'F');
    
20.           output;
    
21.         end;
    
22.   end;
    
23. run;
    
24. /*proc print data=GTL_GS_HeatmapParm; run;*/
    
25. 
    
26. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
27. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180308' image_dpi=300;*/
    
28. ods listing gpath='C:\Users\Administrator\Desktop\test20180308';
    
29. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
30. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
31. /*图片名称picture，输出类型pdf*/
    
32. ods graphics off;
    
33. 
    
34. /*Range Attributes Map – Discrete Colors*/
    
35. proc template;
    
36. define statgraph Fig_7_3_1;
    
37. begingraph;
    
38.    rangeattrmap name='Resp';
    
39.       range -2 - -1 / rangecolor=darkred;
    
40.       range -1 - 0 / rangecolor=red;
    
41.       range 0 - 1 / rangecolor=white;
    
42.       range 1 - 2 / rangecolor=blue;
    
43.    endrangeattrmap;
    
44.    rangeattrvar attrvar=AttrValue var=value attrmap='Resp';
    
45.    entrytitle "Heat Map with Discrete Range Attribute Map";
    
46.    layout overlay / xaxisopts=(display=(ticks tickvalues))
    
47.          yaxisopts=(display=(ticks tickvalues));
    
48.       heatmapparm x=x y=y colorresponse=AttrValue / name='a';
    
49.       continuouslegend 'a';
    
50.    endlayout;
    
51. endgraph;
    
52. end;
    
53. run;
    
54. proc sgrender data=GTL_GS_HeatmapParm template=Fig_7_3_1;
    
55. run;
    
56. 
    
57. ods _all_ close;
    
58. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码

1. %macro comment();    /*这是一种比较特殊的注释方式*/
   
2.    Picture: Range Attributes Map – Gradient Colors
   
3.             (using template and sgrender procedure)
   
4.    Software： SAS 9.4 TS Level 1M2
   
5.    Name: justsoso
   
6.    Date: 20180308
   
7.    From: Getting Started with the Graph Template Language in SAS:
   
8.          Examples, Tips, and Techniques for Creating Custom Graphs
   
9.    History: 20180308;
   
10. %mend;
    
11. 
    
12. data GTL_GS_HeatmapParm;
    
13.   drop pi val;
    
14.   pi=constant('Pi');
    
15.   do x=1 to 20;
    
16.     do y=1 to 10;
    
17.           value=sin(x*pi/10) + cos(y*pi/10);
    
18.           val=ranuni(2);
    
19.       group=ifc(val < 0.5, 'M', 'F');
    
20.           output;
    
21.         end;
    
22.   end;
    
23. run;
    
24. /*proc print data=GTL_GS_HeatmapParm; run;*/
    
25. 
    
26. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
    
27. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180308' image_dpi=300;*/
    
28. ods listing gpath='C:\Users\Administrator\Desktop\test20180308';
    
29. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
    
30. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
    
31. /*图片名称picture，输出类型pdf*/
    
32. ods graphics off;
    
33. 
    
34. /*Range Attributes Map – Gradient Colors*/
    
35. proc template;
    
36. define statgraph Fig_7_3_2;
    
37. begingraph;
    
38.    rangeattrmap name='Resp';
    
39.       range -2 - -1 / rangecolormodel=(darkred red);
    
40.       range -1 - 0 / rangecolormodel=(red white);
    
41.       range 0 - 1 / rangecolormodel=(white blue);
    
42.       range 1 - 2 / rangecolormodel=(blue darkblue);
    
43.    endrangeattrmap;
    
44.    rangeattrvar attrvar=AttrValue var=value attrmap='Resp';
    
45.    entrytitle "Basic Heat Map";
    
46.    layout overlay / xaxisopts=(display=(ticks tickvalues))
    
47.          yaxisopts=(display=(ticks tickvalues));
    
48.       heatmapparm x=x y=y colorresponse=AttrValue / name='a';
    
49.       continuouslegend 'a';
    
50.    endlayout;
    
51. endgraph;
    
52. end;
    
53. run;
    
54. proc sgrender data=GTL_GS_HeatmapParm template=Fig_7_3_2;
    
55. run;
    
56. 
    
57. ods _all_ close;
    
58. /*ods html;      最好将ods destination恢复至默认的html*/

复制代码