r - ggplot2 - 使用自定义函数的误差线

Question

我希望ggplot2使用通用函数为每个数据点生成错误栏，该函数使用该names函数提取相同的列名。下面是一段演示代码：

plotfn <- function(data, xind, yind, yerr) {
    yerrbar <- aes_string(ymin=names(data)[yind]-names(data)[yerr], ymin=names(data) [yind]+names(data)[yerr])
    p <- ggplot(data, aes_string(x=names(data)[xind], y=names(data)[yind]) + geom_point() + geom_errorbar(yerrbar)
    p
}

errdf <- data.frame('X'=rnorm(100, 2, 3), 'Y'=rnorm(100, 5, 6), 'eY'=rnorm(100))
plotfn(errdf, 1, 2, 3)

运行它会出现以下错误：

Error in names(data)[yind] - names(data)[yerr] : 
  non-numeric argument to binary operator

有什么建议么？谢谢。

Answer 1

您需要传递一个包含-( 'a-b'not 'a'- 'b')的字符串

例如，

ggplot(mtcars,aes_string(y = 'mpg-disp',x = 'am')) + geom_point()

在你的例子中

plotfn <- function(data, xind, yind, yerr) {
  # subset the names now so it is slightly less typing later
  yerr_names <- names(data)[c(yind,yerr)]

  yerrbar <- aes_string(ymin = paste(yerr_names, collapse = '-'), 
                         ymax = paste(yerr_names,collapse='+'))
   p <- ggplot(data, aes_string(x=names(data)[xind], y=names(data)[yind])) + 
     geom_point() + 
     geom_errorbar(mapping = yerrbar)
          p
}

# a slightly smaller, reproducible example
set.seed(1)
errdf <- data.frame('X'=rnorm(10, 2, 3), 'Y'=rnorm(10, 5, 6), 'eY'=rnorm(10))
plotfn(errdf, 1, 2, 3)