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