如何绘制时间序列的互相关矩阵？

Question

我的数据的时间序列表示如下（没有行和列）注释：

      L1 L2 L3 L4
t=1    0  1  1  0
t=2    0  1  1  1
t=3    1  0  1  1
t=4    0  1  1  0

我将其读入 R 为：

timeseries = read.table("./test", header=F)

我正在使用互相关函数绘制 L1 的时间序列

ts.plot(timeseries$V1)

并将互相关函数绘制为：

ccf(timeseries$V1, timeseries$V2)

现在，有人可以告诉我吗如何绘制互相关矩阵来显示该函数 L1-L4 的输出？基本上是这样的（在我的例子中，是一个 4x4 的图矩阵）：

Answer 1

似乎还有另一种简单的方法可以做到这一点！

timeseries = read.table("./test", header=F)
acf(timeseries)

给我一个相关图矩阵。当然，如果需要协方差，还有其他选项可以传递给 acf。

Answer 2

执行此操作的一种简单方法是简单地在绘图设备上创建一个绘图矩阵，并将每个 ccf 绘图逐一放置：

M <- matrix(sample(0:1,40,replace = TRUE),nrow = 10)

par(mfrow= c(4,4))
for (i in 1:4){
    for (j in 1:4){
        ccf(M[,i],M[,j])
    }
}

但是如果您稍等一下，了解时间序列包的人就会发现更亲密地可能会使用一个可以更好地完成此操作的函数。

Answer 3

试试这个，其中 M 是 joran 的帖子：

pnl <- function(x, y = x) { par(new = TRUE); ccf(x, y) }
pairs(as.data.frame(M), upper.panel = pnl, diag.panel = pnl, cex.labels = 1)