跳到主要内容

早稲田大学 創造理工学研究科 経営システム工学専攻 2024年8月実施 概率统计

Author

思齐塾, 祭音Myyura

Description

2次元のデータ (xi,yi)(x_i, y_i) (i=1,2,,n)(i = 1, 2, \dots, n) に対して、次の回帰モデルを仮定する。

yi=βxi2+ϵi,ϵiN(0,σ2)y_i = \beta x_i^2 + \epsilon_i, \quad \epsilon_i \sim N(0, \sigma^2)

このとき, β\beta の最小二乗推定量を導け。

题目描述

对二维数据 (xi,yi)(x_i,y_i)i=1,2,,ni=1,2,\ldots,n),假设回归模型

yi=βxi2+ϵi,ϵiN(0,σ2).y_i=\beta x_i^2+\epsilon_i, \qquad \epsilon_i\sim N(0,\sigma^2).

推导 β\beta 的最小二乘估计量。

Kai

最小二乗推定量(Least Squares Estimator, LSE)を求める。 残差二乗和(Residual Sum of Squares, RSS)を最小化する β\beta を求める。

RSS(β)=i=1n(yiβxi2)2RSS(\beta) = \sum_{i=1}^n (y_i - \beta x_i^2)^2

β\beta で微分して0とおく。

dRSSdβ=i=1n2(yiβxi2)(xi2)=0\frac{dRSS}{d\beta} = \sum_{i=1}^n 2(y_i - \beta x_i^2)(-x_i^2) = 0
i=1n(yiβxi2)xi2=0\sum_{i=1}^n (y_i - \beta x_i^2)x_i^2 = 0
i=1nyixi2βi=1nxi4=0\sum_{i=1}^n y_i x_i^2 - \beta \sum_{i=1}^n x_i^4 = 0
βi=1nxi4=i=1nyixi2\beta \sum_{i=1}^n x_i^4 = \sum_{i=1}^n y_i x_i^2
β^=i=1nyixi2i=1nxi4\hat{\beta} = \frac{\sum_{i=1}^n y_i x_i^2}{\sum_{i=1}^n x_i^4}

ただし、この一意な推定量が存在するためには

i=1nxi4>0\sum_{i=1}^n x_i^4>0

すなわち、少なくとも1つの xix_i が0でないことが必要である。すべての xi=0x_i=0 の場合、 RSS(β)=iyi2RSS(\beta)=\sum_i y_i^2β\beta に依存しないため、任意の β\beta が最小化し、 β\beta は識別できない。

よって, ixi4>0\sum_i x_i^4>0 の条件下で β\beta の最小二乗推定量は

β^=i=1nyixi2i=1nxi4\hat{\beta} = \frac{\sum_{i=1}^n y_i x_i^2}{\sum_{i=1}^n x_i^4}