Several recent studies have used matrix factorization algorithms to assess the hypothesis that behaviors might be produced through the combination of a small number of muscle synergies. Although generally agreeing in their basic conclusions, these studies have used a range of different algorithms, making their interpretation and integration difficult. We therefore compared the performance of these different algorithms on both simulated and experimental data sets. We focused on the ability of these algorithms to identify the set of synergies underlying a data set. All data sets consisted of nonnegative values, reflecting the nonnegative data of muscle activation patterns. We found that the performance of principal component analysis (PCA) was generally lower than that of the other algorithms in identifying muscle synergies. Factor analysis (FA) with varimax rotation was better than PCA, and was generally at the same levels as independent component analysis (ICA) and nonnegative matrix factorization (NMF). ICA performed very well on data sets corrupted by constant variance Gaussian noise, but was impaired on data sets with signal-dependent noise and when synergy activation coefficients were correlated. Nonnegative matrix factorization (NMF) performed similarly to ICA and FA on data sets with signal-dependent noise and was generally robust across data sets. The best algorithms were ICA applied to the subspace defined by PCA (ICAPCA) and a version of probabilistic ICA with nonnegativity constraints (pICA). We also evaluated some commonly used criteria to identify the number of synergies underlying a data set, finding that only likelihood ratios based on factor analysis identified the correct number of synergies for data sets with signal-dependent noise in some cases. We then proposed an ad hoc procedure, finding that it was able to identify the correct number in a larger number of cases. Finally, we applied these methods to an experimentally obtained data set. The best performing algorithms (FA, ICA, NMF, ICAPCA, pICA) identified synergies very similar to one another. Based on these results, we discuss guidelines for using factorization algorithms to analyze muscle activation patterns. More generally, the ability of several algorithms to identify the correct muscle synergies and activation coefficients in simulated data, combined with their consistency when applied to physiological data sets, suggests that the muscle synergies found by a particular algorithm are not an artifact of that algorithm, but reflect basic aspects of the organization of muscle activation patterns underlying behaviors.

How the CNS selects the appropriate muscle patterns to achieve a behavioral goal is an open question. To gain insight into this process, we characterized the spatiotemporal organization of the muscle patterns for fast-reaching movements. We recorded electromyographic activity from up to 19 shoulder and arm muscles during point-to-point movements between a central location and 8 peripheral targets in each of 2 vertical planes. We used an optimization algorithm to identify a set of time-varying muscle synergies, i.e., the coordinated activations of groups of muscles with specific time-varying profiles. For each one of nine subjects, we extracted four or five synergies whose combinations, after scaling in amplitude and shifting in time each synergy independently for each movement condition, explained 73–82% of the data variation. We then tested whether these synergies could reconstruct the muscle patterns for point-to-point movements with different loads or forearm postures and for reversal and via-point movements. We found that reconstruction accuracy remained high, indicating generalization across these conditions. Finally, the synergy amplitude coefficients were directionally tuned according to a cosine function with a preferred direction that showed a smaller variability with changes of load, posture, and endpoint than the preferred direction of individual muscles. Thus the complex spatiotemporal characteristics of the muscles patterns for reaching were captured by the combinations of a small number of components, suggesting that the mechanisms involved in the generation of the muscle patterns exploit this low dimensionality to simplify control.

Mammalian locomotion patterns share common roots.

Selecting the appropriate muscle pattern to achieve a given goal is an extremely complex task because of the dimensionality of the search space and because of the nonlinear and dynamical nature of the transformation between muscle activity and movement. To investigate whether the central nervous system uses a modular architecture to achieve motor coordination we characterized the motor output over a large set of movements. We recorded electromyographic activity from 13 muscles of the hind limb of intact and freely moving frogs during jumping, swimming, and walking in naturalistic conditions. We used multidimensional factorization techniques to extract invariant amplitude and timing relationships among the muscle activations. A decomposition of the instantaneous muscle activations as combinations of nonnegative vectors, synchronous muscle synergies, revealed a spatial organization in the muscle patterns. A decomposition of the same activations as a combination of temporal sequences of nonnegative vectors, time-varying muscle synergies, further uncovered specific characteristics in the muscle activation waveforms. A mixture of synergies shared across behaviors and synergies for specific behaviors captured the invariances across the entire dataset. These results support the hypothesis that the motor controller has a modular organization.