SAS学习笔记（9）：SAS绘图

本篇SAS读书笔记主要介绍SAS绘制统计图形，主要介绍Gplot语句和Ggraph语句。
1. 散点图和曲线图GPLOT过程

例1：绘制学生信息数据的散点图与曲线图。

proc gplot data=mylib.stu;

plot math*num[选项];

symbol color=red i=jion v=none;     /*绘图的颜色，是否连线，符号等属性*/

run;

symbol语句的选项

2. 直方图和扇形图gchart过程

例2：绘制学生数学成绩直方图

proc gchart data=mylib.stu;

     vbar math;    /*说明分析变量用vbar*/

run;

如果把VBAR改成HBAR则条形方向变为横向。用GCHART绘制的直方图和在INSIGHT中绘制的直方图有所不同，它在横轴标的是区间的中点值，而在INSIGHT中横轴标的是区间的端点值。

例3：分别绘制男女学生数学成绩直方图

proc gchart data=mylib.stu;

     vbar math /group=sex;    

run;

例4：利用pie绘制表示频数的扇形图

proc gchart data=mylib.stu;

     pie sex /type=percent;    /*后面的type属性表示以百分数显示*/

run;

例5：用block语句绘制三维直方图

proc gchart data=mylib.stu;

   block math /group=sex;

run;

GCHART过程可以使用的图形关键字及其所绘制的图形见下：

3. 三维曲面和等高线图

  假设对一个二元函数z=f(x，y)，我们有了x取等间隔值、y取等间隔值时z的值，这时我们可以用G3D过程绘制曲面图形，用GCONTOUR绘制曲面的等高线图。

例6：绘制一个二维正态分布曲面的图形，假设(X，Y)服从联合正态分布，其均值都是0，方差都为1，相关系数为r。

首先建立数据集：

data dnorm2;

r=0.5;

det=1-r*r;

do x=-3 to 3 by 0.3;

   do y=-3 to 3 by 0.3;

       z=1/(2*3.1415926*det)*exp(-0.5/det*(x*x+y*y-2*r*x*y));

         output;

      end;

end;

keep x y z;

run;

然后，我们就可以用G3D过程来绘制曲面图

proc g3d data=dnorm2;

   plot x*y=z;

   run;

用GCONTOUR过程可以绘制曲面对应的等高线图

proc gcontour data=dnorm2;

   plot x*y=z /nolegend autolabel;

   run;