In [23]:

    

from scipy import stats
# Graficando Beta
a, b = 5, 1 # parametros de forma.
def f1(a,b):
    beta = stats.beta(a, b)
    x = np.linspace(beta.ppf(0.01),
                    beta.ppf(0.99), 100)
    fp = beta.pdf(x) # Función de Probabilidad
    plt.plot(x, fp)
f1(0.5,0.5)
f1(5,1)
f1(1,3)
f1(2,2)
f1(2,5)
plt.title('Distribucion Beta')
plt.ylabel('probabilidad')
plt.xlabel('valores')
plt.xlim(0, 1)
plt.ylim(0, 2.5)
plt.show()


    

In [25]:

    

def random_beta(a,b):
    rand_beta = np.random.beta(a,b,size=100)
    return rand_beta
first_rand = random_beta(0.5,0.5)
print(first_rand.mean())


    

    




0.490882546775

In [ ]: