始める無料で始める

事後分布の更新

前の演習で有効性率の事後分布を推定できました。よくできました! ただし、データ数が少ないため、この分布はかなり広く、薬の品質に関する不確実性が大きいことを示しています。幸い、薬の試験は続いており、新たに12人の患者に投与したところ、10人が治癒しました。これらの新しいデータで事後分布を更新しましょう!

ベイズ的アプローチなら簡単です。先ほどと同様にグリッド近似を実行しますが、異なる事前分布を使います。前の演習で得た事後分布(有効性率に関する私たちの知識を表すもの)を新しい事前分布として利用できます。そのうえで新しいデータに対する尤度を計算し、新しい事後分布を得ます!

前の演習で作成した DataFrame df はワークスペースに用意されています。また、scipy.stats から binom がインポート済みです。

この演習はコースの一部です

Pythonで学ぶBayesianデータ分析

コースを見る

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Assign old posterior to new prior and calculate likelihood
df["new_prior"] = ____
df["new_likelihood"] = ____(df["num_patients_cured"], ____, df["efficacy_rate"])
コードを編集して実行