我想使用显示 kruskal-wallis 检验的 p 值的条形图来可视化我的数据。
问题是我想比较显示的数据的子集,即我想比较 vs=0 的汽车的比例,具体取决于汽车的档位数。
这是我第一次使用 ggpubr(),我不知道如何为此目的使用 stat_compare_means()(如果它真的可能的话??)。
谢谢你的帮助!
#Load packages
library(dplyr)
library(ggpubr)
#Prep data
mtcars2 <- mtcars %>% select(gear,vs) %>%
group_by(gear,vs) %>%
summarise (n = n()) %>%
mutate(proportion = round((n / sum(n)*100),1))
mtcars2$gear <- as.factor(mtcars2$gear)
mtcars2$vs<- as.factor(mtcars2$vs)
#Make plot
ggbarplot(mtcars2, "gear", "proportion",
fill = "vs", color = "vs", palette = "Paired",
label = FALSE,
position = position_dodge(0.9)) + stat_compare_means()