2010-02-16 90 views
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):对象 '县' 未找到

有人知道如何解决这一问题?

回答

3

我不知道是什么原因造成的错误,但我能想出是修复有这样一个数据帧来代替for循环:

date.df<-data.frame(d=dates,t=text) 
p <- p + geom_text(aes(x=d,label=t),y=0, 
        data = date.df, 
        size = 3, hjust = 1, vjust = 0) 
p<-p+geom_vline(aes(xintercept=d),data=date.df,alpha=.4) 
+0

因此,技巧是使用第二个数据框...谢谢! – diegovalle 2010-02-16 13:44:12

相关问题