3

我遇到了一个奇怪的问题（我是 R 新手）。我尝试创建如下函数：

library(ggplot2)
median_confidence_interval <- function(x) {
  quart_list<-c()
  return_data<-data.frame(lower_ci=0,median=0,upper_ci=0)
  for(i in 1:1000){ 
    y<-x[as.integer(runif(length(x), min = 1, max = length(x) + 1))]
    median<-median(y)
    quart_list=c(quart_list,median)
  }
  return_data$median<-median(quart_list)
  return_data$lower_ci<-quantile(quart_list,probs=0.025)
  return_data$upper_ci<-quantile(quart_list,probs=0.975)
  p <- ggplot()
  p <- p + geom_density(aes(x=x)) + geom_density(aes(x=quart_list)) 
  p <- p + geom_vline(aes(xintercept = return_data$median, color='red'))
  p <- p + geom_vline(aes(xintercept = return_data$lower_ci, color='blue'))
  p <- p + geom_vline(aes(xintercept = return_data$upper_ci, color='green')) + coord_cartesian(xlim = c(min(x),max(x)))
  png("density_confidence_internal.png")
  plot(p)
  dev.off()
  return_data
}

在这段代码中，我只是想创建一个图并保存它。虽然我能够在函数外部独立执行这些语句，但不能在函数内部执行。该函数编译时没有错误，但在运行该函数时它显示“quart_list”未找到。如果quart_list和return_data存在于工作区中，那么我可以执行该函数并获得结果。当我清除工作区并执行该功能时，我在运行（而不是编译）时遇到了同样的错误。

另一个问题是，当我调用函数时median_confidence_interval(x)，它希望我只提供 'x' 作为参数，而不需要像median_confidence_interval(possum$earconch). 为什么会这样？

有人能指出我的方向吗？

4

1 回答 1

0

物体的环境评价ggplot有点神秘。但是，如果您记得ggplot希望data.frames传递给data参数并且aes应该是列中的值，data.frame您通常会避免出现问题。

为了调试这样的事情，我发现在函数中插入打印语句以理清我得到了多远是很有帮助的。（见注释行）

相应地调整你的功能给出：

median_confidence_interval <- function(x) {
  quart_list<-c()
  return_data<-data.frame(lower_ci=0,median=0,upper_ci=0)
  for(i in 1:1000){ 
    y<-x[as.integer(runif(length(x), min = 1, max = length(x) + 1))]
    median<-median(y)
    # print ('inside for loop')
    quart_list=c(quart_list,median)
  }

  # print('past for loop')
  return_data$median<-median(quart_list)
  return_data$lower_ci<-quantile(quart_list,probs=0.025)
  return_data$upper_ci<-quantile(quart_list,probs=0.975)

  # print('start of ggplot code')
  foo=data.frame(q=quart_list, x=x)
  p <- ggplot()
  p <- p + geom_density(data=foo, aes(x=x)) + geom_density(data=foo, aes(x=q)) 
  # print('past first quart_list reference in ggplot')
  p <- p + geom_vline(data=return_data, aes(xintercept = median, color='red'))
  p <- p + geom_vline(data=return_data, aes(xintercept = lower_ci, color='blue'))
  p <- p + geom_vline(data=return_data, aes(xintercept = upper_ci, color='green')) + coord_cartesian(xlim = c(min(x), max(x)))
  png("/tmp/density_confidence_internal.png")
  plot(p)
  dev.off()
  return_data
}

另外，我认为@DWin 在他的评论中有一个很好的观点！

于 2013-09-11T02:28:47.750 回答