楼主: l1i2n3i4n5g
48359 125

[学习分享] 如何用SAS画统计图，这是我的学习分享！   [推广有奖]

31楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-1 15:24:55
Vector Plot (using template and sgrender procedure)
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Vector 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;

  11. /*建立数据集GTL_GS_vector,用于后续画图*/
  12. data GTL_GS_vector;
  13. format label 3.0;
  14. do i=1 to 10;
  15.    x=ranuni(2); xo=x*(1+(ranuni(2)-0.5));
  16.    y=ranuni(2); yo=y*(1+(ranuni(2)-0.5));
  17.    if (ranuni(2) < 0.5) then Type='A';
  18.    else Type='B';
  19.    label=(2*x+3*y);
  20.    output;
  21. end;
  22. run;

  24. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
  25. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180229' image_dpi=100;*/
  26. ods listing gpath='C:\Users\Administrator\Desktop\test20180229';
  27. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
  28. /*ods graphics on / imagename="picture" outputfmt=pdf width=4in height=2.5in;*/
  29. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
  30. /*图片名称picture，输出类型pdf*/
  31. ods graphics off;

  33. /*Vector Plot*/
  34. proc template;
  35. define statgraph Fig_4_2_9;
  36. begingraph;
  37. entrytitle 'Random Vectors';
  38. layout overlay;
  39.    vectorplot x=x y=y xorigin=xo yorigin=yo / group=type name='a'
  40.       arrowheadshape=barbed datalabel=label
  41.       lineattrs=(thickness=2);
  42.    discretelegend 'a' / sortorder=ASCENDINGFORMATTED
  43.       title='Type:' location=inside halign=right valign=top;
  44. endlayout;
  45. endgraph;
  46. end;
  47. run;
  48. proc sgrender data=GTL_GS_vector template=Fig_4_2_9;
  49. run;

  51. ods _all_ close;
  52. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
32楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-1 15:37:45
Bar Chart
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Bar Chart (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;

  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;

  20. /*Bar Chart*/
  21. proc template;
  22. define statgraph Fig_4_3_1;
  23. begingraph;
  24. entrytitle 'Mileage by Origin and Type';
  25. layout overlay / xaxisopts=(display=(ticks tickvalues))
  26.    yaxisopts=(griddisplay=on);
  27.       barchart x=origin y=mpg_city / group=type stat=mean
  28.       groupdisplay=cluster dataskin=gloss name='a';
  29.    discretelegend 'a' / title='Type:';
  30. endlayout;
  31. endgraph;
  32. end;
  33. run;
  34. proc sgrender data=sashelp.cars(where=(type not in ('Hybrid' 'Truck')))
  35. template=Fig_4_3_1;
  36. run;

  38. ods _all_ close;
  39. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
33楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-1 16:03:12
Axis Table
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Axis Table (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;

  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;

  20. /*Axis Table*/
  21. proc template; /*--SAS 9.4--*/
  22. define statgraph Fig_4_3_2;
  23. begingraph;
  24. entrytitle 'Mileage by Origin and Type';
  25. layout overlay / xaxisopts=(griddisplay=on)
  26.    yaxisopts=(display=(ticks tickvalues)
  27.       discreteopts=(colorbands=even));
  28.    barchart category=type response=mpg_city / orient=horizontal
  29.       groupdisplay=cluster stat=mean dataskin=gloss;
  30.    innermargin / align=right;
  31.    axistable y=type value=mpg_city / class=origin display=(label)
  32.       labelposition=max stat=mean;
  33.    endinnermargin;
  34. endlayout;
  35. endgraph;
  36. end;
  37. run;
  38. proc sgrender data=sashelp.cars template=Fig_4_3_2;
  39. format mpg_city 4.1; run;

  41. ods _all_ close;
  42. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
34楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-1 16:12:31
Line Chart
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Line Chart (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;

  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;

  20. /*Line Chart*/
  21. proc template; /*--SAS 9.4--*/
  22. define statgraph Fig_4_3_3;
  23. begingraph;
  24. entrytitle 'Mileage by Origin and Type';
  25. layout overlay / xaxisopts=(display=(ticks tickvalues)) yaxisopts=(griddisplay=on);
  26.    linechart category=type response=mpg_city / group=origin   stat=mean
  27.       lineattrs=(thickness=5 pattern=solid) datatransparency=0.4   name='a';
  28.    discretelegend 'a' / title='Origin:';
  29. endlayout;
  30. endgraph;
  31. end;
  32. run;
  33. proc sgrender data=sashelp.cars template=Fig_4_3_3;
  34. run;

  36. ods _all_ close;
  37. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
35楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-2 11:15:25
Pie Chart
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Pie Chart (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;

  11. /*建立数据集GTL_GS_Sedans,用于后续画图*/
  12. proc sort data=sashelp.cars out=CarsCyOrigin;
  13.   by origin;
  14.   run;
  15. data GTL_GS_Sedans;
  16.   retain Sedans 0 Rest 0;
  17.   format Pct percent.;
  18.   keep Origin Type Count Pct;
  19.   set CarsCyOrigin;
  20.   by origin;
  21.   if first.origin then do; Sedans=0; Rest=0; end;
  22.   if type eq 'Sedan' then Sedans+1;
  23.   else  Rest+1;
  24.   if last.origin then do;
  25.    Type='Sedans'; Count=Sedans; Pct=count/(sedans+rest); output;
  26.         Type='Rest';   Count=Rest;   Pct=count/(sedans+rest); output;
  27.   end;
  28. run;

  30. ods _all_ close;     /*关闭所有ods destinations，避免不必要的输出，节省计算机资源*/
  31. /*ods listing gpath='C:\Users\Administrator\Desktop\test20180229' image_dpi=100;*/
  32. ods listing gpath='C:\Users\Administrator\Desktop\test20180229';
  33. /*打开ods listing这一destination，指定图片保存路径，★★★★★需提前建立*/
  34. /*ods graphics on / imagename="picture" outputfmt=pdf width=4in height=2.5in;*/
  35. /*ods graphics on / imagename="picture" outputfmt=pdf;*/
  36. /*图片名称picture，输出类型pdf*/
  37. ods graphics off;

  39. /*Pie Chart*/
  40. proc template;
  41. define statgraph Fig_4_3_4;
  42. begingraph;
  43. entrytitle 'Share of Sedans among Cars in USA';
  44. layout region;/*饼图要用layout region！*/
  45.    piechart category=type response=pct / dataskin=sheen
  46.       datalabelcontent=(category response) name='a' start=340;
  47. endlayout;
  48. endgraph;
  49. end;
  50. run;
  51. proc sgrender data=GTL_GS_Sedans(where=(origin='USA')) template=Fig_4_3_4;
  52. run;

  54. ods _all_ close;
  55. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
36楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-2 11:18:22
Histogram
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Histogram (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;

  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;

  20. /*Histogram*/
  21. proc template;
  22. define statgraph Fig_4_4_1;
  23. begingraph;
  24. entrytitle 'Distribution of Mileage';
  25. layout overlay;
  26.    histogram mpg_city;
  27. endlayout;
  28. endgraph;
  29. end;
  30. run;
  31. proc sgrender data=sashelp.cars(where=(type ne 'Hybrid'))
  32. template=Fig_4_4_1;
  33. run;

  35. ods _all_ close;
  36. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
37楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-2 11:20:15
Density Plots
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Density Plots (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;

  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;

  20. /*Density Plots*/
  21. proc template;
  22. define statgraph Fig_4_4_2;
  23. begingraph;
  24. entrytitle 'Distribution of Mileage';
  25. layout overlay;
  26.    densityplot mpg_city / name='n' lineattrs=graphfit
  27.       legendlabel='Normal';
  28.    densityplot mpg_city / kernel() name='k' lineattrs=graphfit2
  29.       legendlabel='Kernel';
  30.    discretelegend 'n' 'k' / location=inside across=1
  31.       halign=right valign=top;
  32. endlayout;
  33. endgraph;
  34. end;
  35. run;
  36. proc sgrender data=sashelp.cars(where=(type ne 'Hybrid'))
  37. template=Fig_4_4_2;
  38. run;

  40. ods _all_ close;
  41. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
38楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-2 11:22:10
Box Plot – Horizontal - Discrete
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Box Plot – Horizontal - Discrete (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;

  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;

  20. /*Box Plot – Horizontal - Discrete*/
  21. proc template;
  22. define statgraph Fig_4_4_3;
  23. begingraph;
  24. entrytitle 'Cholesterol by Death Cause';
  25. layout overlay / yaxisopts=(display=(ticks tickvalues))
  26.       xaxisopts=(griddisplay=on);
  27.    sboxplot y=cholesterol x=deathcause / orient=horizontal;
  28. endlayout;
  29. endgraph;
  30. end;
  31. run;
  32. proc sgrender data=sashelp.heart template=Fig_4_4_3;
  33. run;

  35. ods _all_ close;
  36. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
39楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-2 11:27:43
Box Plot – Vertical – Interval with Groups
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Box Plot – Vertical – Interval with Groups (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;

  11. /*建立数据集GTL_GS_IntervalBoxGroup,用于后续画图*/
  12. data GTL_GS_IntervalBoxGroup;
  13.   format Date date6.;
  14.   drop i;
  15.   do Date='01Jan2009'd, '15Jan2009'd,'15Mar2009'd,
  16.           '01May2009'd, '01Aug2009'd;
  17.     do i=1 to 10;
  18.           Response=rannorm(2); Drug='A'; output;
  19.           Response=ranuni(2);  Drug='A'; output;
  20.           Response=rannorm(2); Drug='B'; output;
  21.           Response=ranuni(2);  Drug='B'; output;
  22.         end;
  23.   end;
  24. run;

  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;

  35. /*Box Plot – Vertical – Interval with Groups*/
  36. proc template; /*--SAS 9.4--*/
  37. define statgraph Fig_4_4_4;
  38. begingraph / attrpriority=color;
  39. entrytitle 'Response by Time and Treatment';
  40. layout overlay / yaxisopts=(griddisplay=on)
  41.    xaxisopts=(type=time timeopts=(interval=month)
  42.       display=(ticks tickvalues));
  43.    boxplot x=date y=response / group=drug groupdisplay=cluster
  44.       name='a';
  45.    discretelegend 'a' / title='Treatment:' location=inside
  46.       halign=right valign=bottom;
  47. endlayout;
  48. endgraph;
  49. end;
  50. run;
  51. proc sgrender data=GTL_GS_IntervalBoxGroup template=Fig_4_4_4;
  52. run;

  54. ods _all_ close;
  55. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
40楼
l1i2n3i4n5g(未真实交易用户) 在职认证  发表于 2018-3-2 11:42:15
Regression Plot
  1. %macro comment();    /*这是一种比较特殊的注释方式*/
  2.    Picture: Regression 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;

  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;

  20. /*Regression Plot*/
  21. proc template;
  22. define statgraph Fig_4_5_1;
  23. begingraph;
  24. entrytitle 'Fit Plot for Weight by Height';
  25. layout overlay / yaxisopts=(griddisplay=on)
  26.    xaxisopts=(griddisplay=on);
  27.    modelband 'Reg' / name='band' legendlabel='95% Confidence';
  28.    scatterplot x=height y=weight;
  29.    regressionplot x=height y=weight / cli='Reg' name='Reg';
  30.    discretelegend 'band' / location=inside halign=right
  31.       valign=bottom;
  32. endlayout;
  33. endgraph;
  34. end;
  35. run;
  36. proc sgrender data=sashelp.class template=Fig_4_5_1; run;

  38. ods _all_ close;
  39. /*ods html;      最好将ods destination恢复至默认的html*/
复制代码
返回列表
12345678910... 13下一页
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002号-2 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2025-12-6 06:16