一日坊主

雰囲気でやっている

2020-07-01から1ヶ月間の記事一覧

PyMCによるベイズ推論

本記事では以下のレポジトリの Ch1_Introduction_PyMC3.ipynb の内容を一部抜粋して実行してみる. github.com まず,データをダウンロードして表示する. import numpy as np import matplotlib.pyplot as plt import os # 'data' ディレクトリを作成する m…

最尤推定によるガンマ分布のフィッティング

scipyを利用してデータのサンプル点をガンマ分布にフィッティングさせる. ガンマ分布の確率密度関数は以下で表される*1. import matplotlib.pyplot as plt from scipy.stats import gamma # ガンマ分布に従う乱数を1,000個生成 a, b = 3, 2 seed = 1 data …