Я хотел бы визуализировать свои данные, используя гистограмму с отображенным значением 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()