example to show convergence on global mean using runningstats

conv_stats.py

import numpy as np
import matplotlib.pylab as plt

# generate some data
x = np.random.rand(1000,1)

# get the class
import RunningStats

# get an instance of the class
rs = RunningStats.RunningStats()

# convergence on global mean
convav = []
for k in x:
   rs.Push(k[0])
   convav.append(rs.Mean())

plt.plot(x); plt.plot(convav,'r'); plt.show()