如何对数据框的所有列进行单向方差分析？答案

【问题标题】：How to do an one way anova with all columns of data fame?如何对数据框的所有列进行单向方差分析？
【发布时间】：2018-05-21 19:04:43
【问题描述】：

我想用我在不同站点的丰富数据进行方差分析，我想知道编写脚本的正确方法是什么，因为所有带有 aov 函数的示例都与一个公式有关，但我不确定这个公式考虑了所有比较它们的平均值的列。

以下数据是我用来尝试进行测试的方式的一个示例。但是列中没有分类变量，因此，我确信这种方式是错误的。

感谢您的帮助

set.seed(200)
D <- data.frame(a=sample(15),b=sample(15), c=sample(15))

A<-aov(a~c, data = D)

【问题讨论】：

您可以发布示例数据吗？请使用dput(df) 的输出编辑问题。或者，如果 dput(head(df, 20)) 的输出太大。（df 是您的数据集的名称。）

标签： r anova

【解决方案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)

【讨论】：