一次回归正交设计

《一次回归正交设计》由会员分享，可在线阅读，更多相关《一次回归正交设计（32页珍藏版）》请在装配图网上搜索。

1、第五讲回归设计及统计分析设目标性状y与z1、z2zm等因素有关，我们可以应用回归分析的方法建立y与诸因素的回归方程，以此对y进行预测和控制，或筛选y的最优指标。z1、z2zm构成一个因子空间，每一组z1、z2zm值对应一个y值。如何在因子空间中选择最适当的试验点，以最少的试验点寻求y的最优区域，这就要将回归分析与正交设计结合起来应用，称为回归正交设计。按回归模型的次数，回归正交设计又分为一次回归正交设计和二次回归正交设计。一、一次回归正交设计一次回归正交设计主要是应用2水平正交表进行设计，其设计和分析步骤如下。1确定试验因素的变化范围例如研究m个栽培因素z1、z2zm与作物产量y的数量关系，首

2、先需确定各个栽培因素的变化范围。设因素zj的变化区间为(z1j，z2j)，则z1j和z2j分别为因素zj的下水平和上水平。那么 为因素zj的零水平。 为因素zj的变化区间。2.对各因素的水平编码编码就是对各个因素的取值作如下线性变换： 式中xj为编码值。如：这样就建立了zj与xj的一一对应关系：下水平 z1j x1j (-1)零水平 z0j x0j (0 )上水平 z0j x0j (+1)通过上面的编码可知，当zj在区间(z1j，z2j)变化时，它的编码值xj就在区间(-1，+1)内变化。多个因素的编码工作可在因素水平编码表(表1)上进行。表1 因素水平编码表zj因 素Z1Z2Zm下水平Z11

3、Z12Z1m零水平Z01Z02Z0m上水平Z21Z22Z2m变化间距j12m对因素的水平进行编码后，y对z1、z2zm的回归问题就转化为对x1、x2xm的回归问题。在z1、z2zm因子空间选择试验点的问题就转化为x1、x2xm为坐标轴的编码空间选择试验点。在二次回归设计中也要进行因素的编码工作。3.选择合适的二水平正交表常用的二水平正交表有L4(23)、L8(27)、L12(211)、L16(215)等。选用哪一种二水平正交表要依据因素个数及需要研究的交互作用而定。正交表确定以后，把表中的“2”改为“-1”。这样正交表中的“+1”“-1”既表示因素的不同水平，也表示xj的取值。表2列举了经代换

4、后的几张常用二水平正交表。表2 常用二水平正交表L4(23)试验号x1x 3x3111121-1-13-11-14-1-11L8(27)试验号x1x2x3x1 x2x1 x3x2 x3x1 x2 x311111111211-11-1-1-131-11-11-1-141-1-1-1-1115-111-1-11-16-11-1-11-117-1-111-1-118-1-1-1111-1L12(211)试验号x1x2x3x4x5x6x7x8x9x10x11111111111111211111-1-1-1-1-1-1311-1-1-1111-1-1-141-11-1-11-1-111-151-1-11

5、-1-11-11-1161-1-1-11-1-11-1117-11-1-111-1-11-118-11-11-1-1-1111-19-111-1-1-11-1-11110-1-1-11111-1-11-111-1-11-11-1111-1-112-1-111-11-11-1-11如设计一个3因素试验，可选用L8(27)正交表，表中x1、x2、x3分别代表z1、z2、z3的编码值。若因素间有互作存在，在回归中可用非线性项x1 x2、x1 x3、x2 x3等表示。每种交互作用占改造后二水平正交表的1列，该列的取值可由某两列上元素对应相乘得到。如表2 L8(27)中x1 x2列的元素是由x1与 x2

6、列上的对应元素相乘而得。用二水平正交表设计的这种试验具有正交性。若以表示在第试验中第j个变量的编码值，于是在试验计划中有任一列的和任两列的内积具有以上两个性质的设计为正交设计。4.建立回归方程对于3因素试验，若考虑因素间的交互作用，则回归方程为例如用L8(27)正交表设计该试验，那么它的结构矩阵为信息矩阵(系数矩阵)为相关矩阵为常数项矩阵为为试验结果，于是可算出回归系数矩阵那么各类回归系数即由下式算出 回归系数的具体计算可在正交表上进行(表3).表中bj为各回归系数，Qj为偏回归平方和。从而建立回归方程。表3 3因素一次回归正交设计计算表试验号x0x1x2x3x1 x2x1 x3x2 x3试验

7、结果11111111y12111-11-1-1y2311-11-11-1y3411-1-1-1-11y451-111-1-11y561-11-1-11-1y671-1-111-1-1y781-1-1-1111y8Bj8888888B0/8B1/8B2/8B3/8B12/8B13/8B23/8-从以上计算可看出，各变量的偏回归平方和，与偏回归系数bj的平方成正比。bj的绝对值越大，Qj也越大。这就意味着，在利用正交表所得到的回归方程中，每一个回归系数bj的绝对值大小，反映了对应变量xj对y作用的大小。这是因为经过无量纲编码后，所以变量的取值都是1和-1，它们在所研究的区域内取值是平等的，且不受单

8、位的影响，因此所求回归系数bj直接反映了因素zj作用的大小，回归系数的符号反映因素作用的性质。在要求不太高的情况下，一次回归正交设计可省略方差分析，直接把回归系数与零相差不大的因素从回归方程中剔除，不需重新计算其它回归系数，剔除因素对结果的影响可并入试验误差。但对精度要求较高的试验，应继续进行回归关系的显著性测验。5.回归方程及回归系数的显著性测验一次回归正交设计的方差分析如表4。表4 一次回归正交设计的方差分析表变异来源自由度平方和均方F值回 归离回归总x11xm1x1 x21xm-1 xm1对回归方程的显著性假设测验可通过表4中的F测验进行。但这种测验只是说明m个变量对试验结果的影响是显著

9、的，而在研究区域内回归方程与实测值的拟合情况，即采用一次回归模型是不是最合适，从以上测验中没有得到这方面的信息。为了了解回归方程的拟合情况，需在零水平()安排一些重复试验，如在安排p次重复试验所得试验结果为，其平均数为，则 由此可用t测验对回归方程中的拟合情况进行测验： 若求得，则认为b0与无显著差异，一次回归方程与实测值拟合得较好。反之，就认为用一次回归来描述问题不够确切，可考虑建立高次方程。拟合度测验也可通过F测验进行： 若F1F0.05，则说明回归方程的拟合度较好。一般先进行拟合测验，在方程拟合情况较好的情况下，再进行回归方程的显著性测验。故将拟合度测验的F值记作F1，而回归显著性测验的

10、F值记作F2。回归系数的显著性测验可由表4的F测验进行，也可通过t检验给出。 服从的t分布。研究氮、磷、钾施用量对大豆籽粒产量的影响，试作一次回归正交设计，并对试验结果进行分析。1.确定各试验因素水平并进行编码首先各因素的上、下水平，将3个因素的变化范围分别定在(2,6)、(3,9)、(3,9)，单位为斤/亩。然后计算各因素的零水平和变化间距，列出因素水平编码表(表5)。表5因素水平编码表因素变化间距水平编码-10+1z1氮(尿素)2斤/亩246z2磷(三料磷)3斤/亩369z3钾(硫酸钾)3斤/亩3692.制定实施方案选择L8(27)正交表，实施方案如表6。表6 氮、磷、钾肥料三因素试验实施

11、方案试验号试验设计矩阵 实施方案x1x2x3尿素(斤/亩)三料磷(斤/亩)硫酸钾(斤/亩)1111699211-169331-1163941-1-16335-1112996-11-12937-1-112398-1-1-12339000466100004661100046612000466为对回归方程进行拟合度检验，增设4个零水平试验。经实施后，试验结果及结构矩阵列于表7。为计算常数项b0，在试验结构矩阵中添上x0列，取值皆为1。表7 试验结构及数据分析试验号x0x1x2x3x1 x2x1 x3x2 x3y(公斤/亩)11111111155.32111-11-1-1180.7311-11-11-

12、1180.7411-1-1-1-11178.351-111-1-11121.761-11-1-11-1153.671-1-111-1-1112.081-1-1-1111116.191000000158.3101000000163.3111000000166.0121000000150.81836.8191.624.2-59.0-70.213-55.612888888153.123.953.025-7.375-8.7751.625-6.95-4588.8273.205435.125616.00521.125386.423.计算回归系数，建立回归方程产量结果列在表7的最后一列。计算可在表7上进行。

13、各项数值的计算过程如下： 由以上计算得如下回归方程：4.回归方程的假设测验先计算各项平方和与自由度 首先对回归方程的拟合度进行测验，可用F测验和t测验两种方法。(1)F测验(2)t测验以上两种测验都说明建立的回归方程与实际情况吻合较好，可以用一元回归描述。进一步测验回归方程的显著性，在方差分析表(表8)中进行，测验表明所得回归方程达0.01极显著水平。表8方差分析变异来源自由度平方和均方F值Fo.o1回归66120.71020.1213.0810.67离回归5389.8977.98总116510.59回归系数的显著性测验(1)t测验计算各回归系数的t值如下:查表得(2) F测验查表得与t测验结果相同。将不显著的变量x2和x1 x3从回归方程中剔除，则回归方程为作业：有一3因素一次回归正交设计试验，其试验结构矩阵及结果列于下表，试建立回归方程。试验号111111112.942111-11-1-13.48311-11-11-13.49411-1-1-1-113.9551-111-1-113.4061-11-1-11-14.0971-1-111-1-13.8181-1-1-11114.79910000004.171010000004.091110000004.38答案：