計算幾個時間序列

Question

之間的關係。如果我有一個矩陣：

data = rand(365,5); 

什麼是計算每列，其餘列的平均值之間的相關性的最合適的方式。例如，對於第一列：

R = nonzeros(tril(corrcoef(data(:,1),mean(data(:,2:end)')'),-1)); 

我怎麼能重複這個過程，使我有5個相關值，即每個系列？

編輯：

感謝您的意見。這也可以在一行中完成：

R = arrayfun(@(x)nonzeros(tril(corrcoef(data(:,x),... 
    mean(data(:,setdiff(1:size(data,2),x))')'),-1)),1:size(data,2)); 

對於那些希望避免循環的人。雖然在這種情況下，由於其可讀性，下面所示的方法更好。

Answer 1

for i=1:5 
    x = data(:,i); 
    y = mean(data(:,(1:5) ~= i)')'; 
    R(i) = nonzeros(tril(corrcoef(x,y),-1)); 
end 

Answer 2

稍微簡化版本：

R = zeros(1,5); 
for i=1:5 
    x = data(:,i); 
    y = mean(data(:,(1:5)~=i), 2); 
    R(i) = corr(x,y); 
end