高级生物统计（7）

《高级生物统计（7）》由会员分享，可在线阅读，更多相关《高级生物统计（7）（101页珍藏版）》请在金锄头文库上搜索。

1、第二节 回归旋转设计,二次回归正交设计具有试验规模小、计算简便和避免了回归系数间的相关性等优点。然而，它与一般回归分析一样, 试验点在因子空间的位置不同（各因素所取水平不同），对应的各个预测值的方差也就不相同，致使设计在各个方向上不能提供等精度的估计,因此不能对不同试验点预测值之间进行直接比较，不易寻找最优区域。 为克服这一缺点，本节介绍回归旋转设计(regression-rotable design)。,一、旋转性、旋转设计与旋转性条件,（一）旋转性 设有回归方程 可用预测值方差来评价其“精确度(precision)”，即 可见， 与试验点 在空间的位置有关，且D(bj)与cov(bi,bj

2、)和有关，从而与结构矩阵有关。,设用一次回归正交设计所求得的回归方程为,一次回归正交设计的系数矩阵A和相关矩阵C为,于是，一次回归正交设计所得回归系数的方差和协方差为：,所得一次回归方程预测值的方差为：,(4-27)式表明，位于同一球面上的点的预测值 的方差是相等的。这个性质称为旋转性(rotability)。当利用具有旋转性的回归方程进行预测时，对于同一球面上的点可直接比较其预测值的好坏，从而容易找出预测值相对较优的区域。,误差方差,是m维编码空间内的一个球面，球心在原点，半径为,凡与试验中心点距离相等的球面上各点回归方程预测值 的方差相等的回归设计称为旋转设计(rotation desig

3、n)。 显然，一次回归正交设计具有旋转性。一般二次回归正交组合设计不具旋转性。,（二）旋转设计,回归旋转设计，一方面基本保留回归正交设计的优点：试验次数较少，计算简便，且部分地消除了回归系数间的相关性；另一方面能使二次回归设计具有旋转性。,通过回归的旋转设计使所得回归方程具有了旋转性，这样既有助于克服多元线性回归及二次回归正交组合设计中回归预测值 的方差依赖于试验点在因子空间的位置这个缺点，又可以简单地用的大小表示回归预测值误差的大小，小（试验点距离中心近）误差小，大（距离中心远）误差大。,（三）旋转性条件,下面我们从3个自变量的二次回归方程着手来说明这个问题。,对于m=3，二次回归数据结构式

4、是,共有 个待估计参数。,结构矩阵,对应的信息矩阵 A,由此可见，在三元二次回归中，信息矩阵（系数矩阵）A中元素的一般形式是： 其中指数a1，a2，a3分别可取，等非负整数，但是这些指数的和不能超过，即,例如，当a1=a2=a3=0时，（4-28）就是矩阵A的第行第列上的元素N。仔细地观察，还可把系数矩阵的元素分为两类：一类元素，它的所有指数a1，a2，a3都是偶数或零；另一类元素，它的所有指数a1，a2，a3中至少有一个奇数。,在一般的m元d次回归中，共有 项，对应的信息矩阵A是 阶对称方阵，A的元素的一般形式是: 其中指数a1,a2,am分别可取,2d等非负整数，且满足0 a1+a2+am

5、 2d。A的元素亦可类似地分为两类。在旋转设计中，对这两类元素的值的要求，归纳成著名的G. E. P. Box旋转定理。,定理(旋转性条件),m元d次回归设计满足旋转性的充要条件是其对应的信息矩阵的元素 其中指数a1,a2,am是如上所述的非负整数，为试验次数，a= a1+a2+am ，a为待定参数，它的下标a一定是偶数，特别是0。,二、一次回归旋转设计,对于 结构矩阵X和信息矩阵A,根据（4-29）式，计算得,矩阵A中其它元素为0。故，三元一次回归旋转设计信息矩阵A为：,当2=1时， 表明，一次回归正交设计，也就是2=1时的一次回归旋转设计。,三、二次回归旋转设计,(一) 二次回归旋转设计条

6、件 当m=3，d=2信息矩阵见前，其中,所以三元二次回归旋转设计的信息矩阵的元素有如下特点（三元二次回归旋转性条件）：,A的其它元素皆为，于是,0 1 2 3 12 13 23 11 22 33 2、4待定。,一般，m元二次回归旋转设计的信息矩阵的元素有如下特点（m元二次回归旋转性条件）：,这时，二次旋转设计信息矩阵有如下形式（其中空白处为零）：,经计算，有：,其中 。由此可见，要使|A|0，即矩阵A为非退化的，必须要有： (4-31)式称为m元二次旋转设计的非退化条件(non-degenerative conditions)。,1、非退化性的实现,可以证明只要N个试验点至少分布在两个半径不等

7、的球面上就可以满足非退化条件（4-31）。最简单的情况是把N个试验点分布在两个或三个球面上。,2、旋转性的实现,在组合设计中，N个试验点 分布在个球面上: mc个点分布在半径为 的球面上；m=2m个点分布在半径为 的球面上；m0个点集中在半径 的球面上。因而组合设计总是满足非退化条件（4-31）的。通过调整星号臂的值可以使组合设计满足旋转性条件（4-30）。,在组合设计中,而偶次方元素 都不为零，为了满足旋转性条件(4-30)，即满足 就行了。,在mc=2m（全因子实施）的情况下,要使 成立，即,对于mc=2m-1和mc=2m-2(全面试验的1/2或1/4实施),由(4-32)式或（4-33）

8、式算出的常用的和2值列于表4-21。,(二) 二次回归旋转组合设计中m0的选择,二次回归旋转组合设计具有同一球面预测值 的方差相等的优点，但b0与bjj，bii与bjj间存在相关性。 通过适当地选取m0，能使二次旋转组合设计具有更好的统计性质。,1、使二次回归旋转组合设计具有正交性,（1）使二次回归旋转组合设计具有几乎正交性 所谓几乎正交性（near orthogonalty）是除b0与bjj相关外，其他回归系数之间不存在相关，此处，就是要解决bii与bjj之间的相关问题。由保证cov(bii,bjj)来确定N，进而确定m0。 表4-22提供了进行二次回归几乎正交旋转组合设计的各种参数。,例如

9、对m=3的情况，可选m0=9来进行二次回归几乎正交旋转组合设计。三因素二次回归几乎正交旋转组合设计结构矩阵如表4-23所示。,(2) 使几乎正交成为完全正交,使几乎正交成为完全正交，即解决b0与bjj相关的问题。对平方项列的元素施行中心化变换，即令 便可消除b0与bjj之间的相关，进而获得二次回归正交旋转组合设计结构矩阵。,例如，对于m=3，d=2,因为 所以 三因素二次回归正交旋转组合设计结构矩阵如表4-24所示。,2、使二次回归旋转组合设计具有通用性,二次回归旋转组合设计，具有同一球面上各试验点的预测值 的方差相等的优点，但它还存在不同半径球面上各试验点的预测值 的方差不等的缺点。为了解决

10、这一问题，提出了回归设计的通用性(generality)问题。所谓通用性是指各试验点与中心的距离在01范围内，其预测值 的方差基本相等。具有通用性的设计称为通用设计(common design)，亦称等精度设计(equal precision design)。同时具有旋转性与通用性的组合设计称为通用旋转组合设计(common rotation design)。,表4-25提供了进行二次回归通用旋转组合设计的各种参数。 此处的m0比二次回归几乎正交旋转组合设计m0的小。 表4-26是三因素二次回归通用旋转组合设计结构矩阵。,在二次回归通用旋转组合设计中，常数项b0与平方项回归系数bjj，平方项回

11、归系数bii与bjj（ij）间还存在着相关，所以说通用旋转组合设计是损失了部分正交性而达到了单位球内基本一致精度的要求。,综上所述，m元二次回归旋转组合设计包括：二次回归几乎正交旋转组合设计，二次回归正交旋转组合设计和二次回归通用旋转组合设计。组合设计保证了非退化性，的选择保证了旋转性，m0的选择保证了通用性或几乎正交性，将平方项 列的元素中心化保证了几乎正交成为完全正交。,简而言之，旋转组合设计使预测值 有等距等方差的性质；正交旋转组合设计消除了各项回归系数间的相关性；通用旋转组合设计以损失部分正交性为代价保证了预测值 在单位球内基本等方差。,四、二次回归旋转组合设计的统计分析,二次回归旋转

12、组合设计的统计分析方法，因该设计是具有正交性、几乎正交性还是通用性而有所不同。 下面先介绍分析的基本步骤与有关计算公式，而后通过实例介绍各种分析方法。,(一) 基本公式与步骤,1、 因素水平编码 m个因素：Z1、Z2、Zm； Z2j、Z1j分别为Zj的上、下水平 (j=1, 2, , m)； Z0j=(Z2j+Z1j)/2； 可利用方法 I 或方法 II 编码。,2、列出试验实施方案,编码因素的试验方案即试验设计包含N=mc+m+m0 个试验点。将编码因素试验方案的编码水平换成相应的实际水平即得试验实施方案。试验结果记为y(=1,2, ,N)。,3、 计算回归系数,二次回归正交旋转组合设计回归

13、系数的计算同二次回归正交组合设计； 二次回归通用与几乎正交旋转组合设计回归系数的计算如下：,对于二次回归组合设计，(XX)-1可表示为,对于二次回归通用（G0，E0）与几乎正交（G0，E0）旋转组合设计则有：,其中 即 其中，K、E、F、G、e-1、mc-1的值可查表4-27。横线上方数值用于通用旋转组合设计，横线下方数值用于几乎正交旋转组合设计。,4、回归方程显著性检验,对于正交旋转组合设计,对于通用、几乎正交旋转组合设计,失拟性检验,5、回归系数显著性检验,二次回归正交旋转组合设计的各偏回归系数的显著性检验与二次回归正交组合设计的方法相同，即： 各偏回归系数的显著性检验也可用 t 检验（见

14、本节【例44】）。,对于通用、几乎正交旋转组合设计，回归方程中各偏回归系数的显著性检验可用 t 检验（或 F 检验）。此时，因各偏回归系数（包括常数项b0）的方差分别为：,或者,当FLf不显著时 当FLf显著时,【例43】,研究锌肥（Zn，硫酸锌）、氮肥（N，硝酸铵）、磷肥（P，三料磷）对玉米产量的影响。采用二次回归几乎正交旋转组合设计建立因素与产量指标间的二次回归方程。 根据专业要求，Zn、N、P三种肥料的上、下水平分别取：4.0、1.0；55.0、14.0；17.5、4.5。,（1）因素水平编码,因为影响产量指标的因素有3个，即m=3，采用二次回归几乎正交旋转组合设计，查表4-22，得=1

15、.682，2=2.828 。利用方法对因素水平编码，见表4-28。,（2）列出试验实施方案,由表4-22查得，mc=8，2m=6，m0=9，N=mc+2m+m0=23 。试验设计与实施方案见表4-29，结构矩阵与试验结果计算见表4-30。,（3）求出常数项b0及各偏回归系数,先求出B0，Bj，Bij，Bjj (i,j=1,2,3)，见表430最后一行， 以及T=Bjj=32352.588。 查表4-27（横线下方）得,于是,得回归方程为,（4）回归方程及回归系数的显著性检验, 回归方程显著性检验,方差分析见表4-31。, 回归系数显著性检验 t检验,因为 所以,将不显著的一次项x1、x3，乘积项x1x2剔除得：,将上述方程中编码变量 xj 用实际变量 Zj (xj=(Zj-Z0j)/j)表示，则,【例44】,对【例43】采用二次回归正交旋转组合设计建立因素与产量指标间的二次回归方程。,此时，因素水平编码、试验实施方案与【例43】相同。关于试验结果的分析，只需将表4-30中平方项列的元素中心化为： 其余计算与二次回归正交设计相同，结构矩阵与试验结果计算见表4-32。