K indicatorを計算してみる

この論文のK値を計算してみる
www.medrxiv.org
K値の定義

陽性者数は以下のデータを用いた
GitHub - kaz-ogiwara/covid19: 新型コロナウイルス感染症（COVID-19）の国内における状況を厚生労働省の報道発表資料からビジュアルにまとめた。
CSV
Total PCR Positive Cases in Japan (as of 2020.05.06) · GitHub

コード

clear; close;
A=readtable('PCR.csv');

subplot(211);
plot(A.date, A.positive);
title('Total PCR Positive');
xlabel('Date');
ylabel('Number');

subplot(212);
N=A.positive;
Na=N(8:end);   % current
Nb=N(1:end-7); % 7 days ago
K=1-Nb./Na;
t=A.date(8:end);
plot(t,K)
title('K-indicator');
xlabel('Date');
ylabel('K');

ちょうど丘のプラトーの辺りで緊急事態宣言(4/7)が出されてるようにも見えるけど、どうなんだろう。