我有以下数据集:
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值的ANOVA矩阵(如下所示),测试“测量”超过“分钟”:
Minute0 Minute1 Minute2 Minute3 Minute4 ... Minute10
Minute0
Minute1
Minute2
Minute3
Minute4
...
Minute10
任何想法? GGally允许这个?
答案 0 :(得分:1)
您可以使用在其中指定响应变量的pairwise.t.test {stats}函数,然后使用包含要比较的因子级别的分类说明变量,用逗号分隔
pairwise.t.test(Measurement,Minutes)