MikuHatsune https://blog.hatena.ne.jp/MikuHatsune/ 驚異のアニヲタ社会復帰の予備 https://mikuhatsune.hatenadiary.com/ R Rpackage Rを使いこなす 医学 感染症学 数学いらずの医科統計学 数理モデル rstan 統計 社会医学 新型コロナウイルス 国内感染の状況 こちらの2021年1月13日時点の各年代別感染者数と死亡者数のデータがあったので単純に推定してみる。その1 drcパッケージにあるLL.4関数は、4パラメータモデルのロジスティック曲線に当てはめる。 これで簡単にやるその2 rstan を使う。年齢を経るごとに死亡率が上がるとして、orderedで昇順の値のサンプリングをして、inv_logitにより0~1 の値になるようにする。すなわちこれが死亡率になり、やっていることはロジスティック曲線への変換である。死亡率は以下のようになる(%) age [1,] 5 0.007447185 [2,] 15 0.007… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fmikuhatsune.hatenadiary.com%2Fentry%2F2021%2F01%2F18%2F233158" title="新型肺炎の年齢階級別死亡率 - 驚異のアニヲタ社会復帰の予備" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/M/MikuHatsune/20210118/20210118232728.png Hatena Blog https://hatena.blog 2021-01-18 23:31:58 rich https://mikuhatsune.hatenadiary.com/entry/2021/01/18/233158 1.0 100%