r语言 - 分别获得每个唯一因素的ANOVA汇总统计

Question

嗨，我有一个大型数据集，其中包含不同地点每个季节多个物种的同位素值，我从中进行了方差分析:

Anova <- Isotopes %>% 
group_by(Species) %>% 
do(model = aov(d15N~Season+Location, data = Isotopes))
Anova$model

我想分别看到每个独特物种的测试总结。我怎样才能得到这些数据?

提前谢谢你。

Answer 1

如果我理解了你的问题，你可以使用lapply:

# Reproduce a similar dataset
Anova <- iris %>% 
mutate(new_var = sample(LETTERS[1:3], size = nrow(.), replace = T)) %>%
group_by(Species) %>% 
do(model = aov(Sepal.Length~Petal.Width+new_var, data = .))
summ <- lapply(Anova$model,summary)
names(summ) <- Anova$Species
summ

$setosa
Df Sum Sq Mean Sq F值Pr(>F)
花瓣。宽度1 0.471 0.4709 3.919 0.0538 .
new_var 2 0.090 0.0450 0.374 0.6898
残差46 5.527 0.1202
—显著性代码:0 '' 0.001 '' 0.01 '' 0.05 '。' 0.1 ' ' 1

$versicolor
Df Sum Sq Mean Sq F value Pr(>F)
花瓣。宽度1 3.899 3.899 22.432 2.12e-05 ***
new_var 2 1.162 0.581 3.344 0.0441 *
残差46 7.994 0.174
—显著性。代码:0 '' 0.001 '' 0.01 '' 0.05 '。' 0.1 ' ' 1

$virginica
Df Sum Sq Mean Sq F值Pr(>F)
花瓣。宽度1 1.566 1.5656 4.049 0.0501 .
new_var 2 0.459 0.2295 0.593 0.5566
残差46 17.788 0.3867
—显著性代码:0 '' 0.001 '' 0.01 '' 0.05 '。' 0.1 ' ' 1