20/07/06 14:26:08 RBXc/6V3.net
>>638
>変数は実効再生産数ただひとつだろw
その分布をRとStanを使って求めている。
再生産数の事前分布は西浦のソースではRt ~ normal(2.4 ,2)
平均2.4 標準偏差2の正規分布と設定しているが、
この95%CIは
> hdi(qnorm,0.95,mean=2.4,sd=2)[1:2]
lower upper
-1.519928 6.319928
なので、0も含んでいるし、弱情報事前分布としては適切
感染してから他人にうつすまでの日数:Generatin timeは発症間隔Serial Intervalを使っている。
Stanでの西浦モデルではinfecterとinfecteeが発症するまでの期間serial intervalの分布に
## Serial interval [Nishiura et al 2020 - only certain cases]
param1_SI = 2.305,
param2_SI = 5.452,
// serial interval
vector[K] gt = pweibull(param1_SI, param2_SI, K);
面倒なのが、
「感染してから発表の数字に載るまでの日数の分布」