delay = zeros(5);
for i = 1:4
    for j = (i+1):5
        [corr,lag] = xcorr(data(:,i),data(:,j));
        [max_value,max_index] = max(corr);
        delay(i,j) = -1*lag(max_index);
        delay(j,i) = lag(max_index);
    end
end