演習問題の答え     Last modified: Aug 19, 2015

 偏回帰係数が $0$ であるという帰無仮説の検定に用いる $t$ 値は,それぞれ $27.05056$,$26.53645$ となり( 自由度 $7$ の $t$ 分布に従う ),帰無仮説は棄却される。
 また,定数項が $0$ の検定の $t$ 値は $2.73685$ で,同じく帰無仮説は棄却される( 有意水準 $5\%$ )。


R を用いると以下のようになる。

> Y.hat <- X %*% b + b0 # 予測値
> Se <- sum((Y-Y.hat)^2) # 誤差変動
> St <- Syy # 従属変数の変動 
> Sr <- St-Se # 回帰で説明される変動
> MSe <- Se / (n-p-1) # 誤差の平均平方
> s <- solve(Sxx) # 逆行列
> SE.b <- sqrt(diag(s) * MSe) # 偏回帰係数の標準誤差
> t.b <- abs(b) / SE.b # H0:偏回帰係数=0 の検定統計量
> df <- n-p-1 # t 分布の自由度
> p.b <- pt(t.b, df, lower.tail=FALSE)*2 # P 値(両側検定)
> data.frame(b, SE=SE.b, t=t.b, p=p.b)
           b          SE        t            p
X1 0.2046172 0.007564251 27.05056 2.419227e-08
X2 0.2866338 0.010801511 26.53645 2.763814e-08

> ## 定数項について
> Mean.X <- colMeans(X) # 独立変数の平均値ベクトル
> SE.b0 <- sqrt((1/n + Mean.X %*% s %*% Mean.X)*MSe) # 定数項の標準誤差
> t.b0 <- abs(b0) / SE.b0 # H0:定数項=0 の検定統計量
> p.b0 <- pt(t.b0, df, lower.tail=FALSE)*2 # P 値(両側検定)
> c(constant=b0, SE=SE.b0, t=t.b0, p=p.b0)
  constant         SE          t          p 
0.14917563 0.05450634 2.73684913 0.02904999


・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI