ベータ分布のサンプル

beta_dist.py

# -*- coding: utf-8 -*-

import matplotlib
colors = matplotlib.rcParams['axes.color_cycle']

from matplotlib import pyplot as plt
from scipy.stats import beta
import numpy as np

x = np.linspace(0,1, num=200)
for a, b, c in zip([0.5, 0.5, 1, 2, 3], [0.5, 1, 3, 2, 5], colors):
    plt.plot(x, beta.pdf(x, a, b), lw=2,
             c=c, label=r"$\alpha = {0:.1f},\ \beta={1:.1f}$".format(a, b))

plt.legend(loc=0)
plt.ylabel("$pdf$")
plt.xlabel("$x$")
plt.show()