如何从变量传递ggplot2美学?

时间:2010-09-17 17:37:38

标签: r ggplot2

我无法将存储在变量中的POSIXct作为geom_rect的xmin / xmax传递。我试图构建一个独立的例子,而不是轻视我正在尝试做的事情......

这个想法是采用ggplot2绘图对象,其x是POSIXt,并在特定范围内“放大”。缩放位于前80%,整个系列位于底部20%,并指示顶部放大的部分。

我的问题是我似乎无法将xmin / xmax传递给geom_rect - 我尝试过的每件事(除了手工组装绘图而不是函数)都会给我一个不同的错误。我尝试使用aes(),aes_string(),作为参数而不是美学传递,只传递字符串等。

以下示例告诉我:

Error in eval(expr, envir, enclos) : object 'lims' not found

我认为我的问题是,当美学得到处理时,我用来设置美学的变量不在范围内,但我无法弄清楚如何去做。帮助

library(ggplot2)

subplot <- function(x, y) viewport(layout.pos.col=x, layout.pos.row=y)
vplayout <- function(x, y) {
  grid.newpage()
  pushViewport(viewport(layout=grid.layout(y,x)))
}

anm_zoom <- function(limits, p) {

  lims <- as.POSIXct(limits)
  limlab <- paste(lims, collapse=" to ")

  top <- p + scale_x_datetime(limlab, limits=lims, expand=c(0,0))

  bottom <- p;
  bottom <- bottom + opts(title="")
  bottom <- bottom + opts(legend.position="none")
  bottom <- bottom + opts(axis.title.y=theme_blank())
  bottom <- bottom + scale_x_datetime("", expand=c(0,0))
  bottom <- bottom + geom_rect(aes(xmin=lims[1], xmax=lims[2]),
 ymin=-Inf, ymax=Inf, fill="grey80", alpha=0.01)

  ## Render the plots
  vplayout(1,5)
  print(top, vp=subplot(1,c(1,2,3,4)))
  print(bottom, vp=subplot(1,5))
}


pdate <- seq.POSIXt(from=as.POSIXct("2010-09-09 00:00"),
   to=as.POSIXct("2010-09-10 23:59"), by="2 mins")
var1 <- rnorm(length(pdate))
var2 <- rnorm(length(pdate))
df1 <- data.frame(pdate, var1, var2)

dm  <- melt(df1, id="pdate")

p <- ggplot(dm) + aes(x=pdate, y=value) + stat_summary(fun.y="sum", geom="line")

anm_zoom(c("2010-09-09 12:15", "2010-09-09 12:30"), p)

3 个答案:

答案 0 :(得分:8)

嗯,我认为你需要一个新的aes函数,有点像aes(因为它不会尝试解析它的参数),有点像aes_string(因为它在本地环境中立即评估其参数:)

aes_now <- function(...) {
  structure(list(...),  class = "uneval")
}

然后

bottom <- bottom + geom_rect(aes_now(xmin=lims[1], xmax=lims[2]),
 ymin=-Inf, ymax=Inf, fill="grey80", alpha=0.01)

给你你想要的东西。

答案 1 :(得分:1)

来自 hadley 的回答

重写

由于 hadley 更新了较新版本ggplot2, 处理函数ggplot()的非std评估的更直观的方法是使用aes_q()如下:

  xminName  <- substitute(lims[1]); xmaxName  <- substitute(lims[2])
  bottom <- bottom + 
    geom_rect(aes_q(xmin=xminName, xmax=xmaxName),
              ymin=-Inf, ymax=Inf, fill="grey80", alpha=0.01)

答案 2 :(得分:0)

你只需要更改函数参数限制的名称,因为我认为它正在创建一个范围冲突。我刚刚将它更改为limits1,并且还将代码中的第一行更改为lims = as.POSIXct(limits1)并且它完美地运行。看看!!