MathLabs
定理証明済み

通常最小二乗推定量

内容

すべての直線 y=b0+b1xy=b_0+b_1x の中で、残差平方和 ∑i=1n(yi−b0−b1xi)2\sum_{i=1}^n (y_i-b_0-b_1x_i)^2 を最小にする選択は b1=β^1=SxySxxb_1=\hat\beta_1=\dfrac{S_{xy}}{S_{xx}} および b0=β^0=yˉ−β^1xˉb_0=\hat\beta_0=\bar y-\hat\beta_1\bar x である。ここで Sxy=∑i=1n(xi−xˉ)(yi−yˉ)S_{xy}=\sum_{i=1}^n(x_i-\bar x)(y_i-\bar y)、Sxx=∑i=1n(xi−xˉ)2S_{xx}=\sum_{i=1}^n(x_i-\bar x)^2。

なぜ正しいのか?

残差平方和は b0b_0 と b1b_1 に関して滑らかな(二次で凸な)関数であるから、その最小値はちょうど両方の偏微分が0になる点で見つかる——お椀の底を、あらゆる方向の傾きを0にすることで見つけるのと同じ発想である。

証明の概略

Q(b0,b1)=∑i=1n(yi−b0−b1xi)2Q(b_0,b_1)=\sum_{i=1}^n (y_i-b_0-b_1x_i)^2 とおく。b0b_0 に関する偏微分を0とすると:∂Q∂b0=−2∑i=1n(yi−b0−b1xi)=0\dfrac{\partial Q}{\partial b_0}=-2\sum_{i=1}^n(y_i-b_0-b_1x_i)=0、これは第一正規方程式 ∑i=1nyi=nb0+b1∑i=1nxi\sum_{i=1}^n y_i = nb_0+b_1\sum_{i=1}^n x_i、すなわち b0=yˉ−b1xˉb_0=\bar y-b_1\bar x に簡約される。

b1b_1 に関する偏微分を0とすると:∂Q∂b1=−2∑i=1nxi(yi−b0−b1xi)=0\dfrac{\partial Q}{\partial b_1}=-2\sum_{i=1}^n x_i(y_i-b_0-b_1x_i)=0、第二正規方程式 ∑i=1nxiyi=b0∑i=1nxi+b1∑i=1nxi2\sum_{i=1}^n x_iy_i = b_0\sum_{i=1}^n x_i+b_1\sum_{i=1}^n x_i^2。

第一式の b0=yˉ−b1xˉb_0=\bar y-b_1\bar x を第二式に代入すると:∑i=1nxiyi=(yˉ−b1xˉ)nxˉ+b1∑i=1nxi2=nxˉyˉ+b1(∑i=1nxi2−nxˉ2)\sum_{i=1}^n x_iy_i = (\bar y-b_1\bar x)n\bar x+b_1\sum_{i=1}^n x_i^2 = n\bar x\bar y+b_1\left(\sum_{i=1}^n x_i^2-n\bar x^2\right)。整理して b1b_1 を分離する:b1(∑i=1nxi2−nxˉ2)=∑i=1nxiyi−nxˉyˉb_1\left(\sum_{i=1}^n x_i^2-n\bar x^2\right)=\sum_{i=1}^n x_iy_i-n\bar x\bar y。

直接の代数展開により ∑i=1nxi2−nxˉ2=∑i=1n(xi−xˉ)2=Sxx\sum_{i=1}^n x_i^2-n\bar x^2=\sum_{i=1}^n(x_i-\bar x)^2=S_{xx}、∑i=1nxiyi−nxˉyˉ=∑i=1n(xi−xˉ)(yi−yˉ)=Sxy\sum_{i=1}^n x_iy_i-n\bar x\bar y=\sum_{i=1}^n(x_i-\bar x)(y_i-\bar y)=S_{xy} となるので、b1=Sxy/Sxxb_1=S_{xy}/S_{xx}。これを戻すと b0=yˉ−b1xˉb_0=\bar y-b_1\bar x が得られる。QQ は ∣b0∣,∣b1∣→∞|b_0|,|b_1|\to\infty で際限なく増加する二次形式の平方和であるから、この唯一の停留点が大域的最小値である。

この定理を使うトピック

ステップごとの証明

この定理のステップごとの証明はまだありません。