ggplot2 / ggpubr визуализирует уровни значимости для глобальных kruskal-wallis подмножеств сгруппированных данных - PullRequest
0 голосов
/ 06 июня 2018

Я хотел бы визуализировать свои данные, используя гистограмму с отображенным значением p теста Крускала-Уоллиса.

Проблема в том, что я хотел бы сравнить подмножества отображаемых данных, т.е. я хотел бы сравнить долю автомобилей, для которых vs = 0, в зависимости от количества передач, которые есть в автомобиле.

Я впервые использую ggpubr (), и я не могу понять, как с этой целью работать stat_compare_means () (если это действительно возможно ??).

Спасибо за любую помощь!

enter image description here

#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() 
...