Python練習（不偏分散の続き）

import random
import numpy as np

def print_vars(population, ddof=0, samples=100, sample_size=100):
    random.seed(314)
    print('ddof =', ddof)
    print('samples =', samples)
    print('sample_size =', sample_size)
    print('population_size = ', len(population))
    vs = []
    for k in range(0, samples):
        l = random.sample(population, sample_size)
        vs.append(np.var(l, ddof=ddof))
    a = np.mean(vs)
    b = np.var(population, ddof=ddof)
    print('(a) mean of var = ', a)
    print('(b) var         = ', b)
    print('(a) - (b)       = ', a - b)
    print()

np.random.seed(314159)
na = np.random.rand(10000) * 100
print(na)
print()
s = list(na)
print_vars(s, ddof=0)
print_vars(s, ddof=1)