感染ダイナミックス（2020/5/17)

このサイトはコロナ禍が始まって最初にまとめました。ただし、図1と図2は随時更新しています。また感染制御についてもまとめようとしています。

1 はじめに
新型コロナウィルスに関する検討結果を公開しています。「なぜ外出自粛が必要かを説得力をもって説明することは難しい」という問題意識から検討を始めました（ツイート_R020330、ツイート_R020402）。また「８割削減」に示される数値の根拠に感心が集まっています。筆者の専門は感染学ではなく、システム制御工学ですが、一定の説明はできると思っています。

ポイントは「感染ダイナミックス」への理解です。ここでは次の文献を参照します。
[1] 稲葉寿編著：「感染症の数理モデル」,培風館（2008）
[2] 牧野淳一郎：「3.11以後の科学リテラシーno.89」,科学, vol.90, no.5, 岩波書店(2020)
なお、感染ダイナミックスの数理モデルは、次の文献が先駆けとなっています。
Kermack, W. O. and McKendrick, A. G.（1927）. Contributions to the mathematical theory of epidemics – I, Proceedings of the Royal Society Series A, 115, 700–721

2 SIRモデル
2.1 基本再生産数の概念
感染ダイナミックスの数理モデルのうち最も基本的なものはSIRモデルとして知られています。ここでは文献[1]に従って、その概要を述べます。

いま当該地域の人口 $N$ は次のような構成であるとします。

$(1)\quad S(t)+I(t)+R(t)=N$

ただし、時刻 $t$ において

● $S(t)$ ：感受性(susceptible)状態の人口
● $I(t)$ ：感染性(infectious)状態の人口
● $R(t)$ ：回復・隔離または免疫状態(recovered/removed/immune)状態の人口

以下では $S$ 、 $I$ 、 $R$ をそれぞれ未感染者数、感染者数、回復者数と呼びます。いくつかの仮定のもと、感染ダイナミックスの支配方程式として、次のSIRモデルを考えます。

$(2)\quad\left\{\begin{array}{l} \frac{dS(t)}{dt}=-\beta S(t)I(t) \\ \frac{dI(t)}{dt}=-\gamma I(t)+\beta S(t)I(t)\\ \frac{dR(t)}{dt}=\gamma I(t) \end{array}\right.$

ここで、 $\beta$ は感染率、 $\gamma$ は回復率と呼ばれています。

いま時刻 $t$ における感染状況を３項組 $(S(t),I(t),R(t))$ で表します。感染前の自明な平衡状態は

$(3)\quad (S,I,R)=(N,0,0)$

ですが、これが初期時刻 $t=0$ において、感染者がごく少数 $I_0$ だけ現れ

$(4)\quad (S_0,I_0,R_0)=(N-I_0,I_0,0)$

となったとします。このときの関心事は、感染流行が始まるか、それとも感染収束するかどうかです。これは(2)の第２式において、感染初期では $S(t)\simeq S_0$ として

$(5)\quad \frac{dI(t)}{dt}=\underbrace{(\beta S_0-\gamma)}_{\lambda_0}I(t), \quad I(0)=I_0$

の解

$(6)\quad I(t)=\exp(\lambda_0 t)I_0$

の振舞いを調べることになります。明らかに、初期成長率 $\lambda_0$ の正負によって、感染流行か感染収束かが決まります。これを

$(7)\quad \lambda_0=\beta S_0-\gamma=\gamma(\frac{\beta S_0}{\gamma}-1)$

と変形して、ここで現れた定数

$(8)\quad R_0=\frac{\beta S_0}{\gamma}$

を基本再生産数と名付け

● $R_0>1$ のとき感染流行
● $R_0<1$ のとき感染収束

となると判断します。この方法は閾値原理と呼ばれています。

さて、基本再生産数は次のように表されます。

$(9)\quad R_0=\beta S_0\underbrace{\frac{1}{\gamma}}_{T}=1+\underbrace{(\beta S_0-\gamma)}_{\lambda_0}\underbrace{\frac{1}{\gamma}}_{T}$

ここで、 $T$ は感染性時間（感染力をもつ時間）を表すと解釈できます。実際、感染力がなくなったと仮定すると( $\beta=0$ )、(2)の第２式は

$(10)\quad\frac{dI(t)}{dt}=-\gamma I(t), \quad I(0)=I_0$

となり、この解は

$(11)\quad I(t)=\exp(-\gamma t)I_0\rightarrow0\quad(t\rightarrow\infty)$

と表され、 $T=1/\gamma$ の5～10倍の時間が経てば、感染者数 $I_0$ は十分零になるからです。

そこで(9)において、 $\beta S_0$ は１人が単位時間に２次感染させる人数、これに感染性時間 $T$ をかけるので、基本再生産数は感染者１人が２次感染させる人数を表すと言えます。

その算出には初期成長率 $\lambda_0$ が必要ですが、次式を用います。

$(12)\quad \frac{I(t+\Delta t)}{I(t)}=\frac{\exp(\lambda_0 (t+\Delta t))I_0}{\exp(\lambda_0 t)I_0}=\exp(\lambda_0 \Delta t) \Rightarrow \lambda_0=\ln(\frac{I(t+\Delta t)}{I(t)})/\Delta t$

2.2 基本再生産数の推定例
専門家による基本再生産数の推定法は、感染学上の知見に基づく統計推測手法によるものと思われますが、詳細は明らかではありません。そこで、次のような基本再生産数の推定手法を考えてみます。

基本再生産数 $R_0$ の推定手法１：
$1^\circ$ 当該地域の感染者数データの移動平均をとる
$2^\circ$ 次式を用いて時刻 $t$ における（初期）成長率 $\lambda_0$ を計算する

$(13)\quad\lambda_0(t)=\ln(\frac{I(t)}{I(t-\Delta t)})/\Delta t$

$3^\circ$ 次式を用いて時刻 $t$ における基本再生産数 $R_0$ を計算する

$(14)\quad R_0(t)=1+\lambda_0(t) T$

図1に全国、東京都、福岡県、NYの感染数を7日間移動平均したデータを示します。図2に日ごとに基本再生産数 $R_0$ を推定した例を示します。ただし、 $2^\circ$ で $\Delta t=7$ としています。 $3^\circ$ で回復率は $\gamma=0.04,0.08,0.12,0.16$ の４通りを試し、 $R_0<0$ となる場合は0としています。

図1 感染者数の実データ（移動平均後）
国内データソース：NHK 特設サイト新型コロナウイルス
NYデータソ－ス：New York State Department of Health COVID-19 Tracker

図2 基本再生産数 $R_0$ の推定例

図2の妥当性を検討するために、「新型コロナウイルス感染症対策の状況分析・提言」（2020年5月1日）に記載されている次の公表値に注目します。

●全国の実効再生産数：　 $R_0(3/25)=2$ , 　 $R_0(4/10)=0.7$
●東京都の実効再生産数： $R_0(3/14)=2.6$ , $R_0(4/10)=0.5$

ここで基本再生産数の代わりに実効再⽣産数という術語が用いられています。ここでは両者の理論上の違いまでには踏み込まず、実効再⽣産数は基本再生産数 $R_0$ に実質的に等しいとして話を進めます。上記の値は上記資料に掲載されている次の図版から読み取ることができます。

図3 専門家会議資料における実効再生産数の推定

図3の左側は全国について、右側は東京都について、感染者数と実効再生産数が示されています。下側に示されている棒グラフは上側の感染者数の移動平均をとったものと思われ、これに基づいて実効再生産数が計算されています。ただし、移動平均後のデータは２週間ほど前倒しされています（ピーク値と $R_0$ が１を切る時点が少しずれています）。実際、全国（東京都）の感染者数のピークは4/11（4/9）ですが、移動平均後は3/28（3/27）にピークが現れています。図2において全国の基本再生産数が0.7まで下がるのは早くて4/20当たり、東京都の基本再生産数が0.5まで下がるのは早くて4/23当たりですから、ほぼ２週間遅れとなっています。また、回復率 $\gamma$ の値に依らず、感染者数がピークを取る時点で $R_0$ が１を切り、収束に向かっている傾向が捉えられています。

問題はどの回復率 $\gamma$ が妥当かですが、文献でよく見るのは $\gamma=0.04$ です。しかしこれでは直近で $R_0$ が負となりました（ $\gamma=0.12$ , $0.16$ では常に正でした）。回復率の逆数は感染性時間 $T=1/\gamma$ ですから、 $\gamma=0.04$ , $0.08$ , $0.12$ , $0.16$ に対応して、 $T=25$ , $12.5$ , $8.3$ , $6.25$ です。指数関数的に減衰する場合は、初期値の10%まで減衰する期間は $2.3T$ 、2%まで減衰する期間は $3.9T$ ですから、たとえば、 $\gamma=0.16$ , $T=6.25$ でもそうおかしくはないかと思います。

3 感染シミュレーション
3.1 無次元化SIRモデル
前章では、SIRモデルのうち第２式のみを用いて議論しました。ここでは３本の微分法方程式を連立させて、 $S$ 、 $I$ 、 $R$ の時系列のシミュレーションを行います。

$S$ 、 $I$ 、 $R$ の間には(1)の拘束があるので、実質的には２本の微分方程式を連立させて解きます。(2)の第２式と(1)からSを消去して、第３式を合わせると次式を得ます（第１式はこれらから導出できます）。

$(15)\quad\left\{\begin{array}{l} \frac{dI(t)}{dt}=-\gamma I(t)+\beta \underbrace{(N-I(t)-R(t))}_{S(t)}I(t)\\ \frac{dR(t)}{dt}=\gamma I(t) \end{array}\right.$

ここで３つのパラメータ $N$ 、 $\beta$ 、 $\gamma$ がありますが、無次元化という操作を行って、１つのパラメータにまとめます。これは文献[2]を参考にしています。

まず、人口 $N$ への依存を除くために、各変数の割合（人口比）、すなわち、感染者数の割合、回復者数の割合、未感染者数の割合を、それぞれ

$(16)\quad i(t)=I(t)/N, r(t)=R(t)/N, s(t)=S(t)/N$

のように小文字で表します。(1)より次式が成り立ちます。

$(17)\quad \underbrace{\frac{I(t)}{N}}_{i(t)}+\underbrace{\frac{R(t)}{N}}_{r(t)}+\underbrace{\frac{S(t)}{N}}_{s(t)}=1$

これを用いて、(15)から次式を得ます。

$(18)\quad\left\{\begin{array}{l} \underbrace{\frac{d\frac{I(t)}{N}}{dt}}_{\frac{di(t)}{dt}}=-\gamma \underbrace{\frac{I(t)}{N}}_{i(t)}+\beta N(\underbrace{\frac{N}{N}}_{1}-\underbrace{\frac{I(t)}{N}}_{i(t)}-\underbrace{\frac{R(t)}{N}}_{r(t)})\underbrace{\frac{I}{N}}_{i(t)}\\ \underbrace{\frac{d\frac{R(t)}{N}}{dt}}_{\frac{dr(t)}{dt}}=\gamma \underbrace{\frac{I(t)}{N}}_{i(t)} \end{array}\right.$

次に、感染性時間 $T=1/\gamma$ が単位時間となるように、次のように変形します。

$(19)\quad\left\{\begin{array}{l} \underbrace{\frac{di(t)}{d\frac{t}{T}}}_{\frac{di(t)}{d\tau}}=-i(t)+\underbrace{\frac{\beta N}{\gamma}}_{\bar{R}_0} \underbrace{(1-i(t)-r(t))}_{s(t)} i(t)\\ \underbrace{\frac{dr(t)}{d\frac{t}{T}}}_{\frac{dr(t)}{d\tau}}=i(t) \end{array}\right.$

すなわち、 $\tau=t/T$ とおいて、次式を得ます。

$(20)\quad\left\{\begin{array}{l} {\frac{di(\tau)}{d\tau}}=-i(\tau)+\bar{R}_0 s(\tau)i(\tau)=\underbrace{(-1+\bar{R}_0 s(\tau))}_{\lambda(\tau)}i(\tau)\\ {\frac{dr(\tau)}{d\tau}}=i(\tau) \end{array}\right.$

ここで、次の定数が定義されていることに注意してください。

$(21)\quad \bar{R}_0=\frac{\beta N}{\gamma}>\frac{\beta S_0}{\gamma}$

すなわち、 $\bar{R}_0$ は基本再生産数 $R_0$ の上界を与えています。ただし、感染初期の実データでは $S_0\simeq N$ が成り立ち、 $\bar{R}_0\simeq R_0$ です。以下では、 $R_0$ は $\bar{R}_0$ の意味で使用しています。

注1：(20)における $i(\tau)$ , $s(\tau)$ , $r(\tau)$ , $\lambda(t)$ はそれぞれ $i(\tau T)$ , $s(\tau T)$ , $r(\tau T)$ , $\lambda(\tau T)$ として、離散時間の差分方程式として表すべきかもしれませんが、以下では上式を用います。

3-2 感染シミュレーション
SIRモデル(20)の解の振舞いは、数値計算を行って容易にチェックできます。図3は $i(0)=0.001$ のとき、基本再生産数 $R_0$ の値によって感染ダイナミックスがどのように変わるかを示しています。青線が未感染者数 $S$ 、赤線が感染者数 $I$ 、緑線が回復者数（＝累積感染者数） $R$ の割合です。これから、基本再生産数 $R_0$ の値によって、最終的にどれくらいの割合で感染者が出るか、どれくらいの収束期間がかかるかを読み取ることができます。ただし、どの図も $R_0$ が一定の場合であることに注意してください。これらの図から、もし $R_0$ の値を順次下げていくことができれば、最終的な感染者数の割合は小さくなることが分かります。

図4 基本再生産数 $R_0$ に依存する感染ダイナミックス

注2：図4において、時間軸の単位は感染性時間 $T=1/\gamma$ です。したがって、収束期間を見積もるためには、回復率 $\gamma$ の推定が必要になります。

図5は、いま基本再生産数は $R_0=2.5$ であるとし、ある時点からこれを $(1-q)R_0$ に低減したとき、感染者数がどのように変化するかを示しています( $i(0)=0.001$ の場合です)。 $q=0.8$ のときが8割低減、 $q=0.2$ のときが2割低減です。

図5 基本再生産数 $R_0=2.5$ を $(1-q)R_0$ に低減したときの効果

注3：「８割削減の妥当性」についてですが、これはツイート_R020405が基礎となっているようです。ここでは基本再生産数 $R_0=2.5$ の場合であることが表明されていますので図5が参考になります。これは「基本再生産数を８割減することが必要で、２割減では効果がない」という主張をほぼ裏付けていると言えます。

2-2 感染初期における予測
実データによる感染予測ではSIRモデル(20)を解くわけですが、その際 $R_0$ を推定する必要があります。そのために次の問題設定を考えます。

基本再生産数 $R_0$ の推定手法２：
いま人口を $N$ 、感染者数の時系列を $I(t_1),\cdots,I(t_n)$ とします。またパラメータ $R_0$ を含むSIRモデル(20)の解の $t_1,\cdots,t_n$ における値を $i(t_1),\cdots,i(t_n)$ とします。このとき評価関数

$(22)\quad J(R_0)=(I(t_1)/N-i(t_1))^2+\cdots+(I(t_n)/N-i(t_n))^2$

を最小化するようなパラメータ $R_0$ の値 $R_0^*$ を、数値最適化の手法を用いて求めます。

この手法を用いて、全国の1/24～5/5の新規感染者数の一部のデータ（1/24～4/7、赤色〇）をもとに、基本再生産数を推定したところ $R_0^*$ =1.08となりました。図6は、感染者数の実績値（〇）と、SIRモデルによる計算値（青色実線）を重ねてプロットしたものです。感染初期の感染者数の指数関数的増加傾向がよく捉えられています。
図6 日本における感染者数の推移

首相は4/7日夜の記者会見で『このペースで感染拡大が続けば（感染者が）２週間後には１万人、１カ月後には８万人を越える』と指摘しました。図7は図6の感染者数の累積値を示しています。起点となる日を4/7（75日目）としますと、２週間後（88日目）より早く１万人を超えています（実際に１万人を超えたのは4/18（86日目）でした）。また８万人に到達するのは112日目あたりで、１か月後(105日目)より少し遅くなっています。
図7 日本における累積感染者数の推移

注4：上の問題を解くとき、少数の感染者が出た日を第１日として、 $t_1=1$ 、以後 $t_2=2,t_3=3\cdots$ とします。一方、SIRモデル(20)の解も $t_1=1,t_2=2,\cdots$ において求めます。この方法で感染初期の指数関数的増大の様子をよく捉えることができます。一方、4/11放映のNHKスペシャルより、基本再生産数 $R_0$ は大体1.5程度であることが表明されています。上で求めた $R_0^*$ の値1.08は公表値1.5より低めに出ています。そこで、日ごとの実データに回帰をとる場合 $T=1$ と考えられるので、(9)より $初期成長率\lambda_0$ を求めると

$(23)\quad \underbrace{1.08}_{R_0^*}=1+\lambda_0\times\underbrace{1}_{T}\quad\Rightarrow\quad \lambda_0=0.08}$

となります。これをたとえば感染性時間 $T=1/\gamma=1/0.16=6.25$ を用いて

$(24)\quad R_0=1+\underbrace{0.08}_{\lambda_0}\times\underbrace{6.25}_{T}=1.5}$

のように修正すると、公表値と一致します。

5 おわりに
「なぜ外出自粛は必要か」の問いには、感染爆発を防ぐためには基本再生産数を小さく抑え込むことが必要で、そのための間接的手段として外出自粛やロックダウンが考えられるからと答えられます。感染ダイナミックスの簡易モデル（SIRモデル）は発散から収束に切替る振舞いをうまくモデル化しています。ただ、基本再生産数というたった１つのパラメータで、すべての状況を反映させようとしているところには自ずと限界があります。それでも感染ダイナミクスの本質を捉えていることは間違いないと思います。

「８割削減の妥当性」についてですが、注3で述べたように、基本再生産数 $R_0=2.5$ の場合、図5から「基本再生産数を８割削減することが必要で、２割削減では効果がない」という主張は理解できます。ところが、4/11放映のNHKスペシャルより、基本再生産数 $R_0$ は大体1.5程度であることが表明されています。この場合、図5は違ったものになります。また、5/1の専門家会議では図3のように実効再生産数は4/10の時点で1を切っていることが表明されています。いつの時点で、どのくらいの基本再生産数であるかは、感染対策を議論する上で極めて重要です。もし4/7の時点で基本再生産数が１を切っているのなら、非常事態宣言は医療崩壊を防ぐ意味合いが強かったのかもしれません。

感染流行を制御する決め手は基本再生産数 $R_0$ をどう調整していくかにかかっています。これをSIRモデル(18)の非線形最適制御問題または線形パラメータ変動制御問題として定式化し、その解として望ましい $R_0$ の時系列を求めることが考えられます。もちろん、これが得られたとしても、どのような行動変容に結び付ければよいかは不明です。ただ、厳しい外出自粛をいつまでも続けるわけにはいきません。現在収束過程にあると思われる感染者数データから $R_0$ の時系列を求め、比較検討しておくことは、第２波への対応や出口戦略の検討に役立つと思っております。年内には検討結果をご報告したいと思います。

補遺
5/5に小田垣先生が新型コロナウイルスの蔓延に関する一考察を発表されました。現在指摘されている「PCR検査不足」の主張の妥当性を議論する際のベースを与えることのできるきわめて興味深い論文です。以下では、本稿の枠組みの中でフォローアップさせていただきたいと思います。

当該論文では、感染性人口 $I(t)$ を

● $H(t)$ ：未検査感染性人口（市中感染者数）　⇒　筆者の方で記号をIからHに変更
● $Q(t)$ ：検査後隔離状態にある感染性人口（隔離感染者数）

に分けて、次のSIQRモデルが提案されています。

$(25)\quad\left\{\begin{array}{l} \frac{dS(t)}{dt}=-\beta S(t)H(t) \\ \frac{dH(t)}{dt}=-\gamma H(t)-pH(t)+(1-q)\beta S(t)H(t)\\ \frac{dQ(t)}{dt}=-\gamma Q(t)+pH(t)+q\beta S(t)H(t)\\ \frac{dR(t)}{dt}=\gamma (H(t)+Q(t)) \end{array}\right.$

これをSIRモデル

$(2)\quad\left\{\begin{array}{l} \frac{dS(t)}{dt}=-\beta S(t)I(t) \\ \frac{dI(t)}{dt}=-\gamma I(t)+\beta S(t)I(t)\\ \frac{dR(t)}{dt}=\gamma I(t) \end{array}\right.$

と比較します。それぞれの第 $i$ 式を、SIRi、SIQRiと参照します。

●SIQR4はSIR3に $I=H+Q$ を代入しただけです。
●SIQR1はSIR1の $I=H+Q$ のうち隔離感染者数 $Q$ の影響はないとしたものです。
●SIQR2とSIQR3は、SIR2の $I=H+Q$ を分けたものに相当します。SIR2の $\beta S(t)I(t)$ については

$(26)\quad \beta S(t)I(t)=(1-q)\beta S(t)I(t)+q\beta S(t)I(t)$

と分解し、 $I=H+Q$ のうち隔離感染者数 $Q$ の影響はないとしています。ここでパラメータ $q$ は感染率 $\beta$ に対する削減率とみなせます。またSIQR2とSIQR3には、 $pH(t)$ が差し引きしてあります。この $pH(t)$ は検査後陽性判定を受け隔離される人口で、日々発表される感染者数とされています。ここでパラメータ $p$ は市中感染者に対する隔離率とみなせます。

このとき、SIRモデルの場合の(5)と同様に

$(27)\quad \frac{dH(t)}{dt}=\underbrace{((1-q)\beta S_0-\gamma-p)}_{\lambda_0(q,p)}H(t), \quad H(0)=H_0$

が成り立ち、初期成長率 $\lambda_0(q,p)$ は

$(28)\quad \lambda_0(q,p)=(\gamma+p)((1-q)\underbrace{\frac{\beta S_0}{\gamma+p}}_{\tilde{R}_0}-1)$

と表すことができます。これをSIRモデルの場合の

$(7)\quad \lambda_0=\gamma(\underbrace{\frac{\beta S_0}{\gamma}}_{R_0}-1)$

と比較してみますと、 $\gamma$ が $\gamma+p$ に代わっており、基本再生産数が小さくなっています。したがって、SIQRモデルによれば、隔離率 $p$ を上げれば基本再生産数が小さくなる可能性を理解できます。

(27)の解は

$(29)\quad H(t)=\exp(\lambda_0(q,p)t)H_0\rightarrow0\quad(t\rightarrow\infty)$

と表され、初期値 $H_0$ の10%まで下がる時間は( $\lambda_0(q,p)<0$ に注意して)

$(30)\quad \exp(\lambda_0(q,p)t)=0.1\quad\Rightarrow\quad t=\frac{\ln(0.1)}{|\lambda_0(q,p)|}=\frac{2.3}{|\lambda_0(q,p)|}$

で求められます。

当該論文の重要な貢献は、初期成長率を

$(31)\quad \lambda_0=\beta S_0-\gamma\quad\Rightarrow\quad \lambda_0(q,p)=(1-q)\beta S_0-(\gamma+p)$

のように、隔離率 $p$ と削減率 $q$ で変化させることができることを明らかにしていることです。そして、当該論文では、全国の感染者数データについて、次の結論を示しています。

Case 1) 現状： $\lambda_0(0.6,p)=-0.075$ 　⇒　収束期間31日
Case 2) ８割自粛： $\lambda_0(0.8,p)=-0.101$ 　⇒　収束期間23日
Case 3) ロックダウン： $\lambda_0(1,p)=-0.126$ 　⇒　収束期間18日
Case 4) 検査４倍増： $\lambda_0(0,4p)=-0.288$ 　⇒　収束期間8日
Case 5) ５割自粛＋検査２倍増： $\lambda_0(0.5,2p)=-0.159$ 　⇒　収束期間14日

これらは次を(31)と(30)に代入して確認できます。

$(32)\quad \beta S_0=0.126$ 、 $\gamma=0.03$ 、 $p=0.096$

たとえば、Case 1)については

$(33)\quad \underbrace{0.4}_{1-q}\times \underbrace{0.126}_{\beta S_0}-\underbrace{0.096}_{p}-\underbrace{0.03}_{\gamma}=\underbrace{-0.075}_{\lambda_0}\quad\Rightarrow\quad \underbrace{2.3/0.075}_{2.3/|\lambda_0|}=31$

当該論文では、基本再生産数を

感染頂上期間： $\tilde{R}_0=\frac{\beta S_0}{\gamma+p}=\frac{0.126}{0.03+0.096}=1$
感染減少期間： $\tilde{R}_0=\frac{(1-q)\beta S_0}{\gamma+p}=\frac{0.051}{0.03+0.096}=0.405$

のように計算しています。これから現状の削減率は $q=1-0.051/0.126=0.595$ であることが確かめられます。また、回復率としては $\gamma=0.03$ 、感染性期間が $T=1/\gamma=33$ が採用されています。さらに、感染減少期間に入る前の感染頂上期間では $\tilde{R}_0=1$ となるように、隔離率 $p=0.096$ 、 $1/p=10.4$ が定められています。このように、感染者数データから隔離率 $p$ を見積もることができ、 $1/p$ 倍して市中感染者数を推定できるという意味で大変興味深いと思います。

注5：本稿の計算では、基本再生産数は

感染増大期間： $R_0=\frac{\beta S_0}{\gamma}=\frac{0.128}{0.04}=3.191$
感染減少期間： $R_0=\frac{\beta S_0}{\gamma}=\frac{0.022}{0.04}=0.545$

となります。そこで当該論文と同様の計算をしてみますと

感染頂上期間： $\tilde{R}_0=\frac{\beta S_0}{\gamma+p}=\frac{0.128}{0.04+0.088}=1$
感染減少期間： $\tilde{R}_0=\frac{(1-q)\beta S_0}{\gamma+p}=\frac{0.022}{0.04+0.088}=0.545$

となります。これから削減率 $q=1-0.022/0.128=0.829$ 、隔離率 $p=0.088$ 、 $1/p=11.4$ が分かります。これらを用いると次の計算結果を得ます。

Case 2) ６割自粛： $\lambda_0(0.6,p)=-0.077$ 　⇒　収束期間30日
Case 1) 現状： $\lambda_0(0.829,p)=-0.106$ 　⇒　収束期間22日
Case 3) ロックダウン： $\lambda_0(1,p)=-0.128$ 　⇒　収束期間18日
Case 4) 検査４倍増： $\lambda_0(0,4p)=-0.263$ 　⇒　収束期間9日
Case 5) ５割自粛＋検査２倍増： $\lambda_0(0.5,2p)=-0.151$ 　⇒　収束期間15日

本稿のデータ分析では、すでに８割強の自粛が行われている結果になっていますが、PCR検査を増やすことの効果を確認できたと言えます。ただし、偽陽性と偽陰性の可能性を勘案して、陽性の出る可能性の高い場合に検査を行う必要があることは述べるまでもありません。

謝辞
最後に、本サイトでまとめた感染シミュレーションは、３月末からわたって検討してきました。その間多くの知人に補遺を除く本サイトを見ていただき、コメントをいただき、その都度見直すことができて、大変感謝申し上げております。本来ならばお名前を挙げて謝意を表すべきところですが、まだ最終原稿となっていない今の段階では、文責を筆者のみに帰すため、遠慮させていただきます。

5/17/2020 梶原宏之　kajiwara@cacsd2.sakura.ne.jp

追記
有料サイトのため、全文を読むことができませんが、次の発表がありました。 $\gamma$ の値は0.04（感染性期間25日）は長すぎるかもしれません。

「コロナ感染は発症前から発症5日が最多、6日以後はほとんど感染しない！」
台湾では2020年5月14日現在、新型コロナウイルス感染症（COVID-19）患者440人、死者7人で感染制御に成功しています。その台湾衛生福利部疾病管制署（台湾CDC）からJAMA Internal Medicine（2020年5月1日オンライン版）に「台湾のCOVID-19感染ダイナミクス」の論文が出ました。100人の感染確定者と2,761人の接触者について追跡（contact tracing）を行い、時期による感染力の違いが明らかになりました。最重要点は次の3つです。接触者追跡は発症4日前までさかのぼれ、COVID-19は発症前から発症5日までに感染リスクがあり、6日以降はない！COVID-19の二次感染までの期間（serial interval）は４～５日、SARSは8.4日。
西伊豆健育会病院病院長　仲田和正
2020年05月14日 16:50

制御系CAD

制御を意識したモノつくりを目指して

感染ダイナミックス（2020/5/17)