超過直方圖繪製垂直位數線

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))