我使用此代码执行单向方差分析（ANOVA）
anova_models <- lapply(mock_data[sapply(mock_data, is.numeric)], function(x,f) aov(x~f), f=mock_data$weather)
citation for the code。要计算Fisher LSD测试，并尝试使用代码
library(agricolae)
lsd.test <- LSD.test(anova_models, weather, p.adj = "fdr", group = FALSE)
但是它不起作用。 您能否帮助我如何对数据框执行单向ANOVA和Fisher LSD。 Here is the data file.大多数实例是针对一个数字列和一个分类列的。多个数字列和一个分类列的处理方式是什么？