为线和置信区间创建 ggplot2 geom

Question

我想设计一个几何图形来绘制一条带有置信区间的线。 这将基于的数据框包含以下内容：

1. x 值
2. 主线在每个 x 值处的 y 值
3. 每个 x 的 y=f(x) 的标准误差

例如，

xvals <- seq(0,2*pi,length=100)
df <- data.frame(x=xvals, y=sin(xvals), se=.25)
head(df)
     x           y   se
1 0.00 0.000000000 0.25
2 0.01 0.009999833 0.25
3 0.02 0.019998667 0.25
4 0.03 0.029995500 0.25
5 0.04 0.039989334 0.25
6 0.05 0.049979169 0.25

我按照此处列出的指南编写了以下 geom 函数：

geom_myci <- function(yvar, sevar) {
  list(geom_ribbon(mapping=aes_q(ymin=substitute(yvar-1.96*sevar),
                                 ymax=substitute(yvar+1.96*sevar)),
                   colour="lightgrey", fill="lightgrey"),
       geom_line(mapping=aes_q(y=substitute(yvar)), lwd=2),
       theme_bw())
}

然后可以使用以下方法执行：

ggplot(df, aes(x,y)) + geom_myci(y,se)

这很好用，我唯一不喜欢的是我让用户输入 y 变量两次。 在geom函数中，有没有办法知道已经映射到“y”的变量？

Answer 1

您可以通过aes_string substitute se -Variable，如下所示 - 查看http://adv-r.had.co.nz/Functions.html以了解substitute和延迟评估的工作原理

geom_myci <- function(se_var = se) {
  se_var <- as.character(substitute(se_var))
  list(geom_ribbon(aes_string(ymin = sprintf("y - 1.96 * %s", se_var),
                              ymax = sprintf("y + 1.96 * %s", se_var)),
                   colour="lightgrey", fill="lightgrey"),
       geom_line(lwd=2),
       theme_bw())
}

as.character(substitute(se_var))使se_var成为刺痛。 在这个例子中"se" 。 现在您可以使用它通过sprintf("y - 1.96 * %s", se_var) ，结果是
"y - 1.96 * se" - 我们需要的字符串。

绘制它的结果（按要求）在

ggplot(df, aes(x,y)) + geom_myci()