success_probability <- function(P) {
# 前3次成功概率 50%
p1 <- 0.5
n1 <- 3
# 后4次成功概率 25%
p2 <- 0.25
n2 <- 4
# 计算所有可能成功次数的概率
total_prob <- 0
for (k1 in 0:n1) { # 遍历前3次的成功次数
for (k2 in 0:n2) { # 遍历后4次的成功次数
k_total <- k1 + k2 # 总成功次数
if (k_total >= P) { # 只有满足至少P次成功才累加
total_prob <- total_prob +
dbinom(k1, n1, p1) * dbinom(k2, n2, p2)
}
}
}
return(total_prob)
}
# 示例:计算至少3次成功的概率
P <- 3
result <- success_probability(P)
cat("至少", P, "次成功的概率:", result, "\n")