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