使用xyplot绘制大量时间序列

Question

这是我的数据类型的最小示例 争取策划：

这些曲线来自两个过程。

library("lattice")
x0<-matrix(NA,350,76)
for(i in 1:150)     x0[i,]<-arima.sim(list(order=c(1,0,0),ar=0.9),n=76)
for(i in 151:350)   x0[i,]<-arima.sim(list(order=c(1,0,0),ar=-0.9),n=76)

我想将它们绘制成由两个盒子组成的格子中的线条图。盒子位于上方 将包含前150条曲线（橙色），并显示下面的框  接下来的200条曲线（应为蓝色）。我不需要 标签或图例。我试图使用man-page上显示的示例：

aa<-t(x0)
colnames(aa)<-c(paste0("v0_",1:150),paste0("v1_",1:200))
aa<-as.ts(aa)
xyplot(aa,screens=list(v0_="0","1"),col=list(v0_="orange",v1_="blue"),auto.key=FALSE)

但不知怎的，它不起作用。

Answer 1

这样做没有其他因素（但是agstudy的解决方案不是像这样的黑客攻击）：

# This is equivalent to your for-loops, use whatever you prefer
x0 <- do.call(rbind, lapply(1:350, function(i) {
  arima.sim(list(order=c(1,0,0), ar=ifelse(i <= 150, 0.9, -0.9)), n=76)
}))

plotStuff <- function(indices, ...) {
  plot.new()
  plot.window(xlim=c(1, ncol(x0)), ylim=range(x0[indices,]))
  box()
  for (i in indices)
    lines(x0[i,], ...)
}

par(mfrow=c(2,1), mar=rep(1,4)) # two rows, reduced margin
plotStuff(1:150,   col="orange")
plotStuff(151:350, col="blue")

Answer 2

您应该将数据放在长格式中，如下所示：

      Var1   Var2     value group
1        1   v0_1 2.0696016    v0
2        2   v0_1 1.3954414    v0
      ..... ..........
26599   75 v1_200 0.3488131    v1
26600   76 v1_200 0.2957114    v1

例如使用reshape2：

library(reshape2)
aa.m <- melt(aa)
aa.m$group <- gsub('^(v[0-9])_(.*)','\\1',aa.m$Var2)
xyplot(value~Var1|group,data=aa.m,type='l',groups=group)