【Pythonで学ぶデータ分析】相関があるかどうかをベイズ統計で調べる 〜 年齢と原付事故死傷者数に関係はあるのか？

図2　事前分布と尤度関数の考え方相関係数ρの分布は、台を−1〜1に拡張したベータ分布（拡張ベータ分布）とする。図中のExtBetaという関数はPyMCに存在しないので、ベータ分布を変換して使う。尤度関数は多変量正規分布（この場合は2変数）。多変量正規分布を定義するためのMvNormalクラスには、各変数の平均μ0,μ1と、分散・共分散行列Σを指定する。Σは合計を表す記号ではなく、分散・共分散行列を表す変数の名前であることに注意。