Often you will have a set of data and need to compute error bars. We have the following statistical package for you: stats.py that will help you to do this. To use it, you will need the numpy package installed.

If you are going to be doing the QMC tutorials, download this package now and familiarize yourself with its use. As an example to using this package, imagine you have a set of data contained in a numpy array, like myData=numpy.array([0.2,1.2,0.04,0.9,1.2]). You might want to calculate the mean, variance, error bars, and autocorrelation time for this data. Use the following code to do this:

import stats
myStats= stats.Stats(myData)
print "Mean: ",myStats[0]
print "Variance: ",myStats[1]
print "Error: ",myStats[2]
print "Autocorrelation time: ",myStats[3]