<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Hum. Neurosci.</journal-id>
<journal-title>Frontiers in Human Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Hum. Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1662-5161</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnhum.2017.00298</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Switching Adaptability in Human-Inspired Sidesteps: A Minimal Model</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Fujii</surname> <given-names>Keisuke</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/192029/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Yoshihara</surname> <given-names>Yuki</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/431898/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Tanabe</surname> <given-names>Hiroko</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/431866/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Yamamoto</surname> <given-names>Yuji</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/314710/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Structured Learning Team, Center for Advanced Intelligence Project, Institute of Physical and Chemical Research (RIKEN)</institution> <country>Suita, Japan</country></aff>
<aff id="aff2"><sup>2</sup><institution>Intelligence Mobility Group, Institutes of Innovation for Future Society, Nagoya University</institution> <country>Nagoya, Japan</country></aff>
<aff id="aff3"><sup>3</sup><institution>Graduate School of Arts and Sciences, University of Tokyo</institution> <country>Tokyo, Japan</country></aff>
<aff id="aff4"><sup>4</sup><institution>Research Center of Health Physical Fitness and Sports, Nagoya University</institution> <country>Nagoya, Japan</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Mikhail Lebedev, Duke University, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Rahul Goel, University of Houston, United States; Dai Owaki, Tohoku University, Japan; Kazuki Nakada, Hiroshima City University, Japan</p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x0002A;Correspondence: Keisuke Fujii <email>keisuke198619&#x00040;gmail.com</email></p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>07</day>
<month>06</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection">
<year>2017</year>
</pub-date>
<volume>11</volume>
<elocation-id>298</elocation-id>
<history>
<date date-type="received">
<day>15</day>
<month>03</month>
<year>2017</year>
</date>
<date date-type="accepted">
<day>22</day>
<month>05</month>
<year>2017</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2017 Fujii, Yoshihara, Tanabe and Yamamoto.</copyright-statement>
<copyright-year>2017</copyright-year>
<copyright-holder>Fujii, Yoshihara, Tanabe and Yamamoto</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Humans can adapt to abruptly changing situations by coordinating redundant components, even in bipedality. Conventional adaptability has been reproduced by various computational approaches, such as optimal control, neural oscillator, and reinforcement learning; however, the adaptability in bipedal locomotion necessary for biological and social activities, such as unpredicted direction change in chase-and-escape, is unknown due to the dynamically unstable multi-link closed-loop system. Here we propose a switching adaptation model for performing bipedal locomotion by improving autonomous distributed control, where autonomous actuators interact without central control and switch the roles for propulsion, balancing, and leg swing. Our switching mobility model achieved direction change at any time using only three actuators, although it showed higher motor costs than comparable models without direction change. Our method of evaluating such adaptation at any time should be utilized as a prerequisite for understanding universal motor control. The proposed algorithm may simply explain and predict the adaptation mechanism in human bipedality to coordinate the actuator functions within and between limbs.</p>
</abstract>
<kwd-group>
<kwd>sensory-motor system</kwd>
<kwd>multi-link system</kwd>
<kwd>closed-loop system</kwd>
<kwd>autonomous distributed control</kwd>
<kwd>flexible bipedal locomotion</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="7"/>
<ref-count count="41"/>
<page-count count="11"/>
<word-count count="7989"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>We can adaptively operate our bipedal body by cooperating with others in an emergency (Hutchins, <xref ref-type="bibr" rid="B18">1995</xref>; Fujii et al., <xref ref-type="bibr" rid="B11">2016</xref>) and sometimes competing with others (Yamamoto et al., <xref ref-type="bibr" rid="B39">2013</xref>; Fujii et al., <xref ref-type="bibr" rid="B10">2015b</xref>). Current technology can succeed in reproducing such real-time adaptation in video game tasks (Mnih et al., <xref ref-type="bibr" rid="B27">2015</xref>) and overcoming unpredicted deficits (Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>; Cully et al., <xref ref-type="bibr" rid="B6">2015</xref>), although such adaptation is limited to a certain part of the agent&#x00027;s body. However, with regard to bipedal locomotion, which is more dynamically unstable than that of more than four-legged species (Golubitsky et al., <xref ref-type="bibr" rid="B14">1999</xref>), researchers have not paid attention to the adaptive movements when motor commands change suddenly in response to a change in the situation, such as chase-and-escape behavior (Kamimura and Ohira, <xref ref-type="bibr" rid="B19">2010</xref>; Fujii et al., <xref ref-type="bibr" rid="B10">2015b</xref>), which have been acquired in over the course of evolution as biological (Carvalho et al., <xref ref-type="bibr" rid="B4">2012</xref>) and social (Helbing et al., <xref ref-type="bibr" rid="B17">2000</xref>) features essential for life activities. For example, it is considered that a sudden intentional direction change opposite to the original direction, such as in interpersonal sports (Fujii et al., <xref ref-type="bibr" rid="B10">2015b</xref>), is quite difficult to achieve and thus has been ignored in the fields of robotic engineering (Koolen et al., <xref ref-type="bibr" rid="B21">2016</xref>; Kuindersma et al., <xref ref-type="bibr" rid="B22">2016</xref>) and computational neuroscience (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>; Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>), with the focus primarily being placed on bipedal adaptation to external disturbances. Here, we refer to this as switching adaptation in bipedal locomotion because both motor commands (i.e., situation or task) and motor system requirements will switch in this case.</p>
<p>Although most previous studies on human motor control were based on the optimal control theory (Todorov and Jordan, <xref ref-type="bibr" rid="B36">2002</xref>; Scott, <xref ref-type="bibr" rid="B29">2004</xref>), which is considered to be physiologically related to the cerebellum (Shadmehr and Krakauer, <xref ref-type="bibr" rid="B30">2008</xref>), this theory cannot necessarily apply to universal motor control. The theory focuses on optimizing the system based on the centralized invariant cost functions, such as the deviation of target trajectory (Uno et al., <xref ref-type="bibr" rid="B37">1989</xref>) or motor cost as muscle activity (Anderson and Pandy, <xref ref-type="bibr" rid="B1">2003</xref>), such as in arm movement. However, an unstable multi-link closed-loop system with large inertia and a narrow base of support in abruptly changing situations, such as switching adaptation in bipedal locomotion, is difficult to control optimally. This is because it cannot determine the optimal target trajectory due to the large control component with physiological constraints (e.g., joints and muscles), the observation component with cognitive constraints (e.g., ground and opponent) and the context (e.g., the predicted optimal strategy could be defeated by the opponent&#x00027;s counter-attack; Fujii et al., <xref ref-type="bibr" rid="B10">2015b</xref>). Thus, switching adaptation in bipedal locomotion, which is difficult to control even in current robotics (Koolen et al., <xref ref-type="bibr" rid="B21">2016</xref>; Kuindersma et al., <xref ref-type="bibr" rid="B22">2016</xref>), is an excellent example to shed more light on the mystery of universal motor control.</p>
<p>Human bipedality, which is considered to be the result of adaptations to environmental variabilities (Carvalho et al., <xref ref-type="bibr" rid="B4">2012</xref>), is one of the controversial problems to control. While the efficiency of bipedal locomotion in the optimal control theory (Srinivasan and Ruina, <xref ref-type="bibr" rid="B33">2006</xref>) was explained by the dynamics only in the ground phase, neural oscillator control (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>) that is physiologically located in the spinal central pattern generator (Grillner, <xref ref-type="bibr" rid="B15">1985</xref>; Dimitrijevic et al., <xref ref-type="bibr" rid="B7">1998</xref>) can reproduce the whole of aperiodic adaptive bipedal locomotion in a self-organized manner rather than explicitly calculating the target trajectory or joint torques. However, the oscillator system is considered to be limited in cyclic movement with adaptation only to external disturbances (Thelen et al., <xref ref-type="bibr" rid="B35">1987</xref>; Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>). For example, active adaptation to a changing situation will result in excessive deviation from the aperiodic locomotion generated by the oscillator (e.g., in the opposite direction) because the motor command itself changes drastically. In recent years, using a physiological reflex model, the diversity of walking including a direction change of 50&#x000B0; was reproduced (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>), but in situations such as escape or pursuit, robustly faster direction change at any time (Fujii et al., <xref ref-type="bibr" rid="B9">2015a</xref>) is needed. Furthermore, it is unknown which factors make such adaptive bipedal locomotion difficult because previous locomotion models (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>; Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>) including multiple neural oscillators, peripheral reflexes and multi-link body dynamics were implemented in a complicated manner, whereas as far as the passive walk, the previous model simply accomplished it (McGeer, <xref ref-type="bibr" rid="B25">1990</xref>). Therefore, as a prerequisite for such adaptability, it is important to examine a minimal control model that achieves direction change at any time in the opposite direction with a small number of components and a simple algorithm, and to establish a methodology for evaluating it.</p>
<p>Distributed autonomous control, in which autonomous components implicitly function as a whole by interacting with each other without central control, such as in multi-agent (Couzin et al., <xref ref-type="bibr" rid="B5">2002</xref>) or multi-link (Watanabe et al., <xref ref-type="bibr" rid="B38">2012</xref>) biological systems, is applicable to real-time adaptation to the rapid impairment of components (Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>). This control system is biologically plausible than explicit simulation because the system can perform self-modeling (Bongard et al., <xref ref-type="bibr" rid="B2">2006</xref>) to adapt to the situation beyond its framework. The differences and advantages of the distributed autonomous control compared with the neural oscillator control are that the local components autonomously set the local target and have flexibility in the rule-based interaction among components. Among the autonomous system, self-repairing robots (Bongard et al., <xref ref-type="bibr" rid="B2">2006</xref>; Cully et al., <xref ref-type="bibr" rid="B6">2015</xref>) are remarkable, but the switching adaptation task in this study requires more improvisational adaptation (e.g., within 1 s). The mobility control (Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>) based on the design of autonomous systems, in which an autonomous mobile component moves prior to an immobile component, can execute arm reaching movement when confronting a real-time deficit of the component with improvisational adaptation. We thus assumed that mobility control can be a key factor in the switching adaptability with a minimal algorithm due to the real-time adaptability without the explicit control of the components. However, in bipedal locomotion, in addition to the control of the center of mass in locomotion (equal to endpoint control in arm movement), balance and leg swing control are necessary and often conflict, so not only the operation of equivalent rules for each component but also the switching of rules according to the situation should be important.</p>
<p>In this paper, we adopted switching autonomous system, which extended (i.e., incomplete) distributed autonomous control scheme, because the current task can be accomplished by solving multiple conflicting functions. For example, it would be more advantageous for multiple actuators to switch roles to maintain balance by the leading leg and to move the center of mass by the trailing leg (Yamashita et al., <xref ref-type="bibr" rid="B40">2013</xref>). In neurophysiology, this mechanism may be related to postural control in the reticulospinal tracts found in cats (Mori et al., <xref ref-type="bibr" rid="B28">1998</xref>), but its mechanisms of interaction and switching the function of actuators (i.e., muscles) have remained unknown. We therefore implemented an adaptive bipedal model into role-switching for propulsion, balance, and leg swing control with switching mobility control. The objective of this study is to propose a new control algorithm and evaluation methodology of a switching adaptive model for performing bipedal locomotion as a prerequisite for universal motor control.</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>Materials and methods</title>
<sec>
<title>Model overview</title>
<p>In this study, we constructed a three-mass model as a toy model (i.e., a minimally redundant model) of a sidestep locomotor system (Figure <xref ref-type="fig" rid="F1">1A</xref>). The three masses were linked with three actuators, springs, and dampers that represent the legs and inter-leg (<italic>i</italic> &#x0003D; 1&#x02013;3: inter-leg, right leg and left leg, respectively). For simplicity of the spatially symmetrical configuration of three actuators, inter-leg actuator 1 was modeled as a hip abductor and adductor muscles to swing the legs. Passive parameters are partly based on the human-like model in a previous study (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>), as shown in Table <xref ref-type="supplementary-material" rid="SM1">S1</xref>. In this model, the segments were stretchable, but if a leg exceeded a certain length (1.1 times its initial length), we increased the elastic coefficients (Table <xref ref-type="supplementary-material" rid="SM1">S1</xref>). We also increased the leg elasticity in the foot contact phase compared with that in the flight phase (Table <xref ref-type="supplementary-material" rid="SM1">S1</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Three-mass locomotion model. <bold>(A)</bold> The three masses (trunk and feet) were linked with three actuators, springs, and dampers that represent the hips (actuator 1) and legs (actuators 2 and 3). For simplicity of spatially symmetrical configuration in the three actuators, inter-leg actuator 1 was modeled as a hip abductor and adductor muscles to swing the legs. These segments were stretchable, but if a leg exceeded a certain length, we increased the elasticity (Table <xref ref-type="supplementary-material" rid="SM1">S1</xref>). The horizontal desired velocity is given as a single global command from the experimenter. Unidirectional arrows at the three masses and bidirectional arrows at the three actuators illustrate local desired velocities of autonomous actuators and actuator forces, respectively. <bold>(B)</bold> Local desired velocity without balancing (the same as in a previous study; Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>). For propulsion rightward, the inter-leg (actuator 1) and trailing leg (actuator 3) are designed to be lengthened and the leading leg (actuator 2) is designed to be shortened. <bold>(C)</bold> During grounding of the right leading leg, with a risk of falling forward, the inter-leg actuator and the leading leg are designed to be lengthened. <bold>(D)</bold> During grounding of only the left trailing leg, with a risk of raising the leading leg too high, the inter-leg actuator is shortened. Details in a coefficient <italic>a</italic> within brackets and unit vector <italic>e</italic><sub><italic>xi</italic></sub> are given in Section Materials and Methods.</p></caption>
<graphic xlink:href="fnhum-11-00298-g0001.tif"/>
</fig>
<p>The model can perform 2D lateral locomotion by sending appropriate commands to the actuators (Figure <xref ref-type="fig" rid="F1">1A</xref>) according to the following equation of motion:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>m</mml:mi></mml:mstyle><mml:mover><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>x</mml:mi></mml:mstyle><mml:mrow><mml:mn>..</mml:mn></mml:mrow></mml:mover><mml:mo>=</mml:mo><mml:mstyle displaystyle='true'><mml:msubsup><mml:mo>&#x02211;</mml:mo><mml:mi>i</mml:mi><mml:mn>3</mml:mn></mml:msubsup></mml:mstyle><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mi>a</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>m</mml:mi><mml:mi>g</mml:mi></mml:mstyle><mml:mo>+</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mi>p</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula>
<p>where <italic>x</italic> is a position vector of the three mass points, <italic>F</italic><sub><italic>ai</italic></sub> is an active force vector generated from the three actuators, <italic><bold>g</bold></italic> is a gravity acceleration vector, and <italic>F</italic><sub><italic>passive</italic></sub> is a passive force vector including viscoelasticity of the leg, its extension limit, and auxiliary action in the trunk (Table <xref ref-type="supplementary-material" rid="SM1">S1</xref>). The last auxiliary viscoelasticity prevents falling if the horizontal distance between the trunk and either leg is within 0.15 m. This value of 0.15 m is heuristically determined based on the trade-off of falling and propulsion in observation. We improved distributed autonomous control (Figures <xref ref-type="fig" rid="F1">1B&#x02013;D</xref>), which is based on the rule that the velocity commands are determined from the instantaneous &#x0201C;mobility&#x0201D; of each actuator in real time. This rule will be given in an autonomous decentralized form, which is explained in the paragraphs below.</p>
</sec>
<sec>
<title>Switching mobility control</title>
<p>For switching mobility control, here we consider the velocity command for actuator <italic>i</italic>. Command using positional information is not appropriate in this study because the calculation of the precise target trajectory is not needed. Velocity sensing and command may be reasonable such as due to the utilization of visual optical flow in a self-driven agent. It is assumed that the sensory (i.e., proprioceptive) information of the system including actuator lengths and angles and these derivative values was used. In this section, we consider the following two steps to construct the model: (i) First, the mobility index was defined as the difference between each local actuator&#x00027;s desired velocity (<inline-formula><mml:math id="M2"><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>) and the actual velocity <italic>v</italic><sub><italic>i</italic></sub>. (ii) Based on the mobility index, global desired velocity <italic>v</italic><sub><italic>d</italic></sub> was allocated preferentially to mobile actuators and control input in each actuator was determined.</p>
<p>In the first step, actuator <italic>i</italic> divides <italic>v</italic><sub><italic>d</italic></sub> into two components: a local vector <inline-formula><mml:math id="M3"><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and a residual vector <inline-formula><mml:math id="M4"><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>. The former is the component of <italic>v</italic><sub><italic>d</italic></sub> that actuator <italic>i</italic> could generate through its own stretching and shortening and the latter is the component that actuator <italic>i</italic> is incapable of generating in the current leg posture (Figure <xref ref-type="fig" rid="F1">1B</xref>):</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M5"><mml:mrow><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mi>l</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mtext>&#x02009;</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>e</mml:mi></mml:mstyle><mml:mrow><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x000B7;</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mi>d</mml:mi></mml:msub><mml:mo stretchy='false'>)</mml:mo><mml:mtext>&#x02009;</mml:mtext><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>e</mml:mi></mml:mstyle><mml:mrow><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula>
<disp-formula id="E3"><label>(3)</label><mml:math id="M6"><mml:mrow><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mi>r</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mi>d</mml:mi></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mi>l</mml:mi></mml:msubsup></mml:mrow></mml:math></disp-formula>
<p><italic><bold>e</bold></italic><sub><italic>xi</italic></sub> is a unit vector to produce the force in actuator <italic>i</italic>. <italic><bold>e</bold></italic><sub><italic>xi</italic></sub> in the inter-leg is defined as the unit vector from the trailing leg to the leading leg. Switching coefficient <italic>a</italic><sub><italic>i</italic></sub> is basically 1 for propulsion but switches for balancing and leg swing based on the related segment sensory information (the schematics are shown in Figure <xref ref-type="supplementary-material" rid="SM1">S1</xref>). A notable difference from a previous robot arm study (Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>) is that the local desired velocity is modified by coefficient <italic>a</italic><sub><italic>i</italic></sub> in the situations because of a temporal constraint to apply the force to the ground. When <italic>i</italic> is 2 or 3 (i.e., right or left leg) in the flight phase (Figure <xref ref-type="supplementary-material" rid="SM1">S1</xref> right), <italic>a</italic><sub><italic>i</italic></sub> was set to 1/2 because the leg cannot apply force to the ground and the contribution of the leg to trunk movement halved (Figure <xref ref-type="fig" rid="F1">1C</xref>: <inline-formula><mml:math id="M7"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mn>3</mml:mn></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, Figure <xref ref-type="fig" rid="F1">1D</xref>: <inline-formula><mml:math id="M8"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>). Additionally, when the trunk approaches the anterior leg within 0.2 m, <italic>a</italic><sub><italic>i</italic></sub> was set to &#x02212;1 to prevent falling, regardless of being in the flight or supported phase (Figure <xref ref-type="fig" rid="F1">1C</xref>: <inline-formula><mml:math id="M9"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>). When the posterior leg in the flight phase extended over its natural length, <italic>a</italic><sub><italic>i</italic></sub> was also set to &#x02212;1 to attract the leg to the trunk as the swing for the next step (Figure <xref ref-type="fig" rid="F1">1C</xref>: <inline-formula><mml:math id="M10"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mn>3</mml:mn></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>). These two corrections reflect the balance and leg swing, those conflict with propulsion, respectively.</p>
<p>When <italic>i</italic> is 1 (i.e., inter-leg: Figure <xref ref-type="supplementary-material" rid="SM1">S1</xref> left), <italic>a</italic><sub><italic>i</italic></sub> depended on the phase of both legs. In the double support phase, <italic>a</italic><sub><italic>i</italic></sub> was set to 0 because of a lack of contribution to trunk velocity. In the double flight phase, <italic>a</italic><sub><italic>i</italic></sub> was set to 1/2 in the same manner for both legs. Additionally, when the posterior leg length was over 0.6 times the natural length in the double flight phase or the anterior leg support phase, <italic>a</italic><sub><italic>i</italic></sub> was set to &#x02212;1 to attract the posterior leg to the trunk. This value of 0.6 was heuristically determined based on the following observation: if it is too large, the model sometimes cannot perform the leg swing, and if it is too small, it cannot move in the desired direction. Because of the dependence on kinematic sensory information of other segments, this system is not purely autonomous. However, this switching system contributed to achieving the task by resolving the trade-off between propulsion and balance.</p>
<p>For the adaptation under various environmental conditions, the mobility measure <italic>k</italic><sub><italic>i</italic></sub> must evaluate the instantaneous mobility of each joint appropriately, which requires calculation of kinematic and dynamic mobilities (Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>). The kinematic mobility is the ability of actuator <italic>i</italic> to move the trunk, and represented by the absolute value of the local desired velocity, i.e., geometric state of the actuator. Dynamic mobility is the same ability as determined by the dynamic properties of the actuator, which is represented by the difference between the actual state and the geometric state of the actuator. The mobility measure <italic>k</italic><sub><italic>i</italic></sub> is defined as the corrected ratio of dynamic mobility to kinematic mobility:</p>
<disp-formula id="E4"><label>(4)</label><mml:math id="M11"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mtext>exp</mml:mtext><mml:mo stretchy='false'>[</mml:mo><mml:mo>&#x02212;</mml:mo><mml:mn>4</mml:mn><mml:mo stretchy='false'>(</mml:mo><mml:mi>l</mml:mi><mml:mi>n</mml:mi><mml:mn>2</mml:mn><mml:mo stretchy='false'>)</mml:mo><mml:mo stretchy='false'>(</mml:mo><mml:mo>&#x02225;</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mi>l</mml:mi></mml:msubsup><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mi>i</mml:mi></mml:msub><mml:msup><mml:mo>&#x02225;</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003B5;</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo stretchy='false'>)</mml:mo><mml:mo>/</mml:mo><mml:mo stretchy='false'>(</mml:mo><mml:mo>&#x02225;</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mi>l</mml:mi></mml:msubsup><mml:mtext>&#x000A0;</mml:mtext><mml:msup><mml:mo>&#x02225;</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003B5;</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo stretchy='false'>)</mml:mo><mml:mo stretchy='false'>]</mml:mo></mml:mrow></mml:math></disp-formula>
<p>where <italic><bold>v</bold></italic><sub><italic>i</italic></sub> is the velocity produced at the trunk generated by actuator <italic>i</italic>, and &#x003B5;<sub>1</sub> and &#x003B5;<sub>2</sub> are small values (&#x003B5;<sub>1</sub>: 10<sup>&#x02212;10</sup>, &#x003B5;<sub>2</sub>: 10<sup>&#x02212;4</sup>) to avoid dividing by zero. The denominator and numerator are related to the kinematic and dynamic mobility of actuator <italic>i</italic>, respectively. The mobility <italic>k</italic><sub><italic>i</italic></sub> is supposed to take a value of 0 in an immobile actuator, and 1 in a mobile actuator.</p>
<p>Next, by using <italic>k</italic><sub><italic>i</italic></sub>, <inline-formula><mml:math id="M12"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, and <inline-formula><mml:math id="M13"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> (Equations 2&#x02013;4), we intended to design a real-time controller that would make the most mobile actuator work dominantly, and make the other actuators work cooperatively in order to satisfy <italic>v</italic><sub><italic>d</italic></sub>. Actuator <italic>i</italic> basically tries to move according to its own local vector <inline-formula><mml:math id="M14"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, and require the other actuators to create its residual vector <inline-formula><mml:math id="M15"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>. The required velocity from actuator <italic>j</italic> to actuator <italic>i</italic>, <inline-formula><mml:math id="M16"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, is defined as a projection of <inline-formula><mml:math id="M17"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> to <italic><bold>e</bold></italic><sub><italic>xi</italic></sub>:</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M18"><mml:mrow><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>e</mml:mi></mml:mstyle><mml:mrow><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x000B7;</mml:mo><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>j</mml:mi></mml:mrow><mml:mi>r</mml:mi></mml:msubsup><mml:mo stretchy='false'>)</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>e</mml:mi></mml:mstyle><mml:mrow><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula>
<p>The mobility <italic>k</italic><sub><italic>i</italic></sub> of each actuator then determines how the actuators interact with each other. We express this as follows:</p>
<disp-formula id="E6"><label>(6)</label><mml:math id="M19"><mml:mrow><mml:msub><mml:mrow><mml:mover accent='true'><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mo stretchy='true'>&#x002DC;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mo>&#x0220F;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mtext>&#x02009;</mml:mtext><mml:mo>&#x02260;</mml:mo><mml:mtext>&#x02009;</mml:mtext><mml:mi>i</mml:mi></mml:mrow><mml:mi>m</mml:mi></mml:msubsup><mml:mo stretchy='false'>(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>k</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo stretchy='false'>)</mml:mo><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mi>l</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mtext>&#x02009;</mml:mtext><mml:mo>&#x02260;</mml:mo><mml:mtext>&#x02009;</mml:mtext><mml:mi>i</mml:mi></mml:mrow><mml:mi>m</mml:mi></mml:msubsup><mml:msub><mml:mi>k</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:msubsup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:math></disp-formula>
<p><inline-formula><mml:math id="M20"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> is the velocity command for actuator <italic>i</italic>. The first term functions as an inhibitory interaction from actuator <italic>j</italic>, which prevents actuator <italic>i</italic> from moving according to its own local velocity <inline-formula><mml:math id="M21"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>. In contrast, the second term functions as an excitatory interaction make actuator <italic>i</italic> work cooperatively and generate the residual velocity of actuator <italic>j</italic>. In this study, we considered the minimal model which fully connected among three components (i.e., the nearest neighbors equal to the full connections). In more biological model, note that the number of the connection will increase and we should examine the connection configuration, such as based on the nearest neighbors. The velocity command, <inline-formula><mml:math id="M22"><mml:msubsup><mml:mrow><mml:mi mathvariant='bold'>v</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>, is transformed to torque as follows:</p>
<disp-formula id="E7"><label>(7)</label><mml:math id="M23"><mml:mrow><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mi>a</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mrow><mml:mover accent='true'><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mo stretchy='true'>&#x002DC;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>v</mml:mi></mml:mstyle><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy='false'>)</mml:mo></mml:mrow></mml:math></disp-formula>
<p>where <italic>G</italic><sub><italic>i</italic></sub> is the proportional gain of actuator <italic>i</italic>. We heuristically set it to 3,000 kg/s to perform the task.</p>
</sec>
<sec>
<title>Simulation and statistics</title>
<p>Initial horizontal and vertical positions of the trunk were set as 0 and 0.92 m, respectively. Initially, three masses were kept motionless in an equilateral triangular posture and double support stance. The time step in the simulation was set to 10<sup>&#x02212;5</sup> s. To examine the parameter sensitivity, we ideally should use the human parameter for the verification. Although we used the skeletal parameters based on the human parameters (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>), the passive joint viscoelasticity should be approximated when we expressed it as a scalar value (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>), because muscle-tendon complex has many components with various viscoelastic properties. Based on the assumption that we can learn the (sub) optimal parameters in various motor tasks, we approximated the passive joint viscoelasticity with the invariant scalar value. Instead, we examined the following three parameters from as broad a range as possible: the elasticity of legs, the elasticity of ground, and the proportional gain of actuators. We simulated the models to multiply each parameter by 10<sup>&#x02212;1</sup>, 10<sup>&#x02212;1/2</sup>, 1, 10<sup>1/2</sup>, and 10. The details and the results were given by Figure <xref ref-type="supplementary-material" rid="SM1">S3</xref>. In short, in the case of multiplying 10<sup>&#x02212;1/2</sup> or 10<sup>1/2</sup> by the each original parameter, the model had sensitivity enough to accomplish the direction change at any time.</p>
<p>To quantify the switching adaptability performance, the reaching time was calculated as the time interval from the direction change command to the movement at 2 m displacement after the direction change. To compare the switching mobility model with the conventional models, we reproduced the two model simulations in forward walking models (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>; Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>). In the neural oscillator model<sup>13</sup>, we used the cited parameters, motion equations and set the constant input as 6. The time step in the simulation was set to 10<sup>&#x02212;6</sup> s. In the reflex control model<sup>15</sup>, we used the freely available MATLAB code, and set the type of model to normal walk and the time step of extraction to 10<sup>&#x02212;3</sup> s. Other parameters were set as the defaults, including 1.3 m/s as the initial horizontal moving velocity.</p>
<p>For bivariate correlations, we used Pearson&#x00027;s correlation coefficient. For comparing the reaction time between during two different phases, we used the unpaired <italic>t</italic>-test. Both statistics are described with the corresponding degrees of freedom (denoted by a subscript). For all the statistical calculations, <italic>p</italic> &#x0003C; 0.05 was considered significant. All simulations and statistical analyses were performed using MATLAB 2016a Statistics and Machine Learning Toolbox (The MathWorks, Inc., MA, USA).</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec>
<title>Bipedal locomotion with switching mobility control</title>
<p>We first set the target speed to 2 m/s and simulated straightforward locomotion without direction change (Video <xref ref-type="supplementary-material" rid="SM2">S1</xref>). Figure <xref ref-type="fig" rid="F2">2A</xref> shows the time series of the target and actual velocities of the trunk mass. Similar to actual human bipedal locomotion (Bruijn et al., <xref ref-type="bibr" rid="B3">2013</xref>), the target speed was not always achieved because the acceleration of the body can occur only at the moment when the foot is grounded. The trunk velocity is mostly obtained at the time of contact of the trailing leg (Figures <xref ref-type="fig" rid="F2">2A,B</xref>) because the leading (right) leg was not designed to overtake the trailing leg in the sidestep (Yamashita et al., <xref ref-type="bibr" rid="B40">2013</xref>). In particular, in this model, alternate grounding of the leading and the trailing feet did not always occur (Figure <xref ref-type="fig" rid="F2">2B</xref>, Video <xref ref-type="supplementary-material" rid="SM2">S1</xref>). Thus, the locomotion of our model was neither strictly walking nor running. As is the case for skipping with the repeated same foot contact (Minetti, <xref ref-type="bibr" rid="B26">1998</xref>) and galloping with both the double support and flight phase (Yamashita et al., <xref ref-type="bibr" rid="B40">2013</xref>), in the model, the same foot sometimes repeatedly contacted and showed a double support and flight phase. Our model temporally changed the four gaits characteristics, and thus cannot be categorized. The longer term (60 s) characteristics are shown in Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref>. The mobility index of the three actuators (Figure <xref ref-type="fig" rid="F2">2C</xref>) alternately increased and decreased to play their roles as determined by switching coefficients (Figure <xref ref-type="fig" rid="F2">2D</xref>). Mobility seemed to increase in the flight phase in both leg actuators (2 and 3) and at the phase with either leg grounded in inter-leg actuator 1 (Figure <xref ref-type="fig" rid="F2">2C</xref>). The switching coefficients seemed to switch appropriately to propulsion (Figure <xref ref-type="fig" rid="F2">2D</xref> red, e.g., grounded in actuators 2 and 3), balance (blue, e.g., grounded in actuators 2 and 1; in flight in actuator 3), and swing (orange or light blue, in flight in actuator 2 or 3) separately in each actuator (Figure <xref ref-type="fig" rid="F2">2D</xref>: the algorithm is given in Section Materials and Methods).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Straight lateral locomotion and a new method of evaluating switching adaptability. <bold>(A)</bold> Time series of the target speed (black line) and actual speed (green line) of the trunk mass (top component). Target speed was set to 2 m/s and actual speed fluctuated. <bold>(B)</bold> Contact with the ground of the right leading foot (red) and the left following foot (blue). The left and right feet were not always grounded alternately. Detailed processes are shown in Video <xref ref-type="supplementary-material" rid="SM2">S1</xref>. <bold>(C)</bold> Mobility index <italic>k</italic><sub><italic>i</italic></sub> of inter-leg (I, black, <italic>i</italic> &#x0003D; 1), right (R, red, <italic>i</italic> &#x0003D; 2) and left (L, blue, <italic>i</italic> &#x0003D; 3) legs in the switching mobility control algorithm (a darker color indicates greater mobility). Mobility index of the three actuators alternately increased and decreased to play their roles as determined by the switching coefficients <bold>(D)</bold>. The mobility seemed to increase at the flight phase in both leg actuators and at the grounded phase of either leg in the inter-leg actuator. <bold>(D)</bold> Switching coefficient <italic>a</italic><sub><italic>i</italic></sub> of actuator <italic>i</italic>&#x00027;s desired velocity. I, R, and L are the same as in <bold>(C)</bold>. Red (<italic>a</italic><sub><italic>i</italic></sub> &#x0003D; 1) and blue (<italic>a</italic><sub><italic>i</italic></sub> &#x0003D; &#x02212;1) show propulsion and balance during the grounded phase, respectively. Orange (<italic>a</italic><sub><italic>i</italic></sub> &#x0003D; 1/2) and light blue (<italic>a</italic><sub><italic>i</italic></sub> &#x0003D; &#x02212;1/2) show leg swing for propulsion and balance during the flight phase, respectively. Green (<italic>a</italic><sub><italic>i</italic></sub> &#x0003D; 0) is neutral (i.e., zero velocity command) for the grounded state of either or both legs. The switching coefficients seemed to appropriately switch to propulsion (e.g., grounded in actuators 2 and 3), balance (e.g., grounded in actuators 2 and 1; in flight in actuator 3) and swing (in flight in actuator 2 or 3) separately in each actuator. <bold>(E)</bold> Reaching time (vertical axis) toward 2 m after the direction change command with a 0.03-s interval (horizontal axis). Time series corresponds to the timing of the direction change, which shows high variability (reaching time: 2.734 &#x000B1; 0.567 s). The reaching time and its variation are the model performance and a new method of evaluating the switching adaptation in bipedal locomotion.</p></caption>
<graphic xlink:href="fnhum-11-00298-g0002.tif"/>
</fig>
</sec>
<sec>
<title>Direction change at various timings</title>
<p>During the sidestep in Figure <xref ref-type="fig" rid="F2">2</xref>, we switched the target speed to &#x02212;2 m/s at various timings and moved the trunk to 2 m in the opposite direction from that moment. The vertical axis in Figure <xref ref-type="fig" rid="F2">2E</xref> shows the reaching time toward 2 m after the direction change command with a 0.03-s interval (horizontal axis). First, the switching mobility model achieved direction change at any time using only the three actuators. The reaching time and its variation are the model performance (mean and standard deviation of reaching time: 2.734 &#x000B1; 0.567 s) and a new evaluation method of the switching adaptation in bipedal locomotion. The results showed the reaching time increased during the trailing leg stance compared with the other timings (2.940 &#x000B1; 0.644 s vs. 2.604 &#x000B1; 0.474 s, <italic>t</italic><sub>99</sub> &#x0003D; 3.01, <italic>p</italic> &#x0003D; 0.0033). Figure <xref ref-type="fig" rid="F3">3</xref> shows examples of the faster direction change (reaching time: 1.478 s) after 1.861 s from the start of the simulation (Figure <xref ref-type="fig" rid="F3">3B</xref>, Video <xref ref-type="supplementary-material" rid="SM3">S2</xref>) and in the delayed direction change (reaching time: 4.280 s) after 2.161 s (Figure <xref ref-type="fig" rid="F3">3B</xref>, Video <xref ref-type="supplementary-material" rid="SM4">S3</xref>). The faster trial involved a change in direction to switch the mobility index and the switching coefficient and to include fewer steps in a shorter cycle compared with the slower trial. As a coarse grained explanation at the direction change timing, the trailing leg stance increased the reaching time because the trailing leg propelled the body (before the direction change) and then will make the body difficult to change the inverse direction.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Examples of trials in fast and slow reaching upon direction change. Similar to Figure <xref ref-type="fig" rid="F2">2</xref>, this figure shows the actual trunk (green) and target (black) velocities, and the left and right foot contacts in the faster direction change after 1.861 s from the start time <bold>(A)</bold> and in the delayed direction change after 2.161 s <bold>(B)</bold>. The upper right and lower right stick pictures are kinematic postures at the direction change command. The moment of the direction change command and at reaching 2 m (simulation end) are indicated by black dotted lines.</p></caption>
<graphic xlink:href="fnhum-11-00298-g0003.tif"/>
</fig>
<p>As fundamental kinematic characteristics to investigate the fine-grained fluctuation of the reaction time, we examined the relationship of step numbers and foot height with the direction change performance (Figure <xref ref-type="fig" rid="F4">4</xref>). The reaching time was significantly increased with a greater number of steps for both leading and trailing feet (Figure <xref ref-type="fig" rid="F4">4A</xref>, leading: <italic>r</italic><sub>99</sub> &#x0003D; 0.542, <italic>p</italic> &#x0003D; 4.7 &#x000D7; 10<sup>&#x02212;9</sup>, trailing: <italic>r</italic><sub>99</sub> &#x0003D; 0.509, <italic>p</italic> &#x0003D; 5.4 &#x000D7; 10<sup>&#x02212;8</sup>). It was also significantly increased with maximum foot height (Figure <xref ref-type="fig" rid="F4">4B</xref>) for the leading foot (<italic>r</italic><sub>99</sub> &#x0003D; 0.283, <italic>p</italic> &#x0003D; 4.1 &#x000D7; 10<sup>&#x02212;3</sup>), but not that for the trailing foot (<italic>r</italic><sub>99</sub> &#x0003D; 0.03, <italic>p</italic> &#x0003D; 0.976). These results suggest that the faster direction change was derived from the movement with less motor cost, estimated by fewer steps and a smaller leading foot height. However, an underlying cause of difference between the trials in the faster and the slower reaching time was difficult to explain directly because the behaviors were interrelated and generated from the closed-loop structure. This may be generated from a subtle dynamic state difference and the subsequent accumulation of integration error in the non-integrable system.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Relationships between kinematic characteristics and performance. Performance was evaluated as the reaction time after the direction change command (Figure <xref ref-type="fig" rid="F2">2E</xref>). Kinematic characteristics were quantified as <bold>(A)</bold> number of steps and <bold>(B)</bold> maximum vertical foot height of left foot (blue triangle, leading foot) and right foot (red circle, trailing foot) after the direction change. Note that all of the data including foot definition were for after the direction change. <bold>(A)</bold> Reaching time was significantly increased with the number of steps for both leading and trailing feet. <bold>(B)</bold> Reaching time was significantly increased with maximum vertical foot height for the leading foot (blue triangle, left foot), but not that for the trailing foot (red circle, right foot).</p></caption>
<graphic xlink:href="fnhum-11-00298-g0004.tif"/>
</fig>
</sec>
<sec>
<title>Comparison with conventional models</title>
<p>To reveal the difference in motor output in the different architectures, we reproduced two previous forward walk models with a neural oscillator (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>) and reflex control (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>; details are given in Section Materials and Methods). Note that because there is no sidestep bipedal model, except for ours, detailed comparisons between the proposed and previous models are impossible (furthermore, these models have parameter sensitivity, so we cannot match the locomotion velocities). We thus focused on the fundamental locomotion characteristics which differed greatly beyond the mere specific parameters of the models. We first compared the amplitude and variance of the inter-step interval for both legs with a large number of steps (90 stable steps from the 11th to the 100th step) to quantify the stability in locomotion (Figures <xref ref-type="fig" rid="F5">5A&#x02013;C</xref>). The switching mobility model had much larger variance in step interval, despite the lower horizontal moving velocity (leading leg: 0.298 &#x000B1; 0.247 s, trailing leg: 0.257 &#x000B1; 0.190 s at a mean velocity of 0.575 m/s) compared with the forward walking models (neural oscillator model: 1.119 &#x000B1; 0.001 s at a mean velocity of 1.544 m/s, reflex model: 1.229 &#x000B1; 0.004 s at 1.200 m/s; in both models, the intervals for the left and right legs were the same). The step interval variance in the switching mobility model was also much larger than the measured human sidestep walking data (Yamashita et al., <xref ref-type="bibr" rid="B40">2013</xref>; &#x0007E;1.0 &#x000B1; 0.1 s at 1.3 m/s). In addition, from a visual analysis, in the switching mobility model, sidestepping was performed in a less efficient way with higher foot raising (Figure <xref ref-type="fig" rid="F5">5D</xref>) or a repeated grounded phase in the same leg (Figure <xref ref-type="fig" rid="F2">2C</xref>, Video <xref ref-type="supplementary-material" rid="SM2">S1</xref>) at the expense of switching adaptability. We calculated the peak foot height frequency (Figures <xref ref-type="fig" rid="F5">5D&#x02013;F</xref>) and showed that the leading leg was sometimes raised higher than in the comparable models, in which the foot height was strictly controlled, especially in the computational reflex model (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Comparison of step frequency with conventional locomotion models. Time series of step interval of trailing (blue) and leading (red) legs in <bold>(A)</bold> the switching mobility model, <bold>(B)</bold> the neural oscillator model (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>) and <bold>(C)</bold> the reflex control model (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>). Left and center panels show the first 10 steps magnified and 90 subsequent steps, respectively. Note that the models for comparison were forward locomotion models, but the mean moving velocities were largely different among them (the switching mobility model: 0.575 m/s, neural oscillator model: 1.544 m/s, reflex model: 1.200 m/s). The step interval variability in the switching mobility model was much larger than in the conventional locomotion models. We also included signature leg raising in the switching mobility model and show a histogram of peak foot height frequency in <bold>(D)</bold> the switching mobility model, <bold>(E)</bold> the neural oscillator model and <bold>(F)</bold> the reflex control model. The horizontal axis represents the vertical foot height normalized by the trunk height. Comparable models (especially the computational reflex model; Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>) strictly controlled the foot height, but the leading leg in the switching mobility model was sometimes raised higher than in the comparable models.</p></caption>
<graphic xlink:href="fnhum-11-00298-g0005.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>In this study, we constructed a minimal distributed autonomous model achieving bipedal change in direction at any time with only three actuators, but without accurate features representing the whole human body structure, such as a central pattern generator and a lower limb joint, as reproduced in the previous comparable models (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>; Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>), that is, without explicit optimization and cyclic movement in neural oscillators. Simultaneously, we also developed a new method for evaluating the ability to perform robustly faster direction change during bipedal locomotion (Figure <xref ref-type="fig" rid="F2">2E</xref>). In previous studies of adaptive bipedal locomotion in robotic engineering or neurophysiology, such as when encountering an obstacle in the environment or an external disturbance (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>; Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>; Koolen et al., <xref ref-type="bibr" rid="B21">2016</xref>; Kuindersma et al., <xref ref-type="bibr" rid="B22">2016</xref>), it was not considered whether the bipedal model can achieve direction change at any time, which is not explicitly implemented (i.e., unpredictable for the model). Thus, to discuss this control problem, it is necessary to reconstruct the frameworks from the viewpoints of engineering control and neurophysiology, as explained below. The expansion of the frameworks should be necessary for understanding universal motor control.</p>
<p>From the viewpoint of engineering control, the switching mobility model showed switching adaptability at the expense of efficiency because it is difficult for the bipedal locomotion model to satisfy the criteria of both efficiency and adaptability. Previous research using an inverted pendulum locomotion model (Srinivasan and Ruina, <xref ref-type="bibr" rid="B33">2006</xref>) explained the energy efficiency in various locomotion patterns based on the optimal control theory. With regard to efficiency, the switching mobility model showed highly costly movement in which a period of repeated grounding of the same leg occurred (Figure <xref ref-type="fig" rid="F3">3</xref>), in contrast to that in the conventional models (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>; Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>), even without a change in direction. Furthermore, at the time of changing direction, it sometimes takes a long time to adapt to the command to change direction (Figure <xref ref-type="fig" rid="F3">3B</xref>) because of a greater number of steps (Figure <xref ref-type="fig" rid="F4">4A</xref>) and a higher leading leg raise (Figure <xref ref-type="fig" rid="F4">4B</xref>). This leg raise was controlled in the neural oscillator model (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>) and strictly computed in the reflex model (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>) at a lower level than in the switching mobility model without direction change (Figure <xref ref-type="fig" rid="F5">5</xref>). In adaptive bipedal locomotion, it is difficult to calculate the optimal trajectory, so we confirmed the higher motor cost in the switching mobility model, estimated from the much larger variance in step interval and much higher leg raise, than in the conventional models. This is also because there is no mathematical guarantee of efficient movement in mobility control (Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>), which is different from the explicit optimal control or neural oscillator control in cyclic movement. Instead, switching mobility control in the switching mobility model would have an advantage regarding adaptability, even if the movement switches from a cyclic to a non-cyclic pattern and vice versa and has difficulty in control, such as large instability as in bipedal chase-and-escape in interactive sports (Fujii et al., <xref ref-type="bibr" rid="B10">2015b</xref>,<xref ref-type="bibr" rid="B13">c</xref>). We were not able to perform a direct comparison with these control algorithms because there are no models to show the switching adaptability (and performing a sidestep); however, a future model can be compared by the proposed evaluation method in which the model performs direction change at various timings.</p>
<p>Neurophysiologically, our switching mobility control algorithm suggests the presence of reflex-like switching functions of propulsion, balancing, and leg swing within and between limbs to achieve the task. The algorithm does not directly reflect the neural mechanism, but we can consider similarity with human neurophysiology by a process of elimination. The proposed model does not explicitly control the actuator movements like cerebellum (Shadmehr and Krakauer, <xref ref-type="bibr" rid="B30">2008</xref>) and does not directly generate quasiperiodic movements such as using neural oscillators (Taga et al., <xref ref-type="bibr" rid="B34">1991</xref>) physiologically located in spinal central pattern generators (Grillner, <xref ref-type="bibr" rid="B15">1985</xref>; Dimitrijevic et al., <xref ref-type="bibr" rid="B7">1998</xref>). Instead, the previous work suggests that diverse movements can be reproduced by incorporating the multiple reflexes in the spinal feedback circuitry without the central pattern generator (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>). Among these architectures, we can find the similarity between the reflex control and the autonomous distributed control (Yoshihara et al., <xref ref-type="bibr" rid="B41">2007</xref>; Watanabe et al., <xref ref-type="bibr" rid="B38">2012</xref>), which incorporates interaction among components and environments (i.e., feedback of ground reaction force) and can implicitly adapt to a rapid change of command. In other words, the point at which each actuator autonomously sets and executes a target movement according to the situation in the switching mobility model matches that in the physiological reflex mechanism in which the inputs and outputs are locally automated. Our switching distributed autonomous control algorithm provided the relationship among the actuators (i.e., muscles) to switch functions of propulsion, balancing, and leg swing within and between limb levels. In a manner similar to long-latency reflexes, which would possess an internal model of limb dynamics (Kurtzer et al., <xref ref-type="bibr" rid="B23">2008</xref>), the actuator in the switching mobility model rapidly switched functions according to the situation. There is no physiological evidence from the neural circuitry of a response of such inter-limb reflexes to a drastic change of command in the long neural pathway, but we believe that it would be needed for the switching adaptation in bipedal locomotion. The neural mechanism involved would be complicated because of possible involvement with both voluntary and reflex control overlapping in their neural substrates (Kurtzer et al., <xref ref-type="bibr" rid="B23">2008</xref>); however, it may be simply explained by the simple interaction rule of our switching mobility control. This constructive approach can contribute to understanding intelligent motor control including biologically (Carvalho et al., <xref ref-type="bibr" rid="B4">2012</xref>) and socially (Helbing et al., <xref ref-type="bibr" rid="B17">2000</xref>) essential activities, such as escape from enemies, pursuit of prey, and search for food.</p>
<p>However, there are some problems with the above neurophysiological claims. One is that it claims to be based only on the similarity in the architectures without neurophysiological evidences. This is considered as a general problem in finding evidences of long-latency reflexes, which overlapped with voluntary movements in their neural substrates (Kurtzer et al., <xref ref-type="bibr" rid="B23">2008</xref>). Second is validation with real-world human data satisfying both sidestepping (Yamashita et al., <xref ref-type="bibr" rid="B40">2013</xref>) and unpredicted change direction (Fujii et al., <xref ref-type="bibr" rid="B10">2015b</xref>). The latter study showed that the unpredicted competitive situation (i.e., requiring the faster movement) delayed the first step initiation more than 100 ms, suggesting the human can decrease the delay more than the current switching adaptability model. The real human locomotion mechanism includes central pattern generators, peripheral reflexes, and spring-damper system. The spring-damper property, which contributes high speed locomotion according to a quadruped robot study (Kimura et al., <xref ref-type="bibr" rid="B20">2007</xref>), should be further investigated.</p>
<p>Third is the sensitivity of the simulation to the choice of some of the model parameters. Our supplementary results (Figure <xref ref-type="supplementary-material" rid="SM1">S3</xref>) showed that the switching adaptability model had strong sensitivity to the parameters. We suppose the parameterization may be related with inherent adaptation to the individual musculoskeletal system and might be relatively independent of the motor control adaptation. As a further alternative approach, for example, evolutionary algorithm (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>) and reinforcement learning (Lillicrap et al., <xref ref-type="bibr" rid="B24">2015</xref>) efficiently worked in previous bipedal locomotion studies. These are complementary relationship in terms of inter- and intra-generation progress, respectively. The reflex model (Song and Geyer, <xref ref-type="bibr" rid="B32">2015</xref>) optimized the control parameters with the covariance matrix adaptation evolution strategy (Hansen, <xref ref-type="bibr" rid="B16">2006</xref>). Reinforcement learning enables the acquisition of efficient and adaptive locomotion by trial and error on models, like humans actually do, probably in the basal ganglia (Doya, <xref ref-type="bibr" rid="B8">2000</xref>). However, if learning a case of sudden change through evolutionary algorithm and reinforcement learning (also in learning of humans as an experimental condition; Fujii et al., <xref ref-type="bibr" rid="B12">2013</xref>), there is a possibility that it will not mean &#x0201C;unpredictable sudden change&#x0201D; when performing an adaptive bipedal locomotion (Shinya et al., <xref ref-type="bibr" rid="B31">2009</xref>). In other words, it can be considered that the explicit control rule should be difficult to estimate in principle in these approaches (i.e., human bipedal locomotion can also predictively adapt to the situation; Shinya et al., <xref ref-type="bibr" rid="B31">2009</xref>). The evolutionary algorithm and reinforcement learning will also have the advantage to acquire the efficient movement in our less efficient model. Satisfying the requirements of efficiency (e.g., obtained explicitly by using optimal control and implicitly by evolutionary algorithm or reinforcement learning) and switching adaptability proposed by our study will further contribute to solving the general motor control problem.</p>
</sec>
<sec id="s5">
<title>Author contributions</title>
<p>KF conceived the original idea of the model. KF and YYo designed the model. KF, YYo, HT, and YYa analyzed data and wrote the paper.</p>
<sec>
<title>Conflict of interest statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</sec>
</body>
<back>
<ack><p>The authors would like to thank M. Yano of Tohoku University, S. Hagio of Tokyo University and D. Yamashita of Japan Institute of Sports Sciences for useful discussions and S. Song of Carnegie Mellon University for his walking model.</p>
</ack>
<sec sec-type="supplementary-material" id="s6">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="http://journal.frontiersin.org/article/10.3389/fnhum.2017.00298/full#supplementary-material">http://journal.frontiersin.org/article/10.3389/fnhum.2017.00298/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Presentation1.PDF" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Video1.MP4" id="SM2" mimetype="video/mp4" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Video2.MP4" id="SM3" mimetype="video/mp4" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Video3.MP4" id="SM4" mimetype="video/mp4" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>F. C.</given-names></name> <name><surname>Pandy</surname> <given-names>M. G.</given-names></name></person-group> (<year>2003</year>). <article-title>Individual muscle contributions to support in normal walking</article-title>. <source>Gait Posture</source> <volume>17</volume>, <fpage>159</fpage>&#x02013;<lpage>169</lpage>. <pub-id pub-id-type="doi">10.1016/S0966-6362(02)00073-5</pub-id><pub-id pub-id-type="pmid">12633777</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bongard</surname> <given-names>J.</given-names></name> <name><surname>Zykov</surname> <given-names>V.</given-names></name> <name><surname>Lipson</surname> <given-names>H.</given-names></name></person-group> (<year>2006</year>). <article-title>Resilient machines through continuous self-modeling</article-title>. <source>Science</source> <volume>314</volume>, <fpage>1118</fpage>&#x02013;<lpage>1121</lpage>. <pub-id pub-id-type="doi">10.1126/science.1133687</pub-id><pub-id pub-id-type="pmid">17110570</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bruijn</surname> <given-names>S. M.</given-names></name> <name><surname>Meijer</surname> <given-names>O. G.</given-names></name> <name><surname>Beek</surname> <given-names>P. J.</given-names></name> <name><surname>van Die&#x000EB;n</surname> <given-names>J. H.</given-names></name></person-group> (<year>2013</year>). <article-title>Assessing the stability of human locomotion: a review of current measures</article-title>. <source>J. R. Soc. Interface</source> <volume>10</volume>:<fpage>20120999</fpage>. <pub-id pub-id-type="doi">10.1098/rsif.2012.0999</pub-id><pub-id pub-id-type="pmid">23516062</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carvalho</surname> <given-names>S.</given-names></name> <name><surname>Biro</surname> <given-names>D.</given-names></name> <name><surname>Cunha</surname> <given-names>E.</given-names></name> <name><surname>Hockings</surname> <given-names>K.</given-names></name> <name><surname>McGrew</surname> <given-names>W. C.</given-names></name> <name><surname>Richmond</surname> <given-names>B. G.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Chimpanzee carrying behaviour and the origins of human bipedality</article-title>. <source>Curr. Biol.</source> <volume>22</volume>, <fpage>R180</fpage>&#x02013;<lpage>R181</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2012.01.052</pub-id><pub-id pub-id-type="pmid">22440797</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Couzin</surname> <given-names>I. D.</given-names></name> <name><surname>Krause</surname> <given-names>J.</given-names></name> <name><surname>James</surname> <given-names>R.</given-names></name> <name><surname>Ruxton</surname> <given-names>G. D.</given-names></name> <name><surname>Franks</surname> <given-names>N. R.</given-names></name></person-group> (<year>2002</year>). <article-title>Collective memory and spatial sorting in animal groups</article-title>. <source>J. Theor. Biol.</source> <volume>218</volume>, <fpage>1</fpage>&#x02013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1006/jtbi.2002.3065</pub-id><pub-id pub-id-type="pmid">12297066</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cully</surname> <given-names>A.</given-names></name> <name><surname>Clune</surname> <given-names>J.</given-names></name> <name><surname>Tarapore</surname> <given-names>D.</given-names></name> <name><surname>Mouret</surname> <given-names>J.-B.</given-names></name></person-group> (<year>2015</year>). <article-title>Robots that can adapt like animals</article-title>. <source>Nature</source> <volume>521</volume>, <fpage>503</fpage>&#x02013;<lpage>507</lpage>. <pub-id pub-id-type="doi">10.1038/nature14422</pub-id><pub-id pub-id-type="pmid">26017452</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dimitrijevic</surname> <given-names>M. R.</given-names></name> <name><surname>Gerasimenko</surname> <given-names>Y.</given-names></name> <name><surname>Pinter</surname> <given-names>M. M.</given-names></name></person-group> (<year>1998</year>). <article-title>Evidence for a spinal central pattern generator in humans</article-title>. <source>Ann. N. Y. Acad. Sci.</source> <volume>860</volume>, <fpage>360</fpage>&#x02013;<lpage>376</lpage>. <pub-id pub-id-type="pmid">9928325</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Doya</surname> <given-names>K.</given-names></name></person-group> (<year>2000</year>). <article-title>Complementary roles of basal ganglia and cerebellum in learning and motor control</article-title>. <source>Curr. Opin. Neurobiol.</source> <volume>10</volume>, <fpage>732</fpage>&#x02013;<lpage>739</lpage>. <pub-id pub-id-type="doi">10.1016/S0959-4388(00)00153-7</pub-id><pub-id pub-id-type="pmid">11240282</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujii</surname> <given-names>K.</given-names></name> <name><surname>Daichi</surname> <given-names>Y.</given-names></name> <name><surname>Tetsuya</surname> <given-names>K.</given-names></name> <name><surname>Isaka</surname> <given-names>T.</given-names></name> <name><surname>Kouzaki</surname> <given-names>M.</given-names></name></person-group> (<year>2015a</year>). <article-title>Preparatory body state before reacting to an opponent: short-term joint torque fluctuation in real-time competitive sports</article-title>. <source>PLoS ONE</source> <volume>10</volume>:<fpage>e0128571</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0128571</pub-id><pub-id pub-id-type="pmid">26024485</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujii</surname> <given-names>K.</given-names></name> <name><surname>Isaka</surname> <given-names>T.</given-names></name> <name><surname>Kouzaki</surname> <given-names>M.</given-names></name> <name><surname>Yamamoto</surname> <given-names>Y.</given-names></name></person-group> (<year>2015b</year>). <article-title>Mutual and asynchronous anticipation and action in sports as globally competitive and locally coordinative dynamics</article-title>. <source>Sci. Rep.</source> <volume>5</volume>:<fpage>16140</fpage>. <pub-id pub-id-type="doi">10.1038/srep16140</pub-id><pub-id pub-id-type="pmid">26538452</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujii</surname> <given-names>K.</given-names></name> <name><surname>Yokoyama</surname> <given-names>K.</given-names></name> <name><surname>Koyama</surname> <given-names>T.</given-names></name> <name><surname>Rikukawa</surname> <given-names>A.</given-names></name> <name><surname>Yamada</surname> <given-names>H.</given-names></name> <name><surname>Yamamoto</surname> <given-names>Y.</given-names></name></person-group> (<year>2016</year>). <article-title>Resilient help to switch and overlap hierarchical subsystems in a small human group</article-title>. <source>Sci. Rep.</source> <volume>6</volume>:<fpage>23911</fpage>. <pub-id pub-id-type="doi">10.1038/srep23911</pub-id><pub-id pub-id-type="pmid">27045443</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujii</surname> <given-names>K.</given-names></name> <name><surname>Yoshioka</surname> <given-names>S.</given-names></name> <name><surname>Isaka</surname> <given-names>T.</given-names></name> <name><surname>Kouzaki</surname> <given-names>M.</given-names></name></person-group> (<year>2013</year>). <article-title>Unweighted state as a sidestep preparation improve the initiation and reaching performance for basketball players</article-title>. <source>J. Electromyogr. Kinesiol.</source> <volume>23</volume>, <fpage>1467</fpage>&#x02013;<lpage>1473</lpage>. <pub-id pub-id-type="doi">10.1016/j.jelekin.2013.08.001</pub-id><pub-id pub-id-type="pmid">24060389</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujii</surname> <given-names>K.</given-names></name> <name><surname>Yoshioka</surname> <given-names>S.</given-names></name> <name><surname>Isaka</surname> <given-names>T.</given-names></name> <name><surname>Kouzaki</surname> <given-names>M.</given-names></name></person-group> (<year>2015c</year>). <article-title>The preparatory state of ground reaction forces in defending against a dribbler in a basketball 1-on-1 dribble subphase</article-title>. <source>Sports Biomech.</source> <volume>14</volume>, <fpage>28</fpage>&#x02013;<lpage>44</lpage>. <pub-id pub-id-type="doi">10.1080/14763141.2015.1026931</pub-id><pub-id pub-id-type="pmid">25895702</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Golubitsky</surname> <given-names>M.</given-names></name> <name><surname>Stewart</surname> <given-names>I.</given-names></name> <name><surname>Buono</surname> <given-names>P. L.</given-names></name> <name><surname>Collins</surname> <given-names>J. J.</given-names></name></person-group> (<year>1999</year>). <article-title>Symmetry in locomotor central pattern generators and animal gaits</article-title>. <source>Nature</source> <volume>401</volume>, <fpage>693</fpage>&#x02013;<lpage>695</lpage>. <pub-id pub-id-type="doi">10.1038/44416</pub-id><pub-id pub-id-type="pmid">10537106</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grillner</surname> <given-names>S.</given-names></name></person-group> (<year>1985</year>). <article-title>Neurobiological bases of rhythmic motor acts in vertebrates</article-title>. <source>Science</source> <volume>228</volume>, <fpage>143</fpage>&#x02013;<lpage>148</lpage>. <pub-id pub-id-type="doi">10.1126/science.3975635</pub-id><pub-id pub-id-type="pmid">3975635</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hansen</surname> <given-names>N.</given-names></name></person-group> (<year>2006</year>). <article-title>The CMA evolution strategy: a comparing review</article-title>, in <source>Towards a New Evolutionary Computation</source>, eds <person-group person-group-type="editor"><name><surname>Lozano</surname> <given-names>J. A.</given-names></name> <name><surname>Larra&#x000F1;aga</surname> <given-names>P.</given-names></name> <name><surname>Inza</surname> <given-names>I.</given-names></name> <name><surname>Bengoetxea</surname> <given-names>E.</given-names></name></person-group> (<publisher-loc>Berlin; Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>75</fpage>&#x02013;<lpage>102</lpage>.</citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Helbing</surname> <given-names>D.</given-names></name> <name><surname>Farkas</surname> <given-names>I.</given-names></name> <name><surname>Vicsek</surname> <given-names>T.</given-names></name></person-group> (<year>2000</year>). <article-title>Simulating dynamical features of escape panic</article-title>. <source>Nature</source> <volume>407</volume>, <fpage>487</fpage>&#x02013;<lpage>490</lpage>. <pub-id pub-id-type="doi">10.1038/35035023</pub-id><pub-id pub-id-type="pmid">11028994</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hutchins</surname> <given-names>E.</given-names></name></person-group> (<year>1995</year>). <source>Cognition in the Wild</source>. <publisher-loc>Cambridge, MA</publisher-loc>: <publisher-name>MIT Press</publisher-name>.</citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kamimura</surname> <given-names>A.</given-names></name> <name><surname>Ohira</surname> <given-names>T.</given-names></name></person-group> (<year>2010</year>). <article-title>Group chase and escape</article-title>. <source>New J. Phys.</source> <volume>12</volume>:<fpage>053013</fpage>. <pub-id pub-id-type="doi">10.1088/1367-2630/12/5/053013</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kimura</surname> <given-names>H.</given-names></name> <name><surname>Fukuoka</surname> <given-names>Y.</given-names></name> <name><surname>Cohen</surname> <given-names>A. H.</given-names></name></person-group> (<year>2007</year>). <article-title>Biologically inspired adaptive walking of a quadruped robot</article-title>. <source>Philos. Trans. R. Soc. Lond. A Math. Phys. Eng. Sci.</source> <volume>365</volume>, <fpage>153</fpage>&#x02013;<lpage>170</lpage>. <pub-id pub-id-type="doi">10.1098/rsta.2006.1919</pub-id><pub-id pub-id-type="pmid">17148054</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Koolen</surname> <given-names>T.</given-names></name> <name><surname>Bertrand</surname> <given-names>S.</given-names></name> <name><surname>Thomas</surname> <given-names>G.</given-names></name> <name><surname>De Boer</surname> <given-names>T.</given-names></name> <name><surname>Wu</surname> <given-names>T. F.</given-names></name> <name><surname>Smith</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Design of a momentum-based control framework and application to the humanoid robot atlas</article-title>. <source>Int. J. Hum. Rob.</source> <volume>13</volume>:<fpage>1650007</fpage>. <pub-id pub-id-type="doi">10.1142/S0219843616500079</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kuindersma</surname> <given-names>S.</given-names></name> <name><surname>Deits</surname> <given-names>R.</given-names></name> <name><surname>Fallon</surname> <given-names>M.</given-names></name> <name><surname>Valenzuela</surname> <given-names>A.</given-names></name> <name><surname>Dai</surname> <given-names>H. K.</given-names></name> <name><surname>Permenter</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Optimization-based locomotion planning, estimation, and control design for the atlas humanoid robot</article-title>. <source>Auton. Robots</source> <volume>40</volume>, <fpage>429</fpage>&#x02013;<lpage>455</lpage>. <pub-id pub-id-type="doi">10.1007/s10514-015-9479-3</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kurtzer</surname> <given-names>I. L.</given-names></name> <name><surname>Pruszynski</surname> <given-names>J. A.</given-names></name> <name><surname>Scott</surname> <given-names>S. H.</given-names></name></person-group> (<year>2008</year>). <article-title>Long-latency reflexes of the human arm reflect an internal model of limb dynamics</article-title>. <source>Curr. Biol.</source> <volume>18</volume>, <fpage>449</fpage>&#x02013;<lpage>453</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2008.02.053</pub-id><pub-id pub-id-type="pmid">18356051</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lillicrap</surname> <given-names>T. P.</given-names></name> <name><surname>Hunt</surname> <given-names>J. J.</given-names></name> <name><surname>Pritzel</surname> <given-names>A.</given-names></name> <name><surname>Heess</surname> <given-names>N.</given-names></name> <name><surname>Erez</surname> <given-names>T.</given-names></name> <name><surname>Tassa</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Continuous control with deep reinforcement learning</article-title>. <source>arXiv preprint arXiv:1509.02971</source>.</citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>McGeer</surname> <given-names>T.</given-names></name></person-group> (<year>1990</year>). <article-title>Passive dynamic walking</article-title>. <source>Int. J. Rob. Res.</source> <volume>9</volume>, <fpage>62</fpage>&#x02013;<lpage>82</lpage>. <pub-id pub-id-type="doi">10.1177/027836499000900206</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Minetti</surname> <given-names>A. E.</given-names></name></person-group> (<year>1998</year>). <article-title>The biomechanics of skipping gaits: a third locomotion paradigm?</article-title> <source>Proc. Biol. Sci.</source> <volume>265</volume>, <fpage>1227</fpage>&#x02013;<lpage>1235</lpage>. <pub-id pub-id-type="doi">10.1098/rspb.1998.0424</pub-id><pub-id pub-id-type="pmid">9699315</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mnih</surname> <given-names>V.</given-names></name> <name><surname>Kavukcuoglu</surname> <given-names>K.</given-names></name> <name><surname>Silver</surname> <given-names>D.</given-names></name> <name><surname>Rusu</surname> <given-names>A. A.</given-names></name> <name><surname>Veness</surname> <given-names>J.</given-names></name> <name><surname>Bellemare</surname> <given-names>M. G.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Human-level control through deep reinforcement learning</article-title>. <source>Nature</source> <volume>518</volume>, <fpage>529</fpage>&#x02013;<lpage>533</lpage>. <pub-id pub-id-type="doi">10.1038/nature14236</pub-id><pub-id pub-id-type="pmid">25719670</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mori</surname> <given-names>S.</given-names></name> <name><surname>Matsui</surname> <given-names>T.</given-names></name> <name><surname>Kuze</surname> <given-names>B.</given-names></name> <name><surname>Asanome</surname> <given-names>M.</given-names></name> <name><surname>Nakajima</surname> <given-names>K.</given-names></name> <name><surname>Matsuyama</surname> <given-names>K.</given-names></name></person-group> (<year>1998</year>). <article-title>Cerebellar-induced locomotion: reticulospinal control of spinal rhythm generating mechanism in cats</article-title>. <source>Ann. N. Y. Acad. Sci</source>. <volume>860</volume>, <fpage>94</fpage>&#x02013;<lpage>105</lpage>. <pub-id pub-id-type="doi">10.1111/j.1749-6632.1998.tb09041.x</pub-id><pub-id pub-id-type="pmid">9928304</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scott</surname> <given-names>S. H.</given-names></name></person-group> (<year>2004</year>). <article-title>Optimal feedback control and the neural basis of volitional motor control</article-title>. <source>Nat. Rev. Neurosci.</source> <volume>5</volume>, <fpage>534</fpage>&#x02013;<lpage>546</lpage>. <pub-id pub-id-type="doi">10.1038/nrn1427</pub-id><pub-id pub-id-type="pmid">15208695</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shadmehr</surname> <given-names>R.</given-names></name> <name><surname>Krakauer</surname> <given-names>J. W.</given-names></name></person-group> (<year>2008</year>). <article-title>A computational neuroanatomy for motor control</article-title>. <source>Exp. Brain Res.</source> <volume>185</volume>, <fpage>359</fpage>&#x02013;<lpage>381</lpage>. <pub-id pub-id-type="doi">10.1007/s00221-008-1280-5</pub-id><pub-id pub-id-type="pmid">18251019</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shinya</surname> <given-names>M.</given-names></name> <name><surname>Fujii</surname> <given-names>S.</given-names></name> <name><surname>Oda</surname> <given-names>S.</given-names></name></person-group> (<year>2009</year>). <article-title>Corrective postural responses evoked by completely unexpected loss of ground support during human walking</article-title>. <source>Gait Posture</source> <volume>29</volume>, <fpage>483</fpage>&#x02013;<lpage>487</lpage>. <pub-id pub-id-type="doi">10.1016/j.gaitpost.2008.11.009</pub-id><pub-id pub-id-type="pmid">19128973</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>S.</given-names></name> <name><surname>Geyer</surname> <given-names>H.</given-names></name></person-group> (<year>2015</year>). <article-title>A neural circuitry that emphasizes spinal feedback generates diverse behaviours of human locomotion</article-title>. <source>J. Physiol. Lond.</source> <volume>593</volume>, <fpage>3493</fpage>&#x02013;<lpage>3511</lpage>. <pub-id pub-id-type="doi">10.1113/JP270228</pub-id><pub-id pub-id-type="pmid">25920414</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Srinivasan</surname> <given-names>M.</given-names></name> <name><surname>Ruina</surname> <given-names>A.</given-names></name></person-group> (<year>2006</year>). <article-title>Computer optimization of a minimal biped model discovers walking and running</article-title>. <source>Nature</source> <volume>439</volume>, <fpage>72</fpage>&#x02013;<lpage>75</lpage>. <pub-id pub-id-type="doi">10.1038/nature04113</pub-id><pub-id pub-id-type="pmid">16155564</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taga</surname> <given-names>G.</given-names></name> <name><surname>Yamaguchi</surname> <given-names>Y.</given-names></name> <name><surname>Shimizu</surname> <given-names>H.</given-names></name></person-group> (<year>1991</year>). <article-title>Self-organized control of bipedal locomotion by neural oscillators in unpredictable environment</article-title>. <source>Biol. Cybern.</source> <volume>65</volume>, <fpage>147</fpage>&#x02013;<lpage>159</lpage>. <pub-id pub-id-type="doi">10.1007/BF00198086</pub-id><pub-id pub-id-type="pmid">1912008</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thelen</surname> <given-names>E.</given-names></name> <name><surname>Ulrich</surname> <given-names>B. D.</given-names></name> <name><surname>Niles</surname> <given-names>D.</given-names></name></person-group> (<year>1987</year>). <article-title>Bilateral coordination in human infants - stepping on a split-belt treadmill</article-title>. <source>J. Exp. Psychol. Hum. Percept. Perform.</source> <volume>13</volume>, <fpage>405</fpage>&#x02013;<lpage>410</lpage>. <pub-id pub-id-type="doi">10.1037/0096-1523.13.3.405</pub-id><pub-id pub-id-type="pmid">2958589</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Todorov</surname> <given-names>E.</given-names></name> <name><surname>Jordan</surname> <given-names>M. I.</given-names></name></person-group> (<year>2002</year>). <article-title>Optimal feedback control as a theory of motor coordination</article-title>. <source>Nat. Neurosci.</source> <volume>5</volume>, <fpage>1226</fpage>&#x02013;<lpage>1235</lpage>. <pub-id pub-id-type="doi">10.1038/nn963</pub-id><pub-id pub-id-type="pmid">12404008</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Uno</surname> <given-names>Y.</given-names></name> <name><surname>Kawato</surname> <given-names>M.</given-names></name> <name><surname>Suzuki</surname> <given-names>R.</given-names></name></person-group> (<year>1989</year>). <article-title>Formation and control of optimal trajectory in human multijoint arm movement - minimum torque-change model</article-title>. <source>Biol. Cybern.</source> <volume>61</volume>, <fpage>89</fpage>&#x02013;<lpage>101</lpage>. <pub-id pub-id-type="doi">10.1007/BF00204593</pub-id><pub-id pub-id-type="pmid">2742921</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Watanabe</surname> <given-names>W.</given-names></name> <name><surname>Kano</surname> <given-names>T.</given-names></name> <name><surname>Suzuki</surname> <given-names>S.</given-names></name> <name><surname>Ishiguro</surname> <given-names>A.</given-names></name></person-group> (<year>2012</year>). <article-title>A decentralized control scheme for orchestrating versatile arm movements in ophiuroid omnidirectional locomotion</article-title>. <source>J. R. Soc. Interface</source> <volume>9</volume>, <fpage>102</fpage>&#x02013;<lpage>109</lpage>. <pub-id pub-id-type="doi">10.1098/rsif.2011.0317</pub-id><pub-id pub-id-type="pmid">21775323</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yamamoto</surname> <given-names>Y.</given-names></name> <name><surname>Yokoyama</surname> <given-names>K.</given-names></name> <name><surname>Okumura</surname> <given-names>M.</given-names></name> <name><surname>Kijima</surname> <given-names>A.</given-names></name> <name><surname>Kadota</surname> <given-names>K.</given-names></name> <name><surname>Gohara</surname> <given-names>K.</given-names></name></person-group> (<year>2013</year>). <article-title>Joint action syntax in Japanese martial arts</article-title>. <source>PLoS ONE</source> <volume>8</volume>:<fpage>e72436</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0072436</pub-id><pub-id pub-id-type="pmid">24023740</pub-id></citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yamashita</surname> <given-names>D.</given-names></name> <name><surname>Shinya</surname> <given-names>M.</given-names></name> <name><surname>Fujii</surname> <given-names>K.</given-names></name> <name><surname>Oda</surname> <given-names>S.</given-names></name> <name><surname>Kouzaki</surname> <given-names>M.</given-names></name></person-group> (<year>2013</year>). <article-title>Walk-, run- and gallop-like gait patterns in human sideways locomotion</article-title>. <source>J. Electromyogr. Kinesiol.</source> <volume>23</volume>, <fpage>1480</fpage>&#x02013;<lpage>1484</lpage>. <pub-id pub-id-type="doi">10.1016/j.jelekin.2013.08.005</pub-id><pub-id pub-id-type="pmid">24055531</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yoshihara</surname> <given-names>Y.</given-names></name> <name><surname>Tomita</surname> <given-names>N.</given-names></name> <name><surname>Makino</surname> <given-names>Y.</given-names></name> <name><surname>Yano</surname> <given-names>M.</given-names></name></person-group> (<year>2007</year>). <article-title>Autonomous control of reaching movement by &#x0201C;mobility measure&#x0201D;</article-title>. <source>Int. J. Rob. Mech.</source> <volume>19</volume>, <fpage>448</fpage>&#x02013;<lpage>458</lpage>. <pub-id pub-id-type="doi">10.20965/jrm.2007.p0448</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> This work was supported by a Grant-in-Aid for JSPS fellows 26-407 and Exploratory Research 16K12995. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p>
</fn>
</fn-group>
</back>
</article>