Rを使用して臨界t値を計算するにはどうすればよいですか？

これが新しい質問の場合は申し訳ありません。初めて統計を教えようとしています。基本的な手順は下がっていると思いますが、Rで実行するのに苦労しています。

だから、私は形の多重線形回帰における回帰係数の有意性を評価しようとしています

\hat{y} = X \hat{β}

$\hat y = X \hat \beta$

をテストするためのt統計は $H_0: \hat \beta_j = 0, H_a: \hat \beta_j \neq 0$

t_{0} = \frac{{\hat{β}}_{j} - 0}{se ({\hat{β}}_{j})} = \frac{{\hat{β}}_{j}}{\sqrt{{\hat{σ}}^{2} C_{j j}}} = \frac{{\hat{β}}_{j}}{\sqrt{C_{j j} S S_{R e s} / (n - p)}}

$t_0 = \frac{\hat \beta_j - 0}{\text{se}(\hat \beta_j)} = \frac{\hat \beta_j}{\sqrt{\hat \sigma^2 C_{jj}}} = \frac{\hat \beta_j}{\sqrt{C_{jj} SS_{Res}/(n-p)}}$ ここで、は

C_{j j}

$C_{jj}$

j^{t h}

$j^{th}$

対角要素です。

(X^{'} X)^{- 1}

$(X'X)^{-1}$

ここまでは順調ですね。Rの行列演算を使用してこれらすべての値を計算する方法を知っています。しかし、ヌルを拒否するために、この本にはが必要だと書かれています

| t_{0} | > t_{α / 2, n - p}

$|t_0| > t_{\alpha/2,n-p}$

R を使用してこの重要な値をどのように計算できますか？ $t_{\alpha/2,n-p}$ 現在、これらの値を見つける方法を知る唯一の方法は、本の巻末の表を調べることです。もっと良い方法があるはずです。

r statistical-significance multiple-regression

— ジュラ紀
ソース

qt（）関数がこれを行います。

— ゲスト

$t_{\alpha/2, n-p}$ qt(1-alpha/2, n-p)

$95\%$ $df=1, 2, \dots, 10$
> qt(.975, 1:10) [1] 12.706205 4.302653 3.182446 2.776445 2.570582 2.446912 2.364624 2.306004 2.262157 2.228139

— クリストファー・アデン
ソース