Критические значения t в R

Мне нужно определить критические значения t для односторонних хвостов 75% и 99%, для 40 градусов свободы.

Ниже приведен код для двухсторонних значений критического значения на 99%:

qt(0.01, 40)

но как я могу определить для одностороннего критического значения t?

Ответ 1

Код, который вы опубликовали, дает критическое значение для одностороннего теста (см. здесь. Следовательно, ответ на ваш вопрос просто:

abs(qt(0.25, 40)) # 75% confidence, 1 sided (same as qt(0.75, 40))
abs(qt(0.01, 40)) # 99% confidence, 1 sided (same as qt(0.99, 40))

Заметим, что t-распределение является симметричным. Для двухстороннего теста (скажем, с доверием 99%) вы можете использовать критическое значение

abs(qt(0.01/2, 40)) # 99% confidence, 2 sided

Ответ 2

Josh комментарии на месте. Если вы не знакомы с критическими значениями, я предлагаю играть с qt, читая руководство (?qt) в сочетании с поиском таблицы поиска (LINK). Когда я впервые перешел из SPSS в RI, создал функцию, которая сделала критическое значение t довольно простым (я никогда не буду использовать это сейчас, так как занимает слишком много времени и с значениями p, которые обычно предоставляются на выходе, это спорная точка). Вот код для этого:

critical.t <- function(){
    cat("\n","\bEnter Alpha Level","\n")
    alpha<-scan(n=1,what = double(0),quiet=T)
    cat("\n","\b1 Tailed or 2 Tailed:\nEnter either 1 or 2","\n")
    tt <- scan(n=1,what = double(0),quiet=T)
    cat("\n","\bEnter Number of Observations","\n")
    n <- scan(n=1,what = double(0),quiet=T)
    cat("\n\nCritical Value =",qt(1-(alpha/tt), n-2), "\n")
}

critical.t()

Ответ 3

Продолжая ответ @Ryogi выше, вы можете воспользоваться параметром lower.tail следующим образом:

qt(0.25/2, 40, lower.tail = FALSE) # 75% доверия

qt(0.01/2, 40, lower.tail = FALSE) # 99% доверия