中国湖北省武漢市(人口約1100万人)が今月中旬から全市民を対象に実施していたPCR検査(遺伝子検査)がおおむね終了した。
地元紙「湖北日報」(電子版)によると、直近に検査を済ませた人や乳児を除く900万人以上から検体を採取。
武漢では1日に100万人以上の検査が可能で、24日までに約650万人の検査が終わり、無症状感染者が218人見つかった。

事前分布に有病率は一様分布、感度は40-60%、特異度は90−100%のベータ分布を設定してStanでMCMCしてみた。
mean se_mean sd 2.5% 50% 97.5% n_eff Rhat
spc 0.999977 0.000000 0.000008 0.999965 0.999976 0.999994 12593 1.00001
sen 0.452612 0.002116 0.220226 0.073638 0.444200 0.871391 10834 1.00000
prev 0.000041 0.000001 0.000118 0.000001 0.000023 0.000183 8206 1.00004
p 0.000034 0.000000 0.000002 0.000030 0.000034 0.000039 36412 0.99998

spc:特異度、sen:感度  prev:有病率 p:陽性率

900万人での陽性者数の予想

最頻値 306.15

95%信頼区間
lower upper
266.46 346.70

1100万人のうちの有病者数予想

> summary(ms$prev*11e6) ; MODE(ms$prev*11e6)[1] ; hdi(ms$prev*11e6)[1:2]
Min. 1st Qu. Median Mean 3rd Qu. Max.
0 108 251 449 483 137796
x
179.24
lower upper
0.020213 1360.523395