在R中如何使用ggplot绘制正态分布的尾部区域？

Question

我找到了一种“破解”ggplot 的方法，方法是结合两个 geom_area 图来创建带有尾部区域的正态分布：

library(ggplot2)
mean <-  0
standard_deviation <- 1
Zscore <- -1.35

observation = (Zscore*standard_deviation) + mean
(tail_area <- round(pnorm(observation),2))

ggplot(NULL, aes(c(-5,5))) +
    geom_area(stat = "function", fun = dnorm, fill="sky blue", xlim = c(-5, -1.35)) +
    geom_area(stat = "function", fun = dnorm,  xlim = c(-1.35, 5))

是否有使用 ggplot 来创建正态分布并像上面那样突出显示尾部区域的“不那么老套”的方法？

Answer 1

首先，我喜欢你的方法；不确定这是否不那么“hackey”，但这是使用gghighlight的另一个选项

# Generate data (see comment below)
library(dplyr)
df <- data.frame(x = seq(-5, 5, length.out = 100)) %>% mutate(y = dnorm(x))

# (gg)plot and (gg)highlight
library(ggplot2)
library(gghighlight)
ggplot(df, aes(x, y)) + geom_area(fill = "sky blue") + gghighlight(x < -1.35)

据我了解，gghighlight 需要一个 data 参数，因此它不能单独使用 geom_area（意思是：没有 data 但使用 stat = "function"）或使用 stat_function。这就是为什么我首先生成数据df。

---

更新

回应您关于如何“突出显示 1 和 -1 之间的区域”的评论；您可以执行以下操作

ggplot(df, aes(x, y)) + geom_area(fill = "sky blue") + gghighlight(abs(x) < 1)

更新 2

要突出显示区域1.5 < x < 2.5，只需使用条件语句x > 1.5 & x < 2.5

ggplot(df, aes(x, y)) + geom_area(fill = "sky blue") + gghighlight(x > 1.5 & x < 2.5)

---

为了抢占潜在的后续问题：此方法仅适用于连续区域。意思是，我还没有找到在单个 gghighlight 语句中突出显示 x < -2.5 & x > 2.5 的方法。