统计学课件12

1、 本章教学目标l了解方差分析可以解决那些实际问题；l了解应用方差分析的基本条件；l掌握方差分析的基本概念及其分析方法；l正确使用 excel 软件求解单因素和双因素方差分析问题及其运行输出结果分析.1第第6章章 方差分析方差分析 本章主要内容本章主要内容 4.1 方差分析概述 4.2 单因素方差分析 4.3 双因素方差分析 本章重点本章重点：考虑交互作用的双因素方差分析 2 在生产经营管理过程中，通常有很多因素会影响产品的质量、产量、销售量等指标。 如农作物的产量受品种、肥料、气候、雨水、光照、土壤、播种量等众多因素的影响； 产品销售量受品牌、质量、价格、促销手段、竞争产品、顾客偏好、季节、居

2、民收入水平等众多因素的影响； 化工产品的得率受温度、压力、催化剂、原料配比等因素的影响。 因此需要了解： 哪些因素会对所研究的指标产生显著影响； 这些影响因素在什么状况下可以产生最好的结果。 方差分析就是解决这类问题的一种统计分析方法。 36.1 方差分析概述方差分析概述 某大型连锁超市为研究各种促销方式的效果，选择下属 4 个门店，分别采用不同促销方式，对包装食品各进行了4 个月的试验。 试验结果如下：4超市管理部门希望了解：不同促销方式对销售量是否有显著影响？哪种促销方式的效果最好？ 促销方式 与上年同期相比(%) a1 (广告宣传) 104.8 95.5 104.2 103.0 a2 (

3、有奖销售) 112.3 107.1 109.2 99.2 a3 (特价销售) 143.2 150.3 184.7 154.5 a4 (买一送一) 145.6 111.0 139.8 122.7 【案例【案例1 1】哪种促销方式效果最好】哪种促销方式效果最好？ 影响某化工厂化工产品得率的主要因素是反应温度和催化剂种类。 为研究产品的最优生产工艺，在其他条件不变的情况下，选择了四种温度和三种催化剂，在不同温度和催化剂的组合下各做了一次试验，测得结果如下： 化工产品得率试验(得率：%)5催化剂 温度 b1 b2 b3 a1(60 oc) 66 73 70 a2(70 oc) 81 96 53 a3(

4、80 oc) 97 79 66 a4(90 oc) 79 76 88 【案例【案例2 2】如何确定最优生产工艺】如何确定最优生产工艺 案例 2 要研究的问题 温度是否对该产品的得率有显著影响? 若有显著影响，应将温度控制在什么范围内可使得率最高? 催化剂是否对该产品的得率有显著影响? 若有显著影响，哪种催化剂的效果最好？ 温度和催化剂的不同组合是否对产品得率有显著影响？ 如有显著影响，哪种温度和催化剂的组合可使得率最高？6 记 a, b, c 为试验中状态发生变化的因素， 称因素在试验中所取的不同状态为水平水平。 设因素 a 有 a 个水平，记为 a1, a2, , aa；因素 b 有 b个水

5、平，记为 b1, b2, , bb 等。 若试验中只有一个变动的因素，就称为单因素试验； 若有两个变动的因素，就称为双因素试验； 若有两个以上的变动因素，则称为多因素试验。二二.方差分析的基本假设方差分析的基本假设 设因素 a 在水平在水平 ai 下的某项指标为总体下的某项指标为总体 xi，则假定 xi n( i, 2 )， xi 相互独立 7一一. 方差分析的基本概念方差分析的基本概念 就是要检验原假设 h0：1 = 2 = = a 是否成立。 若拒绝 h0，就说明因素 a 对试验结果有显著影响，进一步还应确定使效果达到最佳的水平。 若不能拒绝 h0，则说明因素 a 对该项指标无显著影响，试

6、验结果中的差异主要是由其他未加控制的因素和试验误差所引起的。 虽然可以用两两 t 检验法来检验各 i 间是否存在显著差异，但但 t 检验无法检验多个因素间的交互效应检验无法检验多个因素间的交互效应，而这正是方差分析要解决的主要问题。 8三三. .方差分析的目的方差分析的目的一一.基本概念基本概念 记水平 ai 下的 ni 个试验结果为 xij ，则xij = i + iji = 1, 2, a；j = 1, 2, ni ij n(0, 2 )，且相互独立其中 ij 是由各种无法控制的因素引起的随机误差随机误差。 上式说明，试验结果 xij 受到两方面的影响： 因素 a 的水平 ai 的均值 i

7、 随机误差 ij 96.2 6.2 单因素方差分析单因素方差分析为一般平均。 称 i = i - ； i = 1, 2, , a 10 nniii1称) (为试验总次数iinn 从而要检验的原假设可改写为： h0：1= 2 = = a = 0 为水平 ai 的效应效应，反映了水平 xi 的均值与一般平均的差异。11二二. .方差分析的基本方法方差分析的基本方法方差分析 的基本思路：将因素的不同水平和随机误差对试验结果的影响进行分离，并比较两者中哪一个对试验结果 xij 的影响起主要作用。若因素的不同水平对试验结果 xij 的影响是主要的，就拒绝 h0，说明因素 a 对试验结果有显著影响；若试验

8、结果 xij 中的差异主要是由随机误差引起的，就不能拒绝 h0，说明因素 a 对试验结果无显著影响。为此，需要对总的偏差平方和进行分解。 12称jijiixnx1)1(ijijxnx1.1.总的偏差平方和总的偏差平方和ijijtxxs2)(为总的偏差平方和总的偏差平方和， 它反映了样本数据 xij 间总的差异量的大小。为便于对 st 进行分解，记水平 ai 下的样本均值为 2. 偏差平方和的分解13ijijtxxs2)(ijiijexxs2)(iiiaxxns2)(其中iiiijiiijijiijxxnxxxxxx22)()(2()(iiiijiijxxnxx22)()(ijiiijxxxx2

9、)(aess 反映了各样本(同一水平)内的数据差异，主要是由随机误差所引起的，称为误差平方和误差平方和或组内平方和组内平方和。反映了各样本(不同水平)间数据的差异，主要是由因素a的不同水平效应间的差异引起的，称为因素a的平方和的平方和 或 组间平方和组间平方和。利用 sa 和 se 之比就可以构造出检验 h0 的统计量。 故在给定水平 下，若 f f (a-1, n-a)14)/() 1/(ansasfea可以证明， f (a-1, n-a)3. 检验检验 h0 的统计量的统计量当 h0 为真时，统计量就拒绝 h0，说明各水平 ai 的效应间存在显著差异，或称因素 a 的作用是显著的。由于 s

10、a /(a-1) 和 se /(n-a) 分别是组间数据和组内数据的样本方差，故称这种基于检验样本方差比的方法为方差分析方差分析。 若 f f0.001(a-1, n-a)，称因素 a 的作用极高度显著； 若 f f0.01(a-1, n-a)，称因素 a 的作用高度显著； 若 f0.01 (a-1, n-a) f f0.05(a-1, n-a)，称因素 a 的作用一般显著； 若 f f0.05(a-1, n-a)，则认为因素 a 的作用不显著。 15)/() 1/(ansasea来源平方和自由度均方和f 比显著性asaa-1sa/(a-1)误差sen-ase/(n-a)总和stn-1单因素方

11、差分析表三三. .方差分析表方差分析表 2. 因素因素 a 的各水平间是否都存在显著差异？的各水平间是否都存在显著差异？ 这可以通过对各 i 进行两两 t 检验来确定。 如果存在多个效果最好的水平(它们间并无显著差异)，就可以有多种最优水平可供决策者选择。 16ix 若因素 a 的作用是显著的，接下来还需要确定： 1. 因素因素 a 的哪一水平效果最好？的哪一水平效果最好？四四. .进一步的分析进一步的分析这可以通过求出各来确定， 它是 i 的优良估计。 某大型连锁超市为研究各种促销方式的效果，选择下属4个门店，分别采用不同促销方式，对包装食品各进行了4个月的试验。 试验结果如下：17超市管理

12、部门希望了解：不同促销方式对销售量是否有显著影响？哪种促销方式的效果最好？ 促销方式 与上年同期相比(%) a1(广告宣传) 104.8 95.5 104.2 103.0 a2(有奖销售) 112.3 107.1 109.2 99.2 a3(特价销售) 143.2 150.3 184.7 154.5 a4(买一送一) 145.6 111.0 139.8 122.7 【案例【案例1 1】哪种促销方式效果最好】哪种促销方式效果最好？ 可用 excel 的【工具】“数据分析”“方差分析：单因素方差分析”求解单因素方差分析问题。 案例 1 的方差分析表18 其中：p-value p 值，为检验中达到的

13、显著性水平，其含义与 t 检验中“p(t=t)单尾”相同。 f crit 在水平 （默认0.05）下拒绝域的临界值 f。 p-value = 0.00014 f ( a-1，(a-1)(b-1) ) 时，拒绝 h01； 当 fb f ( b-1，(a-1)(b-1) ) 时，拒绝 h02。 25) 1)(1/() 1/(basbsfebb同样可以证明：当 h01 为真时，统计量) 1)(1/() 1/(basasfeaa f( a-1, (a-1)(b-1) )当 h02 为真时，统计量 f(b-1, ( a-1)(b-1) )4. 检验用的统计量检验用的统计量 无交互作用的双因素方差分析表2

14、6来源 平方和 自由度 均方和 f 比 a sa a-1 sa /(a-1) b sb b-1 sb /(b-1) 误差 se (a-1)(b-1) se /(a-1)(b-1) 总和 st ab-1 ) 1)(1/() 1/(basasea) 1)(1/() 1/(basbseb5. 方差分析表方差分析表 影响某化工厂化工产品得率的主要因素是反应温度和使用的催化剂种类。为研究产品的最优生产工艺条件，在其他条件不变的情况下，选择了四种温度和三种催化剂，在不同温度和催化剂的组合下各做了一次试验，测得结果如下： 化工产品得率试验(%)27催化剂 温度 b1 b2 b3 a1 (60 oc) 66

15、73 70 a2 (70 oc) 81 96 53 a3 (80 oc) 97 79 66 a4 (90 oc) 79 76 88 【案例【案例2 2】如何确定最优生产工艺？】如何确定最优生产工艺？ 可用 excel 【工具】“数据分析”“方差分析：无重复双因素分析”求解无交互作用的双因素方差分析问题。 案例2的方差分析表28来源 平方和 自由度 均方和 f 比 p-value a 250 3 83.3 0.42 0.743 b 360.5 2 180.3 0.92 0.449 误差 1179.5 6 196.6 总和 1790 11 因素 a 的 p-value = 0.743 0.05 因

16、素 b 的 p-valu = 0.49 0.05 故温度和催化剂对该化工产品的得率都无显著影响!? 案例案例 2 分析分析29 显然，很可能是由于未考虑因素间的交互作用，才导致了错误的分析结果。 问题出在哪里？问题出在哪里？催化剂 温度 b1 b2 b3 a1 (60 oc) 66 73 70 a2 (70 oc) 81 96 53 a3 (80 oc) 97 79 66 a4 (90 oc) 79 76 88 以上结论既不符合实际情况， 也违背化学反应的基本常识。由本案例的试验数据可以明显看出，温度和催化剂不同搭配下的得率之间是存在显著差异的。 1. 考虑交互作用时的双因素试验考虑交互作用时

17、的双因素试验 交互效应交互效应 记 a, b 间的交互作用为 ab， , i , j 的定义同前。 由于存在交互作用，因此 ij + i + j 称 ( )ij = ij - - i - j为 ai 与 bj 的交互效应交互效应，它反映了两因素间不同水平的组合对试验结果的影响。 因此，总体 xij 的均值可表示为 ij = + i + j + ( )ij 30二二. .考虑交互作用的双因素方差分析考虑交互作用的双因素方差分析(2)要检验的假设 由于考虑了交互作用，因此要检验的原假设有以下三个： h01：1= 2 = = a = 0 h02：1= 2 = = b = 0 h03：( )ij =

18、0；对一切 i, j 为检验交互作用，就需要在每一 ai bj 水平组合下进行重复试验。 记 xijk 为在 ai bj 组合下的第 k 次试验的结果。 31 其中 se , sa , sb 的含义同前； sab 则主要是由交互效应所引起的差异，称为 ab 间平间平方和方和。 32ijkijktxxs2)( 完全类似地，可以将总的偏差平方和 st 作如下分解：2.2.偏差平方和的分解偏差平方和的分解babaessss2. .1112. .11()()abnijija bijlabijijijsxxxxnxxxx故 h01, h02, h03 的拒绝域分别为： fa f( a-1，ab(n-1)

19、 ) fb f( b-1，ab(n-1) ) fab f( (a-1)(b-1)，ab(n-1) ) 33当 h02为真时，) 1(/) 1/(nabsbsfebb 同样可以证明：当 h01为真时，) 1(/) 1/(nabsasfeaa f( a-1, ab(n-1) ) 当 h03为真时，) 1(/) 1)(1/(nabsbasfebaba f( (a-1)(b-1), ab(n-1) )f( b-1, ab (n-1) )3.3.检验用的统计量检验用的统计量 有交互作用的双因素方差分析表34来源 平方和 自由度 均方和 f 比 a sa a-1 sa /(a-1) b sb b-1 sb

20、 /(b-1) ab sab (a-1)(b-1) sab /(a-1)(b-1) 误差 se ab(n-1) se /ab(n-1) 总和 st abn-1 ) 1(/) 1/(nabsasea) 1(/) 1/(nabsbseb) 1(/) 1)(1/(nabsbaseba4.4.方差分析表方差分析表 影响某化工厂产品得率的主要因素是反应温度和使用的催化剂种类。为研究产品的最优生产工艺条件，在其他条件不变的情况下，选择了四种温度和三种催化剂，在不同温度和催化剂的组合下各做了两次试验，测得结果如下： 化工产品得率试验(%)35催化剂 温度 b1 b2 b3 a1 (60 oc) 66, 58 73, 68 70, 65 a2 (70 oc) 81, 79 96, 97 53, 55 a3 (80 oc) 97, 95 79, 69 66, 56 a4 (90 oc) 79, 71 76, 56 88, 82 【案例【案例2 2】如何确定最优生产工艺？】如何确定最优生产工艺？ 可用 excel 【