方差分析检验.pptx

方差分析检验一、方差分析得概念一、方差分析得概念 方方差差分分析析(ANOVA)又又称称变变异异数数分分析析或或F检检验验,其其目目得得就就是是推推断断两两组组或或多多组组资资料料得得总总体体均均数数就就是是否否相相同同,检检验验两两个个或或多多个个样样本本均均数数得得差差异异就就是是否否有有统统计计学意义。学意义。离均差平方与又称为变异数。离均差平方与又称为变异数。方方差差得得概概念念:方方差差即即为为标标准准差差得得平平方方,又又称称为为平平均均离均差平方与离均差平方与,简称简称“均方均方”。二、方差分析得用途二、方差分析得用途1 1、对对两两组组或或多多组组计计量量资资料料得得均均数数差差异异进进行行显著性检验。显著性检验。2 2、可可用用于于完完全全随随机机设设计计(单单因因素素)、随随机机区区组组设设计计(双双因因素素)、析析因因设设计计、拉拉丁丁方方设设计计与正交设计等资料得统计处理。与正交设计等资料得统计处理。3 3、可可对对两两因因素素间间交交互互作作用用差差异异进进行行显显著著性性检验。检验。4 4、进行方差齐性检验。、进行方差齐性检验。三、方差分析得条件三、方差分析得条件1 1、被比较得资料要有可比性。、被比较得资料要有可比性。2 2、被比较得资料要有实际意义。、被比较得资料要有实际意义。3 3、被比较得资料要呈正态分布。、被比较得资料要呈正态分布。4 4、被比较得资料各组方差齐同。、被比较得资料各组方差齐同。四、方差分析得基本思想四、方差分析得基本思想 例例 1 1 某医院用三种不同疗法治疗同种疾病某医院用三种不同疗法治疗同种疾病,以体温降至正以体温降至正常所需要得天数为指标常所需要得天数为指标,15,15例患者体温降至正常所需要得天数资料例患者体温降至正常所需要得天数资料如下如下:15 15例患者体温降至正常所需要得天数例患者体温降至正常所需要得天数 甲法甲法 乙法乙法 丙法丙法 5 5 7 5 5 7 5 5 9 5 5 9 5 7 9 5 7 9 7 7 9 7 7 9 7 7 9 7 7 9 29 33 43 29 33 43 5 5 5 n=15 5 5 5 n=15 5 5、8 68 6、6 86 8、6 x=7 6 x=7 173 173、0 2290 229、0 3730 373、0 0 15例患者体温降至正常所需要得天数甲法 乙法 丙法 5 5 7 5 5 9 5 7 9 7 7 9 7 7 9 1515例患者体温降至正例患者体温降至正常所需要得天数不尽相同常所需要得天数不尽相同,如果把每个患者所需得如果把每个患者所需得天数与总均数得差异之与天数与总均数得差异之与称为总变异称为总变异,总变异用离总变异用离均差平方与均差平方与(SS)来描述。来描述。答案答案 1 1 造成总变异得造成总变异得原因有两种可能原因有两种可能:1 1、由于个体差异得存在、由于个体差异得存在而引起而引起,即组内变异引起。即组内变异引起。组内变异组内变异,即由于随机误即由于随机误差得原因使得各组内部患者差得原因使得各组内部患者体温降至正常所需要得天数体温降至正常所需要得天数各不相等。各不相等。问题问题 1 1造成总变异得造成总变异得原因就是什么？原因就是什么？15例患者体温降至正常所需要得天数甲法 乙法 丙法 5 5 7 5 5 9 5 7 9 7 7 9 7 7 92 2、由于治疗方法不同而引、由于治疗方法不同而引起起,即组间变异引起。即组间变异引起。组间变异组间变异,即由于各组治即由于各组治疗方法疗方法(处理因素处理因素)不同使得不同使得各组均数大小不等。各组均数大小不等。15例患者体温降至正常所需要得天数甲法 乙法 丙法 5 5 7 5 5 9 5 7 9 7 7 9 7 7 9SS总总=SS组间组间+SS组内组内 v总总=v组间组间+v组内组内总变异与组间变异与组内变异得关系为总变异与组间变异与组内变异得关系为:总变异总变异=组间变异组间变异 +组内变异组内变异10大家应该也有点累了,稍作休息大家有疑问得大家有疑问得大家有疑问得大家有疑问得,可以询问与交流可以询问与交流可以询问与交流可以询问与交流处理因素治疗方法抽样误差个体差异组 间 变 异组 内 变 异总 变 异总变异来源分析总变异来源分析 因组间变异数大小与组数因组间变异数大小与组数(组间自由度组间自由度K-1)有关有关,故用组故用组间变异数除以自由度所得组间均方来表示组间变异间变异数除以自由度所得组间均方来表示组间变异。k=组数 因组内变异数大小与各样本含量大小即组内自由度因组内变异数大小与各样本含量大小即组内自由度(ni 1)有关有关,故用组内变异数除以组内自由度所得组内均方来表示故用组内变异数除以组内自由度所得组内均方来表示组内变异。组内变异。问问题题 2 2 当当处处理理因因素素不不起起作作用用时时,从从理理论论上上来来说说 F F 值就是多大？值就是多大？答答案案 2 2 当当处处理理因因素素不不起起作作用用时时,从从理理论论上上来说来说T=0,T+E=E,MST=0,T+E=E,MS 组间组间=MS=MS组内组内,F=1 F=1。当当处处理理因因素素起起作作用用时时,从从理理论论上上来来说说T0,T0,T T+E E E,E,F F 1,1,多多大大得得F F值值才才认认为为均均数数间间得得差差异异在在统统计计学学上上有有显显著著性性意意义义,可可根根据据自自由由度度查查F F值值表表,确确定定P P值后作出推断。值后作出推断。五、方差分析得方法五、方差分析得方法 根根据据实实验验设设计计方方法法不不同同,完完全全随随机机设设计计(单单因因素素)、随随机机区区组组设设计计(双双因因素素)、析析因因设设计计、拉拉丁丁方方设设计计与与正正交交设设计计等等资资料料应应选选用用不不同同得得方方差差分分析方法。析方法。问题问题 1 1 右边右边 资料就是用何种实验资料就是用何种实验设计方法进行实验所设计方法进行实验所得资料？得资料？答案答案 1 1 用用完全随机设计完全随机设计(单因素单因素)方法进行实验所得资方法进行实验所得资料。料。15例患者体温降至正常所需要得天数 甲法 乙法 丙法 5 5 7 5 5 9 5 7 9 7 7 9 7 7 91、完全随机设计(单因素)资料得方差分析H0:H1:15例患者体温降至正常 所需要得天数 甲法 乙法 丙法 5 5 7 5 5 9 5 7 9 7 7 9 7 7 9 29 31 43 29 31 43 5 5 5 n=15 5 5 5 n=15 5 5、8 68 6、6 86 8、6 x=6 x=7 7 173 173、0 2290 229、0 3730 373完全随机设计多个样本均数比较得方差分析公式 变异来源 SS MS F总 N1组间 k1 组内 Nk 方差分析表变异来源 SS v MS F组 间 20、8 2 10、4 6、5组 内 19、2 12 1、6 确定P值:查F值表,F0、05(2,12)=3、88,F0、01(2,12)=6、93,本例F=6、5,3、886、5P0、01。放置不同时间得血糖浓度 受试者号 放置时间(分)(区组)0 45 90 135 1 5、27 5、27 4、94 4、61 2 5、27 5、22 4、88 4、66 3 5、88 5、83 5、38 5、00 4 5、44 5、38 5、27 5、00 5 5、66 5、44 5、38 4、88 6 6、22 6、22 5、61 5、22 7 5、83 5、72 5、38 4、88 8 5、27 5、11 5、00 4、44 随机区组设计资料随机区组设计资料 方差分析方差分析 研究酵解作用对血糖浓度得影响,从8名健康人中抽取了血液并制备成血滤液,每个受试者得血滤液分成四份,再随机把4份血液分别放置0、45、90、135分钟后测定其血糖浓度,试分析放置不同时间得血糖浓度有无变化。放置不同时间得血糖浓度 受试者号 放置时间(分)(区组)0 45 90 135 1 5、27 5、27 4、94 4、61 2 5、27 5、22 4、88 4、66 3 5、88 5、83 5、38 5、00 4 5、44 5、38 5、27 5、00 5 5、66 5、44 5、38 4、88 6 6、22 6、22 5、61 5、22 7 5、83 5、72 5、38 4、88 8 5、27 5、11 5、00 4、44随机区组设计资料得类型随机区组设计资料得类型 1、单因素多水平多组配对设计、单因素多水平多组配对设计资料。资料。2、多因素单水平多组配对设计、多因素单水平多组配对设计资料资料 3、双因素多水平设计资料。、双因素多水平设计资料。配对配对得方法得方法:1、同体配对研究、同体配对研究 (1)处理前后配对研究处理前后配对研究 (2)同体左右配对研究同体左右配对研究 2、异体配对研究、异体配对研究 (1)动物动物:种系、性别、胎次相种系、性别、胎次相同同,体重体重10%。(2)人群人群:种族、性别、病种相种族、性别、病种相同同,病情或并发症相似病情或并发症相似,年龄年龄5岁。岁。随机区组设计资料变异来源分析随机区组设计资料变异来源分析总变异=处理组间变异+区组间变异+误差变异SS总=SS处理组间+SS区组间+SS误差 不同治疗方法大鼠体重增重变化不同治疗方法大鼠体重增重变化 常规疗法 甲 法 乙 法区组号 治疗前 治疗后 治疗前 治疗后 治疗前 治疗后 1 160 240 160 264 165 290 2 178 253 176 306 172 317 3 156 262 158 258 158 288 4 164 272 164 284 164 298 5 172 265 172 290 174 320 6 168 277 168 286 170 305 7 165 282 168 290 168 310 8 170 290 172 300 168 312 9 164 276 162 280 162 290 10 174 268 176 296 174 300 四种大鼠注射不同剂量雌激素(mg/Kg)后子宫重量如下表,试进行分析。鼠 种 剂 量 0、1 0、2 0、4 0、8 A 98 106 116 145 B 40 42 68 115 C 56 70 111 133 D 38 42 63 87 案例案例 2 2 现有8 名贫血患儿治疗前后血红蛋白含量资料如下:8 名贫血患儿治疗前后血红蛋白含量(g/L)患者编号 治疗前 1个月 3个月 b 1 115 104 120 339 2 107 110 130 347 3 98 109 125 332 4 110 115 130 355 5 90 109 125 324 6 117 115 115 347 7 115 120 135 370 8 115 123 127 365 ni 8 8 8 N=24 867、0 905、0 1007、0 x 108、4 113、1 125、9 94617、0 102657、0 127029、0 1515例患者体温降至正常例患者体温降至正常 所需要得天数所需要得天数 甲法甲法 乙法乙法 丙法丙法 5 5 7 5 5 7 5 5 9 5 5 9 5 7 9 5 7 9 7 7 9 7 7 9 7 7 9 7 7 9 现有现有8 8 名贫血患儿治疗前后血红蛋白含量资料如下名贫血患儿治疗前后血红蛋白含量资料如下:8 8 名贫血患儿治疗前后血红蛋白含量名贫血患儿治疗前后血红蛋白含量(g/L)(g/L)患者编号患者编号 治疗前治疗前 1 1个月个月 3 3个月个月 b b 1 115 104 120 339 1 115 104 120 339 2 107 110 130 347 2 107 110 130 347 3 98 109 125 332 3 98 109 125 332 4 110 115 130 355 4 110 115 130 355 5 90 109 125 324 5 90 109 125 324 6 117 115 115 347 6 117 115 115 347 7 115 120 135 370 7 115 120 135 370 8 115 123 127 3658 115 123 127 365 问问题题 2 2 例例2 2得得总总变变异异来来源源与与例例1 1有有何何异异同同点点？答答案案 2 2 共共同同点点就就是是其其总总变变异异来来源源都都就就是是来来自自于于处处理理因因素素变变异异与与抽抽样样误误差差变变异异,这这不不仅仅就就是是它它们们得得共共同同点点,而而且且就就是是所所有有方方差差分分析析资资料料总变异来源得共同点。总变异来源得共同点。不同点就是例不同点就是例1 1只有只有一个处理因素一个处理因素,例例2 2有两有两个处理因素个处理因素,因此总得因此总得变异来源可分解为变异来源可分解为:SSSS总总=SS=SS因素因素1 1+SS+SS因素因素2 2+SS+SS误差误差例1例2随机区组设计随机区组设计(双因素双因素)资料得方差分析方法为资料得方差分析方法为:HoHo:不同时间不同时间HBHB含量相同含量相同,H H1 1:不同时间不同时间HBHB含量不同或部分不同。含量不同或部分不同。H H0 0:不同得受试者不同得受试者HBHB含量相同含量相同,H H1 1:不同得受试者不同得受试者HBHB含量不同或部分不同。含量不同或部分不同。例例 2 2 现有现有8 8 名贫血患儿治疗名贫血患儿治疗前后血红蛋白含量资料如下前后血红蛋白含量资料如下:8 8 名贫血患儿治疗前后血红蛋白含量名贫血患儿治疗前后血红蛋白含量(g/L)(g/L)编号编号 治疗前治疗前 1 1个月个月 3 3个月个月 b b 1 115 104 120 339 1 115 104 120 339 2 107 110 130 347 2 107 110 130 347 3 98 109 125 332 3 98 109 125 332 4 110 115 130 355 4 110 115 130 355 5 90 109 125 324 5 90 109 125 324 6 117 115 115 347 6 117 115 115 347 7 115 120 135 370 7 115 120 135 370 8 115 123 127 365 8 115 123 127 365 n ni i 8 8 8 N=24 8 8 8 N=24 867 867、0 9050 905、0 10070 1007、0 0 x 108 x 108、4 1134 113、1 1251 125、9 9 94617 102657 127029 94617 102657 127029 配伍组方差分析得计算公式 变异来源 SS MS F总 N1 处理间 k1 配伍间 b1 误差 SS总SS处理SS配伍 总-处理-配伍