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

我的数据有一个timeseries表示，如下所示（没有行和列）注释：

      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矩阵图）：

---

一种简单的方法是在绘图设备上创建一个绘图矩阵，然后将每个

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])
    }
}

---

M试试这个
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)

pnl似乎还有另一种简单的方法

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

给我一个相关图矩阵。当然，如果需要协方差，还可以将其他选项传递给
acf
 应该发送给交叉验证吗？嗯。。我可以问一下为什么吗？这是一个与在有自己语言的
R
中绘图有关的问题。您是否正在寻找一种方法来在单个绘图上绘制每对列的每个
ccf
！确切地我更新了我的问题来展示一个例子。好的，我给出了一个非常基本的方法。但我会耐心等待更多的眼球，因为很可能时间序列包（
zoo
，
xts
）中有一个函数以更漂亮的方式执行类似操作。+1太棒了！非常感谢。我可能要等半天。如果没有别的事，我就接受这个。非常感谢。