11

我是 ggplot2 的新手，它很棒，但我在一件事上遇到了困难。

我已经绘制了许多跨越一年的时间序列。X 轴是从 class 的变量派生的Date。我已经对情节进行了刻面，以便在具有独立 y 轴的列中有 7 个时间序列。该图的重点是比较每个刻面与顶部刻面的相关性。

我想做的最后一件事是将文本（每个方面与第一个方面之间的估计 pearson 相关性）添加到每个方面右上角的图中。

这被证明是非常困难的，因为geom_text()每一位文本都需要 x 和 y 坐标。当 X 轴是日期并且 Y 轴对于每个方面都不同时，如何指定坐标？这是我到目前为止的一些示例数据和代码，因此您可以重现我到目前为止的内容：

library(ggplot2)

date <- rep(as.Date(1:365,origin='2011-1-1'),7)
location <- factor(rep(1:7,365))
product <- rep(letters[1:7], each=365)
value <- c(sample(1:10, size=365, replace=T),sample(1:3, size=365, replace=T),
           sample(10:100, size=365, replace=T), sample(1:50, size=365, replace=T),
           sample(1:20, size=365, replace=T),sample(50:100, size=365, replace=T),
           sample(1:100, size=365, replace=T))
dat<-data.frame(date,location,product,value)

qplot(date, value, data=dat, geom="line", color=location, group=location, 
      main='Time Series Comparison', xlab='Month (2011)',ylab='Value') + 
        facet_grid(product ~ ., scale = "free_y")
4

1 回答 1

11

这不是最整洁的代码，但我认为它有点像你所追求的：

library(plyr)

corr_dat<-ddply(dat, .(product), summarise, value=value)
corr.df<-unstack(corr_dat, value~product)

corr_plot <- data.frame(date=max(dat$date),
                        label=paste0("rho==",round(cor(corr.df)[,1], 2)), 
                        ddply(dat, .(product), summarise, 
                          value=(min(value)+max(value))/2))

ggplot(dat, aes(x=date, y=value, color=location, group=location)) + 
  geom_line()+
  facet_grid(product ~ ., scale = "free_y")+
  geom_text(data=corr_plot, aes(x=date, y=value, label=label), 
            colour="black", inherit.aes=FALSE, parse=TRUE)

用表达式绘图

于 2012-07-24T07:52:17.997 回答