Как передать эстетику ggplot2 из переменной?

У меня проблемы с передачей POSIXct, хранящегося в переменной как xmin / xmax для geom_rect. Я попытался построить отдельный пример, не упрощая то, что я пытаюсь сделать ...

Идея состоит в том, чтобы взять объект графика 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)