経済統計の使い方
統計初心者向けに、データや分析手法の解説をしています。「まとめページ」をご覧くだされば、全体的な内容がわかると思います。
統計学・計量経済学

【計量経済学】回帰係数βの期待値、分散について|丁寧に式を展開

単回帰分析の回帰係数βは誤差の二乗を最小化することで得られます。

その導出も大事ですが、この記事ではその先の期待値や分散の導出方法について解説します。

まず、最初に結論としての式を書きます。その後にできるだけわかりやすく、解説をつけていきます。

経済統計の使い方では、統計データの入手法から分析法まで解説しています。

統計学・計量経済学のまとめ 統計学に関するまとめのサイトです。 教科書 「回帰分析から学ぶ計量経済学-Excelで読み解く経済の仕組み」 回帰分...

回帰係数β

まず結果を見ます。

ˆ𝛽 =𝑛𝑖=1(𝑋𝑖¯𝑋)(𝑌𝑖¯𝑌)𝑛𝑖=1(𝑋𝑖¯𝑋)2 =𝑛𝑖=1(𝑋𝑖¯𝑋)𝑌𝑖𝑛𝑖=1(𝑋𝑖¯𝑋)2

βの確率的表現

まず結果を見ます。係数を誤差項の関数として表したものです。

ˆ𝛽 =𝛽 +𝑛𝑛=1(𝑋𝑖¯𝑋)𝑢𝑖𝑛𝑛=1(𝑋𝑖¯𝑋)2

βの期待値

𝐸(ˆ𝛽) =𝛽 +𝑛𝑛1(𝑋𝑖¯𝑋)𝐸(𝑢𝑖)𝑛𝑛1(𝑋𝑖¯𝑋)2

最小二乗法の仮定により、「𝑢𝑖の期待値はゼロ」なので、以下の結果になります。

𝐸(ˆ𝛽) =𝛽

βの分散

まず結果を見ます。𝜎は誤差項の分散です。

𝑉(ˆ𝛽) =𝜎2𝑛𝑛1(𝑋𝑖¯𝑋)2

回帰係数βの変形

ˆ𝛽 =𝑛𝑖=1(𝑋𝑖¯𝑋)(𝑌𝑖¯𝑌)𝑛𝑖=1(𝑋𝑖¯𝑋)2

=𝑛𝑖=1(𝑋𝑖¯𝑋)𝑌𝑖¯𝑌𝑛𝑖=1(𝑋𝑖¯𝑋)𝑛𝑖=1(𝑋𝑖¯𝑋)2

𝑋𝑖を1からnまで足したものと、平均¯𝑋をn回足したものは同じなので、𝑛𝑖=1(𝑋𝑖 ¯𝑋) =0

=𝑛𝑖=1(𝑋𝑖¯𝑋)𝑌𝑖𝑛𝑖=1(𝑋𝑖¯𝑋)2

βの確率的表現

𝑌𝑖 =𝛼 +𝛽𝑋𝑖 +𝑢𝑖

この式は平均についても成り立つので、以下の式が成り立ちます。

 ¯𝑌 =𝛼 +𝛽¯𝑋 +¯𝑢

基本の式から平均の式を引くと、以下の式になります。

𝑌𝑖 ¯𝑌𝑖 =𝛽(𝑋𝑖 ¯𝑋) +𝑢𝑖 ¯𝑢

ˆ𝛽は以下のように表せます。

ˆ𝛽 =𝑛𝑖=1(𝑋𝑖¯𝑋)(𝛽(𝑋𝑖¯𝑋)+(𝑢𝑖¯𝑢))𝑛𝑖=1(𝑋𝑖¯𝑋)2

=𝛽𝑛𝑖=1(𝑋𝑖¯𝑋)2𝑛𝑖=1(𝑋𝑖¯𝑋)2 +𝑛𝑖=1(𝑋𝑖¯𝑋)(𝑢𝑖¯𝑢)𝑛𝑖=1(𝑋𝑖¯𝑋)2

=𝛽 +𝑛𝑖=1(𝑋𝑖¯𝑋)𝑢𝑖𝑛𝑖=1(𝑋𝑖¯𝑋)2 𝑛𝑖=1(𝑋𝑖¯𝑋)¯𝑢𝑛𝑖=1(𝑋𝑖¯𝑋)2

𝑛𝑖=1(𝑋𝑖 ¯𝑋)¯𝑢 =(𝑋𝑖 ¯𝑋)𝑛𝑖=1¯𝑢で、𝑛𝑖=1¯𝑢は「誤差の平均はゼロ」なので、全体としてもゼロです。

βの期待値の変形

𝐸(ˆ𝛽) =𝛽 +(𝑋𝑖¯𝑋)𝐸(𝑢𝑖)(𝑋𝑖¯𝑋)2

=𝛽 +𝐸(𝑢𝑖)(𝑋𝑖¯𝑋)(𝑋𝑖¯𝑋)2

最小二乗法の仮定により「𝑢𝑖の期待値はゼロ」なので、𝐸(𝑢𝑖) =0になります。

=𝛽

期待値がβになるために使った仮定は「𝑋𝑖ˆ𝑋は非確率変数」「𝑢𝑖の期待値はゼロ」なので、不均一分散がある場合でも、この式は成り立ちます。

βの分散の変形

分散は、変数の平均からの偏差の二乗の平均なので、𝑉(𝑋) =𝐸((𝑋𝐸(𝑋))2)となります。ˆ𝛽𝛽の不偏推定量なので、ˆ𝛽の平均(期待値)は𝛽になります。これらを勘案すると以下の変形ができます。

𝑉(ˆ𝛽) =𝐸(ˆ𝛽𝐸(ˆ𝛽))2 =𝐸(ˆ𝛽𝛽)2

ˆ𝛽𝛽の期待値の式を代入します。

=𝐸⎢ ⎢ ⎢ ⎢⎜ ⎜ ⎜ ⎜ ⎜ ⎜𝛽+𝑛𝑖=1(𝑋𝑖¯𝑋)𝑢𝑖𝑛𝑖=1(𝑋𝑖¯𝑋)2𝛽⎟ ⎟ ⎟ ⎟ ⎟ ⎟2⎥ ⎥ ⎥ ⎥

=𝐸⎢ ⎢ ⎢ ⎢⎜ ⎜ ⎜ ⎜ ⎜ ⎜𝑛𝑖=1(𝑋𝑖¯𝑋)𝑢𝑖𝑛𝑖=1(𝑋𝑖¯𝑋)2⎟ ⎟ ⎟ ⎟ ⎟ ⎟2⎥ ⎥ ⎥ ⎥

=𝐸[(𝑛𝑖=1(𝑋𝑖¯𝑋)𝑢𝑖)2]𝐸[(𝑛𝑖=1(𝑋𝑖¯𝑋)2)2]

 (𝑛𝑖=1(𝑋𝑖¯𝑋)𝑢𝑖)2 に対して、(𝑛𝑖=1𝑋𝑖)2 =𝑛𝑖=1𝑛𝑗=1𝑋𝑖𝑋𝑗 を使います。

1((𝑋𝑖ˆ𝑋)2)2𝐸[𝑛𝑖=1𝑛𝑗=1(𝑋𝑖 ˆ𝑋)𝑢𝑖(𝑋𝑗 ˆ𝑋)𝑢𝑗]

=1((𝑋𝑖ˆ𝑋)2)2𝐸[𝑛𝑖=1𝑛𝑗=1𝑢𝑖𝑢𝑗(𝑋𝑖 ˆ𝑋)(𝑋𝑗 ˆ𝑋)]

=1((𝑋𝑖¯𝑋)2)2(𝑋𝑖¯𝑋)2𝜎2

=1((𝑋𝑖¯𝑋)2)2𝜎2(𝑋𝑖¯𝑋)2

=𝜎2(𝑋𝑖¯𝑋)2

分子の変形について

上記の変形について、詳しく見ていきます。分子の部分は以下の式に期待値をとったものです。2つの部分に分けて期待値をとります。

𝑛𝑖=1𝑛𝑗=1𝑢𝑖𝑢𝑗(𝑋𝑖 ˆ𝑋)(𝑋𝑗 ˆ𝑋)

=𝑖=𝑗(𝑋𝑖¯𝑋)2𝑢2𝑖 +𝑖𝑗(𝑋𝑖 ¯𝑋)(𝑋𝑗 ¯𝑋)𝐸(𝑢𝑖)𝐸(𝑢𝑗)

第1項の期待値

第1項は、iとjが等しい場合で、1からnまでn個あります。

第1項について期待値をとります。「𝑋𝑖¯𝑋は非確率変数」という最小二乗法の仮定から、前半の期待値オぺーレーターを取ることができます。誤差項の二乗の期待値は、誤差項の分散に等しく、𝐸(𝑢2𝑖) =𝜎2です。

𝐸(𝑛𝑖=1(𝑋𝑖¯𝑋)2𝑢2𝑖)

=𝑛𝑖=1(𝑋𝑖¯𝑋)2𝐸(𝑢2𝑖)

=𝑛𝑖=1(𝑋𝑖¯𝑋)2𝜎2

第2項の期待値

次に、第2項について期待値をとります。「𝑋𝑖¯𝑋は非確率変数」という最小二乗法の仮定から、前半の期待値オぺーレーターを取ることができます。

𝐸[𝑖𝑗(𝑋𝑖¯𝑋)(𝑋𝑗¯𝑋)𝑢𝑖𝑢𝑗]

=𝑖𝑗(𝑋𝑖 ¯𝑋)(𝑋𝑗 ¯𝑋)𝐸(𝑢𝑖𝑢𝑗)

異なる誤差どうしは無相関」という最小二乗法の仮定から、𝐸(𝑢𝑖𝑢𝑗) =0になります。

𝑖𝑗(𝑋𝑖 ¯𝑋)(𝑋𝑗 ¯𝑋)𝐸(𝑢𝑖𝑢𝑗) =0

再掲

まとめると以下の式です。

𝐸(𝑛𝑖=1𝑛𝑗=1𝑢𝑖𝑢𝑗(𝑋𝑖 ˆ𝑋)(𝑋𝑗 ˆ𝑋))

=𝑛𝑖=1(𝑋𝑖¯𝑋)2𝜎2

これを元の式から順に書くと以下の通りになります。

1((𝑋𝑖ˆ𝑋)2)2𝐸[𝑛𝑖=1𝑛𝑗=1𝑢𝑖𝑢𝑗(𝑋𝑖 ˆ𝑋)(𝑋𝑗 ˆ𝑋)]

=1((𝑋𝑖¯𝑋)2)2(𝑋𝑖¯𝑋)2𝜎2

=1((𝑋𝑖¯𝑋)2)2𝜎2(𝑋𝑖¯𝑋)2

=𝜎2(𝑋𝑖¯𝑋)2

COMMENT

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA