20/07/07 10:48:38 vMNjk2cp.net
>>557
再生産数の事前分布を
西浦が公開しているStanのソースで
Rt ~ normal(2.4 ,2)
という記述があるから使っているよ。
正規分布にすると負の数も含むから俺ならガンマ分布にするけどね。
形状母数1.44 尺度母数5/3のガンマ分布にしたら平均2.4 標準偏差2で負の数にはならない分布ができる。
一様分布をそのままつかうとStanが収束しないので[0,10]の平均と分散を出して近似した正規分布でやっても
結果はあんまりかわらんね。
西浦モデルの再生算数の事前分布を変化させてグラフにしてみた。
西浦モデルでのデフォルト
URLリンク(i.imgur.com)
事前分布を一様分布(0,10)近似の正規分布で近似させた場合
URLリンク(i.imgur.com)