R language drawing ｜ drawing mixed density function diagram and adding quantile line

brief introduction

There's a need recently , The following figures need to be drawn in the statistics class ：

I mainly use ggridges In bag stat_density_ridges(). Introduction to this bag , Xiaobian has done an issue of content before , so ： . Readers need to read this courseware further post , As well as some Case study .

Load package

library(ggplot2)
library(ggridges)

Generating data sets

Suppose the data comes from a mixed distribution .

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

Plot density function and add quantile lines

#  mapping 
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

Merge two graphs

Use cowplot package , Merge the two figures . Xiaobian has made several issues of introduction to the package , so ： . There are other ways to merge ：.

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

recommend ： You can save the following photos , stay b The station scans the QR code , perhaps b Station search 【 Zhuang Shanshan 】 watch Rmarkdown Series of video tutorials .Rmarkdown Two new videos （ Write wheel eye slide making ） Need documents in video , Reply in official account 【rmarkdown】

Visual tweet recommendation

R Visible ｜ Spatial geographic data visualization （1）

R Visible | use R Tell the person you want

R Visible ｜ Lollipop chart

R Visible ｜ Merge multiple drawings

R Visible | Contour map

R Visible | Bubble chart