2012-07-17 41 views
11

Tôi cần xác định giá trị t quan trọng cho đuôi một chiều 75% và 99%, cho 40 độ tự do.Giá trị t quan trọng trong R

Sau đây là mã cho một 99% giá trị t quan trọng hai mặt:

qt(0.01, 40) 

nhưng làm thế nào tôi có thể xác định cho một giá trị t quan trọng một chiều?

+7

Điều gì khiến bạn nghĩ rằng 'qt (0,01, 40)' là giá trị quan trọng cho thử nghiệm hai mặt? Tôi muốn đề nghị đọc '? Qt', và sau đó suy nghĩ thêm một chút về ý nghĩa của các bài kiểm tra một và hai mặt. –

+1

Thực sự đây là một câu hỏi về sự hiểu biết những giá trị quan trọng thực sự có ý nghĩa gì. Đây là một câu hỏi thống kê và có lẽ nên được di chuyển đến trang web stackexchange thống kê. – Dason

Trả lời

20

Mã bạn gửi cho giá trị quan trọng đối với một bài kiểm tra một chiều (xem here Do đó câu trả lời cho bạn câu hỏi đơn giản là:..

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

Lưu ý rằng phân bố t là đối xứng Đối với một 2 kiểm tra -sided (nói với 99% độ tin cậy), bạn có thể sử dụng các giá trị quan trọng

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

bình luận của Josh là chỗ trên. Nếu bạn không phải là siêu quen thuộc với các giá trị quan trọng tôi muốn đề nghị chơi với qt, đọc hướng dẫn (?qt) kết hợp với việc xem xét một bảng tra cứu (LINK). Khi lần đầu tiên tôi chuyển từ SPSS sang RI tạo ra một hàm làm cho giá trị quan trọng t nhìn lên khá dễ dàng (tôi không bao giờ sử dụng điều này ngay bây giờ vì nó mất quá nhiều thời gian và với các giá trị p thường được cung cấp trong đầu ra)). Dưới đây là đoạn code cho rằng:

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 là mức độ tự do của bạn mà không có tùy chọn để thay đổi nó? Âm thanh như bạn chỉ sử dụng nó cho hồi quy tuyến tính đơn giản? – Dason

+0

Spot on it is Stats I. Nhưng điều đó nghe có vẻ như nơi áp phích cũng vậy. –

3

Mở rộng @Ryogi trả lời ở trên, bạn có thể tận dụng lợi thế của tham số lower.tail như vậy:

qt(0.25/2, 40, lower.tail = FALSE) # 75% độ tin cậy

qt(0.01/2, 40, lower.tail = FALSE) # 99% độ tin cậy

Các vấn đề liên quan