mycolors = [ 0.0 0.0 0.5;
                      0.0 0.0 1.0;
                      0.0 0.5 1.0; 
                      0.0 1.0 1.0; 
                      0.3 0.8 0.5; 
                      0.9 0.9 0.0; 
                      1.0 0.5 0.0; 
                      1.0 0.0 0.0;
                      0.5 0.0 0.0 ];

set(gca, 'ColorOrder', mycolors);
hold all

h = zeros(1,9);

 % Hysteria: The Telegraph
%   time vector is: time1363
%   data vector is: data1363

h(1) = plot(time1363(end-13:end),data1363(end-13:end), 'x-','LineWidth',2);

% Hysteria: The Guardian
%   time vector is: time1362
%   data vector is: data1362

h(2) = plot(time1362(end-13:end),data1362(end-13:end), 'o-','LineWidth',2);

% Hysteria: The Independent
%   time vector is: time1361
%   data vector is: data1361

h(3) = plot(time1361(end-13:end),data1361(end-13:end), 'x-','LineWidth',2);

% Hysteria: The Daily Star
%   time vector is: time1360
%   data vector is: data1360

h(4) = plot(time1360(end-13:end),data1360(end-13:end), 'o-','LineWidth',2);

% Hysteria: Hysteria: BBC News
%   time vector is: time1359
%   data vector is: data1359

h(5) = plot(time1359(end-13:end),data1359(end-13:end), 'o-','LineWidth',2);

% Hysteria: The Mirror
%   time vector is: time1358
%   data vector is: data1358

h(6) = plot(time1358(end-13:end),data1358(end-13:end), 'o-','LineWidth',2);

% Hysteria: The Sun
%   time vector is: time1357
%   data vector is: data1357

h(7) = plot(time1357(end-13:end),data1357(end-13:end), 'o-','LineWidth',2);

% Hysteria: The Daily Mail
%   time vector is: time1356
%   data vector is: data1356

h(8) = plot(time1356(end-13:end),data1356(end-13:end), 'o-','LineWidth',2);

% Hysteria: The Times
%   time vector is: time1427
%   data vector is: data1427

h(9) = plot(time1427(end-13:end),data1427(end-13:end), 'o-','LineWidth',2);




labels = {sprintf('Telegraph (av: %0.2f)', mean(data1363)), ...
sprintf('Guardian (av: %0.2f)', mean(data1362)), ...
sprintf('Independent (av: %0.2f)', mean(data1361)), ...
sprintf('Daily Star (av: %0.2f)', mean(data1360)), ...
sprintf('BBC News (av: %0.2f)', mean(data1359)), ...
sprintf('Mirror (av: %0.2f)', mean(data1358)), ...
sprintf('Sun (av: %0.2f)', mean(data1357)), ...
sprintf('Daily Mail (av: %0.2f)', mean(data1356)), ...
sprintf('Times (av: %0.2f)', mean(data1427)), ...
};

data = {data1363 data1362 data1361 data1360 data1359 data1358 data1357 data1356 data1427};

ymax = max(cellfun(@max, data));
colmeans = cellfun(@mean, data);
[~, I] = sort(colmeans, 'descend');

legend(h(I), labels(I), 'Location', 'NorthWest');

title('The Hysteriometer')
ylabel('''Hysteresis''')
xlabel('Time')

ylim([0 1.4*ymax])

xlim([time1356(end-13) time1356(end)])
datetick('x', 'd')