方差分析4幻灯片资料.ppt

方差分析;1．1.意义 前一章介绍了2个样本均数比较的假设检验方法，但对于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检验。 ;例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名患者与健康人的血磷值大小不等，称这种变异为总变异。可以用总离均差平方和及N来反映 SST= 总自由度 νT=N-1。 ;2组样本均数不等，这种变异称为组间变异，反映了克山病对血磷值的影响和随机误差;2个组各组内部血磷值也不等，这种变异称为组内变异，其大小可用2组组内离均差平方和及各组例数ni来反映. 自由度νE=N-k（k是组数），它反映了随机误差。 ;SST=SSA+SSE νT=νA+νE;4．直观意义;分析思路——;①??? 不受比较的组数限制。 ②??? 可同时分析多个因素的作用。 ③??? 可分析因素间的交互作用。 ;①??? 各样本是相互独立的随机样本 ②??? 各样本来自正态总体 ③??? 各组总体方差相等，即方差齐 ;例 某社区随机抽取了30名糖尿病患者、IGT异常和正常人进行载脂蛋白(mg/dL)测定，结果如下，问三种人的载脂蛋白有无差别？ ;n1=11 n2=9 n3=10 N=30;85.7 ;所有人的载脂蛋白的变异可分解为两部分： 组内变异（SSE），反映载脂蛋白测定值的随机误差和个体差异; 组间变异（SSA），反映随机误差和不同的生理病理状态对载脂蛋白的影响。 ;（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 ;（3）??? 作出推断结论 查F表：F（0.05, 2,27）=3.35, F（0.01, 2,27）=5.49 按α=0.01水平,拒绝H0，接受H1，认为三种人载脂蛋白的总体均数不同(差异极显著)。 ;随机区组设计又称配伍组设计(Random Block Design)。 在农业中如要比较3种化肥的效果，即要比较施用不同化肥对农作物产量的影响，可以抽取若干块试验田，分成3组，每组施用一种肥料，但是地有贫瘠之分，如果某组分得的地都是瘦地，即使所施化肥效果好，但从农作物产量上也反映不出来，此时应将地的肥沃程度考虑进去，以消除田地的影响，这种因素称作区组因素。 ;例3 对小白鼠喂以A、B、C三种不同的营养素，目的是了解不同营养素增重的效果，采用随机区组设计方法，以窝别作为划分区组的特征，以消除遗传因素对体重增长的影响。 现将同品系同体重的24只小白鼠分为8个区组，每个区组3只小白鼠。三周后体重增量结果（克）如下，问小白鼠经三种不同营养素喂养后所增体重有无差别？ ;本例的总变异中除了不同营养素的变异（处理差异）、随机误差外，还存在着区组之间的变异，是由小白鼠的遗传特征不同所致。 xij表示第i个处理组第j个区组样本观察值， 表示第i个处理组的均数(= )， 表示第