標本分布

gistfile1.py

import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
import numpy as np


sample = 10000
mu, std = 0, 10
data = [np.mean(np.random.normal(mu, std, sample)) for x in range(sample)]

n, bins, patches = plt.hist(data, bins=50, normed=True)
line = mlab.normpdf(bins, np.mean(data), np.std(data))
plt.plot(bins, line, 'r--', linewidth=1)

plt.show()