ベータ分布 (numpy.random.beta)

プログラム例

# -*- coding:utf-8 -*-
import numpy as np
import matplotlib.pyplot as plt

# ----- 条件設定 -----
alpha, beta = 2.0, 5.0
N = 10000                # サンプル数
Nb = 50                 # ヒストグラムプロットのビン数
min_x, max_x = 0, 1.5

# ----- ベータ分布関数の乱数の生成 -----
s = np.random.beta(alpha, beta, N)

# ----- プロット作成 -----
plt.rcParams["font.size"] = 18
fig = plt.figure()
dist_plt = fig.add_subplot(1,1,1)
dist_plt.tick_params(axis='both', length=10, which='major')
dist_plt.tick_params(axis='both', length=5,  which='minor')
dist_plt.set_xlim([min_x, max_x])
dist_plt.hist(s, Nb, range=(min_x, max_x), \
            facecolor='b', edgecolor="black")

plt.show()
fig.savefig('beta.pdf', orientation='portrait', \
            transparent=False, bbox_inches=None, frameon=None)
fig.clf()

実行結果