如何用 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 值（如下所示）的方差分析矩阵，该矩阵在 "Minutes" 上测试 "Measurement" :

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

有什么想法吗？ GGally允许这个？

Answer 1

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

pairwise.t.test(Measurement,Minutes)

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

How to produce an ANOVA matrix with R

r

statistics

anova

ggally