不考虑交互作用的正交试验设计041019200448
正交试验设计表-图
选择合适的正交表
根据因素和水平数量,选择合适 的正交表,确保试验具有足够的 代表性和均衡性。
制定试验计划
按照正交表的安排,制定详细的 试验计划,包括试验顺序、条件 控制等。
试验因素与水平的确定
确定关键因素
01
根据研究目的和问题背景,确定对试验结果有显著影响的因素
作为关键因素。
确定因素水平
02
为关键因素选择合适的水平,确保试验结果具有实际意义和应
农业科学中的应用
肥料配比优化
通过正交试验设计,可以找出最佳的肥料配比方案,提高农作物的 产量和品质。
农药喷施方案优化
正交试验设计可用于评估不同农药喷施方案对农作物病虫害防治效 果的影响,从而选择出最佳的喷施方案。
种植技术优化
在农业种植技术中,正交试验设计可用于分析不同种植方式对农作物 生长的影响,从而选择出最佳的种植技术。
科学决策依据
正交试验设计得出的最优 解具有一定的可信度,可 以为科学决策提供有力的 依据。
正交试验设计的步骤与流程
明确试验目的和因素
首先需要明确试验的目的和考察的因素,为 后续的试验设计提供基础。
确定水平
根据实际情况和试验要求,为每个因素选择合 适的水平。
选择合适的正交表
根据因素和水平数量,选择合适的正交表,确保 试验设计的均衡分散和整齐可比。
用价值。
控制其他因素
Байду номын сангаас03
在试验过程中,应尽量控制其他非关键因素的干扰,以提高试
验的准确性和可靠性。
试验数据的收集与分析
记录试验数据
在试验过程中,详细记录每个试验条件下的数据,包括试验结果、 误差等。
数据整理与清洗
(完整版)正交实验设计
正交实验设计当析因设计要求的实验次数太多时,一个非常自然的想法就是从析因设计的水平组合中,选择一部分有代表性水平组合进行试验。
因此就出现了分式析因设计(fractional factorial designs),但是对于试验设计知识较少的实际工作者来说,选择适当的分式析因设计还是比较困难的。
正交试验设计(Orthogonal experimental design)是研究多因素多水平的又一种设计方法,它是根据正交性从全面试验中挑选出部分有代表性的点进行试验,这些有代表性的点具备了“均匀分散,齐整可比”的特点,正交试验设计是分式析因设计的主要方法。
是一种高效率、快速、经济的实验设计方法。
日本著名的统计学家田口玄一将正交试验选择的水平组合列成表格,称为正交表。
例如作一个三因素三水平的实验,按全面实验要求,须进行33=27种组合的实验,且尚未考虑每一组合的重复数。
若按L9(3)3正交表按排实验,只需作9次,按L18(3)7正交表进行18次实验,显然大大减少了工作量。
因而正交实验设计在很多领域的研究中已经得到广泛应用。
1.正交表正交表是一整套规则的设计表格,用。
L为正交表的代号,n为试验的次数,t为水平数,c为列数,也就是可能安排最多的因素个数。
例如L9(34),(表11),它表示需作9次实验,最多可观察4个因素,每个因素均为3水平。
一个正交表中也可以各列的水平数不相等,我们称它为混合型正交表,如L8(4×24) (表12),此表的5列中,有1列为4水平,4列为2水平。
根据正交表的数据结构看出,正交表是一个n行c列的表,其中第j列由数码1,2,… S j组成,这些数码均各出现N/S次,例如表11中,第二列的数码个数为3,S=3 ,即由1、2、3组成,各数码均出现次。
正交表具有以下两项性质:(1)每一列中,不同的数字出现的次数相等。
例如在两水平正交表中,任何一列都有数码“1”与“2”,且任何一列中它们出现的次数是相等的;如在三水平正交表中,任何一列都有“1”、“2”、“3”,且在任一列的出现数均相等。
正交试验设计表
在这 9 个水平组合中, A 因素各水平下包括 了 B 、 C 因素的 3个水平,虽然搭配方式不同, 但B、C皆处于同等地位,当比较A因素不同水 平时,B因素不同水平的效应相互抵消,C因素 不同水平的效应也相互抵消。所以A因素3个水 平间具有综合可比性。同样, B 、 C 因素 3个水 平间亦具有综合可比性。
10-3 因素水平表
试验因素 水平 加水量 加酶量 (mL/100g) (mL/100g) A B 10 50 90 1 4 7 酶解温度 (℃) C 20 35 50 酶解时间 (h) D 1.5 2.5 3.5
1 2 3
(3) 选择合适的正交表
正交表的选择是正交试验设计的首要问题。确定了 因素及其水平后,根据因素、水平及需要考察的交互作 用的多少来选择合适的正交表。正交表的选择原则是在 能够安排下试验因素和交互作用的前提下,尽可能选用 较小的正交表,以减少试验次数。 一般情况下,试验因素的水平数应等于正交表中的 水平数;因素个数(包括交互作用)应不大于正交表的 列数;各因素及交互作用的自由度之和要小于所选正交
个水平为宜。对主要考察的试验因素,可以多取水平,但不宜过
多(≤6),否则试验次数骤增。因素的水平间距,应根据专业
知识和已有的资料,尽可能把水平值取在理想区域。
上一张 下一张 主 页 退 出
对本试验分析,影响山楂液化率的因素很多, 如山楂品种、山楂果肉的破碎度、果肉加水量、原 料pH 值、果胶酶种类、加酶量、酶解温度、酶解 时间等等。经全面考虑,最后确定果肉加水量、加 酶量、酶解温度和酶解时间为本试验的试验因素, 分别记作A、B、C和D,进行四因素正交试验,各 因素均取三个水平,因素水平表见表10-3所示。
例如,要考察增稠剂用量、pH值和杀菌温度对豆奶稳
试验设计类型之无法考察交互作用的多因素设计:交叉设计
试验设计类型之无法考察交互作用的多因素设计:交叉设计胡完;胡良平【摘要】本文目的是介绍一种无法考察交互作用的多因素设计类型,即交叉设计.通过详细介绍与交叉设计类型对应的"主要内容""操作方法"和"设计变形",全面展示了这种设计类型的核心内容、要领和注意事项,为用户正确、合理地选用交叉设计类型创造了有利的条件.%The purpose of this paper was to introduce a kind of design type of multifactor, that was a crossover design, which can't be used to study in the situation with interactions among factors.The core contents, essentials and precautions of the crossover design was demonstrated comprehensively through definitely introducing the main contents, operation methods and deformed designs corresponding with the crossover design.The favorable conditions were created for users to choose the crossover design correctly and reasonably.【期刊名称】《四川精神卫生》【年(卷),期】2017(030)002【总页数】4页(P101-104)【关键词】试验因素;区组因素;样本含量;单个体型;多个体型;交叉设计【作者】胡完;胡良平【作者单位】军事医学科学院生物医学统计学咨询中心,北京 100850;军事医学科学院生物医学统计学咨询中心,北京 100850;世界中医药学会联合会临床科研统计学专业委员会,北京 100029【正文语种】中文【中图分类】R195.1【例1】某麻醉科医生研究催醒宁对氟呱啶的作用,做了一个试验,其格式见表1。
实验设计之正交实验
• 三、有交互作用的正交试验及结果分析 • 例1-4 乙酰胺苯磺化反应试验 • 试验目的:提高乙酰胺苯的产率 因素-水平表
因素
水平 1 2
A 反应温度 ℃ B 反应时间 小时 C 硫酸浓度 x(%) D 操作方法
50 70
1 2
17 27
搅拌 不搅拌
考虑反应温度与反应时间可能会有交互作用,另外,反应温度与 硫酸浓度也可能会有交互作用,即考虑A×B、A×C
135 171 144 45 57 48 12
1 2 3 3 1 2 2 3 1
31 54 38 53 49 42 57 62 64
因素次序 A C B 主 次
• 指标越大越好,应该选取指标最大的水平。 • 本试验应该选取每个因素中k1、k2、k3最大的哪 个水平。即:A3B2C2
•
也可以选取图形中最高的水平点得到最优生产条件:
指标
70 60 50 40
同时可以估计,随着A的增加,指标还有向上的趋势
A1 A2 A3
B1
B2 B3
C1
C2 C3
指标-因素图
• 选取原则: • (1)对主要因素,选使指标最好的那个水平 于是本例中A选A3,C选C2 • (2)对次要因素,以节约方便原则选取水平 本例中B可选B2或者B1 • 于是用A3B2C2、A3B1C2各做一次验证试验,结果如下:
• 1、指标、因素和水平 • 试验需要考虑的结果称为试验指标(简称指标) 定量指标:直接用数量表示; 定性指标:不能用数量表示。定性指标可以按评定 结果打分或者评出等级。定性指标一般要进行定 量化。 • 试验中对试验指标可能有影响的变量简称为因素, 用大写字母A、B、C…表示 • 每个因素的不同状态称为因素的水平。
正交试验设计法 ppt课件
PPT课件
21
正交试验设计法
5.6 计算分析试验结果
5.6.2.3.2 图示说明
为直观起见,用因素的水平变化为横坐标,指标的算数平均 值为纵坐标,画出水平与指标图,如图1:
从图中可明显看出最佳方案应为:A3B2C2。而正交试验选出 的最佳方案为A3B3C2,即第9号方案,显然,正交试验中的9个方 案中没有A3B2C2这一方案,其是否为最佳方案,需要通过正式试 验来验证。
◆正交表中,任意两列,每1行组成1个数字对,有多少行
就有多少个这样的数字对,这些数字对都是完全有序的
◆各种数字出现的次数必须相同,这是正交表必须满足的
的两个特性。
◆ 其他正交表如:L8(27)、 L12(211)、L18(37)、L32
(49)、L25(56)等都满足这两个特性。
PPT课件
12Βιβλιοθήκη 正交试验设计法PPT课件
22
正交试验设计法
5.6.2.3.2 图示说明
PPT课件
23
正交试验设计法
5.7 验证性实验
为了与正交试验选出的最佳方案进行对比,用A3B2C2方案和 A3B3C2方案各做一次验证试验,转化率分别为74%和65%,说明 A3B2C2方案实为最佳方案。
上例表明,最佳方案虽然不在正交试验9个方案中,但通过计 算分析即可准确选出,这充分说明了正交试验法的科学性。
10
正交试验设计法
再以L9(34)为例:
表2 L9(34)表
行
项目
1
2
3
4
1
1
1
1
1
2
1
2
2
2
3
1
3
3
3
不考虑交互作用的正交试验设计041019200448
§6.5 不考虑交互作用的正交试验设计下面我们先来看一种比较简单的情形,即不考虑因子之间交互作用的情形。
不考虑交互作用的正交试验设计和数据处理,可按下列步骤进行:(1)选正交表 。
)(mn r L 选择的原则是: 要等于因子的水平数; 要大于或等于因子的个数; 是试验r m n 次数,要尽可能小。
例如,问题中有4个因子,每个因子都是2个水平。
选正交表时,首先要选 =2 的r 表。
这样的正交表有 。
在 中,=3 ,小于因子的 ),2(),2(),2(15167834L L L )2(34L m 个数4,所以不符合要求。
在 中, ,大于因子的个数 ),2(),2(151678L L ,15,7=m 4,符合要求。
在符合要求的正交表中,还要选试验次数 尽可能小的那一个表。
显然,n 在 中,=8 为最小,所以,我们最后选定正交表 。
)2(78L n )2(78L (2)设计表头。
将各因子安排在正交表的各列上方,每个因子占1列,这称为表头。
在不考虑交互作用的正交试验设计中,表头上的因子可以任意安放。
表头上不放因子的列,称为空白列。
(3)按照设计做试验,取得试验观测值。
正交表的每一行代表一种水平组合,对每一种水平组合做一次试验。
按照第 行的k 水平组合所做的第 次试验,所得到的观测值记为 。
正交表有 行,所以,一共k k X n 要做 次试验,共得到 个试验观测值: 。
n n n X X X ,,,21 (4)在正交表的每一列中,求出与各水平对应的均值,以及这一列的平方和。
设我们考虑的是第 列。
在这一列中,表示水平的数字 ,每一个都重j r ,,2,1 复出现 次。
设与这一列中的数字 对应的那几行的试验观测值之和为 ,与这r n 1j W 1一列中的数字 对应的那几行的试验观测值之和为 , ,与这一列中的数字 2j W 2 r 对应的那几行的试验观测值之和为 。
将 分别除以 ,就得j r W j r j j W W W ,,,21 r n 到与这一列中各水平对应的均值 。
实验设计方法-正交法实验法
实验设计方法—正交实验法概述正交实验法就是利用排列整齐的表-正交表来对试验进行整体设计、综合比拟、统计分析,实现通过少数的实验次数找到较好的生产条件,以到达最高生产工艺效果。
正交表能够在因素变化范围内均衡抽样,使每次试验都具有较强的代表性,由于正交表具备均衡分散的特点,保证了全面实验的某些要求,这些试验往往能够较好或更好的到达实验的目的。
正交实验设计包括两局部内容:第一,是怎样安排实验;第二,是怎样分析实验结果。
正交试验设计法的根本思想正交试验设计法,就是使用已经造好了的表格--正交表--来安排试验并进行数据分析的一种方法。
它简单易行,计算表格化,使用者能够迅速掌握。
下边通过一个例子来说明正交试验设计法的根本想法。
[例1]为提高某化工产品的转化率,选择了三个有关因素进行条件试验,反响温度(A),反响时间(B),用碱量(C),并确定了它们的试验范围:A:80-90℃B:90-150分钟C:5-7%试验目的是搞清楚因子A、B、C对转化率有什么影响,哪些是主要的,哪些是次要的,从而确定最适生产条件,即温度、时间及用碱量各为多少才能使转化率高。
试制定试验方案。
这里,对因子A,在试验范围内选了三个水平;因子B和C也都取三个水平:A:Al=80℃,A2=85℃,A3=90℃B:Bl=90分,B2=120分,B3=150分C:Cl=5%,C2=6%,C3=7%当然,在正交试验设计中,因子可以是定量的,也可以是定性的。
而定量因子各水平间的距离可以相等,也可以不相等。
这个三因子三水平的条件试验,通常有两种试验进行方法:(Ⅰ)取三因子所有水平之间的组合,即AlBlC1,A1BlC2,A1B2C1,……,A3B3C3,共有33=27次试验。
用图表示就是图1 立方体的27个节点。
这种试验法叫做全面试验法。
全面试验对各因子与指标间的关系剖析得比拟清楚。
但试验次数太多。
特别是当因子数目多,每个因子的水平数目也多时。
试验量大得惊人。
16正交试验设计
10
质量管理学
正交试验设计
表 16-2 因子水平表 - 因子 水平 一 A: 充磁量(10-4T) 900 B: 定位角度 (π/180) ( rad) 10 C: 定子线圈匝数(匝) 70
二 1100 11 80
三 1300 12 90
4.选用合适的正交表,进行表头设计,列出试验计划。首 先根据在试验中所考察的因子水平数选择具有该水平数的 一类正交表,再根据因子的个数具体选定一张表。在本例 中所考察的因子是三水平的,因此选用三水平正交表,又 由于现在只考察三个因子,所以选用L9(34)即可。
8
质量管理学
正交试验设计Βιβλιοθήκη 16.2 无交互作用的正交设计
例 16.1 磁鼓电机是彩色录像机磁鼓组建的部件 之一,按质量要求其输出力矩应大于210g.cm。 某生产厂过去这项指标的合格率较低,从而希望 通过试验找出好的条件,以提高磁鼓电机的输出 力矩。
9
质量管理学
正交试验设计
试验的步骤
1.确试验目的。试验前,首先要明确试验目的,即通过试验想解决什 么问题。是为了改进质量,还是为了提高产量,或是为了保护环境, 等等。在本例中试验的目的是提高磁鼓的电机的输出力矩。 2.明确试验指标。试验指标用来判断试验条件的好坏,在本例中直接 用输出力矩作为考察指标,该指标越大表明试验条件越好。 3.确定因子与水平。 在试验前首先要分析影响指标的因子是什么, 每个因子在试验中取那些水平。在本例中,经分析影响输出力矩的 可能因子有三个,它们是 A: 充磁量 B: 定位角度 C: 锭子线圈匝数
4
质量管理学
正交试验设计
试验指标:衡量试验条件好坏的特性(可以是质 量特性也可以是产量特性或其它)称为指标,它 是一个随机变量。为了方便起见,常用x表示。 正交试验设计:来选择最佳的或满意的试验条件, 即通过安排若干个条件进行试验,并利用正交表 的特点进行数据分析的一种常用的试验设计的方 法。
正交试验设计及结果分析
上一张 下一张 主 页
四因素、三水平的试验因素水平表
水平
1 2 3
试 A B
验
因 C
素 D
3
(3) 选择合适的正交表 正交表的选择是正交试验设计的首要问题。确定了因
素及其水平后,根据因素、水平及需要考察的交互作用的
多少来选择合适的正交表。正交表的选择原则是在能够安 排下试验因素和交互作用的前提下,尽可能选用较小的正
3
上一张 下一张 主 页
退 出
以上选择 ,保证了A因素的每个水平与B因素、C因素的
各个水平在试验中各搭配一次 。对于A、B、C 3个因素来 说, 是在27个全面试验点中选择9个试验点 ,仅是全面试 验的三分之一。 从上图中可以看到,9个试验点在选优区中分布是均衡
的,在立方体的每个平面上,都恰是3个试验点;在立方体
1.2 正交试验设计的基本原理
3
上一张 下一张
主 页
正交设计就是从选优区全面试验点(水平组合)中 挑选出有代表性的部分试验点(水平组合)来进行试验。 上图中标有试验号的九个 “(·)”,就是利用正交表
L9(34)从27个试验点中挑选出来的9个试验点。即:
(1)A1B1C1 (4)A1B2C2 (7)A1B3C3 (2)A2B1C2 (5)A2B2C3 (8)A2B3C1 (3)A3B1C3 (6)A3B2C1 (9)A3B3C2
3
上一张 下一张 主 页
退 出
1.3.2.2 代表性
一方面: (1)任一列的各水平都出现,使得部分试验中包括了所 有因素的所有水平; (2)任两列的所有水平组合都出现,使任意两因素间的 试验组合为全面试验。 另一方面: 由于正交表的正交性,正交试验的试验点必然均衡地 分布在全面试验点中,具有很强的代表性。因此,部分试 验寻找的最优条件与全面试验所找的最优条件,应有一致 的趋势。
第四章 正交试验设计
6.60 7.67 8.00
1.40
因素主次
优化方案 Ⅰ Ⅱ Ⅲ 葛根总 黄酮含 量/% k1 k2 k3 R 因素主次 优化方案 18.6 20.2 22.3 6.20 6.73 7.43 1.23 19.3 20.7 21.1 6.43 6.90 7.03 0.6
CAB
C3A2B2或C3A2B3 20.0 20.2 20.9 6.67 6.73 6.97 0.3 ACB A 3 C 3B 3 18.5 20.5 22.1 6.17 6.83 7.37 1.20
17
4.2 无交互作用的正交试验设计
Ⅰ Ⅱ Ⅲ 葛根素 含量/% k1 k2 k3 R 7.2 7.4 8.1 2.40 2.47 2.70 0.3 7.3 8.1 7.3 2.43 2.70 2.43 0.27 7.7 7.1 7.9 2.57 2.37 2.63 0.26 6.8 7.8 8.1 2.27 2.60 2.70 0.43
0.78 0.76 0.74 0.72 0.7 0.68 0.66 0.64 0.62 1 2 3 4 5
0.8
0.75
0.7
0.65
0.6
甲
乙
丙
13
4.2 无交互作用的正交试验设计
4.2.3多指标正交试验的结果分析
在多指标试验中,不同指标的重要程度往往不一样, 各因素对不同指标的影响程度也不完全相同,为了兼顾各 个指标、因素的取优,可以采用以下两种结果分析方法: 综合平衡法和综合评分法。
L16 (4 4 23 )
L16 (8 28 )
L18 (2 37 )
4.1.3 正交表的性质——正交性原理
1、分配均匀性:正交表中因素的每一水平出现的次数 相同。 2、搭配均匀性:正交表中两列有序数对(水平搭配) 出现的次数相同。 正是基于正交性原理,正交试验设计可以用少数次 分布十分均匀的试验来很好地代表全部试验。
微机试验题目-201212
试验一、正交试验设计及结果分析一、正交试验设计及极差分析1、在猕猴桃果茶的研制中,为寻找最好的配方,进行了调配试验,试验考察的指标由感官标准评定,其因素水平见下表。
试设计一个正交试验求出最好的配方组合。
(无交互作用的)三因素三水平试验。
2、在油炸方便面生产中,主要原料质量和主要工艺参数对产品的质量有影响,今欲通过试验确定最佳生产条件。
根据专业知识和实践经验,确定试验因素和水平见表。
选择L9(3)安排试验方案。
二、方差分析1、不考虑交互作用的正交试验方差分析:在双歧杆菌奶研制中,为选择最佳发酵条件,用L8(27)正交表安排了正交试验,试验因素与水平见表1,其方案和结果见表2,试对试验结果进行方差分析。
(两水平)表2 正交试验设计及结果2、考虑交互作用的正交试验设计运动发酵单细胞菌是一种酒精生产菌,为了确定其发酵培养基的最佳配方,进行了四因素三水平正交试验,试验指标为酒精浓度(g/ml ),试验因素水平见下表。
试对试验结果进行方差分析。
1 2 35 15 250 0.5 1.0025 30 355.06.07.0试验二、部分因子试验设计及结果分析1、在探讨竹叶中黄酮化合物的超声波辅助提取条件过程中,影响提取得率的因素很多,其中溶剂浓度、料液比、提取温度、超声波功率和提取时间等是比较常见需要考虑的因子,因子水平表见表1,请采用Design Expert软件设计两水平5因子的部分试验设计,根据试验要求,请选择16组试验作为设计表。
表1 部分因子水平表2、若所选择的16组试验部分因子设计结果见表2,请根据结果找出影响提取得率的较显著因子,请写出方差分析结果和回归模型方程。
表2 部分因子试验设计及其结果试验三、响应面试验设计及结果优化分析采用两因素五水平的中心组合设计研究胰蛋白酶水解螺旋藻蛋白质时,pH值、酶-底物浓度比(E/S)两因素及二因素间交互作用对水解蛋白收率的影响。
因素编码及试验方案和结果见表1和表2。
无交互作用的正交设计实验设计报告
实验报告
(与程序设计有关的课程用)
课程名称:实验设计及优化B
实验题目:无交互作用的正交设计
班级学号:
姓名:
成绩:
沈阳理工大学
2012 年11 月26 日
T3
10.
10. 12.
10. 12.
R
由直观看的好条件是3号实验,最有水平组合为A1B3C3,由极差R的大小得到各因素的主次顺序为A>B>C。
2)直观分析得趋势图如下
1
T
2
T
3
T
二.方差分析结果分析
用方差分析得表如下:
实验号 A B C
转化率y%
1 1 1 1
2 1 2 2
3 1 3 3
4 2 1 2
5 2 2 3
6 2 3 1
7 3 1 3
8 3 2 1
9 3 3 2 水平1 10. 水平2 10. 12. 水平3 10. 12. 误差 SS 91. 14. 11. df 2 2 2 8 2 MS 45.
F 68. 11. P 值
(1)计算离差平方和
①总离差平方和 (1)
其中n 是正交表的行数,正交表的每行确定一个实验处理,每个处理得到一个实验数据,共有n 个实验数据,记为了y 1,y 2,…y n 。
,本例n=9。
y 是n 个实验数据的平均值。
②因素的离差平方和。
因素A 的离差平方和为
(2)。
数理统计课件4-3
较好的生产条件
A3 B2C2
温度900C 加碱量48kg 乙种催化剂
9
是否每个因子都要取使平均收率最高 的水平呢?
通过分析极差寻找主要因子,极差越大 说明该因子对指标影响越大
温度80-900C 收率82-60 极差22 碱量35-55kg 收率75-67 极差 8 催化剂甲乙丙 收率79-65 极差14
1
80 35 甲
51
2
80 48 乙
71
3
80 55 丙
58
4
85 35 乙
82
5
85 48 丙
69
6
85 55 甲
59
7
90 35 丙
77
8
90 48 甲
85 ★
9
90 55 乙
84
收率试验结果计算分析表
L9 (34 )
6
试验号 因素 A 温 B 碱
1
80 35
2
80 48
3
80 55
4
85 35
14
因子 A, B,C在水平 Ai , Bi ,Ci 上的效应分 别为 a i , bi , ci ,试验值用 Y j 表示, 于是有
Y1 a1 b1 c1 1
Y2 a1 b2 c2 2
3
其中 ai 0
i 1
Y3 a1 b3 c3 3
Y4 a2 b1 c2 4 Y5 a2 b2 c3 5
1
一. 无交互作用的正交试验
例1 某化工厂生产的一种产品,收率较低 且不稳定,一般在60%-80% 之间波动.现 要求通过试验,摸清规律,找出好的生产条件, 实现稳产高产(收率=实际产量/理论最大产量)
根据经验认为反应温度的高低, 加碱 量的多少和催化剂的不同是造成收率不稳 定的三个主要原因.
第八章正交试验设计
均值
y1
y2
…
yr 1 , yr 2 ,, yrm
…
yr
记第i 水平下的数据和为Ti,Ti yij ;
j 1
m
记第i水平下的数据均值为 y i ,总均值为 y 。此 时共有n=rm个数据,这n个数据不全相同,它们的 波动(差异)可以用总离差平方和ST去表示
ST ( yij y )2
ST=121492-12002/12=1492, SA=485216/4-12002/12=1304, Se= 1492-1304=188, fT=3×4-1=11 fA=3-1=2 fe=11-2=9
(4)列方差分析表: 本例的方差分析表
来源 因子A 误差e 总计T 偏差平方和 自由度 均方和 F比
~ N (0,1) ,
m
2
SE
~ 2 (n r )
由正态分布性质, yi 与 S E
i 1
r
2 ( y y ) ij i 独立 j 1
1
m ( yi i ) / m ( yi i ) ~ t (n r ). SE MSE /( n r ) 2
A1
A2 A3 A4
利用上表中的数据可算得各平方和:
ST
T2 z 160.7895, 20 i 1 jSA 63.2855, 20 i 1 5 S e ST S A 97.5040.
来源 因子A 平方和 63.258 自由度 3 均方和 21.0952 F比 3.46
的估计是 20.9 4.57
i
各水平均值 i 的区间估计
i 置信度 1 置信区间
yi t (n r ) MSE / m , 1 2 yi t (n r ) MSE / m 1 2
9正交试验设计
Y3
Y4
4
5 6 7 8 9
Y5
Y6 Y7
A3 B1 C 3 A3 B 2 C 1 A3 B 3 C 2
Y8 Y9
假定因素A,B,C没有交互作用。 设因素A在水平 A1 , A2 , A3 上的效应分别为 a1 , a 2 , a 3 因素B在水平 B1 , B 2 , B 3 上的效应分别为 b1 , b 2 , b3 因素C在水平 C 1 , C 2 , C 3 上的效应分别为 c1 , c 2 , c 3
A
K 1 141,
B
K 2 165,
B
K 3 144
B
K 1 1 3 5,
C
k 2 1 7 1,
C
k3 144
C
=>试验值的平均值
k 1 4 1,
A
k 2 4 8,
A
k3 61
A
k 1 47,
B
k 2 55,
B
k 3 48
B
k 1 4 5,
C
k2 57,
5
6 7 8
2
2 2 2
1
1 2 2
2
2 1 1
L ——正交表 2——水平数
8——试验次数 7——表中最多可安排的因素数
L9 3
4
列号 1 1 2 3 4 5 6 7 8 9 1 1 1 2 2 2 3 3 3
正交表
试验号
2
1 2 3 1 2 3 1 2 3
3
1 2 3 2 3 1 3 1 2
B
1 3
K2 ,
C
B
k3
B
1 3
16正交实验设计
B2 (91+94)/2=92.5 (83+88)/2=85.5
显然:A2B1最优,C2最优,因子D可以
任选
1 - 40
质量管理
学实验 MINITAB:工作表
1 - 41
MINITAB:方差分析
命令:统计——方差分析——一般线性模型
质学量实管验理MINITAB:输出结果
1 - 43
MINITAB:方差分析
上面的例子中对因子A与B应该选择A2B2,因 子C可以任选,譬如为节约材料可选择C1
1 - 25
MINITAB:方差分析
命令:统计——方差——一般线性模型——结果
质学量实管验理MINITAB:输出结果
显然:A2B2最优,因子C可以任选
1 - 27
质量管理
学实验 (四)验证试验
对A2B2C1进行三次试验,结 果为:234,240,220,平均值 为231.3此结果是满意的。
可以认为水平的差异确实是有影响的呢?
2. 数据的方差分析 要把引起数据波动的原因进行分解,数据的
波动可以用偏差平方和来表示。
正交表中第j列的偏差平方和的计算公式:
Sj
i
Ti2j T2 n/q n
其中Tij为第j列第i水平的数据和,T为数 据总和,n为正交表的行数,q为该列的水平
数。该列表头是哪个因子,则该Sj即为该因子 的偏差平方和,譬如SA=S1
1 - 21
质量管理
学实验 MINITAB:工作表
1 - 22
MINITAB:方差分析
命令:统计——方差分析——一般线性模型
质学量实管验理MINITAB:输出结果
1 - 24
质量管理 学实验
3. 最佳条件的选择
不考虑交互作用的正交试验设计041019200448
不考虑交互作用的正交试验设计************§6.5 不考虑交互作用的正交试验设计下面我们先来看一种比较简单的情形,即不考虑因子之间交互作用的情形。
不考虑交互作用的正交试验设计和数据处理,可按下列步骤进行:(1)选正交表。
)(mn r L 选择的原则是:要等于因子的水平数;要大于或等于因子的个数;是试验r m n 次数,要尽可能小。
例如,问题中有4个因子,每个因子都是2个水平。
选正交表时,首先要选=2 的r 表。
这样的正交表有。
在中,=3 ,小于因子的 ),2(),2(),2(15167834L L L )2(34L m 个数4,所以不符合要求。
在中,,大于因子的个数 ),2(),2(151678L L ,15,7=m 4,符合要求。
在符合要求的正交表中,还要选试验次数尽可能小的那一个表。
显然,n 在中,=8 为最小,所以,我们最后选定正交表。
)2(78L n )2(78L (2)设计表头。
将各因子安排在正交表的各列上方,每个因子占1列,这称为表头。
在不考虑交互作用的正交试验设计中,表头上的因子可以任意安放。
表头上不放因子的列,称为空白列。
(3)按照设计做试验,取得试验观测值。
正交表的每一行代表一种水平组合,对每一种水平组合做一次试验。
按照第行的k 水平组合所做的第次试验,所得到的观测值记为。
正交表有行,所以,一共k k X n 要做次试验,共得到个试验观测值:。
n n n X X X ,,,21 (4)在正交表的每一列中,求出与各水平对应的均值,以及这一列的平方和。
设我们考虑的是第列。
在这一列中,表示水平的数字,每一个都重j r ,,2,1 复出现次。
设与这一列中的数字对应的那几行的试验观测值之和为,与这r n 1j W 1一列中的数字对应的那几行的试验观测值之和为,,与这一列中的数字 2j W 2 r 对应的那几行的试验观测值之和为。
将分别除以,就得j r W j r j j W W W ,,,21 r n 到与这一列中各水平对应的均值。
如何正确处理无法考察交互作用的三因素或四因素设计定量资料
如何正确处理无法考察交互作用的三因素或四因素设计定量资料[ 09-08-13 17:03:00 ] 作者:高辉, 胡良平,郭编辑:studa20【关键词】统计学; 医学; 统计分析; 交互作用Keywords: statistics; medicine; statistical analysis; interaction1 引言在实际科研工作中,研究者有专业知识为依据可认为多个因素间相互独立(不存在交互作用)或根据研究目的不需要考察多个因素之间的交互作用或受试验条件及样本量等的限制无法考察多因素间的交互作用时,往往需要采用一些比较特殊的试验设计类型。
当试验涉及到3个或4个因素,而由于前述各项原因无法或不必要考察这些因素之间的交互作用时,可根据试验目的和试验因素等选用拉丁方设计、交叉设计、三因素或四因素嵌套设计。
本文着重论述这几种试验设计的特点、区别、适用场合及如何根据资料的类型选用合适的统计分析方法,以期能帮助广大科研工作者提高处理无法考察交互作用的三因素或四因素设计定量资料的统计学水平。
2 3种试验设计的特点及适用场合2.1 拉丁方设计的特点及适用场合拉丁方设计可用来考察一个试验因素和多个区组因素(重要非试验因素)对观察结果的影响,要求这些因素水平数相同且相互之间不存在交互作用或交互作用可以忽略不计。
根据试验涉及的区组因素的个数,可将拉丁方设计分为三因素拉丁方设计或四因素拉丁方设计。
前者即通常所说的拉丁方设计(如无特殊说明,本文后面所讨论的拉丁方设计指三因素拉丁方设计),后者称希腊拉丁方设计或正交拉丁方设计。
2.1.1 拉丁方设计拉丁方设计可安排一个试验因素和2个区组因素,视为随机区组设计的扩展,因为它比随机区组设计多考察一个区组因素对观测结果的影响。
拉丁方设计要求3个因素之间不存在交互作用或交互作用可以忽略不计,要求3个因素水平数相同(以试验因素水平数为基准)。
2.1.2 希腊拉丁方设计希腊拉丁方设计可安排一个试验因素和3个区组因素,要求4个因素之间不存在交互作用或交互作用可以忽略不计,要求4个因素水平数相同(以试验因素水平数为基准)。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
的 各 列 的 平 方 和 。 SS eSS STS SSA 是B 误 差 平 方 和 。 可 以 证 明 ,
分别是表 头为
A, B,
m
SST SS ,即总j 平方和等于各列的平方和之和,所以, SSe 也就是空白列的平方和 j1
之和。
fT n 1 是 总 自 由 度 , f f A r B 1 是 各 因 子 的 自 由 度 ,
;若因子 的
作
A
位置会有一个向右的偏移。
若因子 B 的作用不显著,则
F
B
MS B MSe
~ F( f B, f ) e
B
;若因子 的作用显著,
则
FB 的值会偏大,统计量 FB 的分布,相对于 F( f , f ) 分布B 来e说,峰值的位置会有一
个向右的偏移。
所以,像在方差分析中一样,只要给定显著水平 ,就可以用 F 分布检验因子 A, B, 的作用是否显著。
W
W
到与这一列中各水平对应的均值 X
1j
1j
nr
,
X2j
2j
nr ,
,X
rj
rj n r
。
然后,从
X 1,jX , 2 j , X r j 出 发 , 求 出 这 一 列 的 平 方 和
SS j
n r
r i1
(X
ij
X
)2 ,
其中,X
1 n X 是 总 均 值 。 具 体 算 法 是 , 把
n k1 k
X 1,jX , 2 j , X r j 看作样本, X 就是样本均值, SS 就是j 样本方差再乘以 n (或修
161
正样本方差再乘以 n(r 1) r ) 。
(5)列方差分析表,作显著性检验。
来源 平方和
自由度
均方
A
SS A
f A r 1 MS A SS A f A
B
SS B
19.0
F0.95 (2, 2) 19.0
因为 F B 6.33 19.0 F ( 1f, f ) B,所e以因子 B 作用不显著。
因为 F C 13.00 19.0 F ( 1f, f )C,所e以因子 C 作用也不显著。
(6) 寻找最优水平组合。
对于因子 A ,因为 A 的均1值 X 41 ,1A1 的均值 X 2 48 , A3 的均21
§6.5 不考虑交互作用的正交试验设计
下面我们先来看一种比较简单的情形,即不考虑因子之间交互作用的情形。不考虑交 互作用的正交试验设计和数据处理,可按下列步骤进行:
(1)选正交表 Lnr( m ) 。
选择的原则是: r 要等于因子的水平数; m 要大于或等于因子的个数; n 是试验
次数,要尽可能小。
表头
ABC
列号 试验号
1
2
34
1
1
1
11
观测值(收得率) X k
31
2
1
2
22
54
3
1
3
33
38
4
2
1
23
53
5
2
2
31
49
6
2
3
12
42
7
3
1
32
57
8
3
2
13
62
9
3
3
21
64
X1j
41 47 45 48
X2j
48 55 57 51
X3j
61 48 48 51
X 1 n X 50 n k1 k
试验观测
复出现 n r 次。设与这一列中的数字 1 对应的那几行的试验观测值值之:和X为, W1 j ,与这
一列中的数字 2 对应的那几行的试验观测值之和为 W2 j , ,与X这, 一,列X中的数字 r
n
对应的那几行的试验观测值之和为
Wr
j
。将
W
,
W
,
1
j
, W 分别除以 n r ,就得
2j
rj
W
f e f Tf f A 是B 误差自由度。
MS A SS f A, MAS SS fB , B是各B因子的均方, MS SS fe 是误 e
e
差均方。
可以证明,若因子 A 的作用不显著,则
F
A
MS A MSe
~
F(
f A, f
)e
用显著,则 FA 的值会偏大,统计量 FA 的分布,相对于 F( f , f ) 分布A 来说e ,峰值的
162
(6)寻找最优水平组合。 对每一个因子,在以它为表头的那一列中,比较各水平的均值的大小,可以确定哪一 个水平最优。 由于不考虑交互作用,所以,只要将各因子的最优水平组合起来,就是最优水平组合。 下面看一个实际例子。
例 1 某化工厂为提高产品的收得率(单位:%) ,进行 3 因子 3 水平正交试验。
所取的因子和水平分别为:
因子 A 是反应温度, A 是180C, A 是 852C, A 是 90C3 ;
Hale Waihona Puke 因子 B 是反应时间, B 是190 分钟, B 是 1220 分钟, B3 是 150 分钟 ;
因子 C 是用碱量,C 是15% ,C 是 6%2,C3 是 7% 。
要求进行不考虑交互作用的正交试验设计,检验因子 A, B, C 的作用是否显著(显
值 X 3348 ,其中 X 57 最23 大,所以 C2 是最优水平。
把 3 个因子的最优水平组合起来,就得到最优水平组合 (A , B3, C 2) ,即2 反应温度
为 90C,反应时间为 120 分钟,用碱量为 6% 。
由于因子 B 很不显著,即反应时间的不同对于收得率没有显著的影响,为了节省反 应时间,也可以考虑把反应时间改为 90 分钟,选用水平组合 (A , B 3, C 1) 。 2
值 X 3161 ,其中 X 61 3最1 大,所以 A 是最优水平3 。
对于因子 B ,因为 B 的均1值 X 47 ,1B2 2 的均值 X 55 , B3 的均22
值 X 3248 ,其中 X 55 最22 大,所以 B2 是最优水平。
对于因子 C ,因为 C 的均1值 X 45 ,13C 的均值 X 2 57 , C3 的均23
因子 C 也不十分显著,即用碱量的不同对于收得率也没有太大的影响,如果希望节 省用碱量,还可以考虑把用碱量改为 5%,选用水平组合 (3A , 1B , C1 ) 。
得到最优水平组合后,还可以对它以及在它的附近再做几次试验,看看它是否确实最 优,是否还可以作改进,进一步得到更好的结果。
164
57
F 6.33
117
B
FC 13.00
误差 SSe 18
8222 2 9
总和 SST 984
n 1 8
因为 F A 34.33 19.0 F ( 1f, f ) A,所e以因子 A 作用显著。
F 0.95
分位数
(2, 2)
19.0
F 0.95
(2,
2)
(3)按照设计做试验,取得试验观测值。
正交表的每一行代表一种水平组合,对每一种水平组合做一次试验。按照第 k 行的
水平组合所做的第 k 次试验,所得到的观测值记为 X 。正交k 表有 n 行,所以,一共
12
要做。 n 次
(4)在正交表的每一列中,求出与各水平对应的均值,以及试这验一,列共的平方和。
设我们考虑的是第 j 列。在这一列中,表示水平的数字 1, 2得, 到 ,nr个,每一个都重
f r 1 MS B SS B f B
B
F值
F A
MS
A
MSe
F B
MS B MSe
分位数
F1(f , fA) e
F(f,f)
1 B e
误差 总和
SSe SST
其中, SS T
fe f n 1
T
n
( X X)2
k1
k
MSe SS fee 是SS总A, 平SS方,和B。
SS j
618 114 234 18
n
SST (X X k) k 1
m
2 j1
SS 984
j
163
(5)列方差分析表,作显著性检验。
来源 平方和
自由度
均方
F值
A SS A 618
r 1 2
309 FA 34.33
B SSB 114 C SSC 234
r 1 2 r 1 2
16
,大于因子的个数
4,符合要求。在符合要求的正交表中,还要选试验次数 n 尽可能小的那一个表。显然,
在 L(2 7) 中, =8n为最小,所以,我们最后选定正交表
8
(2)设计表头。
8 7 L (2 )
。
将各因子安排在正交表的各列上方,每个因子占 1 列,这称为表头。在不考虑交互作
用的正交试验设计中,表头上的因子可以任意安放。表头上不放因子的列,称为空白列。
著水平 0.05 ) ,并且找出最优水平组合。
解
(1)选正交表。按照 r 3, m 3, n 尽可能小的原则,选用
(2)设计表头。将因子 A, B, C 依次安排在第 1, 2, 3 列。
L(34) 。
9
(3)按照设计做试验,取得试验观测值。试验得到的观测值见下表。 (4)求各列与各水平对应的均值和各列的平方和。计算结果见下表。
例如,问题中有 4 个因子,每个因子都是 2 个水平。选正交表时,首先要选 r =2 的