方差分析淮安疾控信息网首页淮安疾控,健康热线,12320,疾(免费阅读).ppt

方差分析 analysis of variance 一．方差分析的基本思想 1．1.意义 前一章介绍了两个样本均数比较的假设检验方法，但对于3个、4个、5个均数或更多个的比较，t检验或u检验就无能为力了，或许有人会想起将几个均数两两比较分别得到结论，再将结论综合，其实这种做法是错误的。试想假设检验时通常检验水平α取0.05，亦即弃真概率控制在0.05以内，但将3个均数作两两比较，要作三次比较，可靠度成为 (1-0.05)3=0.857 四均数比较作6次 (1-0.05)6=0.735 五均数比较作10次 (1-0.05)10=0.599 六均数比较作15次 (1-0.05)15=0.463 鉴于以上的原因，对多组均数的比较问题我们采用方差分析（analysis of variance），简称ANOVA。 方差分析采用F检验统计量，也称F检验。 2．基本思想 先讲述几个记号： xij表示第i组第j个样本观察值， 例1 某克山病区测得11例克山病患者与13名健康人的血磷值（mmol/L）如下，问该地急性克山病患者与健康人的血磷值是否不同？ 患者x1：0.84 1.05 1.20 1.20 1.39 1.53 1.67 1.80 1.87 2.07 2.11 健康人x2：0.54 0.64 0.64 0.75 0.76 0.81 1.16 1.20 1.34 1.35 1.48 1.56 1.87 24名患者与健康人的血磷值大小不等，称这种变异为总变异。可以用总离均差平方和 SST= 及N来反映，总自由度 νT=N-1。 2组样本均数也不等，这种变异称为组间变异，反映了克山病对血磷值的影响和随机误差 3．三者关系 SST=SSA+SSE νT=νA+νE 4．直观意义 如果两组样本来自同一总体，即克山病患者与健康人血磷值相同，则理论上F应等于1，因为两种变异都只反映随机误差。由于抽样误差的影响，F值未必是1，但应在1附近。若F较小，我们断定2组均数相同，或者说来自同一总体，F较大，推断不是来自同一总体。 5．优点 ①??? 不受比较的组数限制。 ②??? 可同时分析多个因素的作用。 ③??? 可分析因素间的交互作用。 6．应用条件 ①??? 各样本是相互独立的随机样本 ②??? 各样本来自正态总体 ③??? 各组总体方差相等，即方差齐 二．成组设计的多个样本均数比较 例 某社区随机抽取了30名糖尿病患者、IGT异常和正常人进行载脂蛋白(mg/dL)测定，结果如下，问三种人的载脂蛋白有无差别？ 所有人的载脂蛋白的变异可分解为两部分： 组内变异，反映载脂蛋白测定值的随机误差和个体差异; 组间变异，反映随机误差和不同的生理病理状态对载脂蛋白的影响。 （1）?????? 建立假设和确定检验水平 H0: 三种人载脂蛋白的总体均数相等, μ1=μ2=μ3 H1: 三组总体均数不相等 α=0.05 （2）计算 C=(Σx) 2/N=(3309.5) 2/30=365093 SST=Σx2-C=372974.87-365093=7881.87 SSE=SST-SSA=7881.87-2384.026=5497.84 νT=N-1=29 νA=k-1=2 νE=N-k=30-3=27 MSA=SSA/νA =1192.01 MSE=SSE/νE =203.62 F=MSA/MSE=5.8540 得到方差分析表，查表确定P值： （3）??? 作出推断结论 按α=0.05水平拒绝H0，接受H1，认为三种人载脂蛋白的总体均数不同。 三．随机区组设计的方差分析 随机区组设计又称配伍组设计(Random Block Design)。 在农业中如要比较三种化肥的效果，即要比较施用不同化肥对农作物产量的影响，可以抽取若干块试验田，分成三组，每组施用一种肥料，但是地有贫瘠之分，如果某组分得的地都是瘦地，即使所施化肥效果好，但从农作物产量上也反映不出来，此时应将地的肥沃程度考虑进去，以消除田地的影响，这种因素称作区组因素。 在医学研究中，也存在区组因素，比如动物试验中动物的窝别、年龄、性别、地区的差异等。 例3 对小白鼠喂以A、B、C三种不同的营养素，目的是了解不同营养素增重的效果，采用随机区组设计方法，以窝别作为划分区组的特征，以消除遗传因素对体重增长的影响。现将同品系同体重的24只小白鼠分为8个区组，每个区组3只小白鼠。三周后体重增量结果（克）如下，问小白鼠经三种不同营养素喂养后所增体重有无差别？