小编chr*_*nsn的帖子

我目前正在绘制ggplot中许多回归模型的一些不同的第一个差异分布.为了便于解释差异,我想标记每个分布的2.5%和97.5%百分位数.由于我将做很多图,并且由于数据按二维(模型和类型)分组,我想在ggplot环境中定义和绘制各自的百分位数.使用facets绘制分布使我能够准确到达我想要的位置,除了百分位数.我当然可以手动完成这项工作,但理想情况下我希望找到一个我仍然可以使用的解决方案facet_grid,因为这样可以避免我在试图将不同的情节融合在一起时遇到很多麻烦.

以下是使用模拟数据的示例:

df.example <- data.frame(model = rep(c("a", "b"), length.out = 500), 
                      type = rep(c("t1", "t2", "t2", "t1"), 
                      length.outh = 250), value = rnorm(1000))

 ggplot(df.example, aes(x = value)) +
 facet_grid(type ~ model) +
 geom_density(aes(fill = model, colour = model))

我试图以两种方式添加分位数.第一个产生错误消息:

 ggplot(df.example, aes(x = value)) +
 facet_grid(. ~ model) +
 geom_density(aes(fill = model, colour = model)) +
 geom_vline(aes(x = value), xintercept = quantile(value, probs = c(.025, .975)))

Error in quantile(value, probs = c(0.025, 0.975)) : object 'value' not found

Run Code Online (Sandbox Code Playgroud)

而第二个获取完整变量的分位数而不是子密度.也就是说,绘制的分位数对于所有四个密度是相同的. …