1、SAS 基础,SAS 功能及历史（略） INSIGHT，Analyst，Program 是统计分析的三种方法 逻辑库与 SAS 文件（主要包括数据集和程序，访问形式：逻辑库名.文件名，不加逻辑库名，视为Work） SAS 中变量的名字（最多8个字符，不区分大小写）,数据步和程序步（可能还会包括一些全程语句） SAS 语句书写规则（关键词开始，分号“;”结束）,SAS 编程初步,data exam; /* 考察某些人考试成绩的数据集*/input name $ sex $ math computer;cards; 阿基米德 男 92 98 亚里士多德 男 89 95 居里夫人女 86 90 Ga
2、ga 女 25 82 高斯 男 98 83 ; *以下为过程步 proc print data = exam; proc means data =exam mean;var math computer; run;,指定逻辑库,SAS 编程初步,libname (“” “” );例：libname mylib “d:sasdata”;,显示所有已指定的逻辑库,libname _ALL_ list;,取消逻辑库名,/*取消指定逻辑库名*/ libname clear; /*取消所有指定逻辑库名，不包括系统的四个*/ libname _ALL_ clear ;,数据步,SAS 编程初步,data ;
3、 /*若不指定逻辑库名，缺省建立Work中的临时表*/input $ $ $; /*$为字符型*/cards; d11 d12 d1k d21 d22 d2k ; run;,/*若cards后的数据行有重复的域，可以在input语句后增加行停留符”*/ data ; input $ $ $ ;cards; d11 d12 d1k d21 d22 d2k ; run;,过程步,SAS 编程初步,proc dada=;/;/; run;,/*过程步常用语句之一：var，指定分析变量*/ var ;/*过程步常用语句之二：by，指定分类变量，一般先用sort排序*/ by ;/*过程步常用语句之三：
4、class，指定分类变量，不需要先排序*/ class ;,单因素方差分析,SAS 应用方差分析,/* ANOVA过程*/ proc anova data = ;class ; /*必须，且位于model之前*/model =/; /*主、交互及嵌套效应模型*/means =/; /*计算各水平下因变量均值、标准差，并进行组间多重比较*/ run;*适用条件： *1. 各水平的观测彼此独立 *2. 各水平的观测均为正态分布的样本 *3. 各组方差相等,单因素方差分析,SAS 应用方差分析,/* 例题*/ data exam;input hangye $ tousu ;cards; 零售 57
5、零售 66 零售 49 零售 40 零售 44 旅游 68 旅游 39 旅游 29 旅游 45 旅游 56 航空 31 航空 49 航空 21 航空 34 航空 40 家电 44 家电 51 家电 65 家电 77 家电 58 ; proc anova data = exam;class hangye;model tousu = hangye; run;,线性回归分析,SAS 应用回归分析,/* REG过程*/ proc reg data = /;var ; model =/; /*指定要拟合的回归模型*/print ;plot =/; /*对两个变量绘制散点图，符号缺省即可，如果指定显示符号
6、，则需单引号括起来*/ run;,线性回归分析,SAS 应用回归分析,/* REG过程例*/ data sandian;input y x1 x2 ;cards; 6.2 215 136.5 7.5 250 136.5 4.8 180 136.5 5.1 250 138.5 4.6 180 138.5 4.6 215 138.5 2.8 180 140.5 3.1 215 140.5 4.3 250 140.5 4.9 215 138.5 4.1 215 138.5 ; proc reg data = sandian;var y x1 x2;model y=x1 x2; run;,多项式回归分
7、析REG,SAS 应用回归分析,/* REG过程*/ proc reg data = /;var ; model =/; /*指定要拟合的回归模型*/print ;plot =/; /*对两个变量绘制散点图，符号缺省即可，如果指定显示符号，则需单引号括起来*/ run; /* 关键是数据步，需要增加多项式项 data 新数据集;set 原数据集;新变量=原来变量的多项式表达式; */,多项式回归分析REG,SAS 应用回归分析,/* REG过程例*/ data sandian;input y x1 x2 ;cards; 6.2 215 136.5 7.5 250 136.5 4.8 180 1
8、36.5 5.1 250 138.5 4.6 180 138.5 4.6 215 138.5 2.8 180 140.5 3.1 215 140.5 4.3 250 140.5 4.9 215 138.5 4.1 215 138.5 ; data sandianx;set sandian; x1x2=x1*x2; proc reg data = sandian;var y x1 x2 x1x2;model y=x1 x2 x1x2; run;,多项式回归分析GLM,SAS 应用回归分析,/* GLM过程*/ proc glm data=;class ;model =/选项 run;/* 该过程
9、省去了很多选项和参数， 只把比较有用的列了出来。 */,多项式回归分析GLM,SAS 应用回归分析,/* GLM过程例*/ data exam; input g x y zb; cards; 1 30 10 280 1 25 11 260 1 35 13 330 1 40 14 400 1 45 14 410 2 20 12 270 2 18 11 210 2 25 12 280 2 25 13 300 2 23 13 290 3 40 14 410 3 45 15 420 3 48 16 425 3 50 18 450 3 55 19 470 ; proc glm data=exam; mo
10、del zb=x y x*x x*y y*y; run;,判别分析,SAS 应用判别分析,/* DISCRIM过程*/ proc discrim data=;class ;priors ; /*缺省为0.5*/var ; run;/* DISCRIM过程例，由于没有真实的数据集，该程序只是演示*/ proc discrim data = exam distance listerr; /*distance表明可根据距离判别样品有无差异，listerr列出错误归类的观测结果*/class type;var x1-x4; run;,聚类分析系统聚类,SAS 应用聚类分析,/* CLUSTER过程,聚
11、类分析系统聚类,SAS 应用聚类分析,/* CLUSTER过程,主成分分析,SAS 应用主成分分析,/* PRINCOMP过程*/ proc princomp data= 选项列表; var ; run; /* 选项列表中，out= 输出原始数据集和 主成分得分的数据集；outstat=统计量 的数据集；covariance或者 cov 表示从 协方差阵出发计算，缺省从相关阵出发 计算；n要计算的主成分个数，缺省时 计算全部；standard或者std要求在out= 数据集中把主成分得分标准化，默认为 相应的特征值 */,主成分分析,SAS 应用主成分分析,/* PRINCOMP过程例*/ p
12、roc princomp data=exam n=4 out = ex1 outstat=ex2; var x1-x10; proc print data = ex1; run;/* 该例子没有真实的数据集，所以只是作 为演示用，具体问题具体分析 */,SAS的内容比较多，前面都是针对考试的五个要求而设的，力求简洁实用。 Learn more by yourself,SAS学习,/*或许不能运行，但规则如此*/ data learnSAS; input what you do tomorrow;cards; 1 2 3 4 ; proc work_life_may_be data = learnSAS; run;,
