第四章回归设计与方分析.doc

第四章 回归设计与分析 古典回归分析：只能对已收集到的数据作被动处理；对如何安排试验很少提出合理的要求；对所求的回归方程的精度很少研究。 回归设计与分析：五十年代由正交试验法与回归分析结合而产生；主动地把试验的安排、数据的处理和回归方程的精度统一考虑与研究。目的在于选用较少的适当的试验点，获得较高精度的回归方程。 分类：按类型划分：正交设计、旋转设计、最优设计、均匀设计等 按次数划分：一次设计、二次设计等 第一节 回归正交设计 一、一次回归正交设计 利用二水平正交表如L4(23)、L8(27)、L12(211)等进行设计。 设试验研究m个因素：Z1，Z2，…，Zm；Z1j，Z2j为因素Zj的下水平与上水平。 (一) 对因素的水平进行编码 表4-1 编码表 编码xj Z1 … Zj … Z m 下水平 -1 Z11 … Z1j … Z1m 上水平 +1 Z21 … Z2j … Z2m 零水平 0 … … 变化间隔 Δj … … 编码的几何解释： 对Z1、Z2、…、Zm的回归设计问题转化为对x1、x2、…、xm的回归设计问题。 (一)选用合适的正交表列出编码因素的试验方案(即试验设计)，利用试验结果建立回归方程。 例如有3个因素x1、x2、x3 (已编码)无互作，欲建立 =b1+b1x1+b2x2+b3x3 (b0, b1, b2, b3待计算) 选用 列号 试验号 1 2 3 4 5 6 7 1 1 1 1 1 1 1 1 2 1 1 1 2 2 2 2 3 1 2 2 1 1 2 2 4 1 2 2 2 2 1 1 5 2 1 2 1 2 1 2 6 2 1 2 2 1 2 1 7 2 2 1 1 2 2 1 8 2 2 1 2 1 1 2 并将 中的“1”改为“+1”、“2”改为“-1”，且把x1、x2、x3放在1、2、4列上，可得编码因素的试验方案： 表4-2 试验设计表 试验号 x1 x2 X3 试验结果ya 1 1 1 1 y1 2 1 1 -1 y2 3 1 -1 1 y3 4 1 -1 -1 y4 5 -1 1 1 y5 6 -1 1 -1 y6 7 -1 -1 1 y7 8 -1 -1 -1 y8 试验方案具有正交性：即 ；-1，+1不仅表示因素的状态，还表示变量xj的取值；数据结构式： 结构矩阵： 待估参数 观测值： 列向量 列向量： 则根据最小二乘原理 (XˊX)b=XˊY, b=(XˊX)-1(XˊY) 而 表明b0，b1，b2，b3相互独立， 这是正交设计的优越性之一。 记 则 b0=B0/8， b1=B1/8， b2=B2/8， b3=B3/8。 若x1、x2、x3间有互作，欲建立 =b0+b1 x 1+b2 x2+b3 x3+b12x 1x2+b13x1x3+b23x 2x3 数学模型为：y=β0+β1x1+β2x2+β3x3+β12x1x2+β13x1x3+β23x2x3+ε 此时矩阵形式： Y=Xβ+ε 仍有 (XˊX)b=XˊY， b=(XˊX)-1(XˊY) 表明，b0，b1，…b2 相互独立 记B=XY=(B0 B1 B2 B3 B12 B13 B23)ˊ 则b0=B0/8, b1 =B1/8, b2=B2/8, b3=B3/8, b12=B12/8, B13=B13/8, b23=B23/8 (三) 显著性检验(对于有互作的情况) 1、回归方程显著性检验 划分式： 在一次回归设计下