Implemented Jensen Shannon  divergence; moved all stats stuff out into stats.py