多元线性回归

多元线性回归模型

一、多元线性回归模型的一般形式

设随机变量y与一般变量x1,x2,?,xp的线性回归模型为:

y??0??1x1??2x2????pxp?? 其中:

写成矩阵形式为:y?X???

?1?y1????1y2???y? X?????????y??n??1x11x21?xn1x12x22?xn2???x1p???0???1??????x2p?1?? ???? ???2?

?????????????xnp?????n???p??二、多元线性回归模型的基本假定

1、解释变量x1,x2,?,xp是确定性变量,不是随机变量,且要求

ran(kX)?p?1?n。这里的rank(X)?p?1?n表明设计矩阵X中自变量列之间

不相关,样本容量的个数应大于解释变量的个数,X是一满秩矩阵。

E(?i)?0,i?1,2,?,n????2,i?j2、随机误差项具有0均值和等方差,即:?

cov(?i,?j)??,(i,j?1,2,?,n)??0,i?j?E(?i)?0,即假设观测值没有系统误差,随机误差?i的平均值为0,随机误差?i的协方差为0表明随机误差项在不同的样本点之间是不相关的(在正态假定下即

为独立),不存在序列相关,并且具有相同的精度。

??i~N(0,?2),i?1,2,?,n2?~N(0,?In),3、正态分布的假定条件为:,矩阵表示:???1,?2,??n相互独立由该假定和多元正态分布的性质可知,随机变量y服从n维正态分布,回归模型的期望向量为:E(y)?X?;var(y)??2In 因此有y~N(X?,?2In) 三、多元线性回归方程的解释

对于一般情况含有p个自变量的回归方程E(y)??0??1x1??2x2????pxp的解释,每个回归系数?i表示在回归方程中其他自变量保持不变的情况下,自变量xi每增加一个单位时因变量y的平均增加程度。因此通常把多元线性回归的回归系数称为偏回归系数。下面看个例子,考虑国内生产总值GDP和三次产业增加值的关系,这个问题中GDP=x1?x2?x3是确定性的函数关系,可以看作误差项为

0的特殊回归关系。3个回归系数都是1,对?2解释为第二产业增加值x2每增加1亿元GDP也增加1亿元。假设做GDP对x2的一元线性回归,得到回归方程为

??5289.9?1.8554x2,对这个方程回归系数的解释是第二产业增加值每增加1y亿元GDP增加1.8554亿元。两个回归方程对同样的经济现象给出了不同的解释,问题出在什么地方呢?多元回归系数表示在回归方程中其他自变量保持不变的情况下,相应自变量每增加一个单位时因变量的平均增加速度。因此在用多元回归方程GDP=x1?x2?x3解释?2=1时,一定要强调是在x1和x3保持不变的情况下,

??5289.9?1.8554x2解x2每增加1亿元GDP也增加1亿元。在用一元回归方程y释回归系数时,要强调的是在方程之外的有关变量也相应变化时x2每增加1亿元GDP增加1.8554亿元。GDP增加的1.8554亿元中x2的直接贡献只用1亿元,回归方程外的x1和x3的贡献是0.8554亿元。这里又出现一个问题,为什么回归方程外的x1和x3贡献是0.8554亿元,而不是2亿元呢?可以通过考察数据,x2的增加幅度远大于x1和x3的增加幅度,假如x2增加1亿元,x1和x3相应的增加幅度都达不到1亿元。 四、参数估计

要想用OLSE估计多元线性回归模型的未知数,样本容量必须不少于模型中参数的个数。

在正态假定下,回归参数?的MLE(最大似然估计)与OLSE(最小二乘估计)

2?完全相同,即???(X?X)?1X?y,误差项方差?2的MLE为??L1nSSE?1n(e?e),这

是?2的有偏估计,但它满足一致性,在大样本的情况下,是?2的渐近无偏估计量。

参数估计量的性质:

性质1,??是随机向量y的一个线性变换 性质2,??是?的无偏估计 性质3,D(??)??2(X?X)?1

性质4,高斯-马尔科夫(G-M)定理

(1)c???是c??的无偏估计

(2)c???的方差要小

高斯-马尔科夫定理 在假定E(y)?X?,D(y)??2In时,?的任一线性函数

?,其中c是任一p+1维常数向量,??是?的c??的最小方差线性无偏估计为c??最小二乘估计。

此定理说明了用OLSE估计得到的估计量??是理想的估计量。关于这条性质,

需要注意以下四点:

第一,取常数向量c的第j(j?0,1,?,p)分量为1,其余分量为0,这时

?是?的最小方差线性无偏估计。 G-M定理表明最小二乘估计?jj

第二,可能存在y1,y2,?,yn的非线性函数,作为c??的无偏估计,比最小二

乘估计c???的方差更小。

第三,可能存在c??的有偏估计量,在某种意义(例如均方差最小)下比最

小二乘估计c???更好。

第四,在正态假定下,c???是c??的最小方差无偏估计。

性质5,cov(??,e)?0,在正态假定下??与e不相关等价与??与e独立,从而??与SEE=e?e独立。 性质6,当

?~N(?,?2(X?X)?1)??y~N(X?,?In)时,则?22?SEE?~?(n?p?1)2

五、自变量的显著性

如何剔除多余的不显著的自变量?y对自变量x1,x2,?,xp线性回归的残差平方和为SSE,回归平方和为SSR,在剔除掉xj后,用y对其余的p-1个自变量作回归,所得的残差平方和记为SSE(j),回归平方和为SSR(j),则自变量xj对回归的

贡献为:?SSR(j)?SSR?SSR(j),称为xj的偏回归平方和。由此可以构造偏F统计量:Fj??SSR(j)1SSE(n?p?1),当原假设H0j:?j?0成立时,偏F统计量Fj服从

自由度为(1,n-p-1)的F分布,此F检验与回归系数的t检验是一致的,当从回归方程中剔除变量时,回归平方和减少,残差平方和增加。反之,当往回归方程中引入变量时,回归平方和增加,残差平方和减少,两者的增减量同样相等。 六、关于拟合优度

R2?SSRSST,R2与回归方程中自变量的数目以及样本容量n有关,当样本容量n

与自变量个数接近时,R2易接近1,其中隐含着一些虚假成分。由R2决定模型优劣时还需慎重。 七、中心化和标准化

因为多元回归涉及的数据量很大,就可能由于舍入误差而使计算结果不理想。产生舍入误差有两个主要原因,一是回归分析计算中数据量级有很大差异,比如数据10000与0.1111这样的大小相差悬殊的数据出现在同一个计算中;二是设计矩阵X 的列向量近似线性相关时,X?X为病态矩阵,其逆矩阵(X?X)?1就会产生较大的误差。 1、中心化

多元线性回归模型的一般形式为:y??0??1x1??2x2????pxp??

????x???x?????x ???其经验回归方程为:y01122pp此经验方程进过样本中心(x1,x2,?,xp;y),将坐标原点移至样本中心,即作坐标

??xij?x变换:xijyi??yi?yi?1,2,?,n;j?1,2,?,p上述经验方程即

????x????x??????x?即为中心化经验回归方程。中心化经验????转变为:y01122pp?保持不变,回归方程的常数项为0,而回归系数的最小二乘估计值?因为坐标系j平移变化只改变直线的截距,不改变直线的斜率。

2、标准化回归系数

为了消除量纲不同和数量级的差异所带来的影响,就需要将样本数据作标准化处理,然后用最小二乘法估计未知参数,求得标准化系数。

??样本数据标准化公式:xijxij?xjLjjnijyi??yi?yLyy2i?1,2,?,n;j?1,2,?,p

n其中:Ljj??(xi?1?xj)2,Lyy??(yi?1i?y)

标准化回归系数与最小二乘回归系数之间存在关系式:???j?LjjLyy? ?j?表示在其他变量不变的情况下,自变量x的每单位的绝对普通最小二乘估计?jj??表示自变量x的1%相变化引起的因变量均值的绝对变化量。标准化回归系数?jj对变化(相对于Ljj)引起的因变量均值的相对变化百分数(相对于Lyy)。

标准化回归系数是比较自变量对y影响程度相对重要性的一种较为理想的方法,有了标准化回归系数后,变量的相对重要性就容易进行比较了。但是,仍要注意对回归系数的解释须采取谨慎的态度,这是因为当自变量相关时会影响标准化回归系数的大小。 八、相关阵与偏相关系数 1、样本相关阵

负相关系数R反映了y与一组自变量的相关性,是整体和共性指标,简单相关系数反映的是两个变量见的相关性,是局部和个性指标。在分析问题时,应该本着整体与局部相结合,共性与个性相结合的原则。求出y与每个自变量xi的相关系

?1?r?1yr??r2y数ryi,得到增广的样本相关阵为:~????r?pyry11r21?rp1ry2r121?rp2????ryp??r1p?r2p????1??

2、偏决定系数

在多元线性回归分析中,当其他变量被固定后,给定的任两个变量之间的相关系数,叫偏相关系数。偏相关系数可以度量p+1个变量y,x1,x2,?,xp之中任意两个变量的线性相关程度,而这种相关程度是在固定其余p-1个变量的影响下的线性相关。偏决定系数测量在回归方程中已包含若干个自变量时,再引入某一个新的自变量时,y的剩余变差的相对减少量,它衡量某个自变量对y的变差减少的边际贡献。

(1)两个自变量的偏决定系数

二元线性回归模型为:yi??0??1xi1??2xi2??i,i?1,2,?,n

记SSE(x2)是模型中只含有自变量x2时y的残差平方和,SSE(x1,x2)是模型中同时含有自变量x1和x2时y的残差平方和。模型中已含有x2时,再加入x1使y的剩余变差的相对减小量为:ry21,2?时,y与x1的偏决定系数。 (2)一般情况

在模型中含有x2,?,xp时,y与x1的偏决定系数为:

ry1;2,?,p?2SSE(x2)?SSE(x1,x2)SSE(x2)此时模型中已含有x2SSE(x2,?,xp)?SSE(x1,x2,?,xp)SSE(x2,?,xp),偏决定系数与回归系数显著性检验

的偏F值是等价的。 3、偏相关系数

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4