population = [139212, 70192, 50000, 21000, 16000, 5000, 2000, 500, 600, 100, 10, 5, 5, 5, 5]

sample = [139212, 70192, 21000, 2000, 100, 5, 5]

>>> import numpy as np
>>> from scipy import stats
>>> population = [139212, 70192, 50000, 21000, 16000, 5000, 2000, 500, 600, 100, 10, 5, 5, 5, 5]
>>> sample = [139212, 70192, 21000, 2000, 100, 5, 5]
>>> np.mean(population)
20308.933333333334
>>> np.mean(sample)
33216.285714285717
>>> stats.sem(population)
10082.215934371847
>>> stats.sem(sample)
20135.645011620902

>>> import numpy as np
>>> from scipy import stats
>>> population = [139212, 70192, 50000, 21000, 16000, 5000, 2000, 500, 600, 100, 10, 5, 5, 5, 5]
>>> sample = [139212, 70192, 21000, 2000, 100, 5, 5]

>>> stats.norm.interval(0.9, loc=sample[0])
(139210.35514637304, 139213.64485362696)

>>> stats.norm.interval(0.8, loc=sample[0])
(139210.71844843446, 139213.28155156554)

>>> stats.norm.interval(0.7, loc=sample[0])
(139210.9635666105, 139213.0364333895)