Krytyczne wartości t w R

Question

Muszę określić krytyczne wartości t dla jednostronnych ogonów 75% i 99%, dla 40 stopni swobody.

Poniżej kod dwustronny 99% wartości krytycznych t:

qt(0.01, 40) 

ale jak mogę określić na jednostronnym krytycznej wartości t?

Answer 1

Kod pan pisał daje wartość krytyczną dla jednostronnego testu (patrz here Stąd odpowiedź na ty pytanie jest prosta:..

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

Zauważ, że t-rozkład jest symetryczny Przez 2 -sided Test (powiedzmy z 99% pewnością) można użyć wartości krytycznej

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

Answer 2

komentarze Josha są na miejscu. Jeśli nie są bardzo znane z wartościami granicznymi Sugerowałbym gry z Qt, czytanie instrukcji (?qt) w połączeniu z patrzeniem tabela wyszukiwania (LINK). Kiedy po raz pierwszy przeniosłem się z SPSS do RI utworzyłem funkcję, która sprawiała, że ​​krytyczna wartość t wyglądała całkiem łatwo (nigdy bym jej nie używał, ponieważ zajmuje to zbyt dużo czasu, a przy wartościach p, które są generalnie dostarczane na wyjściu, jest to kwestia sporna). Oto kod na to:

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

Rozszerzanie @Ryogi odpowiedzieć na powyższe, można skorzystać z parametru lower.tail tak:

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

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