Hi,
Hope this helps
% Time series data
time = datetime({'1-Jan-18', '1-Feb-18', '1-Mar-18', '1-Apr-18', '1-May-18', '1-Sep-18', '1-Oct-18', '1-Nov-18', '1-Dec-18'});
GNSS = [2.5114, 2.5854, 2.5992, 2.6136, 2.6344, 2.6572, 2.6262, 2.6023, 2.5399];
SAS = [2.4904, 2.4996, 2.4970, 2.5007, 2.4989, 2.4979, 2.4976, 2.4979, 2.4967];
HOP = [2.5075, 2.5161, 2.5141, 2.5174, 2.5157, 2.5143, 2.5143, 2.5148, 2.5138];
% Calculate standard deviations
stds = std([GNSS; SAS; HOP], 0, 2);
% Calculate correlation coefficients
corrs = corrcoef([GNSS; SAS; HOP]');
% Create a Taylor plot
figure;
errorbar(stds, corrs(1, :), 'ko');
hold on;
scatter(stds(2:end), corrs(2:end, 1), 'filled', 'MarkerFaceColor', 'b');
% Set plot title and axis labels
title('Taylor Plot');
xlabel('Standard Deviation');
ylabel('Correlation Coefficient');
legend('GNSS', 'SAS', 'HOP');
grid on;
hold off;
