如何在晶格中给定条件下绘制叠加多个时间序列？

Question

假设我有一个数据框 df，如下所示：

f t1 t2 t3
h 1  3  4
h 2  4  3
t 3  4  5
t 5  6  8

f 是一个因子，$t 属性是与时间排序事件相关的数值。 我可以使用 par(new=T) 覆盖时间序列 t1 到 t3 并手动按因子隔离。 但我想知道是否有某种方法可以用格子来做到这一点，其中重叠的时间序列 受因素制约。因此，我们将有两个面板，其中重叠的时间序列对应于条件因子 f。我见过的大多数例子每个因素只使用一个时间序列（向量）。我也考虑过使用平行图，但时间信息丢失了。 我也尝试过类似的方法

xyplot(df$t1+df$t2+df$t3 ~seq(3) | factor(df$f))

，但它丢失了行序列连接。有人知道这是否可能吗？

这是使用非格方法的非常粗略的说明。

x<-matrix(seq(12),4,3)
f<-c('a','a','b','b')
df<-data.frame(f,x)

layout(1:2); yr<-c(0,12); xr<-c(1,3);
plot(as.numeric(df[1,2:4])~seq(3),type='o',ylim=yr,xlim=xr,ylab='A')
par(new=T)
plot(as.numeric(df[2,2:4])~seq(3),type='o',ylim=yr,xlim=xr,ylab='A')

plot(as.numeric(df[3,2:4])~seq(3), type='o',ylim=yr,xlim=xr,ylab='B')
par(new=T)
plot(as.numeric(df[4,2:4])~seq(3),type='o',ylim=yr,xlim=xr,ylab='B')

Answer 1

我添加了一个 ID 变量并与 package:reshape2 融为一体，

 dat
  f t1 t2 t3 ID
1 h  1  3  4  1
2 h  2  4  3  2
3 t  3  4  5  3
4 t  5  6  8  4
datm <- melt(dat, id.vars=c("ID","f"), measure.vars=c("t1", "t2", "t3"))
> datm
   ID f variable value
1   1 h       t1     1
2   2 h       t1     2
3   3 t       t1     3
4   4 t       t1     5
5   1 h       t2     3
6   2 h       t2     4
7   3 t       t2     4
8   4 t       t2     6
9   1 h       t3     4
10  2 h       t3     3
11  3 t       t3     5
12  4 t       t3     8

因为您要求将其“覆盖”，所以我使用组参数来保持 ID 的独立和“|”运算符为您提供“h”和“t”的两个面板：

xyplot(value~variable|f, group=ID, data=datm, type="b")

Answer 2

(1) 这可以使用 xyplot.zoo 紧凑地完成。第一个语句将数据框转换为动物园系列（系列存储在动物园对象的列中），第二个语句对其进行绘制，以便屏幕参数定义每个系列显示在哪个面板中：

library(zoo)
library(lattice)

z <- zoo(t(df[-1]))

xyplot(z, screen = df$f, type = "o")

（2）或者如果需要在 X 轴上显示 df 的列名称，然后将 z 定义如下（然后发出上面的 xyplot 命令）

z <- zoo(t(df[-1])), factor(names(df[-1])))

： >xyplot 使用第一个点中的 z 看起来像这样（第二个点除了 X 轴标签之外是相同的）：

编辑：简化 (2)