Как я могу рассчитать критическое значение т, используя R?

13

Извините, если это новый вопрос; Я пытаюсь научить себя статистике в первый раз. Я думаю, что у меня есть базовая процедура, но я изо всех сил пытаюсь выполнить ее с R.

Итак, я пытаюсь оценить значимость коэффициентов регрессии в множественной линейной регрессии формы

y^=Xβ^

Я думаю, что t-статистика для тестирования определяется какH0:β^j=0,Ha:β^j0

t0=β^j0se(β^j)=β^jσ^2Cjj=β^jCjjSSRes/(np)
где - это диагональная запись в .Cjjjth(XX)1

Все идет нормально. Я знаю, как вычислить все эти значения, используя матричные операции в R. Но чтобы отклонить ноль, книга говорит, что мне нужно

|t0|>tα/2,np

Как я могу вычислить это критическое значение используя R? tα/2,np Прямо сейчас я знаю, как найти эти значения, только посмотрев в таблицу в конце книги. Должен быть лучший способ.

юра
источник
Функция qt () делает это.
Гость

Ответы:

22

Добро пожаловать на сайт Stats Exchange. Вы можете вычислить критическое значение в R делая .tα/2,npqt(1-alpha/2, n-p)

В следующем примере я прошу R дать мне критическое значение для . Результатом является список первых десяти критических значений для t-распределения с заданным уровнем достоверности:95%df=1,2,,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

Кристофер Аден
источник