Hàm dt(),qt(),rt() và pt() trong R là gì?

Những hàm này được sử dụng đối với phân phối Student T test trong R

• dt(): Hàm trả giá trị của hàm măt độ xác suất của phân phối Student T test tại giá trị x và bậc tự do là df
• pt(): Hàm trả giá trị của Hàm xác suất tích lũy của phân phối Student T test tại giá trị x và bậc tự do là df
• qt(): Hàm trả giá trị nghịch đảo của hàm xác suất tích lũy trong Student T test tại giá trị x và bậc tự do là df
• rt(): Hàm trả kết quả là vector của các biến ngẫu nhiên tuân theo phân phối Student T test. Chiều dài của vector là n và bậc tự do là df

R Syntax in R

• dt(x,df)

VÍ dụ: Tìm giá trị của hàm mật độ xác suất tại x = 3 với bậc tự do là 10 
dt(x = 3, df = 10)

## [1] 0.01140055

Tìm xác suất xảy ra tại giá trị x trong hình chuông (hàm pdf)

- Bước 1: Tạo một tập hợp các giá trị x ngẫu nhiên
x <- seq(-3,3, length = 500)  

- Bước 2: Tạo một vector của các giá trị, tại đó cho thấy chiều cao của phân phối xác suất
y <- dt(x = x, df = 20)

- Bước 3: Đồ thị hóa 

plot(x,y, type = "l")
axis(1, at = -3:3, labels = c("-3", "-2", "-1", "0", "1", "2", "3"))

• pt(x, df, lower.tail = FALSE)

Tìm diện tích phần bên trái tại giá trị x = 0.5 và bậc tự do df = 15

pt(0.5, df = 15, lower.tail = TRUE)  # TRUE nghĩa là chấp nhận phần diện tích bên trái tại x = 0.5 và df = 14

## [1] 0.6878349

-rt(n,df)

Tạo 1000 biến ngẫu nhiên bất kỳ tuân theo phân phối Student T test với bậc tự do df = 45

a <- rt(n = 1000, df = 45)

Để xem đồ thị các biến tạo ra như thế nào, ta plot

hist(a, breaks = 50, xlim = c(-5,5))

_ qt(x,df)

qt( 0.9, df = 10)

## [1] 1.372184