<?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. Robot. AI</journal-id>
<journal-title>Frontiers in Robotics and AI</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Robot. AI</abbrev-journal-title>
<issn pub-type="epub">2296-9144</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/frobt.2020.00108</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Robotics and AI</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>A Novel User Control for Lower Extremity Rehabilitation Exoskeletons</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Karunakaran</surname> <given-names>Kiran K.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/967786/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Abbruzzese</surname> <given-names>Kevin</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1056813/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Androwis</surname> <given-names>Ghaith</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/375621/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Foulds</surname> <given-names>Richard A.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/152380/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Kessler Foundation</institution>, <addr-line>West Orange</addr-line>, <addr-line>NJ</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Biomedical Engineering</institution>, <addr-line>New Jersey Institute of Technology</addr-line>, <addr-line>Newark</addr-line>, <addr-line>NJ</addr-line>, <country>United States</country></aff>
<aff id="aff3"><sup>3</sup><institution>Stryker Corporation</institution>, <addr-line>Mahwah</addr-line>, <addr-line>NJ</addr-line>, <country>United States</country></aff>
<aff id="aff4"><sup>4</sup><institution>Really Useful Robots</institution>, <addr-line>LLC</addr-line>, <addr-line>Langhorne</addr-line>, <addr-line>PA</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Carlos A. Cifuentes, Escuela Colombiana de Ingenieria Julio Garavito, Colombia</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Guillermo As&#x000ED;n-Prieto, Consejo Superior de Investigaciones Cient&#x000ED;ficas (CSIC), Spain; Camilla Pierella, &#x000C9;cole Polytechnique F&#x000E9;d&#x000E9;rale de Lausanne, Switzerland</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Kiran K. Karunakaran  <email>kkk7&#x00040;njit.edu</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Biomedical Robotics, a section of the journal Frontiers in Robotics and AI</p></fn></author-notes>
<pub-date pub-type="epub">
<day>08</day>
<month>09</month>
<year>2020</year>
</pub-date>
<pub-date pub-type="collection">
<year>2020</year>
</pub-date>
<volume>7</volume>
<elocation-id>108</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>05</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>13</day>
<month>07</month>
<year>2020</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2020 Karunakaran, Abbruzzese, Androwis and Foulds.</copyright-statement>
<copyright-year>2020</copyright-year>
<copyright-holder>Karunakaran, Abbruzzese, Androwis and Foulds</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) and the copyright owner(s) 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>Lower extremity exoskeletons offer the potential to restore ambulation to individuals with paraplegia due to spinal cord injury. However, they often rely on preprogrammed gait, initiated by switches, sensors, and/or EEG triggers. Users can exercise only limited independent control over the trajectory of the feet, the speed of walking, and the placement of feet to avoid obstacles. In this paper, we introduce and evaluate a novel approach that naturally decodes a neuromuscular surrogate for a user&#x00027;s neutrally planned foot control, uses the exoskeleton&#x00027;s motors to move the user&#x00027;s legs in real-time, and provides sensory feedback to the user allowing real-time sensation and path correction resulting in gait similar to biological ambulation. Users express their desired gait by applying Cartesian forces via their hands to rigid trekking poles that are connected to the exoskeleton feet through multi-axis force sensors. Using admittance control, the forces applied by the hands are converted into desired foot positions, every 10 milliseconds (ms), to which the exoskeleton is moved by its motors. As the trekking poles reflect the resulting foot movement, users receive sensory feedback of foot kinematics and ground contact that allows on-the-fly force corrections to maintain the desired foot behavior. We present preliminary results showing that our novel control can allow users to produce biologically similar exoskeleton gait.</p></abstract>
<kwd-group>
<kwd>lower extremity exoskeletons</kwd>
<kwd>gait</kwd>
<kwd>spinal cord injury</kwd>
<kwd>rehabilitation robotics</kwd>
<kwd>robot control systems</kwd>
</kwd-group>
<counts>
<fig-count count="8"/>
<table-count count="3"/>
<equation-count count="14"/>
<ref-count count="51"/>
<page-count count="14"/>
<word-count count="9572"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Individuals with complete paraplegia due to spinal cord injury (SCI) have impaired motor control and sensory feedback that limits their ability to walk (Shepherd Center). While wheelchairs provide alternative mobility to individuals with paraplegia, they are not a complete substitute for natural ambulation. Current research has addressed this issue with wearable lower extremity exoskeletons (Dollar and Herr, <xref ref-type="bibr" rid="B10">2008</xref>; Contreras-Vidal et al., <xref ref-type="bibr" rid="B8">2016</xref>). The past decade has witnessed a dramatic growth in the study and implementation of such technology, not only for those with paraplegia due to SCI, but also for individuals with cerebral palsy, stroke, traumatic brain injury and multiple sclerosis (Canela et al., <xref ref-type="bibr" rid="B7">2013</xref>; Murray et al., <xref ref-type="bibr" rid="B39">2015</xref>; Kozlowski et al., <xref ref-type="bibr" rid="B29">2017</xref>; Lerner et al., <xref ref-type="bibr" rid="B35">2017</xref>; Patan&#x000E9; et al., <xref ref-type="bibr" rid="B41">2017</xref>; Androwis et al., <xref ref-type="bibr" rid="B3">2019</xref>; Karunakaran et al., <xref ref-type="bibr" rid="B24">2019</xref>). These exoskeletons are mechanically similar, consisting of a set of linkages that parallel the wearer&#x00027;s thighs, calves, and feet, and augmented with actuators to provide alternatives to muscle torque at the joints. Unlike devices developed for military and industrial tasks, most rehabilitation exoskeletons rely on the subject&#x00027;s use of crutches or canes to provide balance, as the devices lack sensory feedback and balance compensation (Contreras-Vidal et al., <xref ref-type="bibr" rid="B8">2016</xref>).</p>
<p>Where current rehabilitation exoskeletons differ significantly is in their detection of user initiation of gait patterns (Dellon and Matsuoka, <xref ref-type="bibr" rid="B9">2007</xref>; Strickland, <xref ref-type="bibr" rid="B47">2012</xref>)<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref><sup>,</sup><xref ref-type="fn" rid="fn0002"><sup>2</sup></xref>. The Ekso (Ekso Bionics) has motors at the hip and knee, with passive springs at the ankles to provide gait only in the sagittal plane (Strickland, <xref ref-type="bibr" rid="B47">2012</xref>)<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref>. The Ekso has two options to initiate the gait cycle. (1) The first allows a clinician to control gait by means of an external switch pad for training and therapy. (2) For independent control, sensors embedded in the device detect changes in the hip position<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref>. A step can be initiated by the user moving the hip forward and laterally or by changing tilt angle and making ground contact with sensors on the crutches<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref>. Goldfarb et al. developed what has been commercialized as the Indego (Parker Hannifin) (Farris, <xref ref-type="bibr" rid="B14">2012</xref>; Farris et al., <xref ref-type="bibr" rid="B15">2012</xref>; Quintero et al., <xref ref-type="bibr" rid="B44">2012</xref>). It also has powered degrees of freedom (DOF) at the hip and knee, and passive ankle support in the sagittal plane. It uses Hall effect sensors, potentiometers, and accelerometers to detect the user&#x00027;s center of pressure (COP) (Farris, <xref ref-type="bibr" rid="B14">2012</xref>; Quintero et al., <xref ref-type="bibr" rid="B44">2012</xref>). When the user leans forward with both crutches touching the ground, the COP shifts in the direction of movement, and the exoskeleton initiates swing of the most rearward leg (Farris, <xref ref-type="bibr" rid="B14">2012</xref>; Quintero et al., <xref ref-type="bibr" rid="B44">2012</xref>). The Rewalk (Rewalk Robotics) also has two active DOFs with the ankle consisting of a simple orthotic joint with limited motion and spring assisted dorsiflexion (Esquenazi et al., <xref ref-type="bibr" rid="B12">2012</xref>). The control system includes a tilt sensor to determine changes in trunk motion and center of gravity. Shifts in the center of gravity initiate the preprogrammed hip and knee displacement in the appropriate leg (Esquenazi et al., <xref ref-type="bibr" rid="B12">2012</xref>). The HAL (Cyberdyne) employs a combination of EMG gait initiation detection with an accelerometer and gyroscope to sense body posture (Lee and Sankai, <xref ref-type="bibr" rid="B33">2002</xref>, <xref ref-type="bibr" rid="B34">2003</xref>; Hayashi et al., <xref ref-type="bibr" rid="B21">2005</xref>). In contrast to the above exoskeletons that have two active degrees of freedom, the REX (REX Bionics) has 5 motors per leg and is the only available assistive exoskeleton to provide movement in the coronal as well as sagittal planes, and to be self-balancing (i.e., no crutches<xref ref-type="fn" rid="fn0002"><sup>2</sup></xref>). It is controlled by a joystick that signals one of 8 discrete directions of ambulation, and has button selection for sitting and rising<xref ref-type="fn" rid="fn0002"><sup>2</sup></xref>.</p>
<p>Electrophysiological signals have also been employed to initiate exoskeleton gait (Kilicarslan et al., <xref ref-type="bibr" rid="B27">2013</xref>; Kwak et al., <xref ref-type="bibr" rid="B30">2015</xref>; Lebedev and Nicolelis, <xref ref-type="bibr" rid="B32">2017</xref>). Contreras-Vidal et al. have demonstrated the use of EEG triggers to select various REX exoskeleton&#x00027;s discrete control commands (Farris et al., <xref ref-type="bibr" rid="B15">2012</xref>). Similarly, other investigators have focused on detecting gait initiation and termination events using EEG, EOG, evoked potentials, and other bioelectric signals (Nicolelis, <xref ref-type="bibr" rid="B40">2003</xref>; Kilicarslan et al., <xref ref-type="bibr" rid="B27">2013</xref>; Kwak et al., <xref ref-type="bibr" rid="B30">2015</xref>).</p>
<p>Beyond the triggering of preplanned gait patterns, extensive research has been pursued on control methods that allow users to plan and execute novel gait patterns similar to those of individuals with no disability. Gancet et al. have tried to interpret EEG signals from the motor cortex to calculate the kinematics of the gait cycle. A dynamic recurrent neural network was used to train the network to detect the gait patterns in the EEG signal (Gancet et al., <xref ref-type="bibr" rid="B18">2012</xref>). Lebedev and Nicolelis (<xref ref-type="bibr" rid="B32">2017</xref>) has also attempted to use BCI to communicate the user&#x00027;s desired gait cycle kinematics. Unfortunately, both groups have reported numerous challenges such as difficulty in identifying the user intention of each joint, removal of mechanical artifacts caused by relative movement of the EEG cap, and physiological artifacts due to muscle activity in the vicinity of the cap. Even with extensive signal processing, they were not able to completely isolate the relevant signals at all time periods. L&#x000F3;pez-Larraz et al. (<xref ref-type="bibr" rid="B36">2016</xref>) have concluded that the current state of the art of non-invasive BCI knowledge is insufficient for precise decoding of neutrally intended leg kinematics.</p>
<p>Similar to our intention of redirecting controllable neuromuscular activities to define real-time novel gait patterns, Durandau et al. (<xref ref-type="bibr" rid="B11">2019</xref>) have explored the use of residual muscle force activity by using detecting EMG signals. These investigators explain that no other exoskeleton has the ability to amplify weak user muscle strength. Another group (Ferris and Lewis, <xref ref-type="bibr" rid="B16">2009</xref>) have considered the use of proportional EMG signals to activate pneumatic muscles that power a lower extremity exoskeleton. Yet another group has recognized the contralateral synchronization of the arms and legs in unimpaired walking and has used the shoulder angles to define unique sets of hip, knee, and ankle angles that are provide to the exoskeleton (Fang et al., <xref ref-type="bibr" rid="B13">2017</xref>).</p>
<p>Similar to those efforts, the long-term objective of our work is to allow users to intuitively express their desired gait kinematics and dynamics using their arms and hands as sensory and motor alternatives to their legs and feet. Our prior lower extremity research (Karunakaran et al., <xref ref-type="bibr" rid="B23">2014</xref>, <xref ref-type="bibr" rid="B25">2017</xref>) evaluated the feedback conditions required by the hand to produce gait kinematics. That study included 18 subjects controlling virtual feet using hand movements to produce gait trajectories in a virtual environment. Our results indicated that users, provided with haptic through a physical link, and visual feedback (both sensations felt by the hands and visual observation), produced hand and virtual foot trajectories similar to biological gait trajectories (Karunakaran et al., <xref ref-type="bibr" rid="B23">2014</xref>, <xref ref-type="bibr" rid="B25">2017</xref>). We showed that for this to be a viable exoskeleton control method, the hands must be haptically connected either contralaterally or ipsilaterally to the feet. This ensures that the hands and feet move in precise time synchrony, and that hands sense the movement of the feet in order to provide the central nervous system with both position and force feedback from the feet (Karunakaran et al., <xref ref-type="bibr" rid="B23">2014</xref>, <xref ref-type="bibr" rid="B25">2017</xref>). This study as well as our other prior work has shown that neural control of arm and finger walking-like movements generates kinematics very similar to biologically intact neurally determined leg movements and foot placement (DeMarco and Foulds, <xref ref-type="bibr" rid="B45">2002</xref>; Birmingham et al., <xref ref-type="bibr" rid="B6">2003</xref>; Karunakaran et al., <xref ref-type="bibr" rid="B23">2014</xref>, <xref ref-type="bibr" rid="B25">2017</xref>).</p>
<p>The goal of this work was to develop and user test an intuitive control mechanism that independently controls both legs, while producing symmetrical gait kinematics using trajectories generated by hand movements in real time. Balance and co-ordination are being addressed in the later phase of the work, and will be presented in a subsequent paper.</p>
<p>We have chosen to allow the user to control movement of each foot with movements of the ipsilateral hand and arm. This may appear counter intuitive since the arms swing contralaterally with respect to leg movements during normal walking. However, we believe that ipsilateral control presents little if any physiological impediment to a successful user/exoskeleton interface, while its advantages are significant.</p>
<p>Although contralateral arm swing is commonly observed, its function in human gait is not entirely understood. Recent studies (Meesen et al., <xref ref-type="bibr" rid="B37">2006</xref>; Meyns et al., <xref ref-type="bibr" rid="B38">2013</xref>) agree there may be an enhancement in balance and stability, however, they also note that walkers can carry objects and make purposeful arm movements without compromising balance.</p>
<p>Contralateral arm movement is the preferred pattern at walking speeds above 0.8 m/s, yet at speeds below 0.8 m/s, individuals will often adopt an ipsilateral pattern or not swing the arms at all (Ford et al., <xref ref-type="bibr" rid="B17">2007</xref>). In studies of walking with constrained or impaired arm movements (Meesen et al., <xref ref-type="bibr" rid="B37">2006</xref>; Ford et al., <xref ref-type="bibr" rid="B17">2007</xref>) there were only minor reductions in walking speed, that could be voluntarily corrected by participants. Studies of ladder climbing (Armstrong et al., <xref ref-type="bibr" rid="B4">2009</xref>) show a preference for ipsilateral hand and arm coordination. Most importantly, in Meesen et al. (<xref ref-type="bibr" rid="B37">2006</xref>) subjects were asked to make four types of simultaneous arm/leg movements: ipsilateral same direction movements (e.g., right arm and right leg raised up and down, or adducted/abducted), opposite direction ipsilateral movements (e.g., right arm up while right leg down, or right arm adducted and right leg abducted), contralateral same direction movements and contralateral opposite direction movements. When examining the quality to interlimb coordination, the investigators found that the ipsilateral, same direction movements were modestly more accurate in absolute position/angle and phase than both contralateral conditions (Meesen et al., <xref ref-type="bibr" rid="B37">2006</xref>). Thus, we are confident that we could employ either ipsilateral or contralateral arm/hand control method.</p>
<p>For our purposes, ipsilateral control of foot movement has several significant advantages. Admittance control requires the user to apply forces to the foot of the exoskeleton. This is most easily accomplished via a rigid trekking pole on the same side as the foot. A pole on the same side of the body facilitates directing the foot to move vertically and horizontally in the sagittal plane as well as controlling ab/adduction in coronal plane. Also, the poles allow the user&#x00027;s to feel the impact of the ground reaction through the ipsilateral hands. This quality of sensory feedback is unavailable in any proposed exoskeleton control method, and has been shown in our earlier work (Karunakaran et al., <xref ref-type="bibr" rid="B25">2017</xref>) to be of vital importance to controlling the movement of the feet.</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>Materials and Methods</title>
<sec>
<title>Apparatus</title>
<p>A 1/2 scale biped robot representing a lower extremity exoskeleton was built to test our control method. Each leg has 2 links, from hip to knee, and from knee to ankle, and a foot scaled to the anthropometry of the human leg.</p>
<p>Each robot leg has 5 DOF (<xref ref-type="fig" rid="F1">Figure 1</xref>). The hip has 2 DOF for flexion/extension and abduction/adduction, the knee has 1 DOF for flexion/extension, and the ankle has 2 DOF for plantarflexion/dorsiflexion and inversion/eversion. Using previously published data (Hamill et al., <xref ref-type="bibr" rid="B22">2013</xref>), maximum thigh angular velocities and angular accelerations were found to be 28 radians/second (rad/s) and 35 radians/second<sup>2</sup>, respectively. Using this angular acceleration with an estimate of the robot leg moment of inertia relative to the hip, the motor torque required to achieve the maximum acceleration was computed to be 4.2 Newton-meter (N-m). The similarly computed maximum knee motor torque is smaller due to the smaller moment of inertia. We selected the Dynamixel MX-106 smart servomotor (Robotis, USA)<xref ref-type="fn" rid="fn0003"><sup>3</sup></xref> as the actuator for all joints, since both its angular velocity, 42 rad/s, resolution of 0.088 degrees and maximum torque, 8.5 N-m, allow the robot to match the physiological leg segment velocities and accelerations.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p><bold>(A)</bold> Mount to attach the foot to ankle motors. <bold>(B)</bold> Foot of the biped with extrusion to mount (1) Optoforce. <bold>(C)</bold> Front view of 10 DOF biped robot designed based on anthropometric data. L1 is the link length between hip and knee, L2 is link length between knee and foot. <bold>(C)</bold> The coordinate system X, Y, and Z used for robot&#x00027;s movement on treadmill.</p></caption>
<graphic xlink:href="frobt-07-00108-g0001.tif"/>
</fig>
<p>Dynamixel servos employ Maxon motors supported by 32-bit internal microcontrollers providing proportional/integrative/derivative (PID) control at 1,000 Hz. All motors are daisy chained by a 3-wire bus on which they are group addressed from MATLAB software at 1 Mbits/second so that the motors are activated simultaneously. A 3 DOF Optoforce sensor<xref ref-type="fn" rid="fn0004"><sup>4</sup></xref> detects the forces exerted by the user on carbon-fiber trekking poles that are rigidly attached to the sensor.</p>
</sec>
<sec>
<title>User Control Algorithm</title>
<p>The control algorithm consists of an outer admittance loop running at 100 Hz and an inner impedance loop running at 1,000 Hz. The Dynamixel motor&#x00027;s internal PID controller serves as the inner loop. Our admittance loop receives the Cartesian forces applied to the trekking poles by the user&#x00027;s hands and generates desired Cartesian kinematics for the end-effectors of the robot every 10 milliseconds (ms). Admittance control offers a very intuitive control mechanism; where the robot end-effector will be directed to move in the Cartesian directions proportional to the force applied by the user (Van Der Linde et al., <xref ref-type="bibr" rid="B48">2002</xref>). The force can be scaled to accommodate the user&#x00027;s needs. It is safe and easy for human interaction (Van Der Linde et al., <xref ref-type="bibr" rid="B48">2002</xref>; Haidegger et al., <xref ref-type="bibr" rid="B20">2009</xref>).</p>
<p>Custom inverse kinematics algorithms transform the Cartesian positions generated by the admittance loop into joint angles used to command the motors. Algorithms are written in MATLAB, with time-dependent functions coded in C to maintain a 100 Hz update rate.</p>
</sec>
<sec>
<title>Sagittal Plane Control</title>
<p>Robot control in the sagittal plane is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. User forces applied in the X and Y-direction are read at 1,000 Hz from two axes of the 3 DOF Optoforce force sensor, and averaged to provide 100 samples/second. For every cycle of the admittance loop these forces, virtual mass, virtual damping, as well as the foot&#x00027;s Cartesian position and velocity are passed to custom-written ordinary differential equation (ODE) that is implemented using the C source code Variable ordinary differential equation (CVode) solver and compiled for use in MATLAB (Van Riel, <xref ref-type="bibr" rid="B49">2012</xref>). The ODE is shown in <xref ref-type="fig" rid="F2">Figure 2</xref> as a double integration that provides the desired positions and velocities in X and Y to be achieved in the next 10 ms. The value assigned to the mass allows scaling of the forces to meet the capabilities of the user&#x00027;s hands. The ODE also incorporates the admittance loop&#x00027;s damping that maintains stability of the system. No admittance stiffness is included in this version of the software, but can be added if required in future situations. The ODE function is coded to solve the second-order differential equation shown below:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msup><mml:mtext class="textrm" mathvariant="normal">X</mml:mtext><mml:mrow><mml:mo>&#x02033;</mml:mo></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext class="textrm" mathvariant="normal">t</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mtext class="textrm" mathvariant="normal">F</mml:mtext><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext class="textrm" mathvariant="normal">X</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">M</mml:mtext></mml:mrow></mml:mfrac><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:mtext class="textrm" mathvariant="normal">B</mml:mtext><mml:mo>*</mml:mo><mml:msup><mml:mtext class="textrm" mathvariant="normal">X</mml:mtext><mml:mrow><mml:mo>&#x02032;</mml:mo></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext class="textrm" mathvariant="normal">t</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">M</mml:mtext></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where, F(X) = force (Newton), M = virtual mass (kilogram), B = virtual damping (Newton-second/meter), X&#x02032; (t) =velocity (meter/second), X&#x02033; (t) = acceleration (meter/second<sup>2</sup>).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>The Admittance Control algorithm for control of biped gait in the sagittal plane. &#x003B8;1 is hip angle, &#x003B8;2 is knee angle, and &#x003B8;3 is the ankle angle computed based on Equations (1-3), respectively. Velocities in X and Y are represented by vx and vy, respectively. Acceleration in X and Y are represented by ax, and ay, respectively. Position in X and Y are represented by px and py, respectively.</p></caption>
<graphic xlink:href="frobt-07-00108-g0002.tif"/>
</fig>
</sec>
<sec>
<title>Inverse Kinematics</title>
<p>Since the robot has revolute joints, inverse kinematics converts the X-Y position of the foot to angles of the hip and knee in the sagittal plane. For simplicity, the angle of the ankle is computed to keep the foot parallel to the floor. The angles are calculated using the law of cosines (Equations 2&#x02013;4). The joint angles are in turn converted to hip, knee and ankle motor units in the sagittal plane. These values are in turn fed to the corresponding motors to generate the required torque to perform the movement intended by the user.</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B8;</mml:mi><mml:mn>2</mml:mn><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>2</mml:mn><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:msup><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mtext>&#x000A0;</mml:mtext><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>L</mml:mi><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:mi>L</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>X</mml:mi><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:mi>Y</mml:mi><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>X</mml:mi><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:mi>Y</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>L</mml:mi><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:mi>L</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E3"><label>(3)</label><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B8;</mml:mi><mml:mn>1</mml:mn><mml:mo>=</mml:mo><mml:mtext>tan</mml:mtext><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>L</mml:mi><mml:mn>2</mml:mn><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:mi>L</mml:mi><mml:mn>2</mml:mn><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:mtext>tan</mml:mtext><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>Y</mml:mi><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>X</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E4"><label>(4)</label><mml:math id="M4"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B8;</mml:mi><mml:mn>3</mml:mn><mml:mo>=</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mn>2</mml:mn><mml:mo>-</mml:mo><mml:mtext>tan</mml:mtext><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>Y</mml:mi><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>X</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mtext>tan</mml:mtext><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>L</mml:mi><mml:mn>2</mml:mn><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:mi>L</mml:mi><mml:mn>2</mml:mn><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where, X1, Y1 is the desired end-effector position, L1 is the link length between hip and knee, L2 is link length between knee and ankle, &#x003B8;1 is hip angle, &#x003B8;2 is knee angle, and &#x003B8;3 is the ankle angle.</p>
<p>The Dynamixel motors have sufficiently fast mechanical and electrical response times (&#x0007E;4 ms) to relocate the robot end-effector to the desired location within the 10 ms loop period<xref ref-type="fn" rid="fn0003"><sup>3</sup></xref>. The new position is read at the beginning of the next admittance cycle and serves as the initial conditions for the ODE. As the user&#x00027;s hand is rigidly connected to the robot foot, the user receives real-time sensation of foot movement. The user can modulate his/her forces to alter the speed of foot movement independently in X and Y, and also in response to external forces that may impede the foot.</p>
</sec>
<sec>
<title>Singularity Check</title>
<p>The algorithm verifies that the predicted X-Y position is within reach of the robot. If this check fails, the next position is set to a location on the boundary of the robot&#x00027;s range on a line between the former and predicted position. This maintains stability and smoothness at the singularity. The robot reaches singularity when the estimated position is outside the range of motion of the robot.</p>
</sec>
<sec>
<title>Coronal Plane Control</title>
<p>Force applied by the user to a trekking pole in the z-direction is similarly converted to rotation of the hip in the coronal plane as shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. The user input is treated as an applied torque, so the ODE implements the following rotational equation of motion:</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M5"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x02033;</mml:mo></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext class="textrm" mathvariant="normal">t</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mtext class="textrm" mathvariant="normal">T</mml:mtext><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>I</mml:mi></mml:mrow></mml:mfrac><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:mtext class="textrm" mathvariant="normal">B*</mml:mtext><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext class="textrm" mathvariant="normal">t</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">I</mml:mtext></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E6"><label>(6)</label><mml:math id="M6"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B8;</mml:mi><mml:mn>4</mml:mn><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mo>&#x0222B;</mml:mo></mml:mstyle><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext class="textrm" mathvariant="normal">t</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E7"><label>(7)</label><mml:math id="M7"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B8;</mml:mi><mml:mn>5</mml:mn><mml:mo>=</mml:mo><mml:mn>180</mml:mn><mml:mo>-</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mn>4</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where, T(&#x003B8;) = user applied torque (Newton-meter), I = virtual moment of inertia (kilogram-meter<sup>2</sup>), B = rotational damping (Newton-second/meter), &#x003B8;&#x02032; (t) = angular velocity (meter/second), &#x003B8;&#x02033; (t) = angular acceleration (meter/second<sup>2</sup>). &#x003B8;4 is hip angle, and &#x003B8;5 is the ankle angle in the coronal plane.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>The Admittance Control algorithm for trekking pole control of coronal plane movement of the robot leg. &#x003B8;4 is hip angle, and &#x003B8;5 is the ankle angle and is computed using Equations (6) and (7). Angular velocity and angular acceleration is denoted by vz and az, respectively.</p></caption>
<graphic xlink:href="frobt-07-00108-g0003.tif"/>
</fig>
<p>The ab/adduction angle is executed by the second hip motor, with the new angle serving as the initial condition for the next admittance cycle, and the movement of the robot leg provides a haptic sensation to the user&#x00027;s hand. Similar to the sagittal plane, the inversion/eversion angle of the foot is computed to keep the foot parallel to the floor.</p>
</sec>
<sec>
<title>Biped Control Strategy for Users</title>
<p>Our biped control strategy allows the user to execute the swing and the stance phases of each robot leg independently by applying ipsilateral hand forces at the top of the trekking poles, as shown in <xref ref-type="fig" rid="F4">Figure 4</xref>. The user forces are converted to desired foot positions by the admittance control software, with servomotor response occurring within 10 ms.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>User controlling the biped by holding the trekking poles on the ipsilateral side using hands. Constant force springs connect the robot to the overhead frame to maintain balance.</p></caption>
<graphic xlink:href="frobt-07-00108-g0004.tif"/>
</fig>
<p>Forward progression in gait requires activation of both legs. The swing leg proceeds from ground contact (toe-off) to swing and ultimately to a subsequent heel contact in a new forward position on the ground. At the same time, the stance foot remains in the same position on the ground, with rotation of ankle and hip of the stance leg allowing forward progression of the torso (i.e., shifting the robot&#x00027;s center of gravity forward).</p>
<p>The swing foot movement is controlled by the user applying a time-varying upward hand force to define the Y-axis trajectory and a forward force to signify the desired X-axis foot trajectory. Simultaneously, the stance side hand applies a rearward force only in the negative X direction. As this foot remains in contact with the ground (while the swing foot is in the air), the admittance software controls the ankle, knee, and hip motors of the stance leg for forward progression.</p>
<p>Haptic feedback is provided through the physical link (trekking pole) between the hands and the feet. During the swing, as the user force is converted to foot movement, the trekking pole follows the foot, allowing the hand to move synchronously with the trekking poles to walk the biped. Near the end of swing, the user applies downward force to the trekking pole to bring the foot to heel strike. At heel strike, the ground precludes further downward movement, and the user&#x00027;s hand feels the ground reaction.</p>
</sec>
<sec>
<title>Evaluation of the Accuracy of the Control</title>
<p>A slow gait-like movement was performed by one healthy female participant (who was a member of the research team) for a period of 60 s in the air to evaluate the accuracy of the admittance control and inverse kinematics algorithms, and the time delay. This work was approved by the New Jersey Institute of Technology&#x00027;s (NJIT) institutional review board (IRB).</p>
<p>The accuracy of the control algorithm for the sagittal and coronal plane was assessed by comparing the desired and the actual Cartesian position of the feet in X, Y, and Z direction. The &#x0201C;desired position&#x0201D; is the position of the foot computed by the algorithm based on user input, and the &#x0201C;actual position&#x0201D; is the position reached by the foot of the biped. The accuracy of the Dymanixel&#x00027;s internal impedance control was evaluated by comparing the desired and actual joint angle. The accuracy of the hip and knee angles in the sagittal plane was evaluated while the participant performed the gait using the hand. The accuracy of the joint angle of the hip in the coronal plane was evaluated while the participant performed adduction/abduction. The &#x0201C;desired joint angles&#x0201D; are the angles computed by the algorithm, and the &#x0201C;actual angles&#x0201D; are the angles reached by each joint motor. The motor angles at the end of each iteration from each motor were converted to joint angles of the hip and knee. These joint angles were compared to the desired joint angles at every given time point. The accuracy in sagittal plane was evaluated using a forward kinematics algorithm developed to obtain the X and Y positions of the foot. This position was in turn compared with the desired X and Y positions (X and Y position computed using the admittance control) to evaluate the accuracy of the inverse kinematics algorithm.</p>
<p>A Pearson&#x00027;s r correlation was performed to quantify the similarity between the actual and desired joint angles (hip flexion/extension, knee flexion/extension, and hip adduction/abduction) and Cartesian positions (X, Y, and Z positions).</p>
<p>The average time delay was computed as the iteration time (time required to move all the joints by the motor from the time the forces were applied by the user).</p>
</sec>
<sec>
<title>Evaluation of Gait Using the Control Mechanism</title>
<p>After validating the accuracy of the experimental apparatus, data were collected for an extended study evaluating the control mechanism using seven na&#x000EF;ve participants who controlled the gait of the biped on a treadmill using control mechanism.</p>
<sec>
<title>Participants</title>
<p>The study included seven na&#x000EF;ve participants (2 male and 5 female) for biped walking and one reference participant (control) from whom we collected data on human gait. All participants were between ages 18 and 35, with fully functional upper and lower extremities. Exclusion criteria included any disability to the upper or lower extremities or inability to perform normal gait. The study was approved by the NJIT IRB and the experiment was performed with the participants&#x00027; consent.</p>
</sec>
<sec>
<title>Experimental Setup</title>
<p>The experimental setup included a Pro-form J6 treadmill around which a custom frame was built using 80&#x02013;20 aluminum to support the biped. The frame allowed the users to have complete view of the treadmill and the robot. An Optitrack Trio motion capture system was used to record the biped gait as well as the reference participant gait. As the treadmill was designed for human use, its lowest speed would not accommodate the small-sized robot, thus its speed was reduced when used for biped walking by adding a power resistor in series with the motor. Optitrack markers were placed on the hip, knee, and ankle of both the legs to track the biped gait and the reference participant gait.</p>
</sec>
<sec>
<title>Biped Walking</title>
<p>The biped was placed on the treadmill and each participant was seated in a comfortable chair behind the treadmill. The participants were instructed to control the gait of the biped during each trial by applying small forces to the pole extending from the sensor on each leg in the direction of the intended movement. The study consisted of eight trials. Each trial lasted 1 min, followed by a 30 s rest. The speed of the treadmill during each trial was varied as shown in <xref ref-type="table" rid="T1">Table 1</xref>. The speed variations for the biped were 0.1 (low), 0.2 (medium), and 0.3 (high) mph.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Speed of treadmill for each trial.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Trial</bold></th>
<th valign="top" align="center"><bold>1</bold></th>
<th valign="top" align="center"><bold>2</bold></th>
<th valign="top" align="center"><bold>3</bold></th>
<th valign="top" align="center"><bold>4</bold></th>
<th valign="top" align="center"><bold>5</bold></th>
<th valign="top" align="center"><bold>6</bold></th>
<th valign="top" align="center"><bold>7</bold></th>
<th valign="top" align="center"><bold>8</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Speed</td>
<td valign="top" align="center">Medium</td>
<td valign="top" align="center">Medium</td>
<td valign="top" align="center">High</td>
<td valign="top" align="center">Low</td>
<td valign="top" align="center">Medium</td>
<td valign="top" align="center">High</td>
<td valign="top" align="center">Low</td>
<td valign="top" align="center">Medium</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The participants performed familiarization sessions before the start of the actual session and those data were not included in the analysis. The first familiarization session was performed without the treadmill for 1 min, where the participants controlled the leg of robot in the air to get accustomed to kinematics of the leg. The second familiarization session included eight trials, where each trail lasted 1 min and was performed on the treadmill with the lowest speed. The third familiarization session included controlling the biped at different speeds for a minute each, as shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<p>The participant from whom the reference gait data were collected for one trial used the same treadmill at a self-selected speed for comfortable walking.</p>
</sec>
<sec>
<title>Data Analysis</title>
<p>Horizontal and vertical trajectories collected at 120 Hz of the ankle, hip and knee were filtered using a 4th order, zero-lag Butterworth low-pass filter. The filtered data were used for further analysis. The data were further divided into gait cycles.</p>
<sec>
<title>Spatial and temporal symmetry</title>
<p>The foot trajectories in the sagittal plane by 7 participants and 1 reference participant were evaluated for the effectiveness of the control using hand trajectories using temporal and spatial symmetry outcomes.</p>
<p>The swing and stance time of each foot during each gait cycle was computed, and equation 8 was used to compute the temporal symmetry (Patterson et al., <xref ref-type="bibr" rid="B42">2008</xref>). Similarly step length and step height of each foot were computed for each gait cycle, and equation 9 and 10 were used to compute spatial and step height symmetry (Patterson et al., <xref ref-type="bibr" rid="B42">2008</xref>). The average temporal, spatial, and step height symmetries were computed for all participants in each trial.</p>
<p>Statistical analysis was performed on all trials across the seven participants. Shapiro-Wilk test (<italic>p</italic> &#x0003E; 0.05) of normality showed that data were normal for spatial and temporal symmetry. Repeated measures analysis of variance (ANOVA) was performed on the spatial and temporal symmetry, respectively to determine the effect of change in speed on the spatial and temporal symmetry, respectively. Further, a Greenhouse- Geisser test was performed, since the data showed significance with Mauchly&#x00027;s test for sphericity. Shapiro-Wilk test (<italic>p</italic> &#x0003C; 0.05) of normality showed that data were not normal for step height symmetry. Hence, Friedman Test was used to determine the difference between different trials.</p>
<disp-formula id="E8"><label>(8a)</label><mml:math id="M8"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>T</mml:mi><mml:mi>e</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi> <mml:mi>s</mml:mi><mml:mi>y</mml:mi><mml:mi>m</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mo>)</mml:mo><mml:mo>/</mml:mo><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mo>)</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mi>O</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>y</mml:mi><mml:mi>m</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E10"><label>(8b)</label><mml:math id="M10"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mo>=</mml:mo><mml:mtext class="textrm" mathvariant="normal"></mml:mtext><mml:mstyle class="math"><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>R</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>y</mml:mi><mml:mi>m</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mi>e</mml:mi><mml:mi>f</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>y</mml:mi><mml:mi>m</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mtext class="textrm" mathvariant="normal"></mml:mtext></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E11"><label>(9)</label><mml:math id="M11"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>y</mml:mi><mml:mi>m</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>R</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>p</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mi>e</mml:mi><mml:mi>f</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>p</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E12"><label>(10)</label><mml:math id="M12"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>S</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>p</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>H</mml:mi><mml:mi>e</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>y</mml:mi><mml:mi>m</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>R</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>p</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>h</mml:mi><mml:mi>e</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mi>e</mml:mi><mml:mi>f</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>p</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>h</mml:mi><mml:mi>e</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
</sec>
<sec>
<title>Duty cycle</title>
<p>The percentage of stance and swing phase for each gait cycle was calculated using equations 11, 12, and 13. The average duty cycle of all gait cycles was computed for all participants in each trial and for the reference participant.</p>
<p>Statistical Analysis was performed on all trials across the seven participants. Shapiro-Wilk test (<italic>p</italic> &#x0003E; 0.05) of normality showed that data were normal. Mixed design analysis of variance (ANOVA) was performed on the stance and swing duty cycle, respectively to determine the effect of change in speed as well difference between left and right leg on the stance and swing duty cycle, respectively. The Greenhouse-Geisser test was used, since the data showed significance with Mauchly&#x00027;s test for sphericity to determine the effect of change in speed.</p>
<disp-formula id="E13"><label>(11)</label><mml:math id="M13"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>D</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi><mml:mi>y</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>c</mml:mi><mml:mi>y</mml:mi><mml:mi>c</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mo>=</mml:mo><mml:mi>S</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>P</mml:mi><mml:mi>h</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mi>S</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>P</mml:mi><mml:mi>h</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E14"><label>(12)</label><mml:math id="M14"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>S</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>P</mml:mi><mml:mi>h</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mi>%</mml:mi><mml:mo>=</mml:mo><mml:mn>100</mml:mn><mml:mtext>&#x000A0;</mml:mtext><mml:mo>*</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mi>S</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>p</mml:mi><mml:mi>h</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mo>/</mml:mo><mml:mi>D</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi><mml:mi>y</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>C</mml:mi><mml:mi>y</mml:mi><mml:mi>c</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E15"><label>(13)</label><mml:math id="M15"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>S</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>P</mml:mi><mml:mi>h</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mi>%</mml:mi><mml:mo>=</mml:mo><mml:mn>100</mml:mn><mml:mtext>&#x000A0;</mml:mtext><mml:mo>*</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mi>S</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>g</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>p</mml:mi><mml:mi>h</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mo>/</mml:mo><mml:mi>D</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi><mml:mi>y</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>C</mml:mi><mml:mi>y</mml:mi><mml:mi>c</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
</sec>
<sec>
<title>Joint angles</title>
<p>The joint angles of hip and knee of both the legs of the biped robot walking for all participants were computed from the filtered Cartesian position from the Optitrack data using inverse kinematics. The joint angles were also computed for the gait of the single reference participant.</p>
<p>We statistically compared the similarity of the hip and knee angles of all seven participants with those of the reference participant by computing their Pearson&#x00027;s r correlation after the data had been time-warped to allow a direct comparison (Kianimajd et al., <xref ref-type="bibr" rid="B26">2017</xref>). We followed the method described in La Scaleia et al. (<xref ref-type="bibr" rid="B31">2014</xref>).</p>
</sec>
<sec>
<title>User feedback</title>
<p>Subjective user feedback was obtained from all 7 participants after completing the study. The following questions elicited user experience of the control technique and the responses to each were averaged across all participants.</p>
<list list-type="order">
<list-item><p>Was it easy, moderate or difficult to use hand movements as control? 0 being easy and 10 being difficult.</p></list-item>
<list-item><p>How tired were your hands after each session? 0 being not tired at all and 10 being very tired.</p></list-item>
<list-item><p>Was using hands to the control the leg intuitive? 0 being least intuitive and 10 being very intuitive.</p></list-item>
<list-item><p>How much force was required to move the leg in the direction intended? 0 being least force and 10 being great force.</p></list-item>
<list-item><p>Did you feel the haptic feedback every time foot made contact with the floor? Yes or No.</p></list-item>
</list>
</sec>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec>
<title>Accuracy of the Foot Positions and Joint Angles</title>
<p><xref ref-type="fig" rid="F5">Figures 5A&#x02013;C</xref> show the actual and desired Cartesian positions of X, Y (sagittal plane), and Z (coronal plane), respectively, for a single evaluator of the control mechanism of the robot. The results exhibit minimal positional lag. This is quantitatively shown in <xref ref-type="table" rid="T2">Table 2</xref>, where the mean positional lag (difference between actual and desired position) is &#x0003C;1 cm in the X, Y, and Z-directions, with Pearson&#x00027;s r showing high correlation between the desired and the actual Cartesian positions in X (right <italic>r</italic> = 0.9954, <italic>p</italic> &#x0003C; 0.05, left <italic>r</italic> = 0.9972, <italic>p</italic> &#x0003C; 0.05), Y (left <italic>r</italic> = 0.9995, <italic>p</italic> &#x0003C; 0.05, right <italic>r</italic> = 0.9986, <italic>p</italic> &#x0003C; 0.05), and Z (left <italic>r</italic> = 0.9979, <italic>p</italic> &#x0003C; 0.05). This shows the accuracy of the control algorithm.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Cartesian values of the desired and actual foot position of <bold>(A)</bold> X, <bold>(B)</bold> Y, and <bold>(C)</bold> Z positions vs. time for the left foot. Desired angles generated by the inverse kinematics and actual angles achieved by the motors of <bold>(D)</bold> hip in sagittal, <bold>(E)</bold> knee in sagittal, and <bold>(F)</bold> hip in coronal plane vs. time for left leg. The flat regions in the knee plot indicate the stance phase of gait. The red lines denote the actual position/angle reached by the robot and the blue lines denote the desired position/angle computed by the algorithm.</p></caption>
<graphic xlink:href="frobt-07-00108-g0005.tif"/>
</fig>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>The table shows Mean &#x000B1; Std.</p></caption>
<table frame="hsides" rules="groups">
<tbody><tr>
<td valign="top" align="left" colspan="2" style="background-color:#bbbdc0"><bold>ERROR IN ANGLE</bold></td>
</tr>
<tr>
<td valign="top" align="left">Right&#x02014;hip flexion/extension (degrees)</td>
<td valign="top" align="center">1.02 &#x000B1; 0.006</td>
</tr>
<tr>
<td valign="top" align="left">Left&#x02014;hip flexion/extension (degrees)</td>
<td valign="top" align="center">0.81 &#x000B1; 0.006</td>
</tr>
<tr>
<td valign="top" align="left">Right&#x02014;knee flexion/extension (degrees)</td>
<td valign="top" align="center">0.91 &#x000B1; 0.011</td>
</tr>
<tr>
<td valign="top" align="left">Left&#x02014;knee flexion/extension (degrees)</td>
<td valign="top" align="center">1.021 &#x000B1; 0.011</td>
</tr>
<tr>
<td valign="top" align="left">Left&#x02014;hip abduction/adduction (degrees)</td>
<td valign="top" align="center">0.55 &#x000B1; 0.003</td>
</tr>
<tr>
<td valign="top" align="left" colspan="2" style="background-color:#bbbdc0"><bold>ERROR IN POSITION</bold></td>
</tr>
<tr>
<td valign="top" align="left">Right&#x02014;X position (cm)</td>
<td valign="top" align="center">0.68 &#x000B1; 0.005</td>
</tr>
<tr>
<td valign="top" align="left">Left&#x02014;X position (cm)</td>
<td valign="top" align="center">0.54 &#x000B1; 0.004</td>
</tr>
<tr>
<td valign="top" align="left">Right&#x02014;Y position (cm)</td>
<td valign="top" align="center">0.10 &#x000B1; 0.0008</td>
</tr>
<tr>
<td valign="top" align="left">Left&#x02014;Y position (cm)</td>
<td valign="top" align="center">0.13 &#x000B1; 0.001</td>
</tr>
<tr>
<td valign="top" align="left">Left&#x02014;Z position (cm)</td>
<td valign="top" align="center">0.23 &#x000B1; 0.002</td>
</tr>
<tr>
<td valign="top" align="left">Time delay (s)</td>
<td valign="top" align="center">0.008 &#x000B1; 0.0002</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Error of the (a) difference between actual and desired position (error in position) in the left and right foot in X and Y direction and in the left foot in Z direction, (b) difference between actual and desired hip and knee angle in sagittal plane for right and left leg and hip angle in coronal plane for left leg (error in angle), and (c) time delay (time for each iteration of the loop) to reach the actual position based on user input force</italic>.</p>
</table-wrap-foot>
</table-wrap>
<p><xref ref-type="fig" rid="F5">Figures 5D&#x02013;F</xref> show the accuracy of the impedance control of the hip, knee (in the sagittal plane), and the hip (in the coronal plane), respectively. The actual joint angles follow the desired joint angles with minimal angular lag. This is quantitatively shown in <xref ref-type="table" rid="T2">Table 2</xref>, where the angular lag (difference between actual and desired angles) in the hip and knee (sagittal plane), and the hip (coronal plane) is &#x0003C;1&#x000B0;, and Pearson&#x00027;s <italic>r</italic> showed that correlation was high between desired and actual joint angles in hip flexion/extension (right = 0.9985, <italic>p</italic> &#x0003C; 0.05, left = 0.9984, <italic>p</italic> &#x0003C; 0.05), knee flexion/extension (left = 0.9983, <italic>p</italic> &#x0003C; 0.05, right = 0.9978, <italic>p</italic> &#x0003C; 0.05), and hip abduction/adduction (left = 0.9982, <italic>p</italic> &#x0003C; 0.05). This shows the accuracy of the impedance control.</p>
<p>Taken together, the results show that the biped&#x00027;s foot position reaches the desired position of the user with minimal lag. The time delay or control loop time (time required to move all the joints by the motor from the time the forces were applied by the user) was &#x0003C;10 ms (<xref ref-type="table" rid="T2">Table 2</xref>). Studies have shown that a control loop frequency of 100 Hz is sufficient for human operators to feel smooth, nearly passive, movements of a robot (Van Der Linde et al., <xref ref-type="bibr" rid="B48">2002</xref>). The maximum error in the Cartesian position is &#x0003C;1 cm, and the error in joint angles is &#x0003C;1&#x000B0; (<xref ref-type="table" rid="T2">Table 2</xref>). These results validate that our experimental control method and robot are appropriate for the multi-participant experiments.</p>
</sec>
<sec>
<title>Spatial and Temporal Symmetry</title>
<p><xref ref-type="fig" rid="F6">Figures 6A&#x02013;C</xref> show the overall temporal, spatial and step height symmetry, respectively, for 7 participants and 1 reference participant. All symmetries are close to 1, irrespective of speed variations, indicating that the users were controlling the biped with a bilaterally symmetrical gait. Repeated measures ANOVA shows no significant difference between trials, indicating that speed did not affect the spatial (<italic>p</italic> &#x0003E; 0.05, <italic>F</italic> = 0.796) and temporal (<italic>p</italic> &#x0003E; 0.05, <italic>F</italic> = 0.424) symmetry. In addition, Cohen&#x00027;s d effect size shows a low effect for both spatial (effect size = 0.21) and temporal (effect size = 0.27) symmetry, again signifying that difference between trials is very low. Friedman&#x00027;s test shows no significant difference between trials, indicating that speed did not affect the step height symmetry (<italic>p</italic> &#x0003E; 0.05, chi-square = 8.095). Kendall&#x00027;s W effect size shows a low effect for step height symmetry (effect size = 0.165), again signifying that difference between trials is very low.</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Mean &#x000B1; std. error of <bold>(A)</bold> spatial symmetry, <bold>(B)</bold> temporal symmetry, and <bold>(C)</bold> step height symmetry, for trials 1 through 8 for biped walking by 7 participants and for 1 trial by the reference participant. <bold>(D)</bold> X and Y positions of one biped foot for one na&#x000EF;ve subject walking the biped on the treadmill at medium speed.</p></caption>
<graphic xlink:href="frobt-07-00108-g0006.tif"/>
</fig>
</sec>
<sec>
<title>Duty Cycle</title>
<p><xref ref-type="fig" rid="F7">Figure 7</xref> shows the average percentage of swing and stance phase in the gait cycle of all 7 participants for trials 1 through 8. In human gait, typical swing phase is &#x0007E;40% and stance phase is 60% of the gait cycle (Winter, <xref ref-type="bibr" rid="B51">2009</xref>). Our robot gait cycle across all trials was slightly over 40% swing and slightly below 60% stance. Mixed design ANOVA shows no significant difference between either trials or legs. This indicates that speed did not affect swing (<italic>p</italic> &#x0003E; 0.05, between trials <italic>F</italic> = 1.177) or stance (<italic>p</italic> &#x0003E; 0.05, between trials <italic>F</italic> = 1.177), and that swing (<italic>p</italic> &#x0003E; 0.05, between legs <italic>F</italic> = 0.022) and stance (<italic>p</italic> &#x0003E; 0.05, between legs <italic>F</italic> = 0.022) were similar in both the legs.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Mean &#x000B1; std. error of duty cycle in <bold>(A)</bold> right leg and <bold>(B)</bold> left leg of 7 participants for trials 1 through 8, and of 1 reference participant (C) during 1 trial. Stance phase is red and swing phase is blue.</p></caption>
<graphic xlink:href="frobt-07-00108-g0007.tif"/>
</fig>
</sec>
<sec>
<title>Comparison of Biped and Human Gait Trajectories in the Sagittal Plane</title>
<p><xref ref-type="fig" rid="F8">Figure 8</xref> shows a visual comparison of the hip and knee angles of a reference participant and the corresponding biped angles. The visual appearance of biped and reference participant hip and knee angles in <xref ref-type="fig" rid="F8">Figure 8</xref> is quite similar, with the human knee slightly flexing during stance. We attribute this to the knee accommodation of foot angle changes at toe-off by the human walker, while the biped maintained the foot parallel to the ground. We observe significant high positive correlations between the joint angles of the biped (for 7 participants) and the reference participant: For the knee, Pearson&#x00027;s mean <italic>r</italic> = 0.7770, std &#x000B1; 0.0953, with <italic>p</italic> &#x0003C; 0.005. For the hip, Pearson&#x00027;s mean <italic>r</italic> = 0.9968, std &#x000B1; 0.0926, with <italic>p</italic> &#x0003C; 0.005. This confirms that biped walking produces human-like knee and hip joint trajectories.</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>Right knee angle of multiple strides of the <bold>(A)</bold> user-controlled biped, and <bold>(B)</bold> reference participant, walking on the treadmill. Right hip angle of multiple strides of the <bold>(C)</bold> user-controlled biped, and <bold>(D)</bold> reference participant, walking on the treadmill.</p></caption>
<graphic xlink:href="frobt-07-00108-g0008.tif"/>
</fig>
</sec>
<sec>
<title>User Feedback</title>
<p>Participant responses show that hand control was quite easy and intuitive (<xref ref-type="table" rid="T3">Table 3</xref>). The participants required little force and did not get tired from using the control method (<xref ref-type="table" rid="T3">Table 3</xref>). Six out of the seven participants responded that they felt the ground impact with their hands.</p>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>The table shows the mean &#x000B1; std. error for user feedback questionnaire.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Question</bold></th>
<th valign="top" align="left"><bold>Ease of using hand to control biped</bold></th>
<th valign="top" align="left"><bold>Tiredness</bold></th>
<th valign="top" align="left"><bold>Intuitiveness of hand control</bold></th>
<th valign="top" align="left"><bold>Force required to move hands</bold></th>
<th valign="top" align="left"><bold>Haptic feedback every time foot made contact to ground</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Options</td>
<td valign="top" align="left">0&#x02014;easy to 10&#x02014;difficult</td>
<td valign="top" align="left">0&#x02014;not tired to 10&#x02014;very tired</td>
<td valign="top" align="left">0&#x02014;least intuitive to 10&#x02014;very intuitive</td>
<td valign="top" align="left">0&#x02014;least force to 10&#x02014;large force</td>
<td valign="top" align="left">Yes or no</td>
</tr>
<tr>
<td valign="top" align="left">Participant responses</td>
<td valign="top" align="left">3.29 &#x000B1; 0.52</td>
<td valign="top" align="left">1.14 &#x000B1; 0.59</td>
<td valign="top" align="left">7.43 &#x000B1; 1.19</td>
<td valign="top" align="left">1.86 &#x000B1; 0.59</td>
<td valign="top" align="left">6 participants answered yes and 1 participant answered no</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>Lower extremity exoskeletons show potential for restoring ambulation in individuals with paraplegia due to spinal cord injury. Currently, many lower extremity exoskeletons produce a preprogrammed gait that can be initiated by the user, but not completely controlled by the user. Thus, users have limited control over their gait (i.e., step length, width, and speed), and little feedback of foot placement and ground contact. In this paper, we proposed a novel control mechanism for lower extremity exoskeletons to address these shortcomings; allowing the users to control their gait in real-time, as well as receiving haptic feedback. The accuracy and efficacy of the control mechanism was evaluated using the following outcome measures: accuracy of the foot position and joint angles; time delay between user-desired and actual robot kinematics; temporal, spatial, and step height symmetry; duty cycle of stance and swing phases; and similarity of biped hip and knee angles to those of human gait.</p>
<p>Our results demonstrate the effectiveness of the control mechanism, which allows the user to interactively control both legs of a &#x000BD; scale robot to produce a gait trajectory closely resembling human gait. As required for smooth real-time human-robot interaction (Van Der Linde et al., <xref ref-type="bibr" rid="B48">2002</xref>), the time delay of the control mechanism is below 10 ms; ensuring smooth and stable movement. The difference between user-desired and actual robot Cartesian position and joint angles were &#x0003C;1 cm and 1&#x000B0;, respectively; showing that real-time accurate trajectories can be obtained over time by our system.</p>
<p>Admittance control provides the user with the following advantages: one-to-one correspondence between hand and foot movement; force amplification; intuitive control (as the user applies desired force in the intended direction of movement); and accurate time-varying trajectory control (the leg moves in the direction of the force) (Glosser and Newman, <xref ref-type="bibr" rid="B19">2002</xref>; Van Der Linde et al., <xref ref-type="bibr" rid="B48">2002</xref>). Research has shown that admittance control is an effective strategy for human-robot interaction, and that it requires a minimum of 100 Hz for acceptable human interaction (Glosser and Newman, <xref ref-type="bibr" rid="B19">2002</xref>; Van Der Linde et al., <xref ref-type="bibr" rid="B48">2002</xref>); which we have achieved.</p>
<p>Human gait is a rhythmic movement that is symmetrical but out-of-phase between the two legs. Hence, an effective exoskeleton control strategy should be able to reproduce the same pattern of movement over a period of time, as well as coordinate inter-limb movement (Vaughan et al., <xref ref-type="bibr" rid="B50">1999</xref>; Pearson, <xref ref-type="bibr" rid="B43">2003</xref>). Inter-limb coordination or symmetry is the ability to maintain temporal and spatial symmetry, which results in healthy gait (Patterson et al., <xref ref-type="bibr" rid="B42">2008</xref>). Temporal and spatial symmetry close to &#x0201C;1,&#x0201D; signifies that both limbs are performing a symmetrical movement (Patterson et al., <xref ref-type="bibr" rid="B42">2008</xref>). Our analysis shows that the participants controlling the biped were able to maintain inter-limb temporal, spatial, and step height symmetry close to 1 across all trials; irrespective of speed variations, and similar to that observed in human gait. Studies have shown that there is a correlation between large deviations in temporal &#x00026; spatial symmetry away from 1 and falls &#x00026; reduced walking speed (Balasubramanian et al., <xref ref-type="bibr" rid="B5">2007</xref>; Patterson et al., <xref ref-type="bibr" rid="B42">2008</xref>). Thus, we can conclude that our biped control mechanism will promote stable and safe human-exoskeleton walking. Our robot gait cycle has an average swing phase of just over 40% and a stance phase of just under 60%, as shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, across all trials. The deviation from the normal human gait cycle is small and likely due to the biped performing a flat-footed gait, as the foot is constrained to remain parallel to the ground, resulting in a shorter stance period. In addition, the percentage of time spent in stance and swing phase was the same between both legs, further indicating consistent temporal characteristics.</p>
<p>The results show that the joint angles produced by the biped (<xref ref-type="fig" rid="F5">Figures 5D&#x02013;F</xref>) are within the range of healthy biological gait in the sagittal plane (Winter, <xref ref-type="bibr" rid="B51">2009</xref>), and that the joint angles of the reference participant and the biped are similar.</p>
<p>The relative horizontal and vertical excursions of biped walking were consistent with our earlier work on the benefits of haptic feedback to the hands for controlling the position of virtual feet (Karunakaran et al., <xref ref-type="bibr" rid="B25">2017</xref>). More importantly, they have the same appearance as similar measures captured from unimpaired biological gait. <xref ref-type="fig" rid="F6">Figure 6D</xref> shows that the right foot data from the biped (plotted as X and Y positions) bears a strong resemblance to unimpaired human foot data (Meesen et al., <xref ref-type="bibr" rid="B37">2006</xref>). The user feedback questionnaire indicated that biped control was easy, intuitive, and required only minimal force.</p>
<p>Current commercial exoskeletons are generally slow, with small step lengths and step heights (Kozlowski et al., <xref ref-type="bibr" rid="B28">2015</xref>). They usually rely on control mechanisms that initiate each step separately, which results in a delay between steps (Dellon and Matsuoka, <xref ref-type="bibr" rid="B9">2007</xref>; Strickland, <xref ref-type="bibr" rid="B47">2012</xref>)<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref><sup>,</sup><xref ref-type="fn" rid="fn0002"><sup>2</sup></xref> and hence a prolonged stance phase (Esquenazi et al., <xref ref-type="bibr" rid="B12">2012</xref>)<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref>. Our control method, on the other hand, provides continuous gait (smooth transition from swing to stance) with no delay between steps. It also provides the user with the ability to vary step length and height, based on the required speed. Thus, it produces a more natural gait in terms of temporal &#x00026; spatial characteristics.</p>
<p>This preliminary work shows the feasibility of our algorithm to control exoskeleton kinematics and produce symmetrical gait patterns, using hands as a controller. However, there are several limitations that will need to be addressed in future work. This paper does not address balance, which will be addressed in a separate future publication. We also acknowledge that the benefits of our control mechanism must be thoroughly evaluated with participants walking independently in a full-scale exoskeleton. A wearable exoskeleton is under development for this purpose (Al Rashdan, <xref ref-type="bibr" rid="B1">2016</xref>; Androwis et al., <xref ref-type="bibr" rid="B2">2017</xref>).</p>
<p>While our work has focused on a specific SCI population of users, we are excited about expansion of the technique to individuals with other disabling conditions, such as people with diplegic cerebral palsy who have sufficient arm control to operate our trekking poles. We also see potential for the application of variations of our control mechanism in individuals with stroke and traumatic brain injury (both for assistance with ambulation, as well as for gait therapy), and multiple sclerosis.</p>
</sec>
<sec sec-type="conclusions" id="s5">
<title>Conclusions</title>
<p>We have developed and tested an admittance-control-based user-robot control strategy that allows the user to control foot trajectories with hand-generated forces and hand-sensed foot kinematics, and thus gait in real-time. This new approach has the potential to be used with a wearable rehabilitation exoskeleton, as a control mechanism that provides users with lower leg disability complete control over gait. When implemented with wearable exoskeletons, our method has the potential to greatly improve community ambulation in individuals with lower extremity paralysis.</p>
</sec>
<sec sec-type="data-availability-statement" id="s6">
<title>Data Availability Statement</title>
<p>The datasets presented in this article are not available because New Jersey Institute of Technology Review Board approval does not include data sharing with external institutions.</p>
</sec>
<sec id="s7">
<title>Ethics Statement</title>
<p>The studies involving human participants were reviewed and approved by New Jersey Institute of Technology. The patients/participants provided their written informed consent to participate in this study.</p>
</sec>
<sec id="s8">
<title>Author Contributions</title>
<p>KK designed the algorithm, evaluated the design, analyzed the data, and drafted and finalized the manuscript. KA assisted with the design. GA reviewed the manuscript. RF assisted with the design, revised, reviewed, and finalized the manuscript. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="s9">
<title>Conflict of Interest</title>
<p>The research for this study was conducted when the authors were employed by NJIT. Currently, KA is now employed by the company Stryker Corporation, and RF consults at Really Useful Robots, LLC, in addition to being Professor Emeritus at NJIT. The companies had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. The remaining 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>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Al Rashdan</surname> <given-names>A.</given-names></name></person-group> (<year>2016</year>). <source>Hand Control of Bipedal Balance in Quiet</source> <publisher-loc>Standing</publisher-loc>: <publisher-name>Implementations for Lower Extremity Exoskeleton</publisher-name>. <publisher-loc>Newark, NJ</publisher-loc>: <publisher-name>New Jersey Institute of Technology</publisher-name>.</citation></ref>
<ref id="B2">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Androwis</surname> <given-names>G. J.</given-names></name> <name><surname>Karunakaran</surname> <given-names>K.</given-names></name> <name><surname>Nunez</surname> <given-names>E.</given-names></name> <name><surname>Michael</surname> <given-names>P.</given-names></name> <name><surname>Yue</surname> <given-names>G.</given-names></name> <name><surname>Foulds</surname> <given-names>R. A.</given-names></name></person-group> (<year>2017</year>). <article-title>Research and development of new generation robotic exoskeleton for over ground walking in individuals with mobility disorders (Novel design and control),</article-title> in <source>2017 International Symposium on Wearable Robotics and Rehabilitation</source> (<publisher-loc>Houston, TX</publisher-loc>: <publisher-name>WeRob</publisher-name>), <fpage>1</fpage>&#x02013;<lpage>2</lpage>. <pub-id pub-id-type="doi">10.1109/WEROB.2017.8383864</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Androwis</surname> <given-names>G. J.</given-names></name> <name><surname>Kwasnica</surname> <given-names>M. A.</given-names></name> <name><surname>Niewrzol</surname> <given-names>P.</given-names></name> <name><surname>Popok</surname> <given-names>P.</given-names></name> <name><surname>Fakhoury</surname> <given-names>F. N.</given-names></name> <name><surname>Sandroff</surname> <given-names>B. M.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Mobility and cognitive improvements resulted from overground robotic exoskeleton gait-training in persons with MS</article-title>. <source>Conf. Proc. IEEE. Eng. Med. Biol. Soc</source>. <volume>2019</volume>, <fpage>4454</fpage>&#x02013;<lpage>4457</lpage>. <pub-id pub-id-type="doi">10.1109/EMBC.2019.8857029</pub-id><pub-id pub-id-type="pmid">31946854</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Armstrong</surname> <given-names>T. J.</given-names></name> <name><surname>Young</surname> <given-names>J.</given-names></name> <name><surname>Woolley</surname> <given-names>C.</given-names></name> <name><surname>Ashton-Miller</surname> <given-names>J.</given-names></name> <name><surname>Kim</surname> <given-names>H.</given-names></name></person-group> (<year>2009</year>). <article-title>Biomechanical aspects of fixed ladder climbing: style, ladder tilt and carrying</article-title>. <source>Proc. Hum. Factors Ergon. Soc.</source> <volume>53</volume>, <fpage>935</fpage>&#x02013;<lpage>939</lpage>. <pub-id pub-id-type="doi">10.1177/154193120905301417</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Balasubramanian</surname> <given-names>C. K.</given-names></name> <name><surname>Bowden</surname> <given-names>M. G.</given-names></name> <name><surname>Neptune</surname> <given-names>R. R.</given-names></name> <name><surname>Kautz</surname> <given-names>S. A.</given-names></name></person-group> (<year>2007</year>). <article-title>Relationship between step length asymmetry and walking performance in subjects with chronic hemiparesis</article-title>. <source>Arch. Phys. Med. Rehabil.</source> <volume>88</volume>, <fpage>43</fpage>&#x02013;<lpage>49</lpage>. <pub-id pub-id-type="doi">10.1016/j.apmr.2006.10.004</pub-id><pub-id pub-id-type="pmid">17207674</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Birmingham</surname> <given-names>C.</given-names></name> <name><surname>Noesner</surname> <given-names>M.</given-names></name> <name><surname>Foulds</surname> <given-names>R.</given-names></name></person-group> (<year>2003</year>). <article-title>Haptic controlled functional electrical stimulation of the lower extremities,</article-title> in <source>Proceedings of the IEEE Annual Northeast Bioengineering Conference, NEBEC</source> (<publisher-loc>Newark, NJ</publisher-loc>), <fpage>287</fpage>&#x02013;<lpage>288</lpage>. <pub-id pub-id-type="doi">10.1109/NEBC.2003.1216107</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Canela</surname> <given-names>M.</given-names></name> <name><surname>del Ama</surname> <given-names>A. J.</given-names></name> <name><surname>Pons</surname> <given-names>J. L.</given-names></name></person-group> (<year>2013</year>). <article-title>Design of a pediatric exoskeleton for the rehabilitation of the physical disabilities caused by cerebral palsy,</article-title> in <source>Biosystems &#x00026; Biorobotic, Vol. 1</source>, eds J. Pons, D. Torricelli, and M. Pajaro (<publisher-loc>Berlin; Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>255</fpage>&#x02013;<lpage>258</lpage>. <pub-id pub-id-type="doi">10.1007/978-3-642-34546-3_40</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Contreras-Vidal</surname> <given-names>J. L.</given-names></name> <name><surname>Bhagat</surname> <given-names>N. A.</given-names></name> <name><surname>Brantley</surname> <given-names>J.</given-names></name> <name><surname>Cruz-Garza</surname> <given-names>J. G.</given-names></name> <name><surname>He</surname> <given-names>Y.</given-names></name> <name><surname>Manley</surname> <given-names>Q.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Powered exoskeletons for bipedal locomotion after spinal cord injury</article-title>. <source>J. Neural Eng.</source> <volume>13</volume>:<fpage>031001</fpage>. <pub-id pub-id-type="doi">10.1088/1741-2560/13/3/031001</pub-id><pub-id pub-id-type="pmid">27064508</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dellon</surname> <given-names>B.</given-names></name> <name><surname>Matsuoka</surname> <given-names>Y.</given-names></name></person-group> (<year>2007</year>). <article-title>Prosthetics, exoskeletons, and rehabilitation [grand challenges of robotics]</article-title>. <source>IEEE Robot. Autom. Mag.</source> <volume>14</volume>, <fpage>30</fpage>&#x02013;<lpage>34</lpage>. <pub-id pub-id-type="doi">10.1109/MRA.2007.339622</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dollar</surname> <given-names>A. M.</given-names></name> <name><surname>Herr</surname> <given-names>H.</given-names></name></person-group> (<year>2008</year>). <article-title>Lower extremity exoskeletons and active orthoses:challenges and state-of-the-art</article-title>. <source>IEEE Trans. Robot.</source> <volume>24</volume>, <fpage>144</fpage>&#x02013;<lpage>158</lpage>. <pub-id pub-id-type="doi">10.1109/TRO.2008.915453</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Durandau</surname> <given-names>G.</given-names></name> <name><surname>Farina</surname> <given-names>D.</given-names></name> <name><surname>As&#x000ED;n-Prieto</surname> <given-names>G.</given-names></name> <name><surname>Dimbwadyo-Terrer</surname> <given-names>I.</given-names></name> <name><surname>Lerma-Lara</surname> <given-names>S.</given-names></name> <name><surname>Pons</surname> <given-names>J. L.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Voluntary control of wearable robotic exoskeletons by patients with paresis via neuromechanical modeling</article-title>. <source>J. Neuroeng. Rehabil.</source> <volume>16</volume>:<fpage>91</fpage>. <pub-id pub-id-type="doi">10.1186/s12984-019-0559-z</pub-id><pub-id pub-id-type="pmid">31315633</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Esquenazi</surname> <given-names>A.</given-names></name> <name><surname>Talaty</surname> <given-names>M.</given-names></name> <name><surname>Packel</surname> <given-names>A.</given-names></name> <name><surname>Saulino</surname> <given-names>M.</given-names></name></person-group> (<year>2012</year>). <article-title>The Rewalk powered exoskeleton to restore ambulatory function to individuals with thoracic-level motor-complete spinal cord injury</article-title>. <source>Am. J. Phys. Med. Rehabil.</source> <volume>9</volume>, <fpage>911</fpage>&#x02013;<lpage>921</lpage>. <pub-id pub-id-type="doi">10.1097/PHM.0b013e318269d9a3</pub-id><pub-id pub-id-type="pmid">23085703</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fang</surname> <given-names>J.</given-names></name> <name><surname>Xie</surname> <given-names>Q.</given-names></name> <name><surname>Yang</surname> <given-names>G. Y.</given-names></name> <name><surname>Xie</surname> <given-names>L.</given-names></name></person-group> (<year>2017</year>). <article-title>Development and feasibility assessment of a rotational orthosis for walking with arm swing</article-title>. <source>Front. Neurosci.</source> <volume>11</volume>:<fpage>32</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2017.00032</pub-id><pub-id pub-id-type="pmid">28203142</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Farris</surname> <given-names>R. J.</given-names></name></person-group> (<year>2012</year>). <source>Design of a Powered Lower-Limb Exoskeleton and Control for Gait Assistance in Paraplegics</source>. <publisher-loc>Nashville, TN</publisher-loc>: <publisher-name>Vanderbilt University</publisher-name>.</citation></ref>
<ref id="B15">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Farris</surname> <given-names>R. J.</given-names></name> <name><surname>Quintero</surname> <given-names>H. A.</given-names></name> <name><surname>Goldfarb</surname> <given-names>M.</given-names></name></person-group> (<year>2012</year>). <article-title>Performance evaluation of a lower limb exoskeleton for stair ascent and descent with Paraplegia,</article-title> in <source>Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBS</source> (<publisher-loc>San Diego, CA</publisher-loc>), <fpage>1908</fpage>&#x02013;<lpage>1911</lpage>. <pub-id pub-id-type="doi">10.1109/EMBC.2012.6346326</pub-id><pub-id pub-id-type="pmid">23366287</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Ferris</surname> <given-names>D. P.</given-names></name> <name><surname>Lewis</surname> <given-names>C. L.</given-names></name></person-group> (<year>2009</year>). <article-title>Robotic lower limb exoskeletons using proportional myoelectric control,</article-title> in <source>Proceedings of the 31st Annual International Conference of the IEEE Engineering in Medicine and Biology Society: Engineering the Future of Biomedicine, EMBC 2009</source> (<publisher-loc>Minneapolis, MN</publisher-loc>), <fpage>2119</fpage>&#x02013;<lpage>2124</lpage>. <pub-id pub-id-type="doi">10.1109/IEMBS.2009.5333984</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ford</surname> <given-names>M. P.</given-names></name> <name><surname>Wagenaar</surname> <given-names>R. C.</given-names></name> <name><surname>Newell</surname> <given-names>K. M.</given-names></name></person-group> (<year>2007</year>). <article-title>Arm constraint and walking in healthy adults</article-title>. <source>Gait Posture</source> <volume>26</volume>, <fpage>135</fpage>&#x02013;<lpage>141</lpage>. <pub-id pub-id-type="doi">10.1016/j.gaitpost.2006.08.008</pub-id><pub-id pub-id-type="pmid">16997561</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Gancet</surname> <given-names>J.</given-names></name> <name><surname>Ilzkovitz</surname> <given-names>M.</given-names></name> <name><surname>Motard</surname> <given-names>E.</given-names></name> <name><surname>Nevatia</surname> <given-names>Y.</given-names></name> <name><surname>Letier</surname> <given-names>P.</given-names></name> <name><surname>de Weerdt</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>MINDWALKER: Going one step further with assistive lower limbs exoskeleton for SCI condition subjects,</article-title> in <source>Proceedings of the IEEE RAS and EMBS International Conference on Biomedical Robotics and Biomechatronics</source> (<publisher-loc>Rome</publisher-loc>), <fpage>1794</fpage>&#x02013;<lpage>1800</lpage>. <pub-id pub-id-type="doi">10.1109/BioRob.2012.6290688</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Glosser</surname> <given-names>G. D.</given-names></name> <name><surname>Newman</surname> <given-names>W. S.</given-names></name></person-group> (<year>2002</year>). <article-title>The implementation of a natural admittance controller on an industrial manipulator,</article-title> in <source>Proceedings of the 1994 IEEE International Conference on Robotics and Automation</source>, <fpage>1209</fpage>&#x02013;<lpage>1215</lpage>.</citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Haidegger</surname> <given-names>T.</given-names></name> <name><surname>Beny&#x000F6;</surname> <given-names>B.</given-names></name> <name><surname>Kov&#x000E1;cs</surname> <given-names>L.</given-names></name> <name><surname>Beny&#x000F6;</surname> <given-names>Z.</given-names></name></person-group> (<year>2009</year>). <article-title>Force sensing and force control for surgical robots</article-title>. <source>IFAC Proc. Vol.</source> <volume>42</volume>, <fpage>401</fpage>&#x02013;<lpage>406</lpage>. <pub-id pub-id-type="doi">10.3182/20090812-3-DK-2006.0035</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hayashi</surname> <given-names>T.</given-names></name> <name><surname>Kawamoto</surname> <given-names>H.</given-names></name> <name><surname>Sankai</surname> <given-names>Y.</given-names></name></person-group> (<year>2005</year>). <article-title>Control method of robot suit HAL working as operator&#x00027;s muscle using biological and dynamical information,</article-title> in <source>2005 IEEE/RSJ International Conference on Intelligent Robots and Systems, IROS</source> (<publisher-loc>Edmonton, Alta, AB</publisher-loc>), <fpage>3063</fpage>&#x02013;<lpage>3068</lpage>. <pub-id pub-id-type="doi">10.1109/IROS.2005.1545505</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hamill</surname> <given-names>J.</given-names></name> <name><surname>Knudsen</surname> <given-names>K. M.</given-names></name> <name><surname>Derrick</surname> <given-names>T. R.</given-names></name></person-group> (<year>2013</year>). <source>Biological Basis of Human Movement, 4th Ed.</source> <publisher-loc>Baltimore, MD</publisher-loc>: <publisher-name>Lippincott Williams &#x00026; Wilkins</publisher-name>.</citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Karunakaran</surname> <given-names>K.</given-names></name> <name><surname>Abbruzzese</surname> <given-names>K.</given-names></name> <name><surname>Xu</surname> <given-names>H.</given-names></name> <name><surname>Ehrenberg</surname> <given-names>N.</given-names></name> <name><surname>Foulds</surname> <given-names>R.</given-names></name></person-group> (<year>2014</year>). <article-title>Haptic proprioception in a virtual locomotor task</article-title>. <source>Conf. IEEE Eng. Med. Biol. Soc. IEEE Eng. Med. Biol. Soc.</source> <volume>2014</volume>, <fpage>3594</fpage>&#x02013;<lpage>3597</lpage>. <pub-id pub-id-type="doi">10.1109/EMBC.2014.6944400</pub-id><pub-id pub-id-type="pmid">25570768</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Karunakaran</surname> <given-names>K.</given-names></name> <name><surname>Ehrenberg</surname> <given-names>N.</given-names></name> <name><surname>Cheng</surname> <given-names>J.</given-names></name> <name><surname>Nolan</surname> <given-names>K. J.</given-names></name></person-group> (<year>2019</year>). <article-title>Effects of robotic exoskeleton gait training on an adolescent with brain injury,</article-title> in <source>IEEE 41h Annual International Conference on Engineering in Medicine and Biology Society (EMBC)</source> (<publisher-loc>Berlin</publisher-loc>), <fpage>4445</fpage>&#x02013;<lpage>4448</lpage>. <pub-id pub-id-type="doi">10.1109/EMBC.2019.8856787</pub-id><pub-id pub-id-type="pmid">31946852</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Karunakaran</surname> <given-names>K. K.</given-names></name> <name><surname>Abbruzzese</surname> <given-names>K. M.</given-names></name> <name><surname>Xu</surname> <given-names>H.</given-names></name> <name><surname>Foulds</surname> <given-names>R. A.</given-names></name></person-group> (<year>2017</year>). <article-title>The importance of haptics in generating exoskeleton gait trajectory using alternate motor inputs</article-title>. <source>IEEE Trans. Neural Syst. Rehabil. Eng.</source> <volume>25</volume>, <fpage>2328</fpage>&#x02013;<lpage>2335</lpage>. <pub-id pub-id-type="doi">10.1109/TNSRE.2017.2726538</pub-id><pub-id pub-id-type="pmid">28715331</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kianimajd</surname> <given-names>A.</given-names></name> <name><surname>Ruano</surname> <given-names>M. G.</given-names></name> <name><surname>Carvalho</surname> <given-names>P.</given-names></name> <name><surname>Henriques</surname> <given-names>J.</given-names></name> <name><surname>Rocha</surname> <given-names>T.</given-names></name> <name><surname>Paredes</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Comparison of different methods of measuring similarity in physiologic time series</article-title>. <source>IFAC-PapersOnLine</source> <volume>50</volume>, <fpage>11005</fpage>&#x02013;<lpage>11010</lpage>. <pub-id pub-id-type="doi">10.1016/j.ifacol.2017.08.2479</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Kilicarslan</surname> <given-names>A.</given-names></name> <name><surname>Prasad</surname> <given-names>S.</given-names></name> <name><surname>Grossman</surname> <given-names>R. G.</given-names></name> <name><surname>Contreras-Vidal</surname> <given-names>J. L.</given-names></name></person-group> (<year>2013</year>). <article-title>High accuracy decoding of user intentions using EEG to control a lower-body exoskeleton,</article-title> in <source>Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBS</source> (<publisher-loc>Osaka</publisher-loc>), <fpage>5606</fpage>&#x02013;<lpage>5609</lpage>. <pub-id pub-id-type="doi">10.1109/EMBC.2013.6610821</pub-id><pub-id pub-id-type="pmid">24111008</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kozlowski</surname> <given-names>A. J.</given-names></name> <name><surname>Bryce</surname> <given-names>T. N.</given-names></name> <name><surname>Dijkers</surname> <given-names>M. P.</given-names></name></person-group> (<year>2015</year>). <article-title>Time and effort required by persons with spinal cord injury to learn to use a powered exoskeleton for assisted walking</article-title>. <source>Top. Spinal Cord Inj. Rehabil.</source> <volume>21</volume>, <fpage>110</fpage>&#x02013;<lpage>121</lpage>. <pub-id pub-id-type="doi">10.1310/sci2102-110</pub-id><pub-id pub-id-type="pmid">26364280</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kozlowski</surname> <given-names>A. J.</given-names></name> <name><surname>Fabian</surname> <given-names>M.</given-names></name> <name><surname>Lad</surname> <given-names>D.</given-names></name> <name><surname>Delgado</surname> <given-names>A. D.</given-names></name></person-group> (<year>2017</year>). <article-title>Feasibility and safety of a powered exoskeleton for assisted walking for persons with multiple sclerosis: a single-group preliminary study</article-title>. <source>Arch. Phys. Med. Rehabil.</source> <volume>98</volume>, <fpage>1300</fpage>&#x02013;<lpage>1307</lpage>. <pub-id pub-id-type="doi">10.1016/j.apmr.2017.02.010</pub-id><pub-id pub-id-type="pmid">28315666</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kwak</surname> <given-names>N. S.</given-names></name> <name><surname>M&#x000FC;ller</surname> <given-names>K. R.</given-names></name> <name><surname>Lee</surname> <given-names>S. W.</given-names></name></person-group> (<year>2015</year>). <article-title>A lower limb exoskeleton control system based on steady state visual evoked potentials</article-title>. <source>J. Neural Eng.</source> <volume>12</volume>:<fpage>056009</fpage>. <pub-id pub-id-type="doi">10.1088/1741-2560/12/5/056009</pub-id><pub-id pub-id-type="pmid">26291321</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>La Scaleia</surname> <given-names>V.</given-names></name> <name><surname>Sylos-Labini</surname> <given-names>F.</given-names></name> <name><surname>Hoellinger</surname> <given-names>T.</given-names></name> <name><surname>Wang</surname> <given-names>L.</given-names></name> <name><surname>Cheron</surname> <given-names>G.</given-names></name> <name><surname>Lacquaniti</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Control of leg movements driven by EMG activity of shoulder muscles</article-title>. <source>Front. Hum. Neurosci.</source> <volume>8</volume>:<fpage>838</fpage>. <pub-id pub-id-type="doi">10.3389/fnhum.2014.00838</pub-id><pub-id pub-id-type="pmid">25368569</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lebedev</surname> <given-names>M. A.</given-names></name> <name><surname>Nicolelis</surname> <given-names>M. A. L.</given-names></name></person-group> (<year>2017</year>). <article-title>Brain-machine interfaces: from basic science to neuroprostheses and neurorehabilitation</article-title>. <source>Physiol. Rev.</source> <volume>97</volume>, <fpage>767</fpage>&#x02013;<lpage>837</lpage>. <pub-id pub-id-type="doi">10.1152/physrev.00027.2016</pub-id><pub-id pub-id-type="pmid">28275048</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S.</given-names></name> <name><surname>Sankai</surname> <given-names>Y.</given-names></name></person-group> (<year>2002</year>). <article-title>Power assist control for leg with HAL-3 based on virtual torque and impedance adjustment</article-title>. <source>IEEE Int. Conf. Syst. Man Cybern.</source> <volume>4</volume>:<fpage>6</fpage>. <pub-id pub-id-type="doi">10.1109/ICSMC.2002.1173329</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S.</given-names></name> <name><surname>Sankai</surname> <given-names>Y.</given-names></name></person-group> (<year>2003</year>). <article-title>Power assist control for walking aid with HAL-3 based on EMG and impedance adjustment around knee joint</article-title>. <source>IEEE/RSJ Int. Conf. Intell. Robot. Syst.</source> <volume>2</volume>, <fpage>1499</fpage>&#x02013;<lpage>1504</lpage>. <pub-id pub-id-type="doi">10.1109/IRDS.2002.1043967</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lerner</surname> <given-names>Z. F.</given-names></name> <name><surname>Damiano</surname> <given-names>D. L.</given-names></name> <name><surname>Park</surname> <given-names>H. S.</given-names></name> <name><surname>Gravunder</surname> <given-names>A. J.</given-names></name> <name><surname>Bulea</surname> <given-names>T. C.</given-names></name></person-group> (<year>2017</year>). <article-title>A robotic exoskeleton for treatment of crouch gait in children with cerebral palsy: design and initial application</article-title>. <source>IEEE Trans. Neural Syst. Rehabil. Eng.</source> <volume>25</volume>, <fpage>650</fpage>&#x02013;<lpage>659</lpage>. <pub-id pub-id-type="doi">10.1109/TNSRE.2016.2595501</pub-id><pub-id pub-id-type="pmid">27479974</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>L&#x000F3;pez-Larraz</surname> <given-names>E.</given-names></name> <name><surname>Trincado-Alonso</surname> <given-names>F.</given-names></name> <name><surname>Rajasekaran</surname> <given-names>V.</given-names></name> <name><surname>P&#x000E9;rez-Nombela</surname> <given-names>S.</given-names></name> <name><surname>del-Ama</surname> <given-names>A. J.</given-names></name> <name><surname>Aranda</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Control of an ambulatory exoskeleton with a brain-machine interface for spinal cord injury gait rehabilitation</article-title>. <source>Front. Neurosci.</source> <volume>10</volume>:<fpage>359</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2016.00359</pub-id><pub-id pub-id-type="pmid">27536214</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Meesen</surname> <given-names>R. L. J.</given-names></name> <name><surname>Wenderoth</surname> <given-names>N.</given-names></name> <name><surname>Temprado</surname> <given-names>J. J.</given-names></name> <name><surname>Summers</surname> <given-names>J. J.</given-names></name> <name><surname>Swinnen</surname> <given-names>S. P.</given-names></name></person-group> (<year>2006</year>). <article-title>The coalition of constraints during coordination of the ipsilateral and heterolateral limbs</article-title>. <source>Exp. Brain Res.</source> <volume>174</volume>, <fpage>367</fpage>&#x02013;<lpage>375</lpage>. <pub-id pub-id-type="doi">10.1007/s00221-006-0471-1</pub-id><pub-id pub-id-type="pmid">16819649</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Meyns</surname> <given-names>P.</given-names></name> <name><surname>Bruijn</surname> <given-names>S. M.</given-names></name> <name><surname>Duysens</surname> <given-names>J.</given-names></name></person-group> (<year>2013</year>). <article-title>The how and why of arm swing during human walking</article-title>. <source>Gait Posture</source> <volume>38</volume>, <fpage>555</fpage>&#x02013;<lpage>562</lpage>. <pub-id pub-id-type="doi">10.1016/j.gaitpost.2013.02.006</pub-id><pub-id pub-id-type="pmid">23489950</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Murray</surname> <given-names>S. A.</given-names></name> <name><surname>Ha</surname> <given-names>K. H.</given-names></name> <name><surname>Hartigan</surname> <given-names>C.</given-names></name> <name><surname>Goldfarb</surname> <given-names>M.</given-names></name></person-group> (<year>2015</year>). <article-title>An assistive control approach for a lower-limb exoskeleton to facilitate recovery of walking following stroke</article-title>. <source>IEEE Trans. Neural Syst. Rehabil. Eng.</source> <volume>23</volume>, <fpage>441</fpage>&#x02013;<lpage>449</lpage>. <pub-id pub-id-type="doi">10.1109/TNSRE.2014.2346193</pub-id><pub-id pub-id-type="pmid">25134084</pub-id></citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nicolelis</surname> <given-names>M. A. L.</given-names></name></person-group> (<year>2003</year>). <article-title>Brain&#x02013;machine interfaces to restore motor function and probe neural circuits</article-title>. <source>Nat. Rev. Neurosci.</source> <volume>4</volume>, <fpage>417</fpage>&#x02013;<lpage>422</lpage>. <pub-id pub-id-type="doi">10.1038/nrn1105</pub-id><pub-id pub-id-type="pmid">12728268</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Patan&#x000E9;</surname> <given-names>F.</given-names></name> <name><surname>Rossi</surname> <given-names>S.</given-names></name> <name><surname>Del Sette</surname> <given-names>F.</given-names></name> <name><surname>Taborri</surname> <given-names>J.</given-names></name> <name><surname>Cappa</surname> <given-names>P.</given-names></name></person-group> (<year>2017</year>). <article-title>WAKE-up exoskeleton to assist children with cerebral palsy: design and preliminary evaluation in level walking</article-title>. <source>IEEE Trans. Neural Syst. Rehabil. Eng.</source> <volume>25</volume>, <fpage>906</fpage>&#x02013;<lpage>916</lpage>. <pub-id pub-id-type="doi">10.1109/TNSRE.2017.2651404</pub-id><pub-id pub-id-type="pmid">28092566</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Patterson</surname> <given-names>K. K.</given-names></name> <name><surname>Parafianowicz</surname> <given-names>I.</given-names></name> <name><surname>Danells</surname> <given-names>C. J.</given-names></name> <name><surname>Closson</surname> <given-names>V.</given-names></name> <name><surname>Verrier</surname> <given-names>M. C.</given-names></name> <name><surname>Richard Staines</surname> <given-names>W.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Gait asymmetry in community-ambulating stroke survivors</article-title>. <source>Arch. Phys. Med. Rehabil.</source> <volume>89</volume>, <fpage>304</fpage>&#x02013;<lpage>310</lpage>. <pub-id pub-id-type="doi">10.1016/j.apmr.2007.08.142</pub-id><pub-id pub-id-type="pmid">18226655</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pearson</surname> <given-names>K. G.</given-names></name></person-group> (<year>2003</year>). <article-title>Common principles of motor control in vertebrates and invertebrates</article-title>. <source>Annu. Rev. Neurosci.</source> <volume>16</volume>, <fpage>265</fpage>&#x02013;<lpage>297</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.ne.16.030193.001405</pub-id><pub-id pub-id-type="pmid">8460894</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quintero</surname> <given-names>H. A.</given-names></name> <name><surname>Farris</surname> <given-names>R. J.</given-names></name> <name><surname>Goldfarb</surname> <given-names>M.</given-names></name></person-group> (<year>2012</year>). <article-title>A method for the autonomous control of lower limb exoskeletons for persons with Paraplegia</article-title>. <source>J. Med. Device.</source> <volume>6</volume>:<fpage>041003</fpage>. <pub-id pub-id-type="doi">10.1115/1.4007181</pub-id><pub-id pub-id-type="pmid">23505407</pub-id></citation></ref>
<ref id="B45">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>DeMarco</surname> <given-names>R.</given-names></name> <name><surname>Foulds</surname> <given-names>R.</given-names></name></person-group> (<year>2002</year>). <article-title>Sensory-motor substitution for improved ambulation.</article-title> in <source>IEEE 28th Annual Northeast Bioengineering Conference</source> (<publisher-loc>Philadelphia, PA</publisher-loc>), 147-148.</citation></ref>
<ref id="B46">
<citation citation-type="web"><person-group person-group-type="author"><collab>Shepherd Center</collab></person-group> (<year>2020</year>). <source>Understanding Spinal Cord Injury</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.spinalinjury101.org/details">https://www.spinalinjury101.org/details</ext-link> (accessed October 01, 2020).</citation></ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Strickland</surname> <given-names>E.</given-names></name></person-group> (<year>2012</year>). <article-title>Good-bye, wheelchair</article-title>. <source>IEEE Spectr</source>. <volume>49</volume>, <fpage>30</fpage>&#x02013;<lpage>32</lpage>. <pub-id pub-id-type="doi">10.1109/MSPEC.2012.6117830</pub-id></citation></ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Van Der Linde</surname> <given-names>R.</given-names></name> <name><surname>Lammertse</surname> <given-names>P.</given-names></name> <name><surname>Frederiksen</surname> <given-names>E. B.</given-names></name></person-group> (<year>2002</year>). <article-title>The HapticMaster, a new high-performance haptic interface,</article-title> in <source>Proceedings of EuroHaptic (Edinburgh)</source>, <fpage>1</fpage>&#x02013;<lpage>5</lpage>.</citation></ref>
<ref id="B49">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Van Riel</surname> <given-names>N.</given-names></name></person-group> (<year>2012</year>). <source>Speeding Up Simulations of ODE Models in Matlab Using CVode and MEX Files</source>. <publisher-loc>Eindhoven, NL</publisher-loc>: <publisher-name>Eindhoven University of Technology</publisher-name>.</citation></ref>
<ref id="B50">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Vaughan</surname> <given-names>C. L.</given-names></name> <name><surname>Davis</surname> <given-names>B. L.</given-names></name> <name><surname>O&#x00027;Connor</surname> <given-names>J. C.</given-names></name></person-group> (<year>1999</year>). <article-title>The three-dimensional and cyclic nature of gait.</article-title> in <source>Dynamics of Human Gait, Vol. 168</source> (<publisher-loc>Cape Town</publisher-loc>: <publisher-name>Kiboho Publishers</publisher-name>), <fpage>16</fpage>&#x02013;<lpage>17</lpage>.</citation></ref>
<ref id="B51">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Winter</surname> <given-names>D. A.</given-names></name></person-group> (<year>2009</year>). <source>Biomechanics and Motor Control of Human Movement.</source> <publisher-loc>Hoboken, NJ</publisher-loc>: <publisher-name>John Wiley &#x00026; Sons</publisher-name>. <pub-id pub-id-type="doi">10.1002/9780470549148</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="fn0001"><p><sup>1</sup><italic>Ekso Bionics</italic>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://eksobionics.com/">https://eksobionics.com/</ext-link></p></fn>
<fn id="fn0002"><p><sup>2</sup><italic>RexBionics</italic>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.rexbionics.com/">https://www.rexbionics.com/</ext-link></p></fn>
<fn id="fn0003"><p><sup>3</sup><italic>Robotis</italic>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://www.robotis.com/xe/dynamixel_en">http://www.robotis.com/xe/dynamixel_en</ext-link></p></fn>
<fn id="fn0004"><p><sup>4</sup><ext-link ext-link-type="uri" xlink:href="https://www.universal-robots.com/">https://www.universal-robots.com/</ext-link></p></fn>
</fn-group>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> This work was supported by Gustavus and Louise Pfeiffer Research Foundation, and the Rehabilitation Engineering Research Center on Wearable Robots (NIDILRR grant 90RE5021).</p>
</fn>
</fn-group>
</back>
</article>