  Extreme-mass-ratio inspirals (EMRIs) will be an important type of
astrophysical source for future space-based gravitational-wave detectors. There
is a trade-off between accuracy and computational speed for the EMRI waveform
templates required in the analysis of data from these detectors. We discuss how
the systematic error incurred by using faster templates may be reduced with
improved models such as augmented kludge waveforms, and marginalised over with
statistical techniques such as Gaussian process regression.
