logo

回帰係数のF検定 📂統計的検定

回帰係数のF検定

仮説検定

線形多重回帰モデルのモデル診断で、残差が線形性、等分散性、独立性、正規性を満たすとする。 nn 個の観測値と pp 個の独立変数を対象とした多重回帰分析における仮説検定は次の通りです。

  • H0H_{0}: β1=β2=⋯=βp=0\beta_{1} = \beta_{2} = \cdots = \beta_{p} = 0 つまり、すべての独立変数が従属変数と相関関係を持たない。
  • H1H_{1}: β1,β2,⋯ ,βp\beta_{1} , \beta_{2} , \cdots , \beta_{p} のうち少なくとも1つは 0 0 ではない。つまり、少なくとも1つの独立変数が意味のある相関関係を持つ。

導出

SST, SSR, SSE: 6. 総平方和(Total Sum of Squares)またはSST(Sum of Squares Total): TSS=SST:=∑i=1n(yi−y‾)2 \text{TSS} =\text{SST} := \sum_{i=1}^{n} ( y_{i} - \overline{y} )^2 7. 回帰平方和(Explained Sum of Squares)またはSSR(Sum of Squares due to Regression): ESS=SSR:=∑i=1n(y^i−y‾)2 \text{ESS} = \text{SSR} := \sum_{i=1}^{n} ( \hat{y}_{i} - \overline{y} )^2 8. 残差平方和(Residual Sum of Squares)またはSSE(Sum of squared Error): RSS=SSE:=∑i=1n(yi−y^i)2 \text{RSS} = \text{SSE} := \sum_{i=1}^{n} ( y_{i} - \hat{y}_{i} )^2

SSTは平均1つを用いるため、自由度が(n−1)(n-1)であり、SSEはpp個の独立変数に基づいて計算されるため、定数項を含む(p+1)(p+1)個の回帰係数を使用し、自由度が(n−(p−1))\left( n-(p-1) \right)である。これらは残差の等分散性と独立性、正規性に従ってσ2\sigma^{2}で分けたとき、それぞれの自由度を持つカイ二乗分布に従う。一方 SST=SSR+SSE  ⟺  SSR=SST−SSE SST = SSR + SSE \iff SSR = SST - SSE であるため、SSRSSRの自由度を持ち、カイ二乗分布χ2(p)\chi^{2} (p)に従う。

F分布の導出: ２つの確率変数 U,VU,Vが独立であり、U∼χ2(r1)U \sim \chi^{2} ( r_{1})、V∼χ2(r2)V \sim \chi^{2} ( r_{2})であるとすると U/r1V/r2∼F(r1,r2) {{ U / r_{1} } \over { V / r_{2} }} \sim F \left( r_{1} , r_{2} \right)

検定統計量 FFを F:=SSR/pSSE/(n−p−1) F := {{ \text{SSR} / p } \over { \text{SSE} / (n-p-1 ) }} として定義すると自由度が(p,n−p−1)(p , n-p-1)のF分布に従う。

方程式で再度表すとF=SSR/pSSE/(n−p−1)∼F(p,n−p−1)\displaystyle F = {{ \text{SSR} / p } \over { \text{SSE} / (n-p-1 ) }} \sim F(p, n-p-1)であり、これを利用して仮説検定を行う。

■

このF検定は回帰係数のt検定が見られるのであれば無意味であり、真の本質はモデルの比較にある。統計学に避けられない「主観性」や「あいまいさ」を払拭して、統計的に意味のある結果を出せるようになることにある。（もちろん、実際の分析ではもっと便利で簡単な統計量を使ってモデルを比較する。）

縮小モデルに対する仮説検定

nn 個の観測値と pp 個の独立変数を対象とした多重回帰分析について、i=0,1,⋯ ,pi=0,1,\cdots,pとする。この回帰モデルを全体モデルfMと呼び、FMからkk個の独立変数を除去したモデルを縮小モデルrMと呼ぶ。

  • H0H_{0}: RMで十分である。つまり、多くの変数を使ってFMを使用する必要はない。
  • H1H_{1}: RMは不十分である。つまり、変数を増やしてでもFMを使用する方が良い。

F=[SSR(RM)−SSR(FM)]/(p+1−k)SSE(FM)/(n−p−1) F = {{ [ \text{SSR} (RM) - \text{SSR} (FM) ] / (p +1 - k) } \over { \text{SSE} ( FM ) / (n-p-1 ) }} は、自由度が(p+1−k,n−p−1)(p + 1 - k , n-p-1)のF分布に従う。有意水準 α\alphaでF≤F(p+1−k,n−p−1;α)F \le F_{ ( p+1-k , n-p-1 ; \alpha ) }であれば、H0H_{0}が採用され、縮小RMが使用できるようになる。

関連項目