Time and frequency domain methods for quantifying common modulation of motor unit firing patterns

BACKGROUND: In investigations of the human motor system, two approaches are generally employed toward the identification of common modulating drives from motor unit recordings. One is a frequency domain method and uses the coherence function to determine the degree of linear correlation between each frequency component of the signals. The other is a time domain method that has been developed to determine the strength of low frequency common modulations between motor unit spike trains, often referred to in the literature as 'common drive'. METHODS: The relationships between these methods are systematically explored using both mathematical and experimental procedures. A mathematical derivation is presented that shows the theoretical relationship between both time and frequency domain techniques. Multiple recordings from concurrent activities of pairs of motor units are studied and linear regressions are performed between time and frequency domain estimates (for different time domain window sizes) to assess their equivalence. RESULTS: Analytically, it may be demonstrated that under the theoretical condition of a narrowband point frequency, the two relations are equivalent. However practical situations deviate from this ideal condition. The correlation between the two techniques varies with time domain moving average window length and for window lengths of 200 ms, 400 ms and 800 ms, the r(2 )regression statistics (p < 0.05) are 0.56, 0.81 and 0.80 respectively. CONCLUSIONS: Although theoretically equivalent and experimentally well correlated there are a number of minor discrepancies between the two techniques that are explored. The time domain technique is preferred for short data segments and is better able to quantify the strength of a broad band drive into a single index. The frequency domain measures are more encompassing, providing a complete description of all oscillatory inputs and are better suited to quantifying narrow ranges of descending input into a single index. In general the physiological question at hand should dictate which technique is best suited.