跳到主要内容

一橋大学 経済学研究科 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[ui∣X1i,X2i]=0,E[ui2∣X1i,X2i]=σ2E[u_i\mid X_{1i},X_{2i}]=0,\qquad E[u_i^2\mid X_{1i},X_{2i}]=\sigma^2

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

(a)​

YiY_i を X1i,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=0 対 H1:β1≠0H_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) の標本行列が列フルランク(D>0D>0)、(b) 以降の母二次モーメント行列 QQ が正定値であると仮定する。(c) ではさらに σ2>0\sigma^2>0 とする。これらは係数の識別と検定に必要である。

(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=S11S22−S122D=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=S22S1Y−S12S2YD,β^2=S11S2Y−S12S1YD.\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}

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

1n∑ixixiT→pQ,1n∑ixiui→dN(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(β^−β)=(1n∑ixixiT)−11n∑ixiui→dN(0,σ2Q−1).\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,σ2bab−c2).\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=1n−2∑iu^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 の下で T→dN(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=∑iX1iYi∑iX1i2=β1+β2∑iX1iX2i∑iX1i2+∑iX1iui∑iX1i2.\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[ui∣X1i,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.

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

Reference​