2012-07-17 18 views
11

40 derecelik serbestlik için% 75 ve% 99 tek taraflı kuyruklar için kritik t değerlerini belirlemem gerekir.Kritik t değerleri R

aşağıdaki iki taraflı 99% kritik t değerleri için kod:

qt(0.01, 40) 

ama tek taraflı bir kritik t değeri için tespit nasıl?

+7

İki taraflı test için 'qt (0.01, 40)' ın kritik değeri olduğunu düşünürsünüz? 'Qt' okumayı ve ardından bir ve iki taraflı testlerin ne anlama geldiğini biraz daha düşünmeyi öneririm. –

+1

Gerçekten de bu kritik değerlerin gerçekte ne anlama geldiğini anlamakla ilgili bir soru. Bu istatistiksel bir sorudur ve muhtemelen stackexchange sitesinin istatistiklerine taşınmalıdır. – Dason

cevap

20

Yayınladığınız kodu (here bkz Dolayısıyla soruya yanıt vermek basitçe tek taraflı test için kritik değer verir: 2 için

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

Not T dağılımı simetrik olduğunu.. taraflı testi

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

Josh'un yorumlar spot kritik değerini kullanabilir (% 99 güvenle söyleyebiliriz). Eğer ben qt ile oynarken öneririm kritik değerleri ile süper aşina değilseniz, manüel okuma (?qt) ile birlikte bir bakış tablosu (LINK). SPSS'den RI'ye ilk kez taşındığımda, kritik t değerini oldukça kolay gösteren bir fonksiyon yarattı (bunu şimdi çok fazla zaman aldığından ve genelde çıktıda verilen p değerleri ile kullanmam)).

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

n-2, serbest kalma dereceniz olarak değiştirmek için bir seçenek yok mu? Basit doğrusal regresyon için kullandığınız gibi mi geliyor? – Dason

+0

Bu noktadakiler İstatistilerdi. Ama bu posterin de olduğu gibi görünüyor. –

3

@Ryogi yukarıda cevap genişletmeniz şöyle lower.tail parametrenin yararlanabilir,:

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

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

İşte bunun için kod
İlgili konular