r - 一个代码中的多个二项式测试

Question

我想使用 binom R 包(binom.test 函数)执行精确的二项式测试，以 4 个不同的概率值获得 4 个不同的成功次数。

我可以单独完成，但我想知道是否可以编写一段代码来仅用一个命令(例如 lapply、for 循环)进行计算。

x <- c(68, 69, 70, 75)  #number of successes
p <- c(1/365, 2/365, 3/365, 4/365)  #probability of success
n <- 265 #number of trials

conf.level = 0.95 和 alternative = "two.side"(因为结果可以是 1 或 0)。

有什么建议吗？

---

我尝试过:

for (i in 1:4) {
     test[i] <- binom.test(x[i], n, p[i], alternative = "two.sided", conf.level = 0.95)
      }

但不起作用。

Answer 1

使用maply:

mapply(binom.test, x, p, n=n, alternative = "two.sided", conf.level = 0.95, SIMPLIFY = FALSE)

maply 只是在其第一个参数中调用该函数，并使用其附加参数和附加命名参数中的所有值。 (请参阅mapply函数的帮助)

如果您只想要结果的一个字段，您可以像这样调用 mapply:

mapply(function(x,p,n, ...){
         binom.test(x, n, p, ...)$p.value
       },
       x, p, n=n, alternative = "two.sided", conf.level = 0.95)