R语言绘图｜绘制混合密度函数图以及添加分位数线

简介

最近有个需求，统计课中需要绘制以下图形：

这里我主要使用 ggridges 包中的 stat_density_ridges()。这个包的介绍，小编以前做过一期内容，可见： 。读者需要进一步阅读课件这篇博文，以及一些案例。

加载包

library(ggplot2)
library(ggridges)

产生数据集

假设数据来源于一个混合分布。

item <- 10000
inds <- rbinom(1, item, 0.5)
x <- c(rnorm(inds, 1, 1), rnorm(item - inds, 8, 1))
data <- data.frame("value" = x, "class" = rep(1, length(x)))

绘制密度函数图并添加分位数线

# 绘图
p1 <- ggplot(data, aes(x = value, y = class, fill = factor(stat(quantile)))) +
  stat_density_ridges(
    geom = "density_ridges_gradient",
    calc_ecdf = TRUE,
    quantiles = c(0.025, 0.975)
  ) +
  scale_fill_manual(
    name = "Probability", values = c("#E2EAF6", "#436FB0", "#E2EAF6")
  ) +
  theme_bw() +
  theme(legend.position = "none", panel.grid = element_blank()) +
  labs(x = "x", y = "Density")
p1

p2 <- ggplot(data, aes(x = value, y = class, fill = factor(stat(quantile)))) +
  stat_density_ridges(
    geom = "density_ridges_gradient",
    calc_ecdf = TRUE,
    quantiles = c(0.005, 0.495, 0.51, 0.99)
  ) +
  scale_fill_manual(
    name = "Probability", values = c("#E2EAF6", "#436FB0", "#E2EAF6", "#436FB0", "#E2EAF6"),
  ) +
  theme_bw() +
  theme(legend.position = "none", panel.grid = element_blank()) +
  labs(x = "x", y = "Density")

p2

合并两图

使用 cowplot 包，将两图进行合并。小编对该包的介绍做过几期，可见： 。 其他合并的方式还有：。

library(cowplot)
# pdf("plot_cow.pdf", width = 8, height = 4)
plot_grid(p1, p2, ncol = 1, nrow = 2)
# dev.off()

推荐： 可以保存以下照片，在 b 站扫该二维码，或者 b 站搜索【庄闪闪】观看 Rmarkdown 系列的视频教程。Rmarkdown 视频新增两节视频（写轮眼幻灯片制作）需要视频内的文档，可在公众号回复【rmarkdown】

可视化推文推荐

R可视乎｜空间地理数据可视化（1）

R可视乎|用R给心仪的对象表白吧

R可视乎｜棒棒糖图

R可视乎｜合并多幅图形

R可视乎|等高线图

R可视乎|气泡图