均匀设计中有重复试验的统计分析.pdfVIP

均匀设计中有重复试验的统计分析 孙尚拱 北京医科大学 摘要：本文首先介绍有重复试验的回归分析理论及公式，再对一个例子重新给 以计算及提出问题；最后笔者对此类数据的处理及SAS统计软件包提出讨论 意见，其中特别对常规的用平均数去取代试验结果的方法及当变量数超过样本 数时用流行的逐步回归选取变量法提出批评。 关键词：加权回归分析；决定系数；逐步回归；SAS统计软件。 一、加权回归分析 设有P个自变量x，，x”．．，昂；因变量为Y，则一般的线性回归是： (1) Y=∥o+∥lxI+∥2x2+…+∥rXp=XI∥+占 叫残差。假设我们有厅z组的均匀设计，第i组的试验条件为x(o=(1，Xil，…xp)’， 它的，2，次重复试验结果为y¨…yiq，则(1)变成 7∥+国 ，21，…，垅 ，=1，…，12i (2) 肋=X(o 记，z=以l+…+，如。，及 y-◇11，…,Yt。一．．弘Jl，…舶“。)7 X’=(x(1)，…，x(1)；…；x(，，，)，…，x(n，)) 其中X及I，分别为，z×0+1)及以×l的阵，它们的前，zt行为第l组试验的 条件及结果，…，最后‰行是第_，7z组试验的条件及结果。这时(2)可就写成 Y=X口+占 (3) x 这里的占是由白组成的n l向量。一般的回归中常假定(2)中的函具有独 立性及N(o，盯2)的分布，其中仃2与试验号无关，称为方差齐性。而在均匀设计 中，由于X布满空间，不同的试验条件往往有很不相同的随机波动，因此方差 齐性的假定不易被满足。因此，我们假定(2)中的国满足独立性及 鼬～N(o，费) (4) 在(4)条件下的回归是最重要也最简单的加权回归分析。令，2×九的权阵为 W=diag(wl，…，W1，…，W胁…，wⅢ)盯2 其中∥的对角线对应于(4)中11个方差的倒数，即141，=l／a；2，而盯2是参数． 91 求b的方法是使下述加权残差平方和踬为最小 巩=(Y—t)’以Y一叠)=∑∑w,fy，，一元)2 (5) f=l，，I 令 72 Z=W1 y， a=W1”X， ，=W‘72占 (6) 这时(3)变为 z=卯+厂 (7) 这里厂是残差向量，它的每个分量的分布为N(O，O-2)，也就是说，虽然(3)是不 等方差，但经(6)变换后的(7)，则是方差齐性的回归分析，且(3)与(7)有相同的 回归系数及相同的残差平方和．这时(5)变为 SSo一(z—z)7(z—z) 幺=W172窜：Qb (8) 于是求(5)的最小值变为求2下的最小值．由～般回归公式得 (Q7e)b=QZ 即b满足 (x7∥∞6=X7WY (9) 由(9)求得6，即得．溉为 SSo=Y’WY一(y7∥∞6 (10) 被称之为回归平方和。但这两项中都舍有常数项60的波动影响．这个影响应 当被消去．方法：如果Y与P个变量全无关系，这时的回归模型即是 Y2风+5 (11)