统计分析软件教程(图表制作).ppt

图表制作与矩阵运算,Page2,第四章,内容与要求,掌握SAS系统的图表制作过程掌握SAS系统的矩阵运算过程,Page3,第四章,Page4,第四章,一统计图,1GPLOT过程GPLOT语句用于表现两个变量之间关系的图形。横、纵坐标分别代表两个变量，每个点对应数据集中的一个或几个观测。散点图语句格式：PROCGPLOT选择项；PLOT图形表达式/选择项；BY变量表；SYMBOLn选择项；,Page5,第四章,语句说明：(1)常用选择项：DATA=SAS数据集UNIFORM要求用BY变量分组打印的散点图的坐标刻度相同。(2)PLOT语句：通过图形表达式和选择项对所绘制的散点图提出具体要求。(3)BY语句：要求PLOT语句把观测按BY变量分组每组数据打印一个散点图，但BY变量必须事先排序。,Page6,第四章,图形表达式,纵坐标*横坐标：指定作为纵、横坐标的变量打印符号用系统设定的缺省值。纵坐标*横坐标打印点字符例如：“Y*X=#”表示纵坐标为变量Y横坐标为变量X在(X,Y)坐标点上打印。纵坐标*横坐标打印点变量例如：“Y*X=Z”当三个变量值分别为X=a,Y=b,Z=F时在图中坐标为(a,b)点上打印字符F。,例1,Page7,第四章,例1,PROCPLOTDATA=DST.BCLASS;PLOTHEIGHT*WEIGHT;RUN;,运行结果,Page8,第四章,返回,Page9,第四章,折线图如果需要生成折线图，则可由SYMBOL语句选项设定Interpol=连线方式i=连线方式。其中：NONE(不连线默认值)、JOIN(直线相连)、SPLINE(光滑曲线连接)、NEEDLE(每个点到横轴画垂直线)。CI=颜色：指明连线的颜色。其中：BLACK(黑)、RED(红)、GREEN(绿)、BLUE(蓝)、CYAN（青）LINE=nL=n：表示线型。0空白线，1实线（默认值），2虚线WIDTH=nW=n：线型的宽度。,例2,例3,例4,Page10,第四章,例2,数据集Data.Airqual记录了美国47个州逐月的两项空气质量指标。它包含如下几个变量。STATE：州名，用两个字符的简称REGION：各州所在的区域，以l一7表示MONTH：月份（l12)AVE_TSP：空气中悬浮颗粒含量的平均值AVE_LEAD：铅含量的平均值以纽约(STATE=NY)州的数据，制作逐月的平均悬浮颗粒含量的散点图。时间为横轴。由于只需要纽约州的数据，只需使用WHERE语句从数据集中取一个子集。并利用SYMBOL语句进行一定的修改。,Page11,第四章,编程,procgplotdata=data.airqual;symbolv=dotcv=redh=1.5;plotave_tsp*month;wherestate=NY;run;,Page12,第四章,运行结果,返回,Page13,第四章,例3,利用数据集Data.Airqual的数据用折线图表示纽约州逐月的平均悬浮颗粒含量。编程：procgplotdata=data.airqual;symboli=jionci=bluev=dotcv=redh=1.2pointlabel;/*选项pointlabel要求标出坐标值*/plotave_tsp*month;wherestate=NY;run;,Page14,第四章,运行结果,返回,Page15,第四章,例4,利用数据集Data.Airqual的数据用光滑曲线表示纽约州逐月的平均悬浮颗粒含量。编程：procgplotdata=data.airqual;symboli=splineci=bluev=pointcv=redw=2l=2pointlabel;plotave_tsp*month;wherestate=NY;run;,Page16,第四章,返回,运行结果,Page17,第四章,多重折线图分析过程中，为了进行对比分析需要一次绘制多条图形，这些图形可以使用BY语句进行分组绘制，也可以在一个GPLOT过程中使用多个POLT语句进行绘制。所生成的图形既可以用不同坐标系表达也可以在一个坐标系中绘制。OVERLAY选项：使用OVERLAY选项是将若干图形在一个坐标系内表现出来。,例5,例6,Page18,第四章,例5,数据集data.gmzb是1989年2003年我国宏观国民经济指标的资料，其中：TIME：年份（19892003）GDP：国内生产总值INC：人口自然增长率POP：年末人口总数AVE：人均国内生产总值要求在一个坐标系内绘制各年的GDP和AVE折线图，时间为横轴。比较两者发展变化趋势。,Page19,第四章,编程,procgplotdata=data.gmzb;symbol1i=splineci=bluev=squarecv=blackw=2l=2;symbol2i=splineci=blackv=diamondcv=redw=1l=1h=1.5;plotgdp*timeave*time/overlaylegend;run;由于需要绘制两条折线，因此用两个SYMBOL语句分别定义。利用OVERLAY选项是这两条折线出现在一个坐标系中。使用LEGEND选项加入图例以说明不同的折线。,Page20,第四章,返回,Page21,第四章,例6,从运行结果上看由于两个指标GDP和AVE的数值差距很大，指标AVE的变化反映的较为粗糙。为了解决这一问题，可以使用PLOT2语句。编程：procgplotdata=data.gmzb;symbol1i=splineci=bluev=squarecv=blackw=2l=2;symbol2i=splineci=blackv=diamondcv=redw=1l=1h=1.5;plotgdp*time/overlaylegend;plot2ave*time/overlaylegend;run;,Page22,第四章,返回,Page23,第四章,2GCHART过程语句格式：PROCGCHART选择项；VBAR绘图变量表/选择项；HBAR绘图变量表/选择项；PIE绘图变量表/选择项；BY变量列表；,Page24,第四章,语句说明：PROCGCHART语句：可跟任意多个绘图语句。常用选择项：DATA=SAS数据集：指出该过程使用的数据集，如省略，用最新建立的SAS数据集。BY语句:语句与PROCGCHART一起使用时，可得到由BY变量定义的几个观测组分别画出的图表。,Page25,第四章,绘图语句：VBAR绘图变量表/选择项：列出要画的垂直条形图的变量。每一个绘图变量对应一张条形图。HBAR绘图变量表/选择项：列出要画的水平条形图的变量。PIE语句：绘图变量表/选择项:绘制饼形图。其选择项与VBAR语句基本相同。,例7,Page26,第四章,例7,数据集data.jszl是某学校具有中级以上职称教师的基本资料，包含年龄(AGE)、性别(SEX)、职称(TITL)、收入(INCOME)等变量。利用该资料绘制垂直条形图、水平条形图和饼形图。,编程,Page27,第四章,编程,procgchartdata=data.jszl;vbartitl;hbartitl;pieage;run;,运行结果,运行结果,运行结果,Page28,第四章,返回,Page29,第四章,返回,Page30,第四章,返回,Page31,第四章,条形图的调整(1)中心值点(MIDPOINT)的调整：MIDPOINTS:中点值表。当绘图变量是连续变量时，用此选项指定变量分组的中点。中点值表可以是一系列中点值，也可以用“起点TO终点BY间隔”的方式来表示。LEVELS：如果绘图变量是连续的，用此选项指定产生N条图形。两者的结果是相同的。,Page32,第四章,(2)变量分组的调整:GROUP：分组变量，要求产生以分组变量值分组的并排图。SUBGROUP:细分变量，要求产生的条形图按细分变量的值再分组，分别表示每个条形图中细分变量的各值所占的比例。,例8,Page33,第四章,例8,利用选择项MIDPOINT、GROUP、SUBGROUP选项绘制水平条形图。其程序为：procgchartdata=data.jszl;hbartitl/group=sex;hbartitl/subgroup=sex;hbarage/midpoints=30to70by10;run;,运行结果,Page34,第四章,运行结果,采用选择项:GROUP,Page35,第四章,对职工按职称和性别分组输出条形图。HBARTITL/GROUP=SEX；表示为作变量TITL的水平条形图，在按变量SEX分组的基础上，按变量TITL对分组数据作图。运行结果：纵坐标为频数，横坐标表现为变量SEX和TITL的分组状况。,运行结果,Page36,第四章,采用选择项:SUBGROUP,运行结果,Page37,第四章,采用SUBGROUP选择项：HBARTITL/SUBGROUP=SEX；表示为作变量TITL水平条形图，按变量SEX对条形图进行分块。运行结果：纵坐标为频数，横坐标表现为变量TITL分组状况，条形图表现为SEX的分块状况。,运行结果,Page38,第四章,返回,采用选择项:MIDPOINTS,Page39,第四章,GCHART过程不仅可以生成平面的图形，还可以绘制三维立体图。其图形语句为：VBAR3D：垂直柱形图；HBAR3D：水平柱形图；BLOCK：平铺柱形图；PIE3D：三维饼形图。,例9,Page40,第四章,例9,利用数据集data.jszl绘制三维图形。编程：procgchartdata=data.jszl;hbar3dtitl;hbar3dtitl/subgroup=sex;pie3dtitl;pie3dtitl/subgroup=sex;blockage/subgroup=titlgroup=sexmidpoints=30to60by10;run;,运行结果,Page41,第四章,三维条形图,运行结果,Page42,第四章,细分的三维条形图,运行结果,Page43,第四章,三维饼形图,运行结果,Page44,第四章,细分的三维饼形图,运行结果,Page45,第四章,细分的三维柱状图,返回,Page46,第四章,二统计表的制作,统计表制作的TABULATE过程：语句格式：PROCTABULATE选择项；CLASS分类变量表；VAR分析变量表；TABLE表达式/选择项；FREQ变量名；WEIGHT变量名；,Page47,第四章,语句说明：PROCTABULATE：制图语句，要求定义两类变量，分类变量和分析变量，分类变量可以是离散或连续变量，分析变量只能是连续变量，常用选择项。（1）DATA=SAS数据集：定义TABULATE过程的数据集，如省略，则使用最新建立的数据集。（2）MISSING:规定TABULATE在分析中包括一个或多个分类变量的缺失值，否则，缺失值不参加运算。（3）FORMAT=格式名：定义表中每个格子的大小。默认值为12.2,Page48,第四章,（4）TABULATE过程的其他说明全类变量ALL：ALL变量代表只有一个值的特殊类变量，当与其他分类变量具有分层关系时，产生小计和总和。百分数PCTN和PCTSUM表达式:(PCTN或PCTSUM)=标识语句功能：分两种类型进行统计计算。PCTN：是频数N的百分数。由观测值分组中的频数与总频数相除。PCTSUM：是SUM的百分数。由分组中的分析变量值与分析变量值的和相除。每个百分数是按定义的分母来计算的。,Page49,第四章,CLASS语句：定义分类变量。VAR语句：定义分析变量。TABLE语句：使用TABULATE过程的基础，可由若干个TABLE语句产生所需要的统计表。注：1、TABLE语句的基本内容：表达式分类变量(来自CLASS语句)或ALL语句分析变量(来自VAR语句)描述统计量格式说明,Page50,第四章,注：2TABLE语句的操作符及作用：逗号“,”:定义表格的维数，分隔页维，行维及列维。星号“*”:把变量分层制表，表示变量之间的包含关系。括号“()”:指定操作顺序控制两个以上操作符的优先级。括号“”:指明分母。空格“”：表示并列关系，在行或列中并行列出。等号“=”:指定所需变量或统计量。,Page51,第四章,注：3表格维数及表达式TABLE语句由一到三维表达式组成，变量之间用逗号隔开，如果定义三维表达式，则最左边维定义为页，中间维定义为行，最右边维定义为列。,例10,Page52,第四章,例10,已知某单位的职工基本情况如表所示。NAMESEXTITLAGEZHOULIMP50LIPINGFAP47WANGMAP32FANGML26YUKEFAP48LIULINGML37ZHANGFL28ZHAOFAP36要求：按职称分类作统计表。,编程,Page53,第四章,计算均值、最大值、最小值的二维表,PROCTABULATE;CLASSTITL;VARAGE;TABLETITLALL,AGE*(NMAXMINMEAN);RUN;,运行结果,Page54,第四章,运行结果,运行结果,Page55,第四章,含ALL变量的二维表,PROCTABULATE;CLASSTITLSEX;TABLETITLALL,SEXALL;RUN;,运行结果,Page56,第四章,运行结果,第四章,运行结果,Page57,第四章,计算各栏目频率的二维表,PROCTABULATE;CLASSTITLSEX;TABLETITLALL,SEX*(NPCTN);RUN;,运行结果,Page58,第四章,运行结果,运行结果,Page59,第四章,列频率表和行频率表,/*列频率*/PROCTABULATE;CLASSTITLSEX;TABLETITL,SEX*(NPCTN)ALL;RUN;/*行频率*/PROCTABULATE;CLASSTITLSEX;TABLETITLALL,SEX*(NPCTN)ALL;RUN;,Page60,第四章,列频率表,Page61,第四章,行频率表,返回,Page62,第四章,三矩阵运算,交互式矩阵语言(InteractiveMatrixLanguage，IML)是SAS系统中专门用于矩阵运算编程的模块，提供了极其丰富的矩阵运算功能。在(IML)中矩阵的SAS表达式为：或者记为：A=a11a12a1n,a21a22a2n,am1am2amn;,Page63,第四章,1矩阵加法例：已知矩阵和。计算C=A+B，D=A+B+1。编程：prociml;a=12,34;b=56,78;c=a+b;d=a+b+1;printa,b,c,d;quit;,Page64,第四章,矩阵加法运行结果,返回,Page65,第四章,2矩阵减法例：已知矩阵和。计算C=A-B，D=A-B+1。编程：prociml;a=12,34;b=56,78;c=a-b;d=a-b+1;printc,d;quit;,Page66,第四章,矩阵减法运行结果,返回,Page67,第四章,3矩阵乘法：矩阵相乘的运算符及法则,Page68,第四章,数乘矩阵例：已知矩阵和。计算C=B*A，D=A*B，E=6#A,F=A#6。编程：prociml;a=123,456;b=6;c=a*b;d=b*a;e=6#a;f=a#6;printc,d,e,f;quit;,Page69,第四章,数乘矩阵运行结果,返回,Page70,第四章,矩阵相乘例：已知矩阵和。计算C=A*B。编程：prociml;a=123,456;b=78,910,1112;c=a*b;printa,b,c;quit;,Page71,第四章,矩阵乘法运行结果,返回,Page72,第四章,矩阵乘方例：已知矩阵，计算C=A*2。编程：prociml;a=13,57;resetprint;c=a*2;quit;,Page73,第四章,矩阵乘方运行结果,返回,Page74,第四章,矩阵的逆例：已知矩阵，计算其逆矩阵C=A*(-1)。编程：prociml;a=22,43;resetprint;c=a*(-1);d=c*a;quit;,Page75,第四章,矩阵的逆运行结果,返回,Page76,第四章,矩阵直接相乘例：已知矩阵和计算C=AB和D=BA。编程：prociml;a=12,45;b=78,910;resetprint;c=ab;d=ba;quit;,Page77,第四章,矩阵直接相乘运行结果,返回,Page78,第四章,对应元素相乘例：已知矩阵和计算C=A#B。编程：prociml;a=12,45;b=78,910;resetprint;c=a#b;quit;,Page79,第四章,对应元素相乘运算结果,返回,Page80,第四章,矩阵元素乘方例：已知矩阵和计算C=A#3,D=A#B。编程：prociml;a=12,45;b=22,43;resetprint;c=a#3;d=a#b;quit;,Page81,第四章,矩阵元素乘方运算结果,返回,Page82,第四章,4矩阵转置:例：已知矩阵，计算其转置C=A。编程：prociml;a=3811,648;resetprint;c=a;quit;,Page83,第四章,矩阵转置运算结果,返回,Page84,第四章,5线性方程的求解:一般线性方程组可以写成矩阵形式，然后利用矩阵运算求其解。在SAS/IML中可