回归的旋转设计PPT课件

.,1,第九章回归的旋转设计,本章内容：1旋转设计的基本原理2二次正交旋转组合设计及其统计分析3通用旋转组合设计及其统计分析,本章学习目的与要求：,1.,2.,3.,.,2,1旋转设计的基本原理,“回归的正交设计”具有试验处理数比较少，计算简便，消除了回归系数之间的相关性等优点。但它也存在一定的缺点，即二次回归预测值的方差随试验点在因子空间的位置不同而呈现较大的差异。由于误差的干扰，就不易根据预测值寻找最优区域。为了克服这个缺点，人们通过进一步研究，提出了回归的旋转设计（whirlydesign）。,1旋转设计的基本原理,1.1回归设计的旋转性,所谓旋转性是指试验因素空间中与试验中心距离相等的球面上各处理组合的预测值的方差具有几乎相等的特性，具有这种性质的回归设计称回归旋转设计。利用具有旋转性的回归方程进行预测时，对于同一球面上的点可直接比较其预测值的好坏，从而找出预测值较优区域。,.,3,如何才能使试验设计具有旋转性呢？这就需要弄清楚旋转性对试验设计有什么要求以及获得旋转性必须满足哪些基本条件。首先必须明确的是：在旋转设计中，试验处理的预测值的方差仅与因素空间中从试验点到试验中心的距离有关而与方向无关，从而克服了通常因为不知道最优点在什么方向的缺陷。,.,4,这里应该解决的是二次回归正交的旋转性问题。下面以试验设计中常用的三元二次回归方程来讨论这个问题。,1旋转设计的基本原理,在3个变量情况下，二次回归模型为：,即,它的结构矩阵为：,.,5,1旋转设计的基本原理,此外，为了使旋转设计成为可能，还必须使信息矩阵A不退化（满秩）。为此，必须有不等式,(1330),式(1330)就是m元二次旋转设计的非退化条件。已经证明，只要使N个试验点不在同一个球面上，就能满足非退化条件。,最简单的情况是把N个试验点分布在2个或3个半径不等的球面上。如m0个点分布在半径为0的球面上（即在中心点重复m0次试验），另外m1Nm0个点均匀分布在半径为(0)的球面上。,.,6,1旋转设计的基本原理,综上所述，为了获得m元二次旋转设计方案，就要求既要满足旋转性条件式(1329)，又要满足非退化条件式(1330)。满足条件式(1329)是旋转设计的必要条件，满足非退化条件式(1330)是使旋转性成为可能的充分条件。两者结合起来才能使旋转性设计得以实现。实际操作上主要借助于组合设计来实现。因为组合设计中N个试验点Nmc+m+m0，分布在3个半径不相等的球面上。即,mc个点分布在半径,的球面上；,m个点分布在半,的球面上；,m0个点分布在半径,的球面上；,因此，采用组合设计选取的试验点，完全能够满足非退化条件式(1330)，即信息矩阵A不会退化。此外，采用组合设计，其信息矩阵A的元素中,而它的偶次方元素,均不等于零，完全符合式(1329)的要求。,.,7,为了获得旋转设计方案，还必须根据旋转性条件式(1329)确定值，,1旋转设计的基本原理,事实上只要,求出值就行了。,在组合设计下，当mc2m（全实施）时，则前式变为,解此方程，即可建立全实施时值的计算式，即,(1331),同理,当,当,当,.,8,1旋转设计的基本原理,表1324二次正交旋转组合设计参数表,为了便于设计，现将m个因素不同实施情况下的值列于表1324。,.,9,2次旋转组合设计具有同一球面预测值的方差相等的优点，但回归统计数的计算较繁琐。如果使它获得正交性就能大大简化计算手续。,1旋转设计的基本原理,1.2正交性的获得,在2次旋转组合计划中，1次项和交互项的回归系数bi和bij仍保持正交，但b0与bij之间，以及bii与bjj之间都存在相关，即不具正交性，它们之间的协方差分别为：,(1332),其中,.,10,1旋转设计的基本原理,同样，对于m元二次旋转组合设计，上式中的mc和m也都是固定的。这样就只能通过调整中心点的试验处理数m0使4/221。由此可见，适当地选取m0，就能使2次旋转组合设计具有一定的正交性。为了方便设计，已将m元不同实施的m0和N列入表1324中。,综上所述，只要对平方项施行中心化变换，并适当调整就能获得二次正交旋转组合设计方案，这方面的计划见表1327和表1328。,对于m个因素的二元旋转组合设计，式（1333）中的m、mc和都是固定的。因此，只有适当地调整N才能使4/221，而试验处理数,Nmc+m+m0,.,11,二次回归旋转组合设计，具有同一球面上各试验点的预测值的方差相等的优点，但它还存在不同半径球面上各试验点的预测值的方差不等的缺点。为了解决这一问题，于是提出了旋转设计的通用性问题。所谓“通用性”，就是试验除了仍保持其旋转性外，还具有各试验点与中心的距离在因子空间编码值区间01的范围内，其预测值的方差基本相等的性质，即同时具有旋转性与通用性。这种设计称为通用旋转组合设计。如何才能满足其通用性呢？,1旋转设计的基本原理,1.3二次旋转组合设计的通用性,首先来看预测值的方差，已知在m个因素情况下，其预测值的方差,(1334),此式是在21的约定下得到的，这种约定并非本质的，只是为了讨论简单起见。由此可知，只有恰当确定4，才能满足通用性的要求。,.,12,1旋转设计的基本原理,(1335),那么，对4有什么要求呢？总的来说，它必须使式中D()在诸i(01)区间的内插点）处的值与1处的值的差的平方和为最小，即：,式中,于是，对于不同的m，均可计算出满足式（1335）的4,.,13,1旋转设计的基本原理,当4确定后，由关系式（见1333）可以计算出不同m的试验处理数N。,当计算结果不是整数时，N可取其最靠近的整数。然后再由m0N-mc-m,计算出不同m值的m0，上述计算结果列于表1325。,表1325二次通用旋转组合设计参数表,.,14,从以上可以看出，正交旋转的好处在于正交性，它是通过增加中心点的试验次数换来的，但有时并不合算。在某些实际问题中，反倒不如选用通用旋转设计。因为通用旋转设计，既能在01的较实用区域使方差D()基本不变，又在一定程度上减少了试验次数。,1旋转设计的基本原理,从上述讨论结果看出，为了满足通用性要求，主要在于确定出适当的m0。因此，只要在中心点安排如表1325所列的m0次试验旋转组合设计便获得通用性。,.,15,2二次正交旋转组合设计及其统计分析,设研究因素为m个，分别以Z1,Z2,Zm,表示。在进行设计时，首先确定每个因素的上、下水平，进而计算零水平，以及变化间距。某因素零水平及变化间距的计算式为,2二次正交旋转组合设计及其统计分析,2.1二次正交旋转设计的一般方法,Z0j(Z1j+Z2j)/2,j(Z2jZ0j)/,式中为待定参数，其值可以从表1324中查出。,对每个因素Zj各水平的取值进行线性变换，以实现其编码,xj(ZjZ0j)/j,这样，就将有单位的自然变量Zj变成了无单位的规范变量xj(j1,2,m)，并可编制出因素水平的编码值表（表1326）。,.,16,表13-26二次正交旋转设计因素水平编码值表,试验因素Z1,Z2,Zm经因素水平编码后，以变量x1,x2,xm表示，选用适当的二水平正交表，即可设计出二次回归正交旋转组合方案。,为了方便设计与统计分析，现将常用的二因素和三因素二次正交旋转组合设计的结构矩阵列于表1327和表1328。,2二次正交旋转组合设计及其统计分析,.,17,表13-27二元二次正交旋转组合设计的结构矩阵,2二次正交旋转组合设计及其统计分析,.,18,表13-28三元二次正交旋转组合设计的结构矩阵,（未完）,2二次正交旋转组合设计及其统计分析,.,19,（续前表13-27）,二次回归正交旋转组合设计试验结果的统计分析，与二次回归正交组合设计试验结果的统计分析方法相似，这里不再赘述。,2二次正交旋转组合设计及其统计分析,.,20,2.2三因素（1/2）实施正交旋转组合设计示例,例133采用三因素二次正交旋转设计组合设计，其试验因素水平编码见表1329。,表13-29试验因素水平编码表,试验结果及统计分析如下：,2二次正交旋转组合设计及其统计分析,.,21,表13-30三因素二次回归正交旋转组合设计结构矩阵与结果计算表,（1）建立回归方程。三因素二次回归正交旋转组合设计结构矩阵与结果计算见表1330。初步得回归方程为：,（未完）,2二次正交旋转组合设计及其统计分析,.,22,初步得回归方程为：,（续前表13-30）,2二次正交旋转组合设计及其统计分析,.,23,表13-31三因素二次回归正交旋转组合设计试验结果方差分析表,（2）回归方程的显著性测验：对所得三元二次回归方程进行方差分析，见表1331。,2二次正交旋转组合设计及其统计分析,.,24,剔除x1,x3,x1x2,x1x3,x2和x3，回归方程变为：,将中心化变换还原为xj2,得：,此时，,2二次正交旋转组合设计及其统计分析,.,25,3通用旋转组合设计及其统计分析,通用旋转组合设计与正交旋转组合设计基本相同，其组合计划中试验处理组合数N，也是由3部分组成，即：,3通用旋转组合设计及其统计分析,3.1通用旋转组合设计的一般方法,Nmc+m+m0,上式中mc和m的数值与正交旋转组合设计完全相同，只是N和m0有所不同，其值可从表1225查出。,现将常用的三因素二次通用旋转组合设计的结构矩阵列于表1332。,.,26,表13-32三元二次通用旋转组合设计的结构矩阵,3通用旋转组合设计及其统计分析,.,27,3通用旋转组合设计及其统计分析,3.2通用旋转组合设计试验结果的统计分析,（1）建立二次回归方程。要建立回归方程，必须计算出回归系数，而回归系数,b(XX)-1(XY),式中：(XX)-1为设计的相关矩阵；(XY)为常数项矩阵B，在通用旋转设计下有：,.,28,3通用旋转组合设计及其统计分析,所以回归系数,(1336),式（1336）中K、E、F、G的值如表1333所示。,表13-33二次通用旋转组合设计K、E、F、G值表,.,29,3通用旋转组合设计及其统计分析,注：令,则,由式（1336）计算出回归系数b，即可建立二次多项式回归方程。,.,30,3通用旋转组合设计及其统计分析,（2）回归方程的显著性检验。,计算平方和及自由度：如果m元二次通用旋转组合设计的N个试验结果以y1,y2,yN表示，则各项平方和及其自由度为：,(1337),在通用旋转组合设计中，一般中心点均需做重复试验。如果重复次数为m0试验结果以y01,y02,y0m0表示，则它们的误差平方和及其自由度为：,(1338),可由误差项与剩余项比较计算失拟平方和及其自由度：,(1339),.,31,3通用旋转组合设计及其统计分析,失拟性检验：失拟性可用统计量,FLfF0.05，表示差异不显著，可直接对回归方程进行显著性检验；如果FLfF0.05，差异显著，则表明存在影响试验结果的其他不可忽略的因素，需要进一步考察其原因，改变二次回归模型。,(1340),回归方程的显著性检验：,(1341),进行显著性检验，如果FRF0.05，则回归关系不显著，说明此回归方程不宜应用；如果FRF0.05和F0.01，则回归关系显著或极显著，表明此回归方程可以应用。,.,32,3通用旋转组合设计及其统计分析,（3）回归系数的显著性检验。当FLf检验结果不显著时，回归方程中各变量作用的大小，可通过t检验来判断。为此，需要计算各回归系数的t值，其计算式为：,(1342),式（1342）中K、mc、F、e已如前述（参见表1334）。,.,33,3通用旋转组合设计及其统计分析,3.3四元二次通用旋转组合示例,例134鸡肉乳酸发酵试验，对鸡肉乳酸发酵的产酸条件进行优化试验，采用二次通用旋转组合设计对盐浓度、糖浓度、发酵温度和发酵时间进行试验，采用四元二次通用旋转组合试验寻求最优发酵条件，试验因素及水平编码见表1334.,表13-34鸡肉乳酸发酵产酸条件的四元二次通用旋转组合设计因素水平表,试验设计方案和试验结果见表1335。,.,34,3通用旋转组合设计及其统计分析,表13-35鸡肉乳酸发酵产酸条件的四元二次通用旋转组合设计方案及结果,（未完）,.,35,3通用旋转组合设计及其统计分析,（续前表13-35）,.,36,3通用旋转组合设计及其统计分析,（1）建立四元二次回归方程。根据计算，可建立四元二次多项式回归方程（计算从略）。,.,37,3通用旋转组合设计及其统计分析,（2）回归方程的显著性检验。对鸡肉乳酸发酵产酸条件数学模型的方差分析见表1336。,.,38,3通用旋转组合设计及其统计分析,从方差分析可以看出，回归达到极显著水平。说明本试验设计及分析效果都很好，各因素间显著与不显著也泾渭分明。因此没有必要做二次回归方差分析，可直接将F1的回归系数去掉而得到含酸量与各因素间的回归方程为：,.,39,SourceDFSeqSSAdjSSAdjMSFPRegression37.7897.7892.59621.080.387Linear37.7897.7892.59621.080.387ResidualError1638.59738.5972.4123Lack-of-Fit1136.05736.0573.27796.450.026PureError52.5402.5400.5079Total1946.385S=1.553R-Sq=16.8%R-Sq(adj)=1.2%,输出结果：线性回归方差分析表,此值很小说明线性回归效果不好,此值小于0.05时表示线性回归模型不正确,此值大于0.05时表示回归的效果不显著,线性回归结果,.,40,SourceDFSeqSSAdjSSAdjMSFPRegression936.46536.4654.05174.080.019Linear37.7897.7892.59622.620.109Square313.38613.3864.46194.500.030Interaction315.29115.2915.09705.140.021ResidualError109.9209.9200.9920Lack-of-Fit5