Kritische t-waarden in R

Ik moet de kritische t-waarden voor eenzijdige staarten van 75% en 99% bepalen, voor 40 vrijheidsgraden.

Het volgende is de code voor een tweezijdige 99% kritische t-waarden:

qt(0.01, 40)

maar hoe kan ik bepalen voor een eenzijdige kritische t-waarde?


Antwoord 1, autoriteit 100%

De code die je hebt gepost geeft de kritieke waarde voor een eenzijdige test (vandaar het antwoord op je vraag is eenvoudig:

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

Merk op dat de t-verdeling symmetrisch is. Voor een 2-zijdige test (zeg met 99% betrouwbaarheid) kunt u de kritische waarde gebruiken

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

Antwoord 2, autoriteit 14%

De opmerkingen van Josh zijn perfect. Als je niet zo bekend bent met kritieke waarden, raad ik aan om met qt te spelen, de handleiding (?qt) te lezen in combinatie met een opzoektabel (LINK). Toen ik voor het eerst van SPSS naar RI ging, creëerde ik een functie waarmee de kritische t-waarde vrij eenvoudig kon worden opgezocht (ik zou dit nu nooit gebruiken omdat het te veel tijd kost en met de p-waarden die over het algemeen in de uitvoer worden gegeven, is het een betwistbaar punt ). Hier is de code daarvoor:

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

Antwoord 3, autoriteit 8%

Als je het antwoord van @Ryogi hierboven uitbreidt, kun je als volgt profiteren van de parameter lower.tail:

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

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

LEAVE A REPLY

Please enter your comment!
Please enter your name here

Other episodes