Bagaimana saya bisa menghitung nilai t kritis menggunakan R?

Maaf jika ini adalah pertanyaan newb; Saya mencoba untuk belajar statistik sendiri untuk pertama kalinya. Saya pikir saya memiliki prosedur dasar, tapi saya berjuang untuk melaksanakannya dengan R.

Jadi, saya mencoba untuk mengevaluasi signifikansi koefisien regresi dalam bentuk regresi linier berganda

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

$\hat y = X \hat \beta$

Saya pikir t-statistik untuk menguji diberikan oleh $H_0: \hat \beta_j = 0, H_a: \hat \beta_j \neq 0$

di manaadalahdiagonal masuk dalam.

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}$

Sejauh ini bagus. Saya tahu bagaimana menghitung semua nilai-nilai ini menggunakan operasi matriks dalam R. Tetapi untuk menolak nol, buku itu mengatakan saya perlu

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

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

Bagaimana saya bisa menghitung nilai kritis ini menggunakan R? $t_{\alpha/2,n-p}$ Saat ini satu-satunya cara saya tahu bagaimana menemukan nilai-nilai ini adalah dengan melihat tabel di belakang buku. Pasti ada cara yang lebih baik.

r statistical-significance multiple-regression

— jurassic
sumber

Fungsi qt () melakukan ini.

— tamu

$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

— Christopher Aden
sumber