Расчет доверительных интервалов медианы в ggplot

Question

У меня есть набор данных, состоящий из виртуальной популяции с несколькими наблюдениями на каждый идентификатор. Это население делится на две группы. Я хотел бы построить график зависимости наблюдений от времени для каждой группы. Я мог успешно вычислить среднее значение и доверительный интервал среднего, используя stat_summary. Однако на втором графике я хотел бы сделать то же самое для медианы. stat_summary(fun.y=median) действительно работает, но я не смог найти решения для расчета правильных доверительных интервалов для медианы. Должен ли я рассчитывать их отдельно?

Вот соответствующая часть кода ggplot:

my.labels <- c("1250 mg/m²\nwith dose\nadjustments",
           "1250 mg/m²\nwithout dose\nadjustments")
ggplot(data, aes(x=TIME, y=P2X))+ 
      stat_summary(geom="ribbon", fun.data=mean_cl_boot, 
    fun.args=list(conf.int=0.95), 
                  aes(fill=factor(GROUP),color=NA), alpha = .55)+
      stat_summary(geom="line", fun.y=mean, aes(linetype=factor(GROUP)))+
      scale_fill_manual("dose",labels=my.labels,values=c("#999999", "lightblue"))+
      scale_linetype_manual("dose",labels=my.labels,values = c("solid","dashed"))+
      scale_color_manual("dose",values=c("black", "black"))+labs(x="Time [cycle]",y="Probability [%]")+
      scale_x_continuous(breaks = c(0,1,2,3,4,5,6))+ theme(legend.title=element_blank())+
      scale_y_continuous(labels = function(x) paste0(x*100))

Заранее спасибо!