% time1427 and data1427 available for use
% for Hysteria: The Times
% 
% time1357 and data1357 available for use
% for Hysteria: The Sun 
% 
% time1358 and data1358 available for use
% for Hysteria: The Mirror
% 
% time1361 and data1361 available for use
% for Hysteria: The Independent
% 
% time1360 and data1360 available for use
% for Hysteria: The Daily Star
% 
% time1356 and data1356 available for use
% for Hysteria: Daily Mail
% 
% time1359 and data1359 available for use
% for Hysteria: BBC News
% 
% Hysteria: The Telegraph
%   time vector is: time1363
%   data vector is: data1363
disp b

paddedData1427 = [nan(length(data1363)-length(data1427), 1); data1427];

disp a

myData = [data1363 data1362 data1361 data1360 data1359 data1358 data1357 data1356 paddedData1427];

disp c

b = sum(~isnan(myData'));
myData(isnan(myData'))=0;
mymean = sum(myData')./b;

[r, c]= size(time1363)
[r, c]= size(mymean)

plot(time1363,mymean, 'b-', 'LineWidth', 2);
hold all

%% Broadsheets

broadsheets = [paddedData1427 data1361 data1359 data1363 data1362]';

b = sum(~isnan(broadsheets))
broadsheets(isnan(broadsheets))=0
myBroadsheetMean = sum(broadsheets)./b

plot(time1363,myBroadsheetMean, ':', 'Color', [58 95 205]/256);

%% Tabloids

tabloids = [data1357 data1358 data1360 data1356]';

b = sum(~isnan(tabloids))
tabloids(isnan(tabloids))=0
myTabloidsMean = sum(tabloids)./b

plot(time1363,myTabloidsMean, 'r:');

plot(time1363, 5*ones(1, length(time1363)), 'r--', 'LineWidth', 2)

datetick

title('Average British Newspaper Hysteria')
legend({'Average Hysteria', 'Average Broadsheet Hysteria', 'Average Tabloid Hysteria', 'Mayan Armageddex Threshold'}, 'Location', 'NorthWest')

ylim([0 5.08])