SAS系统和数据分析三维图形.docVIP

三维图形 SAS系统除了可以绘制二维平面图形外，还可以绘制三维立体图形。使用PROC G3D过程能对在PLOT或SCATTER语句中指定的变量值绘制三维图形。PROC G3D过程产生的三维图形分为两种：三维曲面图和三维散布图。三维曲面图是通过在PLOT语句中指定三个图形变量来产生，三维散布图是在SCATTER语句中必需指定三个图形变量才能产生，同样还可以对三维散布图进行大小、形状、颜色的修改。 PROC G3D过程说明 PROC G3D过程产生三维高分辨率立体图形，一般由下列语句控制： Proc G3D DATA=数据集 ; PLOT Y*X=Z /选项; SCATTER Y*X=Z /选项; TITLEn ‘字符串’; FOOTNOTEn ‘字符串’; By 变量列表 ; Run ; 在PROC G3D过程中至少要有一条PLOT或SCATTER语句，如果两条语句同时存在，将在一页中产生一个三维曲面图，下一页中产生一个三维散布图。语句中的Y*X=Z指定三个图形变量X、Y、Z，三维空间中的图形点由三个变量的值共同决定，X和Y为水平变量，Z为垂直变量。PLOT语句是将X、Y、Z三个变量值确定的空间图形点相互连成曲面，而SCATTER语句是将X、Y、Z三个变量值确定的空间图形点用符号表示或画出图形点向XY平面的垂直线（采用NEEDLE选项）。一个PROC G3D过程中可以有多条SCATTER语句。 PLOT语句的选项 用在PLOT语句中的选项按功能可以分为三类：外观选项、轴选项和描述性选项。主要选项如下： ROTATE=角度列表——指定三维图形Z轴的一个或多个旋转角度，缺省值为70度。如果有多个旋转角度，将画出对应多个不同旋转角度的三维图形。 TILT=角度列表——指定三维图形Y轴的一个或多个倾斜角度，缺省值为70度。如果有多个倾斜角度，将画出对应多个不同倾斜角度的三维图形。 GRID——在所有坐标轴的每一个刻度线上画出网格线。 SIDE——在三维曲面图形中画出侧面墙。 XTICKNUM=n YTICKNUM=n ZTICKNUM=n——指定X、Y、Z坐标轴上的刻度线数目n，缺省值为4。 ZMAX=最大数字值 ZMIN=最小数字值——设定Z轴上最大的数字值和最小的数字值，超出了最大值和最小值的Z变量值部分将被省略。 还有设置图形各个对象颜色的选项，如CBOTTOM、CTOP、CTEXT、CAXIS选项，分别指定顶部、底部、文本和所有轴的颜色。其他还有NOAXES选项表示不画出坐标轴，NOLABEL选项表示不显示轴的刻度值和轴标签。 SCATTER语句的选项 SCATTER语句包含上面PLOT语句的这些选项，还有自己的一些常用选项： NONEEDLE——不画出图形点向XY平面的垂直线。 SHAPE=‘符号名’——指定符号名或字符变量名。符号用于显示散布图的图形点的形状，有BALLOON 、DIAMOND 、PRISM 、CLUB 、FLAG 、PYRAMID、CROSS 、HEART、SPADE 、CUBE 、PILLAR 、SQUARE、CYLINDER 、POINT 、STAR。 SIZE=大小值——指定图形点符号的大小值。 COLOR=颜色名——指定图形点符号的颜色。 曲面图 例如，在前面我们用GCONTOUR过程绘制了表达式Z=2+7X-X2+9Y-Y2二维平面的利润等值线图，现在，我们可以用PROC G3D过程的PLOT语句画出三维立体的利润曲面图。我们仍然采用相同的数据集STUDY.PROFIT来绘制。 程序如下： Goptions reset=global gunit=pct cback=white border htitle=6 htext=3 ftext=swissb colors=(back) ; Libname Study d:\sasdata\mydir ; Proc g3d data=study.profit ; Plot y*x=z / rotate=160 tilt=80 ctop=black caxis=black cbottom=black grid zmax=34 zmin=-16 xticknum=10 yticknum=10 zticknum=11; Title Profit Surface ; Run ; 二维等利润线所构成的一圈圈大小值不同的同心椭圆，是由三维利润曲面中相同的Z值点在XY平面投影所构成的。图形中显示的曲面形状与公式的吻合程度， 与STUDY.PROFIT数据集中的观测条数有关，如果我们在STUDY.PROFIT数据集中生成的X、Y、Z数据点数越多，吻合