Блог

Я хотел бы, чтобы boxplot суммировал распределение некоторых базовых данных, но вместо усов, простирающихся до 5-го и 95-го, я бы хотел, чтобы блоки простирались до 5-го и 95-го.

Стандартный коробочный график с удаленными выбросами и усами:

 library("ggplot2")
p <- ggplot(mtcars, aes(factor(cyl), mpg))
p   geom_boxplot(outlier.shape = NA, coef = 0)
  

Boxplot с усами на 5-м и 95-м процентилях:

 p   stat_summary(geom = "boxplot", 
                 fun.data = function(x) setNames(quantile(x, c(0.05, 0.25, 0.5, 0.75, 0.95)), 
                                                 c("ymin", "lower", "middle", "upper", "ymax")))
  

Но чего я действительно хочу, так это чтобы ячейки (без усов) простирались до 5-го и 95-го, так что комбинация обеих этих модификаций. Есть ли способ указать функцию, генерирующую поле в stat_summary() ?

Это может быть немного «хакерский» способ сделать это, но самым простым способом может быть просто использовать geom_segment для каждого класса цилиндров. Это позволит вам указать ширину boxplot и значения, которые вы хотите, чтобы boxplot достигал. Но затем вы могли бы поиграть с aes() ними, а также добавить срединную линию, используя stat_summary() при желании.

 library(ggplot2)
p <- ggplot(mtcars, aes(factor(cyl), mpg))
p   geom_segment(aes(x = 4, xend = 4, y = quantile(subset(mtcars, mtcars$cyl==4)$mpg,0.95), yend = quantile(subset(mtcars, mtcars$cyl==4)$mpg, 0.05)), color = 'firebrick1', lwd = 28)