武汉大学SAS教案第六讲方差分析.pptxVIP

芍卡方检验;u方差今析 r?i：两个或两个以上组均数的比较;■u方差今析 丄,主效应:由分类变量本身引起，不考虑其他因素的影响, 可以用分类变量本身表示。 . Y=a b c ;其中a be是主效应 ■交叉效应：即交互作用。用*连接两个变量以表示它们之 间的交互作用，检验在A因素的各个水平B因素的效应是 否相同，即一个因素的效应是否依赖于交叉项里其他因 素 . Y=a b c a*b a*c a*b*c ;其中a*b a*c a*b*c 是交互效应 ■把Y称为因变量，abc称为自变量。因变量的个数反映 出是一元或多元方差分析，自变量的个数反映的是方差 分析的因素的个数;f 多组此较的假设检验 ■例：比较5种品牌木板的耐久性。;基本思想： I总变差=组间差别+组内误差;总变差=组I可差别+组内误差;1、单向方差分析;1、单向方差分析 (2)数据集的建立 DATA AN0VA1; INPUT BRAND $ W @@; CARDS; XI 2.3X1 2.1 XI 2.4X1 2.5 X2 2.2 X2 2.3 X2 2.4 X2 2.6 X3 2.2X3 2.0 X3 1.9 X3 2.1 X4 2.4 X4 2.7 X4 2.6 X4 2.7 X5 2.3 X5 2.5 X5 3.3 X5 3.4 RUN;;单向方差分析;1、单向方差分析 （4）结果分析 ?:?输出的第一部分：给出了把数据分组的 变量的名字（BRAND）,此变量所取的水 平、此变量的各个值和数据集中的观测 数。 .:?输出的第二部分：标有BRAND的一行PrF 一列屮给出了比较各组的P值，此例中为 0.0112o得出各个品牌的平均磨损有显 著的差异;于多重比较： ■以上的方差分析只能发现各组间有显著 差异，但无法知道哪些均值不同;菜单实现 ? Solutions/Analysis/Analyst 。打开数据集 ?:? Statistics/ANOVA/One-Way ANOVA w—Dependent brand Independent Means/brands/snk/ok ? OK .：.结论：X3品牌的耐久性较其余品牌的好;程序实现 PROC ANOVA; CLASS BRAND; MODEL W=BRAND; Means brand/snk; RUN; 其中，means语句请求进行组间 均值的两两比较;顼向方差分析的非参数检验 (1) 假定 正态性的假定不满足 (2) 数据集的建立 同上;、单向方差分析的非参数检验;顼向方差分析的非参数检验 （3）用菜单分析 ? Solutions/Analysis/Analyst ?:?打开数据集 ? Statistics/ANOVA/Nonparametric Oneway ANOVA w^Dependent brand Independent ? OK;顼向方差分析的非参数检验 （4）结果分析 看ProbCHISQ右边的数字 0. 0175;、二因素方差分析 例：用4种不同的方法治疗8 名患者，其血浆凝固吋间的 资料如下，试分析影响血浆 凝固的因素;L盘、二因素方差分析 ，;VIEWTABLE : T emp.;产单 ? Solutions/Analysis/Analyst ?:?打开数据集 ?:? Statistics/ANOVA/Factorial ANOVA X 今 Dependent a,b 今 Independent Means/a/snk/ok ? OK;［结果 ■总体有非常显著的差异 . A因素：P0.0025, A因素（治疗方法） 对血浆凝固时间影响很大 . B因素：P0.0001, B因素（不同患者） 血浆凝固时间差异很大 ■组间两两比较结果：第4组与其余3纟卩均 值间差异显著;土工衡数据的方差分析 ［例］健康男子各年龄组淋巴细胞转换率之间 的差异是否显著？（单因素） 10-20岁 58 61 61 62 63 68 70 70 74 78 41-50岁 54 57 57 58 60 60 63 64 66 61-75岁 43 52 55 56 60;VTEWTABLE;Proc glm data=temp.an; Class a; Model x=a; Run;;非平衡数据的方差分析;衡数据的方差分析 ［例］二因素非平衡数据的方差分析;、亿;上虹衡数据的方差分析 用菜单实现： . Solutions/Analysis/Analyst,并打开数据集 ■ Statistics/Anova/Linear Models X- Dependent a,B-class MODEL/STANDARD MODELS/Effects upto 2-way interactions（考虑交互效应） ■ O