回帰係数の確率分布2(線形単回帰)

線形単回帰によって得られる標本回帰式 y ^ = β ^ 0 + β ^ 1 x の係数 β ^ 0 β ^ 1 は,標本ごとに値が変わる. したがって, β ^ 0 β ^ 1 確率変数である.

誤差 ε 正規分布に従うと仮定すると, これらの推定量から真の係数を引き,それぞれの標準偏差をその推定量で置き換えて標準化した統計量

t 0 = β ^ 0 β 0 SE ( β ^ 0 )   ・・・・・・(1)

ただし, SE β ^ 0 β ^ 0 標準誤差(Standard Error)

t 1 = β ^ 1 β 1 SE ( β ^ 1 )   ・・・・・・(2)

ただし, SE β ^ 1 β ^ 1 標準誤差(Standard Error)

は,それぞれ自由度 n 2 t 分布に従う.

■(1),(2)の導出

回帰係数の確率分布(線形単回帰)より

β ^ 0 N β 0 , σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2   ・・・・・・(3)

E β ^ 0 = β 0   ・・・・・・(4)

V β ^ 0 = σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2   ・・・・・・(5)

β ^ 1 N β 1 , σ 2 i = 1 n x i x ¯ 2   ・・・・・・(6)

E β ^ 1 = β 1   ・・・・・・(7)

V β ^ 1 = σ 2 i = 1 n x i x ¯ 2   ・・・・・・(8)

誤差の分散の不偏推定量より

E S S E n 2 = σ 2   ・・・・・・(9)

となる.

(9)より σ 2 不偏推定量 S S E n 2 である.これより, β ^ 0 の分散 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 の不偏推定量は

S S E n 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 = S S E i = 1 n x i 2 n n 2 i = 1 n x i x ¯ 2   ・・・・・・(10)

となる.この不偏推定量の平方根は β ^ 0 標準誤差(Standard Error)と言われ SE β ^ 0 で表される.

SE β ^ 0 = S S E i = 1 n x i 2 n n 2 i = 1 n x i x ¯ 2   ・・・・・・(11)

(4)と(11)より

t 0 = β ^ 0 β 0 SE ( β ^ 0 )

で表される.

同様に, β ^ 1 の分散 σ 2 i = 1 n x i x ¯ 2 不偏推定量

S S E n 2 i = 1 n x i x ¯ 2 = S S E n 2 i = 1 n x i x ¯ 2   ・・・・・・(12)

となる.この不偏推定量の平方根は β ^ 1 標準誤差(Standard Error)と言われ SE β ^ 1 で表される.

SE β ^ 1 = S S E n 2 i = 1 n x i x ¯ 2   ・・・・・・(13)

(7)と(13)より

t 1 = β ^ 1 β 1 SE ( β ^ 1 )

で表される.

t 0 の分布について

t 0 = β ^ 0 β 0 SE β ^ 0 を以下のように式変形する.

t 0 = β ^ 0 β 0 SE β ^ 0 = β ^ 0 β 0 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 SE β ^ 0   ・・・・・・(14)

β ^ 0 β 0 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 は(3)の正規分布標準化したものであるので, N 0 , 1 に従う.  ・・・・・・(A)

(14)に(11)を代入する.

= β ^ 0 β 0 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 S S E i = 1 n x i 2 n n 2 i = 1 n x i x ¯ 2

= β ^ 0 β 0 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 1 S S E σ 2 n 2

S S E σ 2 は自由度 n 2 χ 2 分布に従うことが分かっている(ここを参照).  ・・・・・・(B)

β ^ 0 β 0 σ 2 i = 1 n x i 2 n i = 1 n x i x ¯ 2 = Z 0 S S E σ 2 = W とおくと

= Z 0 W n 2

Z 0 は(A)より N 0 , 1 に従う.

W は(B)より自由度 n 2 χ 2 分布に従う.

そして,標本回帰係数と母回帰係数との差と残差が互いに独立より, Z 0 W は互いに独立である.

以上より, t 0 は自由度 n 2 t 分布に従う.

t 1 の分布について

t 1 = β ^ 1 β 1 SE β ^ 1 を以下のように式変形する.

t 1 = β ^ 1 β 1 SE β ^ 1 = β ^ 1 β 1 σ 2 i = 1 n x i x ¯ 2 σ 2 i = 1 n x i x ¯ 2 SE β ^ 1   ・・・・・・(15)

β ^ 1 β 1 σ 2 i = 1 n x i x ¯ 2 は(6)の正規分布標準化したものであるので, N 0 , 1 に従う.  ・・・・・・(C)

(15)に(13)を代入する.

= β ^ 1 β 1 σ 2 i = 1 n x i x ¯ 2 σ 2 i = 1 n x i x ¯ 2 S S E n 2 i = 1 n x i x ¯ 2

= β ^ 1 β 1 σ 2 i = 1 n x i x ¯ 2 1 S S E σ 2 n 2

S S E σ 2 は自由度 n 2 χ 2 分布に従うことが分かっている(ここを参照).  ・・・・・・(D)

β ^ 1 β 1 σ 2 i = 1 n x i x ¯ 2 = Z 1 S S E σ 2 = W とおくと

= Z 1 W n 2

Z 1 は(C)より N 0 , 1 に従う.

W は(D)より自由度 n 2 χ 2 分布に従う.

そして,標本回帰係数と母回帰係数との差と残差が互いに独立より, Z 1 W は互いに独立である.

以上より, t 1 は自由度 n 2 t 分布に従う.

 

ホーム>>カテゴリー分類>>確率>>統計>>回帰係数の確率分布2(線形単回帰)

 最終更新日: 2026年7月30日