我有一个看起来像这样的数据框。有 2 个独立的组和 5 个不同的变量。
df <- read.table(text="Group var1 var2 var3 var4 var5
1 3 5 7 3 7
1 3 7 5 9 6
1 5 2 6 7 6
1 9 5 7 0 8
1 2 4 5 7 8
1 2 3 1 6 4
2 4 2 7 6 5
2 0 8 3 7 5
2 1 2 3 5 9
2 1 5 3 8 0
2 2 6 9 0 7
2 3 6 7 8 8
2 10 6 3 8 0", header = TRUE)
我正在计算每个变量的显着性,以便使用 T 检验(如下所示)区分两组。但是我想实现置换测试来计算 p 值,因为这是一个非常小的数据集。在 R 中执行此操作的最佳方法是什么?
t(sapply(df[-1], function(x)
unlist(t.test(x~df$Group)[c("p.value")])))