2

binomR package( binom.testfunction) を使用して、4つの異なる確率値で4つの異なる成功数に対して正確な二項検定を実行したいと思います。

私はそれを単独で行うことができますが、1 つのコマンド (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.95and alternative = "two.sided"(結果は 1 または 0 になる可能性があるため) 。

なにか提案を？

私は試した：

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

しかし、うまくいきません。

4

2 に答える 2

1

結果のp値のみに関心がある場合は、次のことができます。

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
test <- numeric(4)

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

test
[1] 6.621447e-111  1.801758e-92  3.467288e-82  2.442975e-81
于 2016-07-28T10:34:32.967 に答える