文档介绍：R语言学****系列28-、,有些随机因素是很难人为控制的,但它们又会对结果产生显著影响。如果忽略这些因素的影响,则有可能得到不正确的结论。这种影响的变量称为协变量(一般是连续变量)。例如,研究3种不同的教学方法的教学效果的好坏。检查教学效果是通过学生的考试成绩来反映的,而学生现在考试成绩是受到他们自身知识基础的影响,在考察的时候必须排除这种影响。协方差分析将那些难以控制的随机变量作为协变量,在分析中将其排除,然后再分析控制变量对于观察变量的影响,从而实现对控制变量效果的准确评价。协方差分析要求协变量应是连续数值型,多个协变量间互相独立,且与控制变量之间没有交互影响。前面单因素方差分析和多因素方差分析中的控制变量都是一些定性变量,而协方差分析中既包含了定性变量(控制变量),又包含了定量变量(协变量)。协方差分析在扣除协变量的影响后再对修正后的主效应进行方差分析,是一种把直线回归或多元线性回归与方差分析结合起来的方法,其中的协变量一般是连续性变量,并假设协变量与因变量间存在线性关系,且这种线性关系在各组一致,即各组协变量与因变量所建立的回归直线基本平行。当有一个协变量时,称为一元协方差分析,当有两个或两个以上的协变量时,称为多元协方差分析。(1)自变量是分类变量,协变量是定距变量,因变量是连续变量;对连续变量或定距变量的协变量的测量不能有误差;(2)协变量与因变量之间的关系是线性关系,可以用协变量和因变量的散点图来检验是否违背这一假设;协变量的回归系数(即各回归线的斜率)是相同的,且不等于0,即各组的回归线是非水平的平行线。否则,就有可能犯第一类错误,即错误地接受虚无假设;(3)自变量与协变量相互独立,若协方差受自变量的影响,那么协方差分析在检验自变量的效应之前对因变量所作的控制调整将是偏倚的,自变量对因变量的间接效应就会被排除;(4)各样本来自具有相同方差σ2的正态分布总体,即要求各组方差齐性。二、=均值+分组变量影响+协变量影响+(1)其中,为所有协变量的平均值。注:在方差分析中,协变量影响是包含在随机误差中的,在协方差分析中需要分离出来。用协变量进行修正,得到修正后的yij(adj)为就可以对yij(adj)做方差分析了。=分组变量离差+协变量离差+随机误差,(1)计算总离差平方和时,记总离差平方和:最终要检验分组自变量对因变量有无显著作用。原假设H0:无显著作用。假设检验是在H0为真条件下进行,可认为ti=0,则按最小二乘法原理线性回归可得到β的估计值记修正的总离差平方和(残差平方和)为Tyy(adj),则,自由度为n-2注:为回归平方和,若(回归线为水平线),表示协变量x对y无作用,用方差分析就可以解决了。(2)计算组内离差平方和时,记组内总离差平方和:根据协方差分析的基本假设:各组内回归系数相等(做协方差分析时需要检验这一点),得到组内回归系数βw的估计值记修正的组内总离差平方和(组内残差平方和)为Eyy(adj),则,自由度为n-k-1其中,为组内回归平方和,当时,组内总离差平方和认为完全是由随机因素引起的,Eyy(adj)就是随机为误差。这里的是的加权平均值。(3)计算分组变量离差平方和Byy(adj),它反映的是各个水平之间的差异。即,分组变量离差=总离差-协变量离差-随机误差。于是,就可以进行组间无差异检验了:,在做协方差分析前,需要依次做两个假设检验:(1)协变量对因变量的影响对与各组来说都是相同的,即各组回归系数相等:;步骤:①先按回归系数相等和不相等分别表示模型并计算出误差平方和其中,.②计算F值若F值小于临界值Fα,则说明各组回归系数无显著差异(相等)。(2),(1)检验数据是否满足假设条件:正态分布性、方差齐性、各分组通过协变量预测因变量的回归斜率相同;(2)检验效应因子的显著性;(3)估计校正的组均值;(4)检验校正的组均值之间的差异。三、R语言实现协方差分析要求数据满足:正态性、方差齐性、各分组通过协变量预测因变量的回归斜率相同。R语言用aov()函数进行协方差分析,基本格式为:aov(formula,data,...)其中,data为数据框;formula为协方差公式形式,形如y~x+A,x为连续型协变量,A为组别因子。例1研究分别接受了3种不同的教学方法的3组学生,在数学成绩上是否有显著差异,数据文件“”。先不考虑数学入学成绩,只以“教学方法”为分组变量,“后测成绩”为因变量进行单因素方差分析:setwd("E:/办公资料/R语言/R语言学