超过直方图绘制垂直位数线

Question

我目前中的R用ggplot生成以下情节：

的数据被存储在单个数据帧具有三列：PDF（y轴中的情节上述），mids（x）和数据集名称。这是从直方图创建的。
我想要做的就是绘制为代表的95位数，就像我手工画下面作为一个例子中，每个数据集彩色编码的垂直线：

我试图用+ geom_line(stat="vline", xintercept="mean")当然，但我正在寻找分位数，而不是中位数，AFAIK ggplot不允许。颜色很好。
我也试过+ stat_quantile(quantiles = 0.95)，但我不确定它究竟做了什么。文档非常稀少。颜色，也很好。

请注意，密度值非常低，低至1e-8。我不知道quantile（）函数是否喜欢。

据我所知，计算直方图的分位数与计算数列表的分位数并不完全相同。我不知道它会有什么帮助，但HistogramTools包中包含用于直方图分位数的函数ApproxQuantile()。

最小工作示例如下。正如你所看到的，我从每个直方图中获得一个数据帧，然后将这些数据帧绑定在一起并绘制出来。

library(ggplot2) 
v <- c(1:30, 2:50, 1:20, 1:5, 1:100, 1, 2, 1, 1:5, 0, 0, 0, 5, 1, 3, 7, 24, 77) 
h <- hist(v, breaks=c(0:100)) 
df1 <- data.frame(h$mids,h$density,rep("dataset1", 100)) 
colnames(df1) <- c('Bin','Pdf','Dataset') 
df2 <- data.frame(h$mids*2,h$density*2,rep("dataset2", 100)) 
colnames(df2) <- c('Bin','Pdf','Dataset') 
df_tot <- rbind(df1, df2) 

ggplot(data=df_tot[which(df_tot$Pdf>0),], aes(x=Bin, y=Pdf, group=Dataset, colour=Dataset)) + 
geom_point(aes(color=Dataset), alpha = 0.7, size=1.5) 

Answer 1

预计算这些值并单独绘制它们看起来像是最简单的选项。与dplyr这样做需要最小的努力：

library(dplyr) 
q.95 <- df_tot %>% 
    group_by(Dataset) %>% 
    summarise(Bin_q.95 = quantile(Bin, 0.95)) 

ggplot(data=df_tot[which(df_tot$Pdf>0),], 
     aes(x=Bin, y=Pdf, group=Dataset, colour=Dataset)) + 
    geom_point(aes(color=Dataset), alpha = 0.7, size=1.5) + 
    geom_vline(data = q.95, aes(xintercept = Bin_q.95, colour = Dataset))