跳到主要内容

一橋大学 経済学研究科 2022年度 統計学・計量経済学 第2題 2

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

以下の (a) から (e) のすべてに答えよ。導出過程は省略しないこと。

{(Yi,X1i,X2i)}i=1n\{(Y_i,X_{1i},X_{2i})\}_{i=1}^n を独立同一分布に従う無作為標本とし、

Yi=β1X1i+β2X2i+ui,Y_i=\beta_1X_{1i}+\beta_2X_{2i}+u_i,
E[uiX1i,X2i]=0,E[ui2X1i,X2i]=σ2E[u_i\mid X_{1i},X_{2i}]=0,\qquad E[u_i^2\mid X_{1i},X_{2i}]=\sigma^2

とする。ただし β20\beta_2\ne0 とする。(b)から(e)では nn\to\infty とし、 大数の法則と中心極限定理に必要な条件を仮定する。

(a)

YiY_iX1i,X2iX_{1i},X_{2i} に回帰して得られる最小二乗推定量 (β^1,β^2)(\hat\beta_1,\hat\beta_2) を導出せよ。

(b)

n(β^1β1)\sqrt n(\hat\beta_1-\beta_1) の漸近分布を導出せよ。

(c)

H0:β1=0H_0:\beta_1=0H1:β10H_1:\beta_1\ne0 を漸近的に有意水準 5%5\% で 検定する手続きを説明せよ。

(d)

E[X1iX2i]=γ0E[X_{1i}X_{2i}]=\gamma\ne0 のとき、X2iX_{2i} を除いて回帰した β1\beta_1 の最小二乗推定量が一致推定量でないことを示せ。

(e)

E[X1iX2i]=0E[X_{1i}X_{2i}]=0 のとき、同じ推定量が一致推定量であることを示せ。

Kai

(a)

Sjk=i=1nXjiXki,SjY=i=1nXjiYi,S_{jk}=\sum_{i=1}^nX_{ji}X_{ki},\qquad S_{jY}=\sum_{i=1}^nX_{ji}Y_i,
D=S11S22S122D=S_{11}S_{22}-S_{12}^2

とおく。正規方程式を解けば

(S11S12S12S22)(β^1β^2)=(S1YS2Y)\begin{pmatrix}S_{11}&S_{12}\\S_{12}&S_{22}\end{pmatrix} \begin{pmatrix}\hat\beta_1\\\hat\beta_2\end{pmatrix} =\begin{pmatrix}S_{1Y}\\S_{2Y}\end{pmatrix}

より

β^1=S22S1YS12S2YD,β^2=S11S2YS12S1YD.\boxed{\hat\beta_1=\frac{S_{22}S_{1Y}-S_{12}S_{2Y}}D},\qquad \boxed{\hat\beta_2=\frac{S_{11}S_{2Y}-S_{12}S_{1Y}}D}.

(b)

xi=(X1i,X2i)T\boldsymbol x_i=(X_{1i},X_{2i})^\mathsf T

Q=E[xixiT]=(accb)Q=E[\boldsymbol x_i\boldsymbol x_i^\mathsf T] =\begin{pmatrix}a&c\\c&b\end{pmatrix}

とおく。大数の法則と中心極限定理より

1nixixiTpQ,1nixiuidN(0,σ2Q).\frac1n\sum_i\boldsymbol x_i\boldsymbol x_i^\mathsf T\xrightarrow{p}Q, \qquad \frac1{\sqrt n}\sum_i\boldsymbol x_iu_i \xrightarrow{d}N(\boldsymbol0,\sigma^2Q).

したがって

n(β^β)=(1nixixiT)11nixiuidN(0,σ2Q1).\sqrt n(\hat{\boldsymbol\beta}-\boldsymbol\beta) =\left(\frac1n\sum_i\boldsymbol x_i\boldsymbol x_i^\mathsf T\right)^{-1} \frac1{\sqrt n}\sum_i\boldsymbol x_iu_i \xrightarrow{d}N(\boldsymbol0,\sigma^2Q^{-1}).

したがって

n(β^1β1)dN(0,σ2babc2).\boxed{ \sqrt n(\hat\beta_1-\beta_1) \xrightarrow{d}N\left(0,\frac{\sigma^2b}{ab-c^2}\right)}.

(c)

u^i=Yiβ^1X1iβ^2X2i\hat u_i=Y_i-\hat\beta_1X_{1i}-\hat\beta_2X_{2i} とし、

σ^2=1n2iu^i2,T=β^1σ^2S22/D\hat\sigma^2=\frac1{n-2}\sum_i\hat u_i^2, \qquad T=\frac{\hat\beta_1}{\sqrt{\hat\sigma^2S_{22}/D}}

とおく。H0H_0 の下で TdN(0,1)T\xrightarrow{d}N(0,1) なので、

T>1.96\boxed{|T|>1.96}

ならば H0H_0 を棄却する。

(d)

X2iX_{2i} を除いた推定量を β~1\widetilde\beta_1 とすると

β~1=iX1iYiiX1i2=β1+β2iX1iX2iiX1i2+iX1iuiiX1i2.\widetilde\beta_1 =\frac{\sum_iX_{1i}Y_i}{\sum_iX_{1i}^2} =\beta_1+\beta_2\frac{\sum_iX_{1i}X_{2i}}{\sum_iX_{1i}^2} +\frac{\sum_iX_{1i}u_i}{\sum_iX_{1i}^2}.

E[X1iui]=E[X1iE[uiX1i,X2i]]=0E[X_{1i}u_i]=E[X_{1i}E[u_i\mid X_{1i},X_{2i}]]=0 より

plimβ~1=β1+β2γE[X1i2]β1.\operatorname{plim}\widetilde\beta_1 =\beta_1+\beta_2\frac{\gamma}{E[X_{1i}^2]}\ne\beta_1.

よって一致推定量でない。

(e)

E[X1iX2i]=0E[X_{1i}X_{2i}]=0 ならば同様に

plimβ~1=β1.\operatorname{plim}\widetilde\beta_1=\beta_1.

したがって β~1pβ1\boxed{\widetilde\beta_1\xrightarrow{p}\beta_1} である。

Reference