stat_compare_means 所选组的方差分析

stat_compare_means anova on selected groups

我想计算一个方差分析,只包括 gr1-gr3,不包括 gr4。这可能使用 stat_compare_means()?

x   value
gr1 3.543
gr1 2.12
gr1 4.56
gr2 3.3
gr2 2.12
gr2 4.1
gr3 3.32
gr3 2.3
gr3 3.3
gr4 3.325
gr4 2.355
gr4 3.34

library(ggpubr)
ggplot(df(aes(x,value)) + geom_boxplot() + stat_compare_means()

这可能会有所帮助。有一个 comparisons 参数可以用于此目的。

ggplot(dta, aes(x, value)) + 
    geom_boxplot() + 
    stat_compare_means(comparisons = list(1:2, 2:3, c(1, 3), c(1, 2, 3)))

编辑:方差分析

ggplot(dta, aes(x, value)) + 
    geom_boxplot() + 
    stat_compare_means(method = "anova", label.y = 4.8, 
                       data = dta %>% filter(x != "gr4"))