2013-08-20 44 views
0

我希望使用通用函數爲ggplot2中的每個數據點生成錯誤欄,該函數使用names函數提取相同的列名。以下是演示代碼:ggplot2 - 使用自定義函數的錯誤條

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 

有什麼建議?謝謝。

回答

7

您將需要通過包含字符串的-'a-b''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) 

enter image description here

相關問題