我想对子组均值进行统计比较,然后用上面的p值在它们之间画一条水平线。我知道我可以使用stat_compare_means在x轴因子级别上执行此操作。但是,我不知道如何比较子组。
在此示例中,我可以比较x轴因数1和2,但我想比较x轴因数1子组4和x轴因数2子组4。
library(ggplot2)
library(ggpubr)
ggplot(mtcars, aes(x=as.factor(carb), y=hp, fill=as.factor(cyl))) +
geom_boxplot(outlier.shape=NA) +
stat_compare_means(comparisons = list(c("1", "2")))