如何繪製時間序列的互相關矩陣？

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的這個函數的輸出嗎？基本上，這樣的事情（在我的情況下，4×4地塊的矩陣）：

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)