r - stat_compare_means 选定组的方差分析

Question

我想计算一个只包括 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()

标签： rggplot2ggpubr

Answer 1

这可能会有所帮助。有一个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"))