统计分析软件教程(图表制作).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 n L n 表示线型 0 空白线 1 实线 默认值 2 虚线WIDTH n W n 线型的宽度 例2 例3 例4 Page10 第四章 例2 数据集Data Airqual记录了美国47个州逐月的两项空气质量指标 它包含如下几个变量 STATE 州名 用两个字符的简称REGION 各州所在的区域 以l一7表示MONTH 月份 l 12 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 年份 1989 2003 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 a11a12 a1n a21a22 a2n am1am2 amn 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 A B和D B A 编程 prociml a 12 45 b 78 910 resetprint c a b d b a 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线性方程的求解 一般线