在ggplot2中具有平均轮廓的刻面时间序列

Question

使用以下模拟时间序列：

n=70
m1 = matrix(rnorm(n), ncol=7)
m2 = matrix(rnorm(n, 0,4), ncol=7)
d = data.frame(rbind(m1,m2), cl=rep(c(1,2), each=5))

（前7列代表时间点，最后一列代表该类）

是否可以使用ggplot2构建包含每个图中平均曲线的分面时间序列？

结果应如下所示：

Answer 1

它可能不是最漂亮的代码，但我相信它可以让你找到你想要的东西，

n=70
m1 = matrix(rnorm(n), ncol=7)
m2 = matrix(rnorm(n, 0,4), ncol=7)
d = data.frame(rbind(m1,m2), cl=rep(c(1,2), each=5))

d <- cbind(paste("d", 1:NROW(d), sep = ""), d)
names(d)[1] <- "id.var"

library(reshape)
longDF <- melt(d, id=c("cl", "id.var"))
library(ggplot2)

p <- ggplot(data = longDF, aes(x = variable, y = value, group = id.var))
p + geom_line() + stat_smooth(aes(group = 1), method = "lm", 
se = FALSE, colour="red") + facet_grid(cl ~ .)

请不要犹豫，改进我的代码。