Мне нужно определить критические значения t для односторонних хвостов 75% и 99%, для 40 градусов свободы.
Ниже приведен код для двухсторонних значений критического значения на 99%:
qt(0.01, 40)
но как я могу определить для одностороннего критического значения t?
Код, который вы опубликовали, дает критическое значение для одностороннего теста (см. здесь. Следовательно, ответ на ваш вопрос просто:
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
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()
Продолжая ответ @Ryogi выше, вы можете воспользоваться параметром lower.tail следующим образом:
qt(0.25/2, 40, lower.tail = FALSE) # 75% доверия
qt(0.01/2, 40, lower.tail = FALSE) # 99% доверия