在R中的数据帧上对prop.test进行矢量化

时间:2017-08-23 00:10:38

标签: r dataframe statistics vectorization statistical-test

这是一个基本问题,我收到以下错误。

  

prop.test中的错误:'x'和'n'必须具有相同的长度

使用此代码

cv_MH$pval <- (prop.test(x = c(cv_MH$search, cv_MH$against), n = c(size, size2)))$p.value

其中sizesize2是大型的常数(> 200,000)。

这是cv_MH的样子

search    against
45        23
384       274
657       883

基本上,我正在尝试在cv_MH中创建另一个计算p值的变量。

感谢。

1 个答案:

答案 0 :(得分:2)

我认为您需要为x中的每个值重复计数(n)。 那怎么样?

cv_MH$pval <- prop.test(x = c(cv_MH$search, cv_MH$against), 
                        n = c(rep(size, length(cv_MH$search)),
                              rep(size2, length(cv_MH$against))))$p.value

x表示成功的次数(感兴趣的事件)和n的事件总数... x应该具有相同的n长度,如错误消息所示