У меня проблемы с передачей 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)