Как построить матрицу взаимной корреляции для временных рядов?

У меня есть представление моих данных в таймсерии в виде следующих аннотаций (без строки и столбца):

      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):