Fractal analysis of muscle activity patterns during locomotion: pitfalls and how to avoid them

Time-dependent physiological data sets are often difficult to interpret objectively. Biosignals such as electromyogram, electroencephalogram, or single-neuron recordings can be interpreted using various linear and nonlinear methods. Each analysis technique aims at the explanation of different data features that might be visible or not to the naked eye. Here, we used linear decomposition based on machine learning to extract motor primitives (the time-dependent coefficients of muscle synergies) from the hindlimb electromyographic activity of mice during normal and mechanically perturbed locomotion. We set out to investigate the effects of calculation parameters and data quality on two nonlinear metrics derived from fractal analysis: the Higuchi's fractal dimension (HFD) and the Hurst exponent (H). Both HFD and H proved to be exceptionally sensitive to changes in motor primitives induced by external perturbations to locomotion. We discuss the potential pitfalls that might arise from fractal analysis by using examples based on surrogate data. We conclude giving some simple, data-driven suggestions to reduce the chance of misinterpretations when metrics such as HFD and H are applied to any biological signal containing elements of periodicity. NEW & NOTEWORTHY Despite the lack of consensus on how to perform fractal analysis of physiological time series, many studies rely on this technique. Here, we shed light on the potential pitfalls of using the Higuchi's fractal dimension and the Hurst exponent. We expose and suggest how to solve the drawbacks of such methods when applied to data from normal and perturbed locomotion by combining in vivo recordings and computational approaches.