4
我正在使用ggplot2来探索不同军事行动对谋杀率的影响。为了显示效果,我在操作发生时绘制了一条垂直线,并在手术前后绘制了一条平滑的谋杀率线。如何从ggplot2中获取geom_vline和facet_wrap以在函数内部工作
我写了一个facet_wrap图来显示一堆县的情况。它的工作原理非常漂亮,但是当转换为函数时,使用局部变量绘制垂直线时会出现错误。
下面是一些示例代码:
drawTS <- function(df, dates, text) {
p <- ggplot(df, aes(date, murders)) +
facet_wrap(~ county, ncol = 1,
scale="free_y") +
scale_x_date() +
geom_smooth(aes(group = group), se = FALSE)
for(i in 1:length(dates)) {
#If it's not a global variable I get an object not found error
temp[i] <<- dates[i]
p <- p + geom_text(aes(x,y), label = text[i],
data = data.frame(x = dates[i], y = -10),
size = 3, hjust = 1, vjust = 0) +
#Here's the problem
geom_vline(xintercept=temp[i], alpha=.4)
}
p
}
library(ggplot2)
df <- data.frame(date = rep(seq(as.Date("2007/1/01"),
length=36, by='1 month'),4),
murders = round(runif(36*4) * 100),
county = rep(rep(factor(1:4),9),each=4),
group = rep(c(rep(1,6), rep(2,12),rep(3,18))), each=4)
dates <- c(as.Date("2007/6/15"), as.Date("2008/6/15"))
temp <- c()
drawTS(df, dates, c("Op 1","Op 2"))
有与全局变量没有错误,但它看起来丑陋。
如果不是的temp[i]
可变我用dates[i]
内geom_vline()
,我得到这样的:在NextMethod( “[”)
错误:对象 '我' 未找到
如果我将变量dates[i]
包装在aes()
中,我得到:
错误的eval(表达式,ENVIR,enclos):对象 '县' 未找到
有人知道如何解决这一问题?
因此,技巧是使用第二个数据框...谢谢! – diegovalle 2010-02-16 13:44:12