新型コロナ・ウィルスの陽性者数と死者数をプロットしてみる

新型コロナ・ウィルスの陽性者数と死者数をプロットしてみます。

陽性者数と死者数のデータは東洋経済のところから貰ってきました。データだけダウンロードする方法が分からなかったので、GitHub から一式落としてきました。

陽性者数は日ごとの値なのに、死者数は累計値になっています。なにか事情があるのかしら…。

MATLABでは次のようなコードでプロットできます。陽性者数は累計を計算し、日ごとの死者数は累計を微分しています。

death = readtable("death_total.csv", 'ReadVariableNames', false);
cases = readtable("pcr_positive_daily.csv", 'ReadVariableNames', false);

death.Properties.VariableNames{1} = 'date';
death.Properties.VariableNames{2} = 'death';

cases.Properties.VariableNames{1} = 'date';
cases.Properties.VariableNames{2} = 'cases';

figure(1)
yyaxis left
plot(cases.date, cumsum(cases.cases)/1e3, '-b', 'LineWidth', 2)
ylabel("PCR 陽性: 累計 (千人)")
hold on
yyaxis right
plot(death.date, death.death /1e3, '-r', 'LineWidth', 2)
hold off
legend("PCR 陽性(累計)", "死者数(累計)", 'Location', 'northwest')
ylabel("死者数: 累計 (千人)")

ax = gca;
ax.FontSize = 14;
ax.FontName = 'Yu Gothic';

figure(2)
yyaxis left
plot(cases.date, cases.cases, '-b', 'LineWidth', 2)
ylim([0 2500])
hold on
ylabel("PCR 陽性者数/日")
yyaxis right
plot(death.date(2:end), diff(death.death), '-r', 'LineWidth', 2)
ylim([0 150])
hold off
legend("PCR 陽性者数/日", "死者数/日", 'Location', 'north')
ylabel("死者数/日")

ax = gca;
ax.FontSize = 14;
ax.FontName = 'Yu Gothic';