如何使用 R 生成方差分析矩阵

Question

我有以下数据集：

Minutes <- rep(0:10, each=5)
Minutes <- factor(Minutes) 
Subject <- rep(1:5,time=11)
Measurement <- rnorm(55, mean=40, sd=7)
df <- data.frame(Subject, Measurement, Minutes)

我想要一个带有 p 值（如下所示）的方差分析矩阵，用于测试"测量"而不是"分钟"：

         Minute0 Minute1 Minute2 Minute3 Minute4 ... Minute10
Minute0
Minute1
Minute2
Minute3
Minute4
...
Minute10

有什么想法吗？允许吗？

Answer 1

您可以使用 pairwise.t.test {stats} 函数，在该函数中指定响应变量，然后使用包含要比较的因子水平的分类解释变量，用逗号分隔

pairwise.t.test(Measurement,Minutes)