我想用不同站点的丰富数据进行ANOVA测试，我想知道编写脚本的正确方法是什么，因为所有带aov函数的例子都与公式有关但我不确定这个公式考虑了比较其均值的所有列。
以下数据是我尝试进行测试的方式的一个示例。但是列中没有分类变量，因此，我确信这种方式是错误的。
感谢您的帮助
set.seed(200)
D <- data.frame(a=sample(15),b=sample(15), c=sample(15))
A<-aov(a~c, data = D)
答案 0 :(得分：1)
首先需要reshape the data.frame from wide to long format。
以下内容使用外部包reshape2来重新整理从长格式到长格式的数据。
molten <- reshape2::melt(D)
head(molten)
model <- lm(value ~ variable, data = molten)
anova(model)