§3.2多元线性回归模型的估计 佔计目标:结构参数β及随机误差项的方差G2 估计方法:OLS、ML或者MM 、普通最小二乘估计 二、最大或然估计 三、矩估计 四、参数估计量的性质 五、样本容量问题 六、估计实例
§3.2 多元线性回归模型的估计 估计方法:OLS、ML或者MM 一、普通最小二乘估计 *二、最大或然估计 *三、矩估计 四、参数估计量的性质 五、样本容量问题 六、估计实例
一、普通最小二乘估计 对于随机抽取的n组观测值(x,X1)1=12,…n,1=012, 如果样本函数的参数估计值已经得到,则有: =B+BX1+B2X2+…+B1X8 根据最小二乘原理,参数估计值应该是下列方程组的解 Q=0 oB Q=0 其中g=∑e2=∑(y-) Q=0 oB Y1-(B+B1X1+B2X2 BLXN) oB
一、普通最小二乘估计 对于随机抽取的n组观测值 Y X i n j k ( i , ji), =1,2, , , = 0,1,2, 如果样本函数的参数估计值已经得到,则有: Yi X i X i ki X Ki ˆ ˆ ˆ ˆ ˆ = 0 + 1 1 + 2 2 ++ i=1,2…n 根据最小二乘原理,参数估计值应该是下列方程组的解 = = = = 0 ˆ 0 ˆ 0 ˆ 0 ˆ 2 1 0 Q Q Q Q k 其中 2 1 1 2 ) ˆ ( = = = = − n i i i n i Q ei Y Y 2 1 0 1 1 2 2 )) ˆ ˆ ˆ ˆ ( ( = = − + + + + n i Yi X i X i k X ki
于是得到关于待估参数估计值的正规方程组: ∑(B6+B1X1+B2X21+…+B4X)=Y Σ(B+B1X1+B2X2+…+BX)X1=∑FX1 Σ(B6+B1X1+B2X2+…+B4X)X21=∑yX2 Σ(B6+B1X1+B2X21+…+BkX)X如=∑yXb 解该(k+1)个方程组成的线性代数方程组,即可得到 (k+1)个待估参数的估计值月,j=012,…k
于是得到关于待估参数估计值的正规方程组: + + + + = + + + + = + + + + = + + + + = i i k ki ki i ki i i i k ki i i i i i k ki i i i i i k ki i X X X X Y X X X X X Y X X X X X Y X X X X Y ) ˆ ˆ ˆ ˆ ( ) ˆ ˆ ˆ ˆ ( ) ˆ ˆ ˆ ˆ ( ) ˆ ˆ ˆ ˆ ( 0 1 1 2 2 0 1 1 2 2 2 2 0 1 1 2 2 1 1 0 1 1 2 2 解该(k+1)个方程组成的线性代数方程组,即可得到 (k+1)个待估参数的估计值 , , , , , j j = 0 1 2 k
正规方程组的矩阵形式 X li ∑XB ∑ X,XB,X ∑X∑XX1…∑X k2 (XX)B=XY 由于ⅹX满秩,故有 B=(XX)X'Y
正规方程组的矩阵形式 = k k kn n n ki ki i ki k i i i ki i ki Y Y Y X X X X X X X X X X X X X X n X X 2 1 1 2 1 1 1 1 2 1 0 2 1 1 2 1 1 1 1 1 1 ˆ ˆ ˆ 即 (XX)β ˆ = XY 由于X’X满秩,故有 β= XX XY −1 ( ) ˆ
将上述过程用矩阵表示如下 寻找一组参数佔计值β,使得残差平方和 0=e=ee=(Y-X B)'(Y-X B) 最小。 即求解方程组:(Y-XB)(Y-XB)=0 (Y'Y-BXY-YXB+BXXB)=0 -( YY-2YXB+BXXB=0 XY+XXB=0 得到 XY=XXB 于是:B=(XX)XY
将上述过程用矩阵表示如下: 即求解方程组: ( ˆ ) ( ˆ ) 0 ˆ Y − Xβ Y − Xβ = β ( ˆ ˆ ˆ ˆ ) 0 ˆ − − + = Y Y βX Y Y Xβ βX Xβ β ( 2 ˆ ˆ ˆ ) 0 ˆ − + = Y Y Y Xβ βX Xβ β − XY + XXβ ˆ = 0 得到: β= XX XY −1 ( ) ˆ X Y X Xβ ˆ = 于是: