重复的随机样本用于R中的t检验

Question

我有两个时期（1999年和2009年）的数据。数据分为两列，一列称为“值”（4.960001,4.847222，...... 3.639985，3.849003），另一列称为“时间”（1999年，1999年，2009年，2009年）：

每年有75,049个值。我想从每年中随机选择30个值（无需替换）并进行配对T检验，然后将结果输入矩阵。我想在同一个数据集上运行该程序1000次以获得T检验值的分布。 关于如何做到这一点的任何想法将不胜感激！感谢。

Answer 1

“....两个栏目，一个名为”值“（4.960001,4.847222，... 3.639985,3.849003），一个名为”时间“（1999年，1999年，...，2009年，2009年）”.. “每年75,049个值。”

听起来可以这样做：

 folded <- matrix(dfrm$values, ncol=2)
 replicate( 1000, { pick <- sample(1:75049, 30)
                    t.test(folded[pick, 1], folded[pick,2], paired=TRUE)$statistic
                   } )