Valores t críticos em R

Preciso de determinar os valores críticos de T para caudas unilaterais de 75% e 99%, para 40 graus de liberdade.

o seguinte código é para um valor t crítico de 99% em duas faces:

qt(0.01, 40)
Mas como posso determinar um valor t crítico unilateral?

 12
r
Author: luciano, 2012-07-17

3 answers

O código que publicou dá o valor crítico para um teste unilateral (daí a resposta à sua pergunta é simplesmente:

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))

Note que a distribuição t é simétrica. Para um teste de 2 lados (digamos com 99% de confiança) pode usar o valor crítico

abs(qt(0.01/2, 40)) # 99% confidence, 2 sided
 21
Author: Ryogi, 2018-04-17 06:58:21
Os comentários do Josh estão certos. Se você não está super familiarizado com os valores críticos, eu sugiro jogar com o qt, ler o manual (?qt) em conjunto com olhar para uma tabela de olhar para cima (LINK). Quando eu primeiro movi de SPSS para R I criou uma função que fez o valor t crítico olhar para cima muito fácil (eu nunca usaria isso agora porque leva muito tempo e com os valores p que são geralmente fornecidos na saída é um ponto discutível). Aqui está o código para isso:
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()
 4
Author: Tyler Rinker, 2012-07-17 15:55:42

Estendendo a resposta @Ryogi acima, você pode tirar proveito do parâmetro lower.tail Assim:

qt(0.25/2, 40, lower.tail = FALSE) # 75% confiança

qt(0.01/2, 40, lower.tail = FALSE) # 99% confiança

 3
Author: sedeh, 2015-05-02 01:08:28