Boxplot показать значение среднего

в этом boxplot мы можем видеть среднее значение, но как мы можем иметь также числовое значение на графике для каждого среднего каждого поля?

 ggplot(data=PlantGrowth, aes(x=group, y=weight, fill=group)) + geom_boxplot() +
     stat_summary(fun.y=mean, colour="darkred", geom="point", 
                           shape=18, size=3,show_guide = FALSE)
45
ggplot2 r

4 ответов

во-первых, вы можете рассчитать группу означает с aggregate:

means <- aggregate(weight ~  group, PlantGrowth, mean)

этот набор данных можно использовать с geom_text:

library(ggplot2)
ggplot(data=PlantGrowth, aes(x=group, y=weight, fill=group)) + geom_boxplot() +
  stat_summary(fun.y=mean, colour="darkred", geom="point", 
               shape=18, size=3,show_guide = FALSE) + 
  geom_text(data = means, aes(label = weight, y = weight + 0.08))

здесь + 0.08 используется для размещения метки над точкой, представляющей среднее значение.

enter image description here

альтернативная версия без ggplot2:

means <- aggregate(weight ~  group, PlantGrowth, mean)

boxplot(weight ~ group, PlantGrowth)
points(1:3, means$weight, col = "red")
text(1:3, means$weight + 0.08, labels = means$weight)

enter image description here

60

Вы можете использовать выходное значение stat_summary()

ggplot(data=PlantGrowth, aes(x=group, y=weight, fill=group)) 
+ geom_boxplot() 
+ stat_summary(fun.y=mean, colour="darkred", geom="point", hape=18, size=3,show_guide = FALSE)
+ stat_summary(fun.y=mean, colour="red", geom="text", show_guide = FALSE, 
               vjust=-0.7, aes( label=round(..y.., digits=1)))
18

вы также можете использовать функцию в stat_summary для вычисления среднего и аргумента hjust для размещения текста, вам нужна дополнительная функция, но без дополнительного фрейма данных:

fun_mean <- function(x){
  return(data.frame(y=mean(x),label=mean(x,na.rm=T)))}


ggplot(PlantGrowth,aes(x=group,y=weight)) +
geom_boxplot(aes(fill=group)) +
stat_summary(fun.y = mean, geom="point",colour="darkred", size=3) +
stat_summary(fun.data = fun_mean, geom="text", vjust=-0.7)

enter image description here

10

в Magrittr способом

Я знаю, что уже есть принятый ответ, но я хотел показать один классный способ сделать это в одной команде с помощью magrittr пакета.

PlantGrowth %$% # open dataset and make colnames accessible with '$'
split(weight,group) %T>% # split by group and side-pipe it into boxplot
boxplot %>% # plot
lapply(mean) %>% # data from split can still be used thanks to side-pipe '%T>%'
unlist %T>% # convert to atomic and side-pipe it to points
points(pch=18)  %>% # add points for means to the boxplot
text(x=.+0.06,labels=.) # use the values to print text

этот код создаст boxplot со средствами, напечатанными в виде точек и значений: boxplot with means

Я разделил команду на несколько строк, чтобы я мог комментировать то, что делает каждая часть, но ее также можно ввести как oneliner. Вы можете узнать больше об этом мой суть.

0