Currently some high-level statistics (mean, std, quantiles, etc) are calculated separately in several different parts of the module. These should be consolidated into high level function(s) that can be called rather than having redundant code.
In response to this comment: #206 (comment)