def plot_scipy_pmf():
"""
https://docs.scipy.org/doc/scipy/reference/generated/
scipy.stats.betabinom.html#scipy.stats.betabinom
"""
bb = BetaBinomial(n=5, alpha=2.3, beta=0.63)
axf = AxesFormatter()
bb.pmf().plot(k=range(0, 5), kind='line', ls='',
marker='o', color='blue', vlines=True, ax=axf.axes)
axf.set_x_lim(-0.5, 4.5)
axf.show()
def plot_wikipedia_cdfs():
"""
https://en.wikipedia.org/wiki/File:Beta-binomial_cdf.png
"""
axf = AxesFormatter(width=12, height=9)
for alpha, beta, color in zip(alpha_wiki, beta_wiki, colors_wiki):
BetaBinomial(n=10, alpha=alpha, beta=beta).cdf().plot(
k=k_wiki, color=color, kind='line', ax=axf.axes,
marker='o'
)
axf.add_legend()
axf.show()