《非负回归系数的线性回归模型的构建》由会员分享，可在线阅读，更多相关《非负回归系数的线性回归模型的构建（8页珍藏版）》请在金锄头文库上搜索。
1、非负回归系数的线性回归模型的构建谢忠秋 （江苏常州，江苏技术师范学院经济管理系，213001） 摘要：本文针对某些线性回归模型负的回归系数不具有实际的物理意义和经济意义的问题， 提出了非负回归系数的线性回归模型构建的新方法。与现有的方法相比较，该方法具有简单 和易操作的特点。在实际中具有一定的应用价值。关键词：非负回归系数 非负回归系数线性回归模型一、引言在许多的实际问题中，必须要求线性回归系数为非负，否则没有实际意义。例如在方开 泰等人研究的配方问题中就是如此。在配方问题中，每个成分的线性回归系数相当于它在总 配方中的比例，若线性回归系数为负就失去了物理意义（1982，1985）。而事实上，
2、在用线 性回归模型反映实际经济问题时，我们也会常常遇到这样的情况，就是线性回归模型的回归 系数从经济意义上进行阐释，应为正数，但根据最小二乘法所确定出的线性回归模型的某一 个或某几个回归系数却偏偏为负数，从而出现了其计算结果与经济分析相互矛盾的情形。对 此种情况，通常的做法就是运用一定的方法将该变量从回归模型中消除。如果说，对于不太 重要的变量，这样做还没有多大影响的话，那么对于一些非常重要的变量，这样做就可能会 产生一些消极的作用，如回归模型失真等等。那么，能否在遵守一定规则的前提下，通过某 种科学的方法使负的回归系数转变为正的回归系数呢？这就是有关非负回归系数的线性回 归模型的构建问题。显
3、然，我们有必要这一问题进行研究，这将有助于人们更好地利用线性 回归模型对实际物理问题和实际经济问题作出合理的阐述和解释。对非负回归系数的线性回归模型的构建问题，Waterman在1974年曾进行过讨论，他建 议用一切可能回归的办法来求最小二乘估计；我国学者方开泰、王东谦、吴国富也在 1982 年对这一问题展开研究，提出了立足于矩阵的消去变换方法；而后方开泰、贺曙东又在1985 年对这一方法作了进一步的改进，使这一方法的应用更具一般性（1）、（2）。显然，他们的贡 献是巨大的。但也应该看到，这些方法在实际应用上还存在着一定的缺陷，最大的缺陷就是 它们的计算过程还相当复杂，而这对于很多的人来说，无
4、疑是一道难以逾越的障碍。本文试 图给出另一种解决方案。该方案对于熟悉统计学和稍有数学基础的人是容易接受和掌握的。二、研究思路和方法设多元线性回归模型为：Y * = b * + b * X + b * X + b * X(2.1)0 1 1 2 2 n n（1） 、方开泰、王东谦、吴国富 一类带约束的回归配方回归计算数学1982 年第1 期 p5769（2）方开泰、贺曙东含有线性约束及非负回归系数的回归模型计算数学1985年第 3期 p237246对于非负回归系数的线性回归模型来说，如果根据最小二乘法所求出的b*、b*、b*1 2 n 俱能在经济意义上解释得通，则意味着上述模型符合要求，直接应用
5、即可。而事实上，在许 多实际的线性回归模型中，当某一个或某几个变量的回归系数表现为负数时，往往表现出与 实际经济意义的解释并不相符的情形。换句话说，只有表现为非负时，才具有经济意义上的 合理性。如果出现这种情况，则意味着该模型不符合非负回归系数线性回归模型的要求。此 时，一个简单而自然的想法就是根据某些规则通过一些变换，将负的回归系数（注意：并不 是所有的负的回归系数，而只是那部分不具有经济意义上的合理性的负的回归系数） 转变 为非负的回归系数。有关变换规则及过程如下：1、对不具有经济意义上合理性的负的回归系数（为简便计，不妨直接称之为负的回归系数，下同）乘以一个负的调整系数，使之变换为非负回
6、归系数。设非负回归系数为 b ，i负的回归系数为b*，负的调整系数为p，则b = p b *， i = 1,2,m。（2.2）i1i1 i2、由于b*变换为b，则必然引起具有经济意义上合理性的回归系数（为简便计，也不ii妨直接称之为正的回归系数，下同）的回归系数的变动。设变动后的正的回归系数为 b ，j变动前的正的回归系数为 b* ， 正的回归系数调整系数为 p ， 则 b = p b* ， j 2 j 2 jj = m = 1， m + 2,n （2.3）3、同样，由于b*变换为b，b*变换为b，也必然会引起截距项的变动。设变动后的 iijj截距项为b，变动前的截距项为b*。004、 为保证
7、p为负，p为正，令= -1，即p = -p o （2.4）12p1225、不管回归系数如何变化，总是要求变化后的剩余变差与变化前的剩余变差相等，即工（Y f ）2 =工（Y Y *）2 =工 e2（2.5）于是，非负回归系数线性回归模型为：+ b X + + b X + b X + + b X22m mm +1m +1n n2.6)b = p b * i = 1,2,mi1 ib = p b * j = m + 1,nj2 jb 0i用最小二乘法，则要求：min( Q) = X (Y b b X b X b X b X b X )2(2.7)01 12 2m mm +1 m + 1n n于是，
8、(2.6)式就可变换为下列模型：min( Q)=工(Y - b - b X - b X -b X - b X - b X )201122m mm +1m +1n n2.8)b - p b * = 0 i = 1,2,，mi1 ib - p b * = 0 j = m + 1, nj2 jb 0i（2.8）式的含义是在（2.2）式、（2.3）式条件下求（2.7）式的最小值。显然，对于（2.8） 式可以用拉格朗日乘数求条件极值的方法求出该模型的解。构成函数：min( Q) = S (Y - b - b X - b X b X - b X b X )201122m mm 十 1 m 十 1nn+ 九
9、(b- p b *) +九(b- p b *) + + 九 (b -p b * )111 1221 2m m 1 m+ 九 (b - p b * ) + 九 (b - p b * ) + + 九(b - p b *)m +1 m +12 m + 1m + 2 m +1 2 m + 2n n 2 n求其对b、b、b的偏导数，并使之为零，得到：0 ijab= (-2)S (Y - b -0bX -1bX22mmbX- -b X )(Xnn=0ab(-2)S (Y - b - b0X11X-2mmm + 1m +1 b X )(X ) + 九=0nna Q=(-2)工(Y - b-b X-b X -
10、bX-b X-bX)=0a b001122mmm + 1m + 1nna q= (-2)S (Y - b-bX-bX -bX-b X-bX)(X ) + X = 0ab101122mmm + 1m + 1nn11a q=(-2)工(Y - b-bX-bX -bX-b X-bX)( X ) + X = 0ab201122mmm + 1m + 1nn22a q-=(-2)工(Y - b-bX-bX-bX-b X-bX)( X ) + X = 0a b01122mmm + 1m + 1nnmmm分别对上述各式进行整理得：SY-nb -bSX -b SX0 1 1 2 2-bSX-b S X = 02
11、.9)SXY-bSX -bSX2-bSX X212 -b S X X -b S X X -b S X X1m1nSX Y-bS20-bS1X X -b S X 21 2 2 -b SXX2mXX2 m + 1 -b S=0XX2nb XX m - = 0 n m2E X Y - b E X - b E X X - b E X X b E X 2 - b E X Xm0m 11 m 22 mmmm+1m+1mEX Y-bm +10EX -bEXX -bEX X11 m+122XXm m+1-b E X2m +1EXYnX X -b E X X1 n 22XXm-b Em+1X X - - b m
12、 +1n显然，要求出b、b、b、b、0bm+2的值，关键在于求出 p 、 p 。12由（2.5）式工（Y即：工 Y 2 - b0 工 Y -bi 工 x 1Y 一 b2 工 X 2 -bm 工X Y -bmX ,Yb E X Y = E e2m + 1n n将（2.2）式和（2.3）式代入上式得：工 e2 = E Y2 A 工 Y n (b* *工 X,Y + b*工 X Y + + b* 工01 1122mY)- p (b* E X Y + + b*Em 2m + 1m + 1nX Y ) (2.12 ) n将（2.4）式代入（2.12）式，则有：工 Y 2 - b0 工 Y - P(b：工
13、 X Y + b；工 X 2Y + + b；工XmY)+ p1(bm*+1X Y + +b*Em + 1n2.13)进一步整理得：E e2 = E Y 2 - b0 E Y - pj(b； E X / + b； E X 2Y + + b： EX Y)- (b*Emm + 1XnY)2.14)X Y) (2.15 )解(2.14)得：工Y将（2.2）式、（2.3）式和（2.15）代入（2.9）得E Y 2 - E e 2 - p (b E Y 2 - E e 2 - p (b * E X Y + b * E X Y + + b * E X Y) - (b * E X Y + + b * E X Y) E Y = nl E X Y + b * E X Y + + b *+ p ( b * E X +1 1 1b*E2X + +b*E X )2mm(b * Em + 1+ + b*E X )nn移项并整理得：E X Y + + b * E X Y)m + 1nn* E Y E X + + b * E Y E X )m + 1m + 1nnn(E Y 2 -E e2)- (E Y)2=p n( b * E
