  We discuss the status of the RBC & UKQCD collaboration's lattice
determination of $\epsilon'$, the measure of direct CP-violation in kaon
decays, focusing in particular on recent developments in statistical techniques
for estimating standard errors and goodness-of-fit metrics for large amounts of
data that have correlations both in the temporal coordinate and also in
molecular dynamics time. A key result is the formulation of a technique for
determining the null distribution of a fit using a bootstrap method that is
free from the usual assumptions of independence, large-$n$ and/or normal data,
that can for instance be applied to compute p-values even for uncorrelated
fits.
