在R中保存qplot外部函数

Question

所以我有一个函数可以绘制估计值及其95％的可信区间，并计算真值在区间内的次数。

使用ggplot2，我绘制间隔，但我想将它保存到我的函数之外的pdf，因为我需要为多个数据集重复这个。

这是我的功能：

require("ggplot2")

CIbias <- function(truevalue, data, cilow, cihigh, PP){ 
  ppn = 1:PP 
limits <- aes(ymax = cihigh, ymin = cilow)

ciplot = qplot(ppn,data)+geom_errorbar(limits, width = 1)+geom_abline(intercept = truevalue,
 slope = 0) 

true = 0  
for(i in 1:PP){
  if(truevalue > cilow[i] & truevalue < cihigh[i]){
  true = true + 1
  } 
}
bias = (true / PP)*100
  return(list(ciplot, bias)) 
}

以下是虹膜数据集的功能示例：

data(iris)

CIbias(6, iris[,2], iris[,2 ] - sd(iris[,2]), iris[,2] + sd(iris[,2]), 150)

Answer 1

使用pdf中的CIbias功能，例如

CIbias <- function(truevalue, data, cilow, cihigh, PP, fname){ 
    ....
    pdf(fname) 
    print(ciplot)
    dev.off()

    return(....)
}