在r中进行大量的2样本比例测试

Question

我正在尝试为比例测试创建一个模拟，以查看多久可以得到一次有意义的结果（即p / value <0.05）。 我假设人口规模不平等，例如

test <- rbinom(100, 5000, 0.05)

ctrl <- rbinom(100, 10000, 0.04)

在这里，我有100个模拟测试，其中测试人口为5000个，对照组为10000个，以及转换率的实际差异（测试为0.05，对照组为0.04）。

我现在想为每个模拟测试运行一个比例测试，并获得p值。 因此，对于第一个测试，我会做：

prop.test(c(test[1], ctrl[2]), c(5000, 10000))$p.value

然后使用for循环遍历此过程，以获取每个模拟测试的p值。

我正在寻找一种更紧凑的方式。 我发现了一个办法做到这一点与地图是1个样品比试验在这里 ，但我坚持如何做到这一点的2样本检验。 看来我应该能够在一行中完成此操作，并且只需获取p值向量作为输出即可。

我尝试了以下

Map(prop.test, c(test,ctrl), c(5000, 10000), alternative="two.sided", correct = F)

但这给了我200个1个样本的比例测试，而不是100个2个样本的比例测试...我有什么办法可以解决？

理想情况下，我将得到一个长度为100的向量，其p.values值使我可以检查哪个部分低于5％，这就是测试功效。

Answer 1

如果我理解正确，则需要每个道具的每个p.value。 测试。 你有尝试过吗？ prop.test（test，ctrl）[[4]]