;
set(gca,'FontSize',fontSize,'fontWeight','normal')
disp('z: tau = ' + string(tau))


%% Time delay and plot x and y variables
tau = mdDelay(data(:,1:2), 'maxLag', 25, 'plottype', 'all');
set(gca,'FontSize',fontSize,'fontWeight','normal')
disp('xy: tau = ' + string(tau))


%% Time delay and plot x and z variables
tau = mdDelay(data(:,[1,3]), 'maxLag', 25, 'plottype', 'all');
set(gca,'FontSize',fontSize,'fontWeight','normal')
disp('xz: tau = ' + string(tau))


%% Time delay and plot y and z variables
tau = 