使用 ggplot2 从函数中绘制 data.frame

Question

我有这个函数来获取 IRT 包 sirt 返回的对象，并为用户可以指定的一组项目绘制项目响应函数：

plotRaschIRF <- function(x,items=NULL,thl=-5,thu=5,thi=.01,D=1.7) {
  if (!class(x)=="rasch.mml") stop("Object must be of class rasch.mml")
  thetas <- seq(thl,thu,thi)
  N <- length(thetas)
  n <- length(x$item$b)
  tmp <- data.frame(item=rep(1:n,each=N),theta=rep(thetas,times=n),b=rep(x$item$b,each=N))
  probs <- exp(D*(tmp[,2]-tmp[,3]))/(1+exp(D*(tmp[,2]-tmp[,3])))
  dat <- data.frame(item=rep(1:n,each=N),theta=rep(thetas,times=n),b=rep(x$item$b,each=N),p=probs)
  #dat$item <- factor(dat$item,levels=1:n,labels=paste0("Item",1:n))
  if (is.null(items)) {
    m <- min(10,n)
    items <- 1:m
    if (10<n) warning("By default, this function will plot only the first 10 items")
  }
  if (length(items)==1) {
    title="Item Response Function"
  } else {
    title="Item Response Functions"
  }
  dat2 <- subset(dat,eval(quote(eval(item,dat) %in% items)))
  dat2$item <- factor(dat2$item,levels=unique(dat2$item),labels=paste0("Item",unique(dat2$item)))
  out <- ggplot(dat2,aes(x=theta,y=p,group=item)) +
   geom_line(aes(color=dat2$item),lwd=1) + guides(col=guide_legend(title="Items")) +
   theme_bw() + ggtitle(title) + xlab(expression(theta)) +
   ylab("Probability") + scale_x_continuous(breaks=seq(thl,thu,1))
  print(out)
}

但它似乎卡在我开始使用 ggplot2 之前的那一行（我将一列 dat2 转换为一个因子）或 ggplotting 本身——不太确定是哪一行。我收到错误消息"Error in eval(expr, envir, enclos) : object 'dat2' not found".

我尝试按照here 的建议阅读this，但要么这是一个不同的问题，要么我就是不明白。当我逐行遍历它时，该功能工作正常。非常感谢任何帮助！

Answer 1

根据您的 cmets，错误几乎肯定在 geom_line(aes(color=dat2$item))。摆脱dat2$，它应该可以正常工作（即geom_line(aes(color=item))）。 aes 中的内容在data 参数（此处为dat2）中进行评估，全局环境作为外壳。值得注意的是，这意味着aes 不能使用函数环境中的内容，除非它是data 的一部分（此处为dat2）。由于dat2 不存在于dat2 中，并且dat2 不存在于全局环境中，因此您会收到该错误。