我是ggplot2的新手并且很棒,但是我遇到了一件事情。
我绘制了一系列跨越一年的时间序列。 X轴派生自类Date
的变量。我已将该情节分面,因此我在具有独立y轴的列中有7个时间序列。该图的重点是比较每个面与顶面的相关性。
我想做的最后一件事是在每个方面的右上角添加文本(每个方面和第一个方面之间的估计的皮尔逊相关性)。
这被证明是非常困难的,因为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")
答案 0 :(得分:10)
这不是最新的代码,但我认为它有点像你所追求的:
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)