Valeurs critiques de t dans R

Question

Je dois déterminer les valeurs critiques de t pour des queues unilatérales de 75% et 99%, pour 40 degrés de liberté.

Le code suivant permet d'obtenir des valeurs t critiques bilatérales à 99 % :

qt(0.01, 40)

Mais comment puis-je déterminer la valeur critique unilatérale de t ?

Answer 1

Le code que vous avez posté donne la valeur critique pour un test unilatéral (la réponse à votre question est donc simple) :

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

Notez que la distribution t est symétrique. Pour un test bilatéral (disons avec un niveau de confiance de 99 %), vous pouvez utiliser la valeur critique suivante

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

Answer 2

Les commentaires de Josh sont tout à fait pertinents. Si vous n'êtes pas très familier avec les valeurs critiques, je vous suggère de jouer avec qt, de lire le manuel ( ?qt ) en conjonction avec la consultation d'une table de recherche ( LIEN ). Lorsque je suis passé de SPSS à R, j'ai créé une fonction qui facilitait la recherche des valeurs t critiques (je ne l'utiliserais plus maintenant, car cela prend trop de temps et, avec les valeurs p qui sont généralement fournies dans les résultats, c'est une question sans intérêt). Voici le code de cette fonction :

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

Answer 3

En prolongeant la réponse de @Ryogi ci-dessus, vous pouvez tirer parti de l'option lower.tail comme suit :

qt(0.25/2, 40, lower.tail = FALSE) # 75% de confiance

qt(0.01/2, 40, lower.tail = FALSE) # 99% de confiance