2020年5月6日水曜日

米国の集団免疫率(3)

米国の集団免疫率(2)からの続き

タイトルはもう変更したほうがいいかもしれない。というのも,抗体検査の結果,ニューヨーク州では12.3%が抗体を持っている(感染済という結果がでているからだ。こちらの結果とはほぼ1桁違うので,我々のモデルの前提や仮定のどれかががまったく間違っているのではないか。しかし,モデルを検討する余力がないので(遠隔授業の準備で手いっぱい),そこは放置したまま,米国の新しいデータに基づいたパラメタ推定を行う。というのも,これまでの6万人から6.5万人という発表に代わり,再び死亡者が10万人を越えるという予想が出ているからだ。前回と同様に,HEMLのCOVID-19 projections のページを見れば,確かに米国全体の死亡数は8月には13.4万人になりそうだとある。

前回同様のSIIDR2モデルで計算する。使用するWHOのデータ(人口1万人当り)はこれ。
- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
ya=[0.014,0.021,0.030,0.038,0.051,0.051,0.051,
0.106,0.107,0.215,0.317,0.462,0.462,0.958,
1.280,1.576,1.929,2.074,2.587,3.136,3.722,
4.268,4.953,5.684,6.483,7.335,8.310,9.327,
10.13,11.03,11.99,12.93,14.00,14.96,15.92,
16.81,17.55,18.33,19.20,20.19,21.10,21.96,
22.80,23.58,24.31,25.19,26.12,27.29,28.28,
29.16,29.85,30.47,31.42,32.39,33.20,34.16]
za=[0.0006,0.0008,0.0009,0.0011,0.0012,0.0012,0.0012,
0.0018,0.0018,0.0030,0.0046,0.0061,0.0061,0.0122,
0.0143,0.0204,0.0268,0.0301,0.0377,0.0506,0.0641,
0.0728,0.0865,0.117,0.146,0.178,0.213,0.254,
0.290,0.329,0.387,0.445,0.504,0.562,0.620,
0.667,0.712,0.785,0.856,0.922,0.984,1.038,
1.089,1.141,1.216,1.284,1.337,1.402,1.456,
1.492,1.532,1.591,1.679,1.742,1.792,1.842]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - -

パラメタはこれ,前回と少し変えている。
#β=0.61,ν=0.12,λ=49,τ=16,α2=5.0/0.20,γ2†=15.0/0.10 

結果はこれ,確かに15〜16万人くらいになりそうだ。
図 米国の感染カーブ(u3=重症感染数,u4=死亡数,u6=新規感染数累計,1万人当)

P. S. 日本でも,神戸の病院の外来患者の約3%に新型コロナウイルス感染症の抗体が検出されたとある。なかなか大きな数字だ。これが1%だとしても,大都市部の数十万人が感染済みということ。

0 件のコメント:

コメントを投稿