R을 사용하여 임계 t 값을 어떻게 계산할 수 있습니까?

이것이 새로운 질문이라면 죄송합니다. 통계를 처음으로 가르치려고합니다. 기본 절차가 중단되었다고 생각하지만 R로 실행하는 데 어려움을 겪고 있습니다.

그래서 저는 선형 선형 회귀 형태의 회귀 계수의 중요성을 평가하려고합니다.

\hat{와이} = 엑스 \hat{β}

$\hat y = X \hat \beta$

테스트를위한 t- 통계량 은 다음과 같이 주어집니다. $H_0: \hat \beta_j = 0, H_a: \hat \beta_j \neq 0$

티_{0} = \frac{{\hat{β}}_{제이} - 0}{se ({\hat{β}}_{제이})} = \frac{{\hat{β}}_{제이}}{\sqrt{{\hat{σ}}^{2} 씨_{제이 제이}}} = \frac{{\hat{β}}_{제이}}{\sqrt{씨_{제이 제이} 에스 {에스}_{아르 자형 이자형 에스} / (엔 - 피)}}

$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에서 행렬 연산을 사용하여 이러한 값을 모두 계산하는 방법을 알고 있습니다. 그러나 null을 거부하려면이 책에서

| 티_{0} | > 티_{α / 2, 엔 - 피}

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

$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

— 크리스토퍼 아덴
소스