This deals with non-parametric stats on metrics that average across trials, leading to sparse data for the purpose of statistics. Requires a function to shuffle (surrogate) data along a particular axis, array of metrics to compute across trials, method to compute the null distribution and evaluate the original metric against it.