A robust, iteratively prewhitened general linear model using multiple short-separation measurements provides the best correction for systemic physiological signals in functional near-infrared spectroscopy.
The use of multiple SS measurements as regressors of no interest implemented in a robust, iteratively prewhitened, general linear model has the best performance of the tested existing methods.
Santosa et al. (Wed,) studied this question.