% Script:   GraphInset_CoV.m
% Author:   Anna Ziff
% Date:     February 15, 2020


% Import data
data = readtable('coronavirus.xlsx','Range','G3:P30');

data{24,9} = NaN;

deaths = data{:,9};
logdeaths = log(deaths);
report = data{:,2};
logreport = log(report);
date = data{:,1};

xlength = length(logreport(report~=24 & report>=8));

X = [ones(xlength,1), logreport(report~=24 & report>=8)];
Y = logdeaths(report~=24 & report>=8);

b = (X'*X)^(-1)*(X'*Y);

yhat = X*b;
yhatexp = exp(yhat);

xaxis = exp(X(:,2));

% Plot
fig = figure;
ax1 = axes(fig, 'Position', [0.1 0.1 0.85 0.85], 'XScale', 'log', 'YScale', 'log') ;
grid on;
hold on;
ax2 = axes(fig, 'Position', [0.18 0.63 0.3 0.3]) ;
grid on;
ylabel(ax1, 'Total Deaths');
xlabel(ax1, 'Time');
ylabel(ax2, 'log(Total Deaths)');
xlabel(ax2, 'Time');
hold on;
loglog(ax1, report,deaths, 'LineWidth', 1.5, 'Color', 'black');
hold on;
scatter(ax1, report, deaths, [], 'black');
hold on;
loglog(ax1, xaxis, yhatexp, 'LineWidth', 1.5, 'Color', 'red');
hold on;
plot(ax2, report, logdeaths, 'LineWidth', 1.5, 'Color', 'black');
hold on;
scatter(ax2, report, logdeaths,[],'black');

saveas(gcf, 'GraphInset_CoV.png');