将多个时间序列数据（长格式）绘制到一个图中？

Question

假设我有 3 个时间序列数据“a”，“b”，“c”，它们每个都有 2 个变量要在 7 天内记录。

示例代码如下：

require(data.table)  
#Create data
DT<-data.table(type=c(rep("a",7),rep("b",7),rep("c",7)), time=1:7,V1=rnorm(7*3,10,5),V2=rnorm(7*3,100,20),key="type")
# plot.zoo
require(zoo)
plot(zoo(DT["a"])[,3:4])

我可以一次绘制一种类型（如上），但我想将所有“a”“b”“c”...绘制到该面板中，其中不同的颜色代表不同类型的时间序列。

所以我正在寻找的是一个有两个面板（用于“V1”和“V2”）的情节，并且在每个面板中，有几行（“a”、“b”、“c”.. .) 用不同的颜色

Answer 1

试试这个动物园经典图形、动物园格子图形和动物园 ggplot2 图形：

library(zoo)
z <- read.zoo(DT, split = 1, index = 2, FUN = identity)
Names <- read.table(text = names(z), sep = ".", col.names = c("screen", "col"))

plot(z, screen = Names$screen, col = Names$col) # also see note 3 below

library(lattice)
xyplot(z, screen = Names$screen, col = Names$col) # also see note 3 below


library(ggplot2) # also see note 4 below
m <- fortify(z, melt = TRUE)
m2 <- transform(m, screen = sub("\\..*", "", Series), col = sub(".*\\.", "", Series))
qplot(Index, Value, data = m2, col = col, geom = "line") + facet_wrap(~ screen)

注意事项

1) 如果我们只是想要单独的面板，那就是 plot(z)、xyplot(z) 和 autoplot(z)。

2) names(z) 和 Names 是：

> names(z)
[1] "V1.a" "V2.a" "V1.b" "V2.b" "V1.c" "V2.c"
> Names
  screen col
1     V1   a
2     V2   a
3     V1   b
4     V2   b
5     V1   c
6     V2   c

3) 我们可以硬编码成这样（在这种情况下我们不需要Names）：

plot(z, screen = 1:2, col = c(1, 1, 2, 2, 3, 3))

xyplot(z, screen = 1:2, col = c(1, 1, 2, 2, 3, 3))

4) 这是不涉及 z 的 ggplot2 替代方案。它使用 data.table 包将DT 转换为长格式，然后执行qplot：

long <- DT[, list(Series = names(.SD), Value = unlist(.SD)), by = list(type, time)]
qplot(time, Value, data = long, col = type, geom = "line") + facet_wrap(~ Series)

Answer 2

以lattice 包为例，你可以这样做：

library(lattice)
xyplot(V1+V2~time,groups=type,data=DT,type='l')

您也可以使用ggplot2：

library(reshape2)
dt.m <- melt(DT,measure.vars=c('V1','V2'))

ggplot(dt.m) +
 geom_line(aes(x=time,y=value,group=type,color=type)) +
  facet_grid(~variable)