In [9]:

    

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

fig, ax = plt.subplots(1, 1)
ax.set_xlim([0,1])
ax.set_ylim([0,3])
a1, b1 = 0.5, 0.5
a2, b2 = 5, 1
a3, b3 = 1, 3
a4, b4 = 2, 2
a5, b5 = 2, 5

mean, var, skew, kurt = beta.stats(a, b, moments='mvsk')
x = np.linspace(beta.ppf(0.01, a1, b1), beta.ppf(0.99, a1, b1), 100)
ax.plot(x, beta.pdf(x, a1, b1),'r-', lw=3, alpha=0.6, label='beta1 pdf')

y = np.linspace(beta.ppf(0.01, a2, b2), beta.ppf(0.99, a2, b2), 100)
ax.plot(y, beta.pdf(y, a2, b2),'b-', lw=3, alpha=0.6, label='beta2 pdf')

z = np.linspace(beta.ppf(0.01, a3, b3), beta.ppf(0.99, a3, b3), 100)
ax.plot(z, beta.pdf(z, a3, b3),'g-', lw=3, alpha=0.6, label='beta3 pdf')

w = np.linspace(beta.ppf(0.01, a4, b4), beta.ppf(0.99, a4, b4), 100)
ax.plot(w, beta.pdf(w, a4, b4),'y-', lw=3, alpha=0.6, label='beta4 pdf')

v = np.linspace(beta.ppf(0.01, a5, b5), beta.ppf(0.99, a5, b5), 100)
ax.plot(v, beta.pdf(v, a5, b5),'p-', lw=5, alpha=0.6, label='beta5 pdf')

plt.show()


    

In [ ]: