二項分布から生成されると仮定される2つのデータがあったときにその差（2つの二項分布のパラメタの差）を推定することを考える。 In [1]: import numpy as np import pandas as pd from pandas import DataFrame, Series from matplotlib import pyplot as p…

前回の記事で二項分布のパラメタ推定をPyMC2で行った。 しかしPyMC2をいまから使っていくのも微妙な気がしてきたので新しいPyMC3で書き直す。 ついでにPyStanでも同じことをして比べてみる。 In [1]: import numpy as np import pandas as pd from pandas im…