def test_RHist():
import numpy as np
import scipy.stats as stats
from bigstats.hist import RHist
xs = stats.norm.rvs(size=10000)
rh = RHist(name='test',decimals=1)
[rh.incr(x) for x in xs]
print('The normal distribution was sampled 10,000 times....')
print('Est. Mean (0): {0}'.format(rh.mean()))
print('Est. Var (1): {0}'.format(rh.var()))
print('Est. Stdev (0.01): {0}'.format(rh.stdev()))
print('N (10000): {0}'.format(rh.n()))
print('** Decimals should be precise to 1 place.**')
fig = rh.plot(fig=None,norm=True)
# TODO: Add a uniform and some asymmetric dist too.
