R - Phân phối nhị thức

Mô hình phân phối nhị thức đề cập đến việc tìm xác suất thành công của một sự kiện chỉ có hai kết quả có thể xảy ra trong một loạt thử nghiệm. Ví dụ, tung đồng xu luôn cho một đầu hoặc một đuôi. Xác suất tìm thấy chính xác 3 đầu khi tung đồng xu liên tục trong 10 lần được ước tính trong phân phối nhị thức.

R có bốn hàm tích hợp để tạo phân phối nhị thức. Chúng được mô tả dưới đây.

dbinom(x, size, prob)
pbinom(x, size, prob)
qbinom(p, size, prob)
rbinom(n, size, prob)

Sau đây là mô tả về các tham số được sử dụng:

x là một vectơ của số.
p là một vectơ xác suất.
n là số quan sát.
size là số lần thử.
prob là xác suất thành công của mỗi lần thử.

dbinom ()

Hàm này cung cấp phân phối mật độ xác suất tại mỗi điểm.

# Create a sample of 50 numbers which are incremented by 1.
x <- seq(0,50,by = 1)

# Create the binomial distribution.
y <- dbinom(x,50,0.5)

# Give the chart file a name.
png(file = "dbinom.png")

# Plot the graph for this sample.
plot(x,y)

# Save the file.
dev.off()

Khi chúng tôi thực thi đoạn mã trên, nó tạo ra kết quả sau:

pbinom ()

Hàm này cung cấp xác suất tích lũy của một sự kiện. Nó là một giá trị duy nhất đại diện cho xác suất.

# Probability of getting 26 or less heads from a 51 tosses of a coin.
x <- pbinom(26,51,0.5)

print(x)

Khi chúng tôi thực thi đoạn mã trên, nó tạo ra kết quả sau:

[1] 0.610116

qbinom ()

Hàm này nhận giá trị xác suất và đưa ra một số có giá trị tích lũy khớp với giá trị xác suất.

# How many heads will have a probability of 0.25 will come out when a coin
# is tossed 51 times.
x <- qbinom(0.25,51,1/2)

print(x)

Khi chúng tôi thực thi đoạn mã trên, nó tạo ra kết quả sau:

[1] 23

rbinom ()

Hàm này tạo ra số lượng giá trị ngẫu nhiên cần thiết của xác suất nhất định từ một mẫu nhất định.

# Find 8 random values from a sample of 150 with probability of 0.4.
x <- rbinom(8,150,.4)

print(x)

Khi chúng tôi thực thi đoạn mã trên, nó tạo ra kết quả sau:

[1] 58 61 59 66 55 60 61 67