def jensen_shannon_divergence(p,q): """Jensen-Shannon distance between distributions p and q""" m = (p+q)/2.0; return stats.entropy(p,m) + stats.entropy(q,m);