用扫帚将 t 检验结果合并到数据框中

【问题标题】：combine t-test results into a dataframe with broom用扫帚将 t 检验结果合并到数据框中
【发布时间】：2015-09-27 17:00:10
【问题描述】：

我正在尝试对两个向量 x 和 y 重复采样 100 次，并在每个采样后计算 t 检验。然后我想使用broom 包将结果组合到一个数据框中。我的尝试如下。然而，我所做的只是对相同的两个向量进行 100 次 t 检验。

library(plyr); library(broom)

x <- rnorm(10000, 3, 3)
y <- rnorm(10000, 5, 3)


x_sample <- sample(x, size = 20, replace = FALSE)
y_sample <- sample(y, size = 20, replace = FALSE)

ldply(1:100, function(x) tidy(t.test(x_sample, y_sample)))

【问题讨论】：

标签： r statistics plyr broom

【解决方案1】：

我们可以使用replicate 和mapply。我们为“x”和“y”复制n 的sample 次，并使用mapply 为相应的list 元素执行t.test，然后使用tidy 输出。

n <- 100
mapply(function(x,y) tidy(t.test(x,y)), 
       replicate(n, sample(x, size=20, replace=FALSE), simplify=FALSE), 
       replicate(n, sample(y, size=20, replace=FALSE), simplify=FALSE))

【讨论】：