"
英国：第1波を20/2/15-20/7/31，第2/3波を20/8/1-とした場合，死亡者数/PCR陽性者数比は，第1波 41189/303181=0.136 → 第2波 26212/1736966=0.015
イタリア：第1波を20/2/15-20/6/30，第2波を20/8/1-とした場合，死亡者数/PCR陽性者数比は，第1波 35286/240599=0.147 → 第2波 33504/1706027=0.020
"
rm(list=ls())
par(mfrow=c(2,2))
k=1e5
de1=41189 # deaths in England at 1st wave
pe1=303181 # positives in England at 1st wave
de2=26212
pe2=1736966
e1=rbeta(k,1+de1,1+pe1-de1) ; hist(e1)
e2=rbeta(k,1+de2,1+pe2-de2) ; hist(e2)
hist(e2-e1)
hist(e2/e1)

di1=35286
pi1=240599
di2=233504
pi2=1706027
i1=rbeta(k,1+di1,1+pi1-di1) ; hist(i1,col=2)
i2=rbeta(k,1+di2,1+pi2-di2) ; hist(i2,col=2)
hist(i2-i1,col=2)
hist(i2/i1,col=2)

layout(1)
BEST::plotPost(i2/i1-e2/e1) # >0
BEST::plotPost((i2/i1)/(e2/e1)) # >1