资源预览内容
第1页 / 共23页
第2页 / 共23页
第3页 / 共23页
第4页 / 共23页
第5页 / 共23页
第6页 / 共23页
第7页 / 共23页
第8页 / 共23页
第9页 / 共23页
第10页 / 共23页
亲,该文档总共23页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
课时授课计划课次序号:15 一、课题:3.2 两因素等重复试验下的方差分析 二、课型:新授课三、目的要求:1.掌握两因素等重复实验下的方差分析理论与方法、模型的建立与显著性检验;2.掌握利用方差分析的SAS过程解决有关实际问题.四、教学重点:方差分析方法的基本理论;利用方差分析的SAS过程解决有关实际应用问题.教学难点:方差分析方法的基本理论;利用方差分析的SAS过程解决有关实际应用问题.五、教学方法及手段:传统教学与上机实验相结合六、参考资料:应用多元统计分析,高惠璇编,北京大学出版社,2005;使用统计方法与SAS系统,高惠璇编,北京大学出版社,2001;多元统计分析(二版),何晓群编,中国人民大学出版社,2008;应用回归分析(二版),何晓群编,中国人民大学出版社,2007;统计建模与R软件,薛毅编著,清华大学出版社,2007.七、作业:3.6 八、授课记录:授课日期班次九、授课效果分析:复习 单因素方差分析1统计模型因变量Y因素,水平,上观测值 , 2显著检验 或 ,拒绝3置信区间 置信度的置信区间, 置信度的置信区间为个的置信度至少的Bonferroni同时置信区间3.2 两因素等重复实验下的方差分析3.2.1 统计模型设影响Y的因素有两个,分别记为A和B,其中A有a个不同水平,B有b个水平在因素A和B的各水平下均做c(c1)次实验,记为水平组合下第k次实验的Y的观测值,则两因素等重复试验下的方差分析数据可表示为表3.7的形式表3.7 两因素等重复方差分析数据 因素B因素A和对于任一水平组合(总体), Y观测值为, 则各样本间是相互独立的样本观察值可看成是来自均值为的总体,即 , 令,为水平组合下Y的随机误差, 则,这样,就是其均值与随机误差迭加而产生的因此,两因素重复试验下方差分析的统计模型: 为便于统计分析,我们需要对水平组合上的样本均值作进一步分解,为此引入如下记号:其中为总平均, ,是因素水平与因素水平在单元上所有观察值的平均,为因素A的水平的效应,为因素B的水平的效应进一步有其中反映了水平组合对Y的效应一般情况, 其差称为与的交互效应因此容易验证: 因此两因素等重复下的方差分析模型等价地改写为如下形式:3.2.2交互效应及因素效应的显著性检验一偏差平方和分解下面先对Y的观测的总平方和进行分解:,观测数据的总(偏差)平方和为 其中因素的平方和由于,为的无偏估计,故度量的各水平效应的估计量的变化因素的平方和由于,为的无偏估计,故度量的各水平效应的估计量的变化交互效应的平方和由于,为的无偏估计,故度量和的交互效应的估计量的变化误差平方和度量了来自各总体的观测值与其样本均值的差异,反映了误差的变化 由于且相互独立,可得分别称为的自由度,称为的自由度,为上述四个自由度的和令 因素的均方, 则 因素的均方, 交互效应均方,误差均方, , 为的无偏估计 二假设检验对两因素的情况,方差分析的主要目的除了考察因素A或B的各水平对因变量Y的影响有无显著差异外,还要考虑A和B之间是否存在交互作用,因为交互作用的存在会直接影响到对A和B影响检验结果的解释涉及如下三个检验问题: 不全相等 不全相等 检验问题也可以改写成: 利用上述结果,构造适当的统计量检验上述假设为的无偏估计,如果假设成立,取值接近,如果假设不成立,则有增大的趋势因此,针对检验分别构造统计量,分别有 如果,各检验统计量的值变大,则拒绝原假设各检验的值分别为 其中为统计量观测值给定显著性水平,如(或),则拒绝(或)否则不能拒绝(或)结果如下:表3.8 两因素(ab)等重复数c试验下的方差分析表形式变异来源source离差平方和SS自由度df均方MSF统计量FP值P因素ASSA因素BSSBAB交互效应SSAB误差SSE总和SST= SSA+ SSB+ SSAB +SSE注意:检验步骤:1)先检验,如不拒绝,即交互作用不显著时,再考察和的效应的显著性因为当()全为0时,即与无交互作用,则在的任何水平上 = ,在的各水平上均相等且完全由在水平和的效应之差确定,因此,检验假设的结论真实的反映了仅由的各水平对的影响是否显著2)如拒绝,交互作用显著时,通过估计和比较因素和各水平组合(A,B)上的均值考察因素的联合影响如果和存在交互作用,则不全为0,对于的两个水平和在的第个水平上的两个组合和下均值差为 因此当()不全为0时,除与有关外,还可能与有关,即处在不同水平,有所不同如图3.1所示,假设A与B均有两个水平和 (a) (b) 图3.1 有交互效应时A的各水平均值在B的不同水平上的差异(a) 在下均值差而认为(即与不全为0),差异主要表现在在水平上的差异;(b) 在下均值差与相反,综合有,主要表现在在和水平上的差异相互抵消,使得综合差异为零因此,在有交互效应时,尤其是交互效应显著,而因素与的效应不显著时,检验每个因素显著差异实际意义不大,应慎重此情况下,要进一步考察各因素对的影响的显著性,只能将一个因素的各个水平逐个给定,在给定的水平上考察另一因素的各水平均值之间的差异来了解该因素对的影响例3.5 某高校为了了解数学专业和计算机科学专业的低年级学生、高年级学生及研究生在人文科学知识方面的差异,从不同专业和不同年级的学生中任选四名学生参加有关考试,其成绩如表3.9所示,假设考试成绩服从两因素的方差分析模型,对其作方差分析表 3.9 人文社科知识的考试成绩 级别专业低年级 高年级 研究生数学计算机81 78 79 78 75 80 78 73 82 80 83 8889 82 77 90 79 80 75 78 93 93 86 95解:因变量Y为成绩,有两个因素,记专业因素为,有两个水平(数学),(计算机);学生级别为因素,有三个水平(低年级),(高年级),(研究生)因此,利用SAS系统proc anova过程作方差分析,程序及结果如下:data examp3_5;input majors $ classes $ grade ; /*输入majors课程,年级classes,成绩grade*/cards;a1 b1 81 a1 b1 78 a1 b1 79 a1 b1 78a1 b2 75 a1 b2 80 a1 b2 78 a1 b2 73a1 b3 82 a1 b3 80 a1 b3 85 a1 b3 88a2 b1 89 a2 b1 82 a2 b1 77 a2 b1 90a2 b2 79 a2 b2 80 a2 b2 75 a2 b2 78a2 b3 93 a2 b3 93 a2 b3 86 a2 b3 95;run;proc anova data=examp3_5; /* 调用方差分析的anova过程 */class majors classes; /* 因素变量名称为majors专业、classes级别 */model grade=majors classes majors*classes; /* 因变量grade,因素变量专业、级别、因素效应 */run;SAS 系统 11:29 Tuesday, October 21, 2008The ANOVA Procedure Class Level Information 因素 水平 因素变量 Class Levels Values 因素变量A majors a=2 a1 a2 因素变量B classes b=3 b1 b2 b3 Number of observations n=abc=24 c=4注:,因素变量为自变量,个数为自由度p-1=ab-1=5先检验,统计量,观测值,检验值 SAS 系统 11:29 Tuesday, October 21, 2008 2 The ANOVA Procedure因变量:成绩 Dependent Variable: grade表 3.10 考试成绩的方差分析结果 Sum of Source DF Squares Mean Square F Value Pr F方差来源 自由度 平方和SS 均方 F值 p值Model(模型) p-1=ab-1=5 =637.0000000 =127.400000 =9.34 0.0002Error(误差) n-p=ab(c-1)=18 =245.5000000 =13.6388889Corrected Total n-1=abc-123 =882.5000000 R-Square Coeff Var Root MSE grade
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号