第九章回归旋转试验设计.doc

1、9回归旋转试验设计本章要点：主要介绍了回归旋转设计的基本原理、实现条件、组合设计的步骤和统计分析方法，并给出二次回归正交旋转试验设计的计算案例。重点：回归正交旋转设计的实现条件、组合设计的方法、方程的建立及显著性检验。难点：回归正交旋转设计正交和旋转的实现条件及其统计分析。9.1回归旋转试验设计的基本原理前面所介绍的“回归正交设计”，具有试验处理数比较少，计算简便、消除回归系数之间的相关性等优点。但它也存在一定的缺点，即二次回归预测值 的方差随试验点在因子空间的位置不同而呈现较大的差异。由于误差的干扰，就不易根据预测值寻找最优区域。为了克服这个缺点，人们通过进一步研究，提出了回归旋转设计（wh

2、irly design）。所谓旋转性是指试验因素空间中与试验中心距离相等的球面上各处理组合的预测值 的方差具有几乎相等的特性，具有这种性质的回归设计称回归旋转设计。这种设计的意义在于可以直接比较各处理组合预测值的好坏，从而找出预测值相对优良的区域。9.1.1回归设计旋转性条件旋转设计包括一次、二次和三次旋转设计，但研究中最常见的设计是二次回归旋转设计。下面以三元二次回归方程来讨论回归正交的旋转性问题。二次正交多项式方程的估计值为： 如果以三因素二次回归正交设计的数学模型为例：因此其信息矩阵为：上述信息矩阵中的各个元素可用一般形式表达为： ，其中x的指数、分别可取0、1、2、3、4等非负整数。但

3、这些指数之和不能超过4。即：例如，当时，就是信息矩阵中的第1行、第1列元素。又如，当时，就是信息矩阵中的第10行、第10列元素 。信息矩阵A中的元素可分为两类：一类是元素的所有指数都是偶数或零，另一类是元素的所有指数中至少有一个是奇数。这一规律在任何因素和次数的回归设计信息矩阵中都可以发现。一般，在个因素次回归方程中，共有 项，对应的信息矩阵是 阶对称方阵，矩阵中的元素一般形式为：其中指数分别可取等非负整数，且满足为了使旋转设计成为可能，元次回归设计必须满足旋转条件和非退化条件，具体要求如下。1) 回归设计旋转条件 在因素次回归旋转设计中，对应的信息矩阵中的元素应有： （9-1）式中为试验次数

4、， 为待定参数，它的下标 ，为偶数，且 。这个条件说明了旋转设计下信息矩阵的具体结构，也是旋转设计的基本要求，故称旋转条件。下面以三因素二次回归旋转设计为例介绍旋转条件的实现。按旋转条件定理，可计算出都为偶数或零时，各元素 值有下列三种： （1） （2）（3） （i，j1，2，3）其他元素皆为0，这时，三因素二次旋转设计的信息矩阵有如下形式：这种旋转设计的信息矩阵虽然具有旋转性，但是否存在其逆矩阵？因为逆矩阵 的存在，回归系数才有唯一解，这是回归系数存在的条件，即非退化条件。2) 回归设计非退化条件我们知道，对信息矩阵取行列式，只要，必然存在。这就是二次旋转设计的非退化条件。如前述的因素二次旋

5、转设计的信息矩阵取行列式，经计算得：由最后结果看出，欲使 ，即保持矩阵是非退化的，必须要有 ，即 （9-2）上式即为二次旋转设计的非退化条件。进行旋转设计时必须使待定参数Q满足此条件。已经证明，只要使个试验点不在同一个球面上，就能满足非退化条件。最简单的情况是把个试验点分布在2个或3个半径不等的球面上。 综上所述，旋转性条件是旋转设计的必要条件，非退化条件是使旋转设计成为可能的充分条件。二者结合起来才能构成一项旋转设计方案。实际操作上主要借助于组合设计来实现。通常的组合设计试验处理由、和三部分组成，它们都有固定的组合搭配，每种组合中三个部分的处理分别分布在三个半径不等的球面上。即个析因点分布在

6、半径为 的球面上；个析因点分布在半径为 的球面上；个析因点分布在半径为 的球面上。因此，采用组合设计选取的试验点，完全能够满足非退化条件（式9-2），即信息矩阵不会退化。此外，采用组合设计，其信息矩阵的元素中而它的偶次方元素均不等于零，完全符合式（91）的要求。 为了获得旋转设计方案，还必须根据旋转条件式（9-1）确定值，事实上只要由 就可以求出值。在组合设计下，当（全实施）时，前式变为 ，解方程得 （9-3）同理，当（1/2实施）时，；当（1/4实施）时，；（1/8实施）时，为了便于设计，现将个因素不同实施情况下的值列于表9-1。表9-1 二次正交旋转组合设计参数表2（全实施）3（全实施）4

7、（全实施）5（全实施）5（1/2实施）6（1/2实施）6（1/4实施）7（1/2实施）7（1/4实施）8（1/2实施）8（1/4实施）8（1/8实施）481632163216643212864324681010121214141616168912171015822133320111623365936593610059177100591.4141.6822.0002.3782.0002.3782.0002.8282.3783.3642.8282.3749.1.2 旋转组合设计正交性的实现 二次旋转组合设计具有同一球面上各试验点的预测值 的方差相等的优点，但回归统计数的计算较繁琐。如果使它获得正交

8、性就能大大简化计算手续。在二次旋转组合设计中，一次项和交互项的回归系数和仍保持正交，但与 之间以及与之间都存在相关，即不具正交性，因此，要使二次回归旋转组合设计获得正交性，就必须消除常数项与平方项之间以及平方项与之间的相关性，方法如下：1）常数项与平方项间相关性的消除：与两者间相关性的消除比较简单，只要对平方项施行中心化变换即可实现，即： （9-4）2）平方项与间相关性的消除： 平方项之间相关性的消除，必须使 或 ，如何才能满足该要求呢？在组合设计中它的值为： （9-5）对于个因素的二次旋转组合设计，式（95）中的、和都是固定的。因此，只有适当地调整才能使 ，而试验处理数中和也是固定的，这样就

9、只能通过调整中心点的试验处理数使 。由此可见，适当地选取，就能使二次旋转组合设计具有一定的正交性。为了方便设计，已将元不同实施的和列入表91中。综上所述，只要对平方项施行中心化变换，并适当调整就能获得二次正交旋转组合设计方案。9.1.3 旋转组合设计的通用性 二次回归旋转组合设计，具有同一球面上各试验点的预测值 的方差相等的优点，但它还存在不同半径球面上各试验点的预测值 的方差不等的缺点。为了解决这一问题，于是提出了旋转设计的通用性问题。所谓“通用性”，就是试验除了仍保持其旋转性外，还具有各试验点与中心的距离在因子空间编码值区间的范围内，其预测值 的方差基本相等的性质，即同时具有旋转性与通用性

10、。这种设计称为通用旋转组合设计。首先来看预测值 的方差，已知在个因素情况下，其预测值 的方差 （9-6）上式是在 的约定下得到的，这种约定并非本质的，只是为了讨论简单起见。由此可知，只有恰当确定4，才能满足通用性的要求。 那么，对4有什么要求呢？总的来说，它必须使式中 在诸 （内插点）处的值与 处的值之差的平方和为最小，即： （9-7）式中：于是，对于不同的，均可计算出满足式（97）的4。当4确定后，由关系式（见95）可以计算出不同的试验处理数。 当计算结果不是整数时，可取与其最靠近的整数。然后再由 计算出不同值的，上述计算结果列于表92。表9-2 二次通用旋转组合设计参数表2345（1/2实

11、施）6（1/2实施）7（1/2实施）8（1/2实施）8（1/4实施）48161632641286446810121416181.4141.6822.0002.0002.3782.8283.3643.8280.810.860.860.890.900.920.930.931320313253921659356769142113从上述讨论结果看出，为了满足通用性要求，主要在于确定出适当的。因此，只要在中心点安排如表92所列的次试验，旋转组合设计便获得通用性。综上可以看出，正交旋转的好处在于正交性，它是通过增加中心点的试验次数换来的，但有时并不合算。在某些实际问题中，反倒不如选用通用旋转设计。因为通用

12、旋转设计，既能在的较实用区域使方差 基本不变，又在一定程度上减少了试验次数。9.2 二次回归正交旋转组合设计及统计分析9.2.1 二次正交旋转组合设计的基本步骤设研究因素为个，以表示。在进行设计时，首先确定每个因素的上、下水平，进而计算零水平，以及变化间距。因素零水平的计算式为 （9-8）变化间距为： （9-9）式中为待定参数，其值可以从表91中查出。对每个因素的各个水平的取值进行线性变换，得到水平的编码为： （9-10）这样，就将有单位的自然变量变成了无单位的规范变量，编码后，试验因素的水平被编为，具体见表9-3表9-3 因素水平的编码表规范变量自然变量上星号臂上水平1零水平0下水平1下星号

13、臂变化间距试验因素经因素水平编码后，以变量表示，选用适当的二水平正交表，即可设计出二次回归正交旋转组合方案。为了方便设计与统计分析，现将常用的二因素和三因素二次正交旋转组合设计的结构矩阵列于表9-4和表9-5。表9-4 二元二次正交旋转组合设计的结构矩阵处理号123411111111111111110.50.50.50.50.50.50.50.5567811111.4141.41400001.4141.41400001.51.50.50.50.50.51.51.5910111213141516111111110000000000000000000000000.50.50.50.50.50.50

14、.50.50.50.50.50.50.50.50.50.51688488表9-5 三元二次正交旋转组合设计结结构矩阵处理号12345678111111111111111111111111111111111111111111111111111111110.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.4060.406910111213141111111.6821.6820000001.6821.6820000001.6821

15、.6820000000000000000002.2342.2340.5940.5940.5940.5940.5940.5942.2342.2340.5940.5940.5940.5940.5940.5942.2342.2341516171819202122231111111110000000000000000000000000000000000000000000000000000000.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.5940.

16、5940.5940.5940.5940.5940.5942313.65813.65813.65888815.88715.88715.887二次回归正交旋转组合设计试验结果的统计分析，与二次回归正交组合设计试验结果的统计分析方法相似，详见第八章相关内容（8.3.4）例91 用木瓜蛋白酶酶解虾蛋白，试应用三元二次回归正交旋转组合设计法研究酶用量、温度、底物浓度三因素对酸溶性肽得率影响方程式，其中酶用量（）的上下水平为6000U/g和3600U/g，温度（）的上下水平为65和55，底物浓度（）的上下水平为5和3。解：1）因子编码查表得1.682， 依据公式，计算得：表9-6 例9-1二次回归旋转组合

17、试验设计因素编码表 酶用量（） 温度（）底物浓度（）600055134800408736007136563605755354.643.430.62）正交组合设计由于因素数3，故可以选择表94进行试验方案设计，并记录试验结果。3）回归方程的建立由组合设计计算表97可得方程的各个回归系数，因此建立回归方程如下：4）回归方程及偏回归系数的显著性检验对该回归方程及其各系数进行显著性检验，具体见方差分析表98。word文档 可自由复制编辑表9-7 例9-1 三因素二次回归正交旋转组合设计计算表试验号12345678910111213141516171819202122231111111111111111

18、11111111111-1-1-1-11.682-1.682000000000000011-1-111-1-1001.682-1.682000000000001-11-11-11-100001.682-1.68200000000011-1-1-1-1110000000000000001-11-1-11-110000000000000001-1-111-1-110000000000000000.4060.4060.4060.4060.4060.4060.4060.4062.2342.234-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.594

19、-0.594-0.594-0.594-0.5940.4060.4060.4060.4060.4060.4060.4060.406-0.594-0.5942.2342.234-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.5940.4060.4060.4060.4060.4060.4060.4060.406-0.594-0.594-0.594-0.5942.2342.234-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.594-0.59429.4330.0132.3831.09

20、30.4529.730.7530.139.4730.932.430.4730.130.436.9531.0934.336.334.634.735.136.433.623750.6932.638713.658316.32471.195219.510913.6583-1.4837-0.10860.161113.65831.60540.11750.18868-3.33-0.41631.73958-0.69-0.08630.05968-1.77-0.22130.391815.8871-2.9935-0.18840.563915.8871-24.2035-1.523536.874815.8871-30.

21、9059-1.945660.1384=24676.9349=175.3925=119.628655.7639方差分析结果见表9-8。表9-8 例9-1方差分析表差异源显著性回归残差19.51090.16110.18861.73950.05960.39180.563936.874860.1384119.628655.763911111111191319.51090.16110.18861.73950.05960.39180.563936.874860.138413.29214.28954.54850.03760.04400.40550.01390.09130.13158.596514.01993.0988*总和175.3925注：方差分析结果表明，所建立的回归方程以及、回归系数达到极显著水