统计分析系统SAS软件专题讲义.ppt

5.4.3 应用实例 方差分析基础 常用试验设计的ANOVA 缺省数据的方差分析 方差分析基础 单向分组资料的ANOVE 1）样本容量相等（model y=trt) 2）样本容量不相等(model y=trt) 3）数据转换（model yp=trt) 单一自由度的比较(model y=trt) 两向分组资料的ANOVA(model y=m h) 系统分组资料的ANOVA(model y=trt pot(trt)) 常用试验设计的ANOVA 完全随机化设计（单向分组资料…） 随机区组设计 1）单因素（model y=blk trt) 2）两因素(model y=blk a b a*b) 3）三因素(model y=blk a b a*b c a*c b*c a*b*c) (model y=blk a|b|c) 拉丁方设计(model y= r c t) 二裂式裂区设计(model y=blk a a*blk b a*b) 多点品种区域试验(model y=r(u) u v u*v) 模型定义中可用|和@n简化表达方法 协方差分析 过程格式 PROC GLM 选项; CLASS 变量表; MODEL 依变量=效应/选项; MEANS 效应/选项; LSMEANS 效应/选项; PROC GLM 语句选项可设定分析数据集等; CLASS 指明分类变量,协方差分析时必须设定,必须在MODEL语句之前; MODEL 定义协方差分析的线性数学模型, model y=a t;一般的协方差分析模型 model y=t x(t);分离斜率模型 model y=t x x*t;共斜率模型 model后SOLUTION的选项给出参数的估计值 MEANS 语句用于计算依变量的平均数,选项用于多重比较; LSMEANS 语句计算效应的最小二乘估计的平均数(LSM) 选项E=效应,设定测验误差项,缺省为试验分析误差. STDERR给出LSM的标准误, TDIEF, PDIEF要求显示测验H0:LSM(i)=LSM(j)的t值和概率值. 成对比较 UNIVARIATE过程 Univariate过程对数值变量给出比较详细的变量分布的描述，其中包括： 变量的极端值。 常用的百分位数，包括四分位数和中位数。 用几个散点图描绘变量的分布。 频数表。 确定数据为正态分布的检验。 Univariate过程常用的选项如下 NOPRINT 禁止统计报告在OUTPUT视窗中输出 PLOT 绘出茎叶图、箱式图和正态概率图 FREQ 给出频数表 NORMAL 对变量进行正态性检验 【关键字】 SAS中用关键字来指定所需要的统计量，事实上结果输出中用的就是各种关键字，常用的关键字有： 基本统计量 N MEAN STD(标准误) CV SUM VAR(方差) RANG 百分位数描述 MIN P1 P5 P10 Q1 MEDIAN Q3 P90 P95 P99 MAX 与假设检验有关的统计量 STDMEAN(标准误) T TTEST过程 组群比较 例 将钩体病人的血清分别用标准株和水生株做凝溶实验，测得稀释倍数如下。问两组的平均效价有无区别。 标准株：100 200 400 400 400 400 800 1600 1600 1600 3200 水生株：100 100 100 200 200 200 200 400 400 统计分析系统SAS软件 敖 雁 Tel:7979358 E-mail:ay@yzcn.net 第五章 方差分析和协方差分析 ANOVA过程 -----速度快，功能有限,平衡资料 GLM过程 -----速度慢，但功能强大（非平衡资料 协方差分析） 其实，这里的速度快慢只是相对而言，SAS的处理速度是首屈一指的。举个例子，一个有6600条记录的，7因素的，交叉设计的方差分析（是不是已经有人喊头痛了？），我先是用SPSS FOR WIN95 7.5来做，运行了大约10分钟才出结果。我又换用SAS FOR WIN95 6.12来做，结果用了――2.47秒！ §5.1 方差分析数据集的建立技巧 5.1.1 方差分析的数据集格式 统计分析所用的数据格式和我们在分析整理资料时所用的格式是不同的: -------数据集中应至少有一个结果变量，用于记录不同处理因素水平下观察值的大小； -------至少有一个处理因素变量，用于记录处理因素的类型及其水平数。 例5.1 某职业病防治院对31名石棉矿工中的石棉肺患者、可疑患者及非患者进行了用力肺活量测定，请给出数据集的结构。 解