试验设计与数据处理-李云雁-全套_第1页
试验设计与数据处理-李云雁-全套_第2页
试验设计与数据处理-李云雁-全套_第3页
试验设计与数据处理-李云雁-全套_第4页
试验设计与数据处理-李云雁-全套_第5页
已阅读5页,还剩441页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、试验设计与数据处理(第二版)Experiment Design and Data Processing0.1试验设计与数据处理的发展概况 20世纪2吟代 英国生物统计学家及数学家费歇 (R. A. Fisher)提出了方差分析 20世纪50年代,日本统计学家田口玄一将试验设计中应用 最广的正交设计表格化数学家华罗庚教授也在国内积极倡导和普及的"优选法” 我国数学家王元和方开泰于1978年首先提出了均匀设计0.2试验设计与数据处理的意义0.2.1试验设计的目的:合理地安排试验,力求用较少的试验次数获得较好结果例:某试验研究了3个影响因素:A: Al, A2, A3B: Bl, B2,

2、B3C: Cl, C2, C3全面试验:27次正交试验:9次0.2.2数据处理的目的 通过误差分析,评判试验数据的可靠性;确定影响试验结果的因素主次,抓住主要矛盾,提高试验效率;确定试验因素与试验结果之间存在的近似函数关系,并 能对试验结果进行预测和优化;试验因素对试验结果的影响规律,为控制试验提供思 路;确定最优试验方案或配方O第1章试验数据的误差分析误差分析(eiTOT analysis):对原始数据的可靠性进 行客观的评定误差(ewor):试验中获得的试验值与它的客观真实 值在数值上的不一致试验结果都具有误差,误差自始至终存在于一切科学实验过程中 客观真实值值1.1真值与平均值1.1.1

3、 真值(true mine)真值:在某一时刻和某一状态下,某量的客观值或实际值 真值一般是未知的相对的意义上来说,真值又是已知的平面三角形三内角之和恒为180。 国家标准样品的标称值 国际上公认的计量值 高精度仪器所测之值多次试验值的平均值1.1.2 平均值(mean)(1) 算术平均值(arithmetic mean)ni=l适合:等精度试验值试验值服从正态分布(2)加权平均值(wei欧ted mean)加权和/7jy “KI吗石+ %龙2+色也一台 / f 応& 妙+ % +叱吕1-乂 W,1=适合不同试验值的精度或可靠性不一致时(3) 对数平均值(logarithmic mean

4、)设两个数:形0,七0,则In石-nx2说明: 若数据的分布具有对数特性,则宜使用对数平均值 对数平均值W算术平均值 如果1/2W场化W2时,可用算术平均值代替(4) 几何平均值(geometric mean)设有n个正试验值:咼,七,竝,则JC ° 兀2当_组试验值取对数后所得数据的分布曲线更加对称 时,宜采用几何平均值。几何平均值W算术平均值(5) 调和平均值(hsinnonic mean)设有n个正试验值:和 花 竝,贝11:1 1 1I石龙2占77常用在涉及到与一些量的倒数有关的场合 调和平均值W几何平均值W算术平均值1.2误差的基本概念1.2.1 绝对误差(absolute

5、error)(1)定义绝对误差=试验值一真值或 Ax=尢一£(2)说明真值未知,绝对误差也未知 可以估计出绝对误差的范围:绝对误差限或绝对误差上界|Ajt| =x - xt xf a x±|Ar|1II max绝对误差估算方法:最小刻度的一半为绝对误差;最小刻度为最大绝对误差;根据仪表精度等级计算:绝对误差=量程X精度等级1.2.2 相对误差(relativeerror)(1)定义:相对误差二绝对误差真值或(2)说明:Er =才xtX - xtXX真值未知,常将乂与试验值或平均值之比作为相对误差:XXEr -Ar¥X可以估计出相对误差的大小范围:ErjtAx&qu

6、ot;相对误差限或相对误差上界 = xi±KI)相对误差常常表示为百分数()或千分数(0)1.2.3 算术平均误差(average discrepancy)定义式:E k - XA =<试验值七与算术平均值示之间的偏差可以反映一组试验数据的误差大小1.2.4 标准误差(standard error)当试验次数无穷大时,总体标准差:试验次数为有限次时,样本标准差:表示试验值的精密度,标准差(,试验数据精密度f1.3试验数据误差的来源及分类1.3.1 随机误差(random error )(1)定义:以不可预知的规律变化着的误差,绝对误差时正时负,时大时小(2) 产生的原因:偶然因

7、素(3) 特点:具有统计规律 小误差比大误差出现机会多 正、负误差出现的次数近似相等 当试验次数足够多时,误差的平均值趋向于零 随机误差不可完全避免的可以通过增加试验次数减小随机误差1.3.2 系统误差(systenurtic error )(1) 定义:一定试验条件下,由某个或某些因素按照某一 确定的规律起作用而形成的误差(2) 产生的原因:多方面(3) 特点:系统误差大小及其符号在同一试验中是恒定的 它不能通过多次试验被发现,也不能通过取多次试验值的平均值而减小 只要对系统误差产生的原因有了充分的认识,才能对它进行校正,或设法消除。1.3.3 过失误差(mistake )(1)定义:一种显

8、然与事实不符的误差(2)产生的原因:实验人员粗心大意造成(3)特点:可以完全避免没有_定的规律14试验数据的精准度1.4.1 精密度(precision)(1) 含义:反映了随机误差大小的程度在一定的试验条件下,多次试验值的彼此符合程度例:甲:11.45, 11.46, 11.45, 11.44乙:11.39, 11.45, 11.48, 11.50(2) 说明: 可以通过增加试验次数而达到提高数据精密度的目的 试验数据的精密度是建立在数据用途基础之上的 试验过程足够精密,则只需少量几次试验就能满足要求(3)精密度判断 极差(range) 标准差(standaTd error)£ 3

9、-才)2J 曰-111 1=1 /=1”7ZF也(兀-劝2 1上1上1177-177-1标准差I ,精密度t 方差(variance)标准差的平方: 样本方差(s2 ) 总体方差(0 2 ) 方差(,精密度t1.4.2 正确度(correctness)(1)含义:反映系统误差的大小(2)正确度与精密度的关系:(c)精密度高并不意味着正确度也高精密度不好,但当试验次数相当多时,有时也会得到好 的正确度1.4.3 准确度(accuracy)(1)含义:反映了系统误差和随机误差的综合表示了试验结果与真值的一致程度(2)三者关系无系统误差的试验EM竽购宜一/ J精密度: 正确度: 准确度:A>B

10、>CA = B = CA>B>C有系统误差的试验精密度:A,> B,> 准确度:Az > Bz >Cz , A' >B, C1.5试验数据误差的统计假设检验1.5.1随机误差的检验1.5.1.1力检验(% 2-test)(1)目的:在试验数据的总体方差/已知的情况下, 对试验数据的随机误差或精密度进行检验0(2)检验步骤: 计算统计量才若试验数据心心心服从正态分布,则?(77 _ 1)/犷=;一b服从自由度为df=n 的旷分布 查临界值加(如a显著性水平一般取0.01或005,表示有显著差异的概率 检验双侧(尾)检验(twodded/ta

11、iled test):2 2 2若窃s < % < z«r?J亍则判断两方差无显著差异,否则有显著差异单侧(尾)检验(ongsided/tailed test):>左侧(尾)检验:若 Z2>Z(l_a)()则判断该方差与原总体方差无显著减小,否则有显著减小 >右侧(尾)检验若 Z2 </«(#)则判断该方差与原总体方差无显著增大,否则有显著增大(3) Excels /2检验中的应用1.5.1.2 F(F-test)(1)目的:对两组具有正态分布的试验数据之间的精密度进行比较(2)检验步骤 计算统计量设有两组试验数据:石£和坷丿&

12、#163;心(2 都月艮从正态分布,样本方差分别为彳和彳,则9 F=服从F分布,第一自由度为df = n、- 第二自由度为必 查临界值给定的显著水平a、必=勺_1 检验双侧(尾)检验(two«sided/tailed test):(1T)T则判断两方差无显著差异,否则有显著差异单侧(尾)检验(oiiesided/tailed test):左侧(尾)检验:若刀 %一a)(劣必)则判断该判断方差1比方差2无显著减小,否则有显著减小右侧(尾)检验若刀巧(勿,必)则判断该方差1比方差2无显著增大,否则有显著增大(3) Excel在F检验中的应用1.5.2系统误差的检验1.5.2.1媳验法(1

13、) 平均值与给定值比较 目的:检验服从正态分布数据的算术平均值是否与给定值 有显著差异 检验步骤:计算统计量:服从自由度莎* = ”-1的t分布(t-distribution)给定值(可以是真值、期望值或标准值)双侧检验:若 kl 4则可判断该平均值与给定值无显著差异,否则就有显著差异 单侧检验左侧检验若/<0则判断该平均值与给定值无显著减小,否则有显著减小右侧检验若/>0则判断该平均值与给定值无显著增大,否则有显著增大(2) 两个平均值的比较目的:判断两组服从正态分布数据的算术平均值有无显著差异 俯算统计量: 两组数据的方差无显著差异时服从自由度df +772 -2的t分布s并标

14、准差:(坷 _1)£ +(2 _1E* 召 +_ 2两组数据的精密度或方差有显著差异时服从t分布,其自由度为:居I恋 I &/2)2(1+1)(z+l)双侧检验:若 W2则可判断两平均值无显著差异,否则就有显著差异 单侧检验左侧检验若/<0则判断该平均值1较平均值2无显著减小,否则有显著减小右侧检验若/>0则判断该平均值1较平均值2无显著增大,否则有显著增大(3) 成对数据的比较目的:试验数据是成对出现,判断两种方法、两种仪器或 两分析人员的测定结果之间是否存在系统误差俯算统计量:服从自由度为df=n 的t分布零或其他指定值工(兀-牙)成对测定值之差的算术平均值:

15、刁=Sdn对试验值之差值的样本标准差:Sd =n_7=177-1nt检验若ki < 4 ,则成对数据之间不存在显著的系统误差,9否则两组数据之间存在显著的系统误差(4) Excel在媳验中的应用1.5.2.2 秩和检验法(rank sum test)(1)目的:两组数据或两种试验方法之间是否存在系统误差、两种方法是否等效等,不要求数据具有正态分布(2)内容:设有两组试验数据,相互独立,m,112分别是两组数据的个数,总假定nDj; 将这个试验数据混在一起,按从小到大的次序排列 每个试验值在序列中的次序叫作该值的秩(rank)将属于第1组数据的秩相加,其和记为RRi第1组数据的秩和(ran

16、k sum)如果两组数据之间无显著差异,则R】就不应该太大或太小根据显著性水平a和!1】,1 检验:旳,可查得Ri的上下限T2和T如果Ri>T2或R1VT1,则认为两组数据有显著差异,另 一组数据有系统误差如果TVRVT2,则两组数据无显著差异,另一组数据 也无系统误差(3)例:设甲、乙两组测定值为:甲:86, 10.0, 9.9, 88, 9.1, 9.1乙:8.7, 84, 9.2, 8.9, 7.4, 80, 7.3, 8.1, 6.8已知甲组数据无系统误差,试用秩和检验法检验乙组测定值是否有系统误差。(a=005)解:(1)排序:秩1234567甲8.6乙6.87.37.48.0

17、8.18.4891011.511.51314158.89.19.19.910.08.78.99.2(2)求秩和RiRi=7+9+11.5+11.5+14+15=68(3)查秩和临界值表对于a=0.05, ni=6, H2=9 得 Tx=33, T2=63,R> T2故:两组数据有显著差异,乙组测定值有系统误差1.5.3异常值的检验可疑数据、离群值、异常值般处理原则为:在试验过程中,若发现异常数据,应停止试验,分析原 因,及时纠正错误试验结束后,在分析试验结果时,如发现异常数据,则应 先找出产生差异的原因,再对其进行取舍在分析试验结果时,如不清楚产生异常值的确切原因,则 应对数据进行统计处

18、理;若数据较少,则可重做一组数据对于舍去的数据,在试验报告中应注明舍去的原因或所选用的统计方法1.5.3.1拉依达JP麵ta )检验法 内容:可疑数据若x 一 才 > 3 s 或 2 sP则应将该试验值剔除。 说明: 3s相当于显著水平a=001,2s相当于显著水平近005 计算平均值及标准偏差s时,应包括可疑值在内可疑数据应逐一检验,不能同时检验多个数据首先检验偏差最大的数剔除一个数后,如果还要检验下一个数,应重新计算平 均值及标准偏差方法简单,无须查表该检验法适用于试验次数较多或要求不高时 3s为界时,要求n>102s为界时,要求n>5 例:有一组分析测试数据:0.128

19、, 0.129, 0.131, 0.133,0.135, 0.138, 0.141, 0.142, 0.145, 0.148,0.167,问其中偏差较大的0167这一数据是否应被舍去? (a=0e01)解:(1)计算元= 0.140 朋= 0.01116(2) 计算偏差-T|=|0.167 -0.140|=0,027(3) 比较3s=3 X 0.01116=0.0335>0.027故按拉依达准则,当a=0.01时,0.167这一可疑值不应舍去(2)格拉布斯(Grubbs)检验法 内容: 可疑数据弓,若I讣丰护则应将该值剔除。G©申) mbbs检验临界值格拉布斯(Grubbs)检

20、验临界值G(a,n)表n显著性水平ct0.050.0250.010.00531.1531.1551.1551.15541.4631.4811.4921.49651.6721.7151.7491.76461.8221.8871.9441.97371.9382.0202.0972.13982.0322.1262.2212.27492.1102.2152.3232.387102.1762.2902.4102.482112.2342.3552.4852.564122.2852.4122.5502.636132.3312.4622.6072.699142.3712.5072.6592.755152.40

21、92.5492.7052.806162.4432.5852.7472.252172.4752.6202.7852.894 说明:计算平均值及标准偏差S时,应包括可疑值在内可疑数据应逐_检验,不能同时检验多个数据首先检验偏差最大的数剔除一个数后,如果还要检验下一个数,应重新计算平均值及标准偏差能适用于试验数据较少时格拉布斯准则也可以用于检验两个数据偏小,或两个数据偏大的情况(3)狄克逊(Dixon)检验法单侧情形将n个试验数据按从小到大的顺序排列:如果有异常值存在,必然出现在两端,即场或竝计算出统计量型2查单侧临界值D,_a(77)检验A检验竝时,当DD3时,可剔除竝双侧情形计算D和D,査双侧临

22、界值检验>当D>D'、D>判断百为异常值> 当 D' > D、D' > ZA_aO) 判断为异常值 说明 适用于试验数据较少时的检验,计算量较小单侧检验时,可疑数据应逐_检验,不能同时检验多个数 据剔除一个数后,如果还要检验下一个数,应重新排序1.6有效数字和试验结果的表示1.6.1 有效数字(significancefigure)能够代表一定物理量的数字有效数字的位数可反映试验或试验仪表的精度数据中小数点的位置不影响有效数字的位数例如:50IIIII,0.050m, 5.0 x 104» m第一个非0数前的数字都不是有效数

23、字,而第一个耳罚数后 的数字都是有效数字例如:29mm和29.00mm第一位数字等于或大于8,则可以多计一位例如:9.991.6.2有效数字的运算(1) 加、减运算:与其中小数点后位数最少的相同(2) 乘、除运算以各乘、除数中有效数字位数最少的为准(3)乘方、开方运算:与其底数的相同:例如:2.42=5.8(4)对数运算:与其真数的相同例如ln684 = 1.92; lgO.00004 = -4(5)在4个以上数的平均值计算中,平均值的有效数字可增加一位(6) 所有取自手册上的数据,其有效数字位数按实际需要 取,但原始数据如有限制,则应服从原始数据。(7) 一些常数的有效数字的位数可以认为是无

24、限制的例如,】周率兀.重力加速度g. . 1/3等(8) 一般在工程计算中,取23位有效数字1.6.3有效数字的修约规则 W4:舍去 M5,且其后跟有非零数字,进1位例如:3.14159 f 3.142=5,其右无数字或皆为0时尾留双”:若所保留的末位数字为奇数则进1若所保留的末位数字为偶数则舍弃例如:3.1415 - 3.14213665 - 1.3661.7误差的传递误差的传递:根据直接测量值的误差来计算间接测量值的 误差1.1.1误差传递基本公式间接测量值卩与直接测量值不之间函数关系:歹=/(兀1,2忑)全微分站 dx2A厂 艺乂+必X+.+巫北dxY dx2 dxH函数或间接测量值的绝

25、对误差为:y二£7=1dx.&相对误差为:Ar v7=5dx. yf 误差传递系数Axz.直接测量值的绝对误差;V 间接测量值的绝对误差或称函数的绝对误差。函数标准误差传递公式:1.7.2常用函数的误差传递公式表141.73误差传递公式的应用(1)根据各分误差的大小,来判断间接测量或函数误差的主要来源:例116(2)选择合适的测量仪器或方法:秩和临界值表口1n20?=0.025a=0.05Tit3Tit2551837193661941204072045224382149234792253255010245626546626522850728563054829613258931

26、6533631033693567统计量D计算公式n检验高端异常值检验低端异常值37厂_巧一巧可一石810厂_巧一巧-12D'=丕一卞巧一1P11 13门兀一兀-2小巧一石£一1 -石14 30第2章 试验数据的表图表示法2.1列表法将试验数据列成表格,将各变量的数值依照一定的形式和顺序一对应起来(1)试验数据表记录表 试验记录和试验数据初步整理的表格 表中数据可分为三类: 原始数据中间数据最终计算结果数据结果表不表表达试验结论应简明扼要表乙2离心泵特性曲线测定实验结果表示表序号流量g訂(甘压头兄Jm轴功率马/W数率和(2)说明:三部分:表名、表头.数据资料必要时,在表格的下方

27、加上表外附加表名应放在表的上方,主要用于说明表的主要内容,为了 引用的方便,还应包含表号表头常放在第一行或第一列,也称为行标题或列标题,它 主要是表示所研究问题的类别名称和指标名称数据资料:表格的主要部分,应根据表头按一定的规律排表外附加通常放在表格的下方,主要是一些不便列在表内 的内容,如指标注释.资料来源.不变的试验数据等(3)注:表格设计应简明合理.层次清晰,以便阅读和使用;数据表的表头要列出变量的名称、符号和单位;要注意有效数字位数;试验数据较大或较小时,要用科学记数法来表示,并记入 表头,注意表头中的与表中的数据应服从下式:数据的实 际值=表中数据;数据表格记录要正规,原始数据要书写

28、得清楚整齐,要记录各种试验条件,并妥为保管。2.2图示法2.2.1常用数据图(1) 线图(line graph/chart)表示因变量随自变量的变化情况 线图分类: 单式线图:表示某一种事物或现象的动态 复式线图:在同一图中表示两种或两种以上事物或现象的动态,可用于不同事物或现象的比较o o o n- o o o Q Q o o0 9 8 7 6 5 4 3 2 11高吸水性树脂保水率与时间和温度的关系流量 qv (L/s)2某离心泵特性曲线(2) XY散点图(scsitterdiagram)表示两个变量间的相互关系 散点图可以看出变量关系的统计规律2. 52.01. 51.00. 50246

29、810(3) 条形图和柱形图湿浸法 碱提法 醇提法超声波法提取方法图4不同提取方法提取率比较12108 6 4、2X用等宽长条的长短或高低来表示数据的大小,以反映各数据点的差异 两个坐标轴的性质不同A数值轴:表示数量性因素或变量 分类轴:表示的是属性因素或非数量性变量分类:单式:复式:超声波法醇提法碱提法湿浸法只涉及一个事物或现象涉及到两个或两个以上的事物或现象|植物2|I植物图5不同提取方法对两种原料有效成分提取率效果比较(4)圆形图和环形图圆形图(circle chart)也称为饼图(pie graph)表示总体中各组成部分所占的比例只适合于包含一个数据系列的情况饼图的总面积看成100%

30、,每8% 1%药品食品圈化妆品其它36。圆心角所对应的面积为1% ,以扇形面积的大小来分别表示各项的比例环形图(ckcular diagram)7全球合成、天然维生素E消费比例比较每一部分的比例用环中的一段表示可显示多个总体各部分所占的相应比例,有利于比较(5) 三角形图(ternary)常用于表示三元混合物各组分含量或浓度之间的关系三角形:等腰RtZJ等边、不等腰Rt等 顶点:纯物质边:二元混合物三角形内:三元混合物A"s)8等腰直角三角形坐标图39等边三角形坐标图(6) 三维表面图(3D surface graph)三元函数巧(X,Y)对应的曲面图,根据曲面图可以看出因 变量Z值

31、随自变量X和Y值的变化情况XI.6O图10三维表面图.4 O-0.4 -06 -D.8-0.0(7) 三维等高线图(contour plot)三维表面图上Z值相等的点连成的曲线在水平面上的投影绘制图形时应注意:(1) 在绘制线图时,要求曲线光滑,并使曲线尽可能通过较多的实验点,或者使曲线以外的点尽可能位于曲线附近,并使曲线两侧的点数大致相等;(2) 定量的坐标轴,其分度不一定自零起;(3) 定量绘制的坐标图,其坐标轴上必须标明该坐标所代表 的变量名称、符号及所用的单位,一般用纵轴代表因变量(4) 坐标轴的分度应与试验数据的有效数字位数相匹配;(5) 图必须有图号和图题(图名),以便于引用,必要

32、时还 应有图注。2.2.2坐标系的选择坐标系(coordinate system)笛卡尔坐标系(又称普通直角坐标系)、半对数坐标系、对数坐标系、极坐标系、概率坐标系、三角形坐标系对数坐标系(semi-logarithmic coordinate system)半对数坐标系双对数坐标系(1)选用坐标系的基本原则: 根据数据间的函数关系 线性函数:普通直角坐标系 麻函数:双对数坐标系指数函数:半对数坐标 根据数据的变化情况两个变量的变化幅度都不大,选用普通直角坐标系;有一个变量的最小值与最大值之间数量级相差太大时,可以选用半对 数坐标;两个变量在数值上均变化了几个数量级,可选用双对数坐标; 在自变

33、量由零开始逐渐增大的初始阶段,当自变量的少许变化引起因变量极大变化时,此时釆用半对数坐标系或双对数坐标系,可使图形轮廓清楚例:X10204060801001000200030004000y2414608010017718118820012普通直角坐标系100013对数坐标系(2)坐标比例尺的确定IIIIIMill在变量讶0/的误差禺丿已知时,比例尺的取法应使试 验"点"的选长为2禺2A/,而且使2厶尸2丿=12 ,若2 丿=2mm,则y轴的庄例尺NyS为: 推荐坐标轴的比例常数M= (1> 2、5) XIO*0 (n为正整2mm2 Ay=(rnrnl y)数),而3、

34、6、7、8等的比例常数绝不可用; 纵横坐标之间的比例不一定取得一致,应根据具体情况选 择,使曲线的坡度介于30。60。之间1=1例2:研究pH值对某溶液吸光度A的影响,已知pH值的测量 误差pH=Ol,吸光度A的测量误差AAuOei。在一定波 长下,测得pH值与吸光度A的关系数据如表所示。试在普通 直角坐标系中画出两者间的去系曲线。PH8.09.010.011.0吸光度A1.341.361.451.36解:ApH = 2AA = 2mmV ApH=O.l, AA=O.O1横轴的比例尺为MpH =2mm 1mmTApH 0.2= 10(刃加/单位/刃值)纵轴的比例尺为亿=筈二誥=】0。(林单位吸

35、光度)pH14坐标比例尺对目8.09.010.011.0pH12.0g形形状的影响23计算机绘图软件在图表绘制中应用2.3.1 Excel在图表绘制中的应用(1) 利用Excel生成图表的基本方法(2) 对数坐标的绘制(3) 双Y轴(X轴)复式线图的绘制(4) 图表的编辑和修改2.3.2 Origin在图形绘制中的应用(1) 简单二维图绘制的基本方法(2) 三角形坐标图的绘制(3) 三维图的绘制表21离心泵特性曲线测定实验的数据记录表流量计读数真空表读数压力表读数/MPq功率表读数/(L/h)/MPaAV附:泵入口管径:mm;泵出口管径: mm;真空表与压力表垂直距离:mm;水温:°

36、C;电动机转速 r/mnio表水的物理性质温度饱和蒸汽压导热系数粘度表面张力几1叫(師府1矿1)“X12"皿咱)crxlO4/(M/酬)00.0061155.11788756.4300.4241061.8801.5712.2600.19920659469.9662.2第3章试验的方差分析方差分析(analysis of wrismce,简称ANOVA)检验试验中有关因素对试验结果影响的显著性试验指标(experimental index)衡量或考核试验效果的参数 因素 (experimental factor)影响试验指标的条件 可控因素(controllsible factor)水

37、平(level of factor)因素的不同状态或内容3.1单因素试验的方差分析(one-way analysis of variance)3.1.1单因素试验方差分析基本问题(1) 目的:检验一个因素对试验结果的影响是否显著性(2)基本命题:设某单因素A有I种水平:AP A2, Ar,在每种水平下 的试验结果服从正态分布在各水平下分别做了勿(z=l, 2,,尸)次试验判断因素A对试验结果是否有显著影响试验次数A12码2花2 J %"i A/ A, xu %3.1.2单因素试验方差分析基本步骤(1)计算平均值 组内平均值:总平均:(2) 计算离差平方和 总离差平方和皿(sum of

38、 squares for total)尸 /_/=i j= 表示了各试验值与总平均值的偏差的平方和 反映了试验结果之间存在的总差异 组间离差平方和 SSA (sum of square for factor A)弟=S -才)2 = 士32z=l y=lz=l 反映了各组内平均值之间的差异程度 由于因素A不同水平的不同作用造成的 组内离差平方和SSe (sum of square for error)r Y_z=l 7=1反映了在各水平内,各试验值之间的差异程度由于随机误差的作用产生三种离差平方和之间关系:(3) 计算自由度(degree of freedom)总自由度:dff=n组间自由度:

39、弧=L 组内自由度:df=nr三者关系:勿=必+必(4)计算平均平方均方=离差平方和除以对应的自由度MSA=SSA!dfA MSeSSJdfeMS、组间均方MSy组内均方/误差的均方(5)冷验组间均方二MS人组内均方顾服从自由度为(dfpdf)的F分布(F distribution)对于给定的显著性水平 从F分布表查得临界值丘(必,勿 如“ >/;«, %),则认为因素A对试验结果有显著衫响否则认为因素A对试验结果没有显著影响(6)方差分析表单因素试验的方差分析表差异源SSdfMSF显著性组间(因素A)SSAr-1MSa=SSa!(r-1)MSA/MSe组内(误差)ssenrM

40、S=SSe!(»r)总和SStzr1若屁 坷01(必,00,称因素A对试验结果有非常显著的影响,用号表示;若均05(必? 09 F 人01(必? 09,则因素A对试验结果有显著的影响,用3"号表示;若F" Fg凤、豕),则因素A对试验结果的影响不显著3.1.3 Excel在单因素试验方差分析中的应用利用Excel “分析工具库”中的“单因素方差分析”工具3.2双因素试验的方差分析讨论两个因素对试验结果影响的显著性,又称,二元方差分析3.2.1双因素无重复试验的方差分析(1)双因素无重复试验Bi Bsn-12 22 Ar片1-r2 TS(2)双因素无重复试验方差分析

41、的基本步骤 计算平均值总平均:TEAX/. - V X . s y=l巧水平时: 计算离差平方和总离差平方和:F s _ 2隅辽£仇-刁=隔+鸥+咫7=1 7=1因素A弓|起离差的平方和:r s_sse = n (形 _ J;/. _ 才丿 + j;)2 z=l 7=1因素B引起离差的平方和:误差平方和: 计算自由度 SSa的自由度:么=r1D SSg的自由度:dfs 1 SSe 的自由度:df=(7-1) (y-1) S&T的自由度:dfr=nl.=rs妇必+%+% 计算均方SSe(尸一 1)(一 1) 胸验MSe血服从自由度为(必o的F分布;人服从自由度为(确,勿的F分布

42、;对于给定的显著性水平a ,查F分布表:Fa (勿,09' Fa (硯,必)若FA>Fa (必,豕),则因素A对试验结果有显著影响,否 则无显著影响;若FB>Fa (必O,则因素B对试验结果有显著影响,否 则无显著影响; 无重复试验双因素方差分析表无重复试验双因素方差分析表差异源dfMSF显著性因素Ar1MSA = S”=丝A M»因素BSSbs1MS启B s-1Fb = MSe误差sse(尸-1)(一1)(尸-1)($-1)总和SST/y13.2.2双因素重复试验的方差分析(1)双因素重复试验方差分析试验表双因素重复试验方差分析试验表因素 BsA】石11 眄 1

43、2“眄元石22,,石2r lP £12,宀lc尢?? X2s2s2'X2sc Ar形1” £12,Xrc形21,巧22 x ,山乜c rsX'1 rsl> > rsc(2)双因素重复试验方差分析的基本步骤计算平均值总平均:1rscrscEESz=l 丿=1任一组合水平z巧)上:X ijk A冰平时: Bj水平时:计算离差平方和 总离差平方和:7- = Z X X 5 - 7)2 = ssA + SSb + SSg + SSe z=l 7=1 力=1因素A引起离差的平方和:鸥=储工(2-孑7=1因素B引起离差的平方和:SS£ = rc (

44、Ly.-)2交互作用AXB引起离差的平方和:7_1r s _ _EaxB (龙沪_龙八+ X)上1 7=1误差平方和:SSe = ±±± (% -乙)2z=l y=l =1计算自由度 SSa的自由度:dfp r1 SSb的自由度:砖s SSAXB的自由度:XB =(r-l)(J-l) SSe的自由度:df=rsc 1) SSp的自由度:dff=nX=rscX妇必+%+必xb+饥(I朋 a pF%【一"V奩5也備44$)胸验=竺MSeMSe若FA>Fa (必则认为因素A对试验结果有显著影响 ,否则无显著影响;若FB>Fa (滋,勿,则认为因素B

45、对试验结果有显著影响 ,否则无显著影响;若FAXB>Fa (必xb,09,则认为交互作用AXB对试验结 果有显著影响,否则无显著影响。重复试验双因素方差分析表有交互作用双因素试验的方差分析表差异源MSF显著性因素观MS.r-1se巴MSg因素B5-1MS更_ 2>_叭交互作用AXBCr-1)£一1)胚"MS,卜口莘1 天差rs (c 1)MS厂汎rs(c 1)总和SSTrsc 13.2.3 Excel在双因素方差分析中的应用(1)双因素无重复试验方差分析利用“分析工具库”中的“无重复双因素方差分析"工具(2)双因素重复试验方差分析利用“分析工具库”中的“重复双因素方差分析“工具试验数据的回归分析4.1基本

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论