Optimization of selective stimulation parameters for multicontact electrodes
 Lee E Fisher^{1}Email author,
 Dustin J Tyler^{1, 2} and
 Ronald J Triolo^{1, 2}
DOI: 10.1186/174300031025
© Fisher et al; licensee BioMed Central Ltd. 2013
Received: 25 April 2012
Accepted: 21 February 2013
Published: 27 February 2013
Abstract
Background
Multicontact stimulating electrodes are gaining acceptance as a means for interfacing with the peripheral nervous system. These electrodes can potentially activate many independent populations of motor units within a single peripheral nerve, but quantifying their recruitment properties and the overlap in stimulation between contacts is difficult and time consuming. Further, current methods for quantifying overlap between contacts are ambiguous and can lead to suboptimal selective stimulation parameters. This study describes a novel method for optimizing stimulation parameters for multicontact peripheral stimulating electrodes to produce strong, selective muscle contractions. The method is tested with fourcontact spiral nervecuff electrodes implanted on bilateral femoral nerves of two individuals with spinal cord injury, but it is designed to be extendable to other electrode technologies with higher densities of contacts.
Methods
To optimize selective stimulation parameters for multicontact electrodes, first, recruitment and overlap are characterized for all contacts within an electrode. Recruitment is measured with the twitch response to single stimulus pulses, and overlap between pairs of contacts is quantified by the deviation in their combined response from linear addition of individual responses. Simple mathematical models are fit to recruitment and overlap data, and a cost function is defined to maximize recruitment and minimize overlap between all contacts.
Results
Results are presented for fourcontact nervecuff electrodes stimulating bilateral femoral nerves of two human subjects with spinal cord injury. Knee extension moments between 11.6 and 43.2 Nm were achieved with selective stimulation through multiple contacts of each nervecuff with less than 10% overlap between pairs of contacts. The overlap in stimulation measured in response to selective stimulation parameters was stable at multiple repeated time points after implantation.
Conclusions
These results suggest that the method described here can provide an automated means of determining stimulus parameters to achieve strong muscle contractions via selective stimulation through multicontact peripheral nerve electrodes.
Keywords
Functional neuromuscular stimulation Selectivity Multicontact electrodesBackground
Multicontact stimulating electrodes have been gaining acceptance as a means for interfacing with peripheral nerves in functional neuromuscular stimulation (FNS) systems [1–3]. These electrodes allow for a high density of contacts to be placed around or in peripheral nerves to independently activate multiple fascicles and motor units.
Independent activation can allow for the control of multiple functions with a single electrode and for recruitment of multiple populations of agonist motor units within a single muscle [4, 5]. Control of multiple functions with a single electrode could reduce the number of implantation sites required to produce a variety of functional joint moments for FNS systems [5]. Control of multiple agonist motor unit populations with a single electrode could allow for better control of joint moment by varying the number of motor units recruited by the electrode [4]. Further, by alternating stimulation of multiple agonist motor unit populations, it may be possible to reduce stimulation duty cycle and prolong the time of muscle contractions before the onset of fatigue [6]. To take advantage of these benefits of multicontact electrodes, it is crucial to be able to measure and minimize overlap in recruitment between the contacts within these electrodes. Overlap in recruitment may limit the options available to control the response to stimulation. Repeated activation of the same motor units by multiple contacts may also contribute to rapid fatigue.
Determining selective stimulation parameters is a complex problem when multiple contacts activate agonist populations of motor units. It is difficult or impossible to separate populations of agonist motor units based on their lines of action or resultant joint moments. Ideally, if the spatial relationship between the electrode and the fascicles within the nerve was known a priori, it would be possible to select stimulation parameters based on computational models and estimates of activation of the neural tissue. Peripheral nerves, however, have a highly branched structure, and there is a high degree of persontoperson variability in the anatomy of those nerves, both in branching structure and in the location of fascicles that innervate specific muscles [7, 8]. Because current stateoftheart imaging technology is insufficient to accurately visualize the fascicular structure of nerves in vivo, it is not possible to know what the spatial relationship is between the contacts within a nervecuff electrode and the fascicles and respective muscles the electrode will activate. Stimulation overlap must, therefore, be inferred from indirect measurements [9]. Electromyogram recordings are commonly used to infer selectivity, but the strong dependence on the spatial relationship between the recording electrodes and the activated motor units limits the practical utility of this method and can lead to over or underestimation of selectivity [4, 10].
Another method for quantifying selectivity takes advantage of the concept that motor units in their absolute refractory period, which lasts between 1.5 and 2.1 ms, will not respond to stimulation [1, 3, 6, 11–15]. Therefore, if a stimulus pulse is applied through one electrode contact within 2.1 ms after a pulse is applied through another electrode contact, motor units activated by the first stimulus will not respond to the second stimulus. This means that, if there is overlap in the stimulation fields of two contacts, the resultant force generated when one contact is activated within 2.1 ms after the other will be less than the linear sum of the individual forces when each contact is activated separately. Conversely, the forces generated by stimulating two completely independent populations of motor units will add linearly, even if one population is stimulated while neurons from the other population are in their absolute refractory period [14]. By varying stimulation parameters including pulse amplitude and pulse width, it may be possible to minimize the deviation from linear addition, and thereby minimize stimulation overlap, while maximizing the magnitude of the force generated by each independent motor unit population to selectively produce strong muscle contractions to lock the knees and improve FNSassisted standing.
While this method provides a useful way to quantify overlap between two contacts, it does not easily scale to larger numbers of contacts. To create a metric of overall overlap for an entire multicontact electrode, others have used this method to quantify overlap between pairs of contacts within the electrode, and then averaged all of those overlaps [3, 6]. While this method provides some insight into the general amount of overlap for a multicontact electrode, it does not provide a clear means of tuning stimulation parameters for each individual contact to reduce overlap while generating functionally useful stimulated joint moments. Furthermore, as the number of contacts within the electrode increases, the number of pairwise combinations of contacts that must be considered for this method increases exponentially. This can quickly lead to impossibly large data sets as the number of contacts increases.
In this report, we propose a method for choosing optimal selective stimulation parameters for multicontact electrodes that minimizes overlap between adjacent contacts while maximizing the joint moment produced by stimulating through each contact. We rely on the method described above to quantify overlap between pairs of contacts, with the addition of a set of mathematical models to reduce the data requirements for characterizing the electrodes, and a cost function that acts to minimize all pairwise overlaps while simultaneously maximizing all joint moments. We show that this method can efficiently characterize overlap and selectivity for multicontact electrodes to produce strong muscle contractions with little or no overlap between stimulated motor unit populations. While the method is designed to accommodate electrodes with high densities of contacts, we test it clinically with the fourcontact spiral nervecuff electrode developed at Case Western Reserve University (CWRU). We demonstrate that the method can select stimulation parameters that generate strong contractions with low overlap between contacts in this electrode, and that the selective responses are stable over months after implantation.
Methods
The process we propose to quantify and optimize selective stimulation for multicontact electrodes consists of four fundamental steps. First, the response to stimulation through the multicontact electrode and the overlap between pairs of contacts are quantified. These responses are the forces generated by muscle twitches, elicited by single stimulus pulses. Muscle twitches are less likely to cause fatigue than tetanic contractions and can be collected more quickly. Next, the relationship between twitch and tetanic responses to stimulation is quantified. This relationship provides a scaling factor so that the twitch responses can be converted to more functionally relevant tetanic responses. Third, mathematical models are fit to the scaled recruitment and overlap data. These models serve the dual purposes of reducing the size of the data set required for optimization and providing a mathematical framework over which optimization can be performed. Finally, the scaled recruitment and overlap models are used as inputs to a cost function that can be minimized to provide optimal selective stimulation parameters.
Subject selection and multicontact electrodes
Informed consent was acquired prior to participation in any experiments, and all experimental protocols were approved by the Institutional Review Board of MetroHealth Medical Center, Cleveland, OH.
Recruitment and overlap characterization
For the first step in optimizing selective stimulation parameters, the response to stimulation and the overlap between pairs of contacts were characterized. With the knee fixed at 20° of flexion and one axis of a 6 degreeoffreedom load cell (JR3, Woodland, CA) aligned with the knee joint center, isometric knee extension moment was recorded in response to stimulus pulses applied to the femoral nerve through each contact of the cuff electrodes. Data were lowpass filtered at 31.25 Hz and sampled at 150 Hz.
To characterize the response to stimulation, pulse width modulated recruitment curves were collected. To characterize overlap between pairs of contacts, a stimulus pulse was applied through one contact, followed by a 2 ms time delay, and then a pulse through a second contact. The pulse widths of all stimuli were varied between 1 and 255 us.
For some multicontact electrodes, it is possible to reduce the size of the data set by eliminating pairwise combinations of contacts that are not adjacent to each other, since elimination of overlap between adjacent contacts will also eliminate overlap between nonadjacent contacts. For example, in the case of an eight contact Flat Interface Nerve Electrode, this would reduce the number of possible pairwise combinations from 28 to 16 [21, 22]. In the case of the fourcontact nervecuff, all contacts are adjacent to each other, so all six pairwise combinations were considered in this study.
Twitch/tetanic relationship
While the twitch response to stimulation can be collected quickly and with less likelihood of causing fatigue, the tetanic response to stimulation is functionally relevant. Studies in both animals and humans have demonstrated that there is a linear relationship between the shape of isometric twitch and tetanic recruitment curves, characterized by a single scalar multiplier [4, 23]. To quantify this scaling factor, twitch and tetanic responses to stimulation were recorded with the knee held in 20° of flexion. The ratio of the maximum twitch and tetanic responses was used as a scaling factor.
Mathematical models of recruitment and overlap
Fitting mathematical models to recruitment and overlap data reduces the size of the data set required for characterizing the electrodes while also providing a framework for optimization of stimulation parameters. To determine the best form of functions for recruitment and overlap, mathematical functions were fit to knee extension moment data in response to stimulation through individual contacts or pairs of contacts within each electrode. For each contact or pair of contacts, 48 data points were collected, with subsets of 32 recruitment or overlap data points fitted to a variety of models and separate subsets of 16 data points used to test for goodnessoffit (GOF). For recruitment data, 1^{st} through 5^{th} order polynomial, sigmoid, Gaussian, and Gompertz functions were tested. For overlap data, which are twodimensional since pulse width can be controlled independently for both contacts in a pairwise combination, 1^{st} through 5^{th} order twodimensional polynomials were tested. To determine GOF, coefficients of determination (R^{2}), and the corrected Akaike Information Criterion (AICc) were calculated for each model. AICc is a measure of how well a model fits a set of data relative to the number of parameters in that model [24]. Models that achieved the best fits were selected for implementation in the optimization described in the following section.
Optimization of selective stimulation parameters
where PW is an Ndimensional vector of pulse widths of stimulus pulses for an Ncontact electrode, O_{ T } quantifies the overlap of all contacts within the electrode, M_{ T } quantifies the joint moment generated by all contacts within the electrode, and ω_{ 0 } and ω_{ 1 } are weighting factors. By minimizing this cost function, M_{ T } will be maximized and O_{ T } will be minimized, producing the strongest possible contractions with the least possible overlap in stimulation between contacts.
where M_{ i } is the moment generated when stimulating through contact i, which is described mathematically by the model function previously fit to recruitment data. The sum of these functions is divided by the sum of the maxima of the functions to normalize the joint moment term. Because of this normalization, the overall joint moment term ranges between 0 and 1, while the summation of individual moments allows stronger joint moments to be weighted more heavily than weaker moments.
where M_{iUj} is the moment generated when stimulating through two contacts with a short time delay, M_{ i∩j } is the overlap between contacts i and j, and M_{ i } and M_{ j } are the moments generated when stimulating through contacts i and j, respectively. By rearranging these terms, the overlap in stimulation can be calculated for a pair of contacts within an electrode.
which ranges between 0 when there is no overlap and 1 when there is 100% overlap. Note also that in this equation, i≠j. During this study, a 2 second interpulse interval was used when stimulus pulses were applied to a pair of contacts within an electrode.
Since both O_{ T } and M_{ T } are normalized, the weighting factors ω_{ 0 } and ω_{ 1 } can be used to emphasize larger joint moments or lower overlap, depending on the particular application. For this study, the terms were weighted equally, i.e., ω_{ 0 } = ω_{ 1 }.
To ensure that either sufficiently large joint moments or sufficiently small overlaps are achieved, a linearly increasing penalty was added to the cost function if joint moment for any contact was less than 5 Nm or overlap was greater than 10% between any two contacts. The slope of these penalties can be tuned to tighten or relax the constraints on minimum joint moment or maximum overlap.
A direct search optimization algorithm (Matlab, Natick, MA) was used to find the minimum of the cost function and the optimal set of pulse widths for selective stimulation.
Stability of overlap
To achieve reliable control in an FNS system, it is important that the response to stimulation is stable over time. Selectivity is a function of both muscle strength and stimulation overlap. While it is expected that muscle strength will change over time as the user exercises and builds muscle mass, if the interface between the electrode and nerve is stable, overlap is less likely to change over time. To test the stability of stimulation, the amount of overlap between pairs of contacts was quantified at multiple time points after implantation. At the first time point, optimal stimulation parameters were chosen by the algorithm described above, and at each subsequent time point, overlap in stimulation between all pairs of contacts was measured using those same stimulation parameters. A onetailed Student's ttest was used to determine if the mean overlap for any pair of contacts was greater than 10%.
Results
Electrode characterization and mathematical models
is the most appropriate model of the overlap in stimulation between two contacts (R^{2} = 0.92 ± 0.11) with the highest average AICc ranking.
Twitch/tetanic relationship
Optimization of selective stimulation parameters
Note that, in Subject 1, only two of four contacts, and in Subject 2, only three of four contacts for either electrode have nonzero stimulus parameters. The results of the optimization consistently demonstrated that removing these one or two contacts from the cost function produced significantly higher joint moments with less overlap than if all contacts were included or other contacts were removed. For contact 3 of the left nervecuff electrode in Subject 2, there was never any measurable motor response to stimulation, indicating that it was probably located over sensory axons or connective tissue rather than fascicles containing motor axons.
Stability of overlap
Discussion
The methods presented in this study are designed to optimize stimulation parameters for multicontact peripheral stimulating electrodes. Results presented here demonstrate that this method can select a set of stimulation parameters that provide strong muscle contractions with low overlap. One of many applications of this technique is in neuroprostheses for standing after paralysis. In the case of these four nervecuff electrodes examined in this study, the optimization determined that it is possible for two or three contacts to produce at least 11.6 Nm, but possibly as much as 43.2 Nm of knee extension moment with less than 10% overlap between pairs of contacts. Biomechanical studies estimate that as much as 27% body weight (BW) of knee extension moment is required to keep the knees locked during standing [25]. The sum of selective moments for any one of the nervecuff electrodes in this study is greater than these requirements. It is also possible to tailor the results of the optimization by adjusting weighting factors in the cost function to favor larger joint moments at the cost of tolerating more overlap, if necessary. The technique can be applied to any number of realworld situations utilizing multicontact or multiple singlecontact electrodes activating synergistic muscle fiber populations, such as diaphragmatic pacing, hand grasping, or stepping. The methodology can readily be adapted to any electrode technology including high density penetrating arrays or multiple musclebased or even surface electrodes.
While the optimization method described here was applied to an electrode with only four contacts, it was designed to be scalable to much higher density electrodes. By focusing only on overlap between adjacent pairs of electrodes, and by fitting mathematical models to both overlap and recruitment, it should be possible to select optimal stimulation parameters for electrodes with many more contacts without becoming prohibitively data intensive. If nonadjacent pairwise combinations of contacts were ignored and a limited data set was used to fit the Gompertz and thirdorder polynomial models, which have three and ten parameters, respectively, the method described above would require approximately 13 hours to completely optimize selective stimulation for a 100contact Utah slant electrode array, as compared to 11.6 days otherwise [3]. This time requirement is a worst case scenario that only applies if every contact within the electrode activates an agonist population of motor units. In reality, it is likely that many contacts could be removed from the optimization because they activate either sensory neurons or nonagonist populations of motor units. Furthermore, the method used here to identify the relationship between twitch and tetanic recruitment was chosen because of its simplicity, but not its efficiency. In fact, of the estimated 13 hours required to optimize a 100contact Utah array, nearly 70% of that time would be devoted to characterizing the relationship between twitch and tetanic responses. Other methods described by Durfee, et al. [23], require significantly less time to characterize the relationship between twitch and tetanic responses, and could reduce the worstcase scenario time to completely optimize a 100contact electrode to approximately six hours.
As presented, the optimization method relies on pulse width modulation, but this is largely a result of the limitations of the implanted stimulator used in this study. While it may be necessary to rely on different mathematical models of recruitment and overlap data, the method should be applicable to pulse amplitude modulation or even modulation of both pulse width and amplitude.
It should be noted that in the case of all four nervecuff electrodes presented here, the optimization produced better results if either one or two contacts were eliminated from the optimization. In the case of Subject 1, the two contacts used in the optimization sit opposite one another around the circumference of the nerve, so it is reasonable to expect that they would have less overlap with one another than contacts that are directly adjacent. In the case of Subject 2, one of the contacts that was excluded never demonstrated any motor response to stimulation, which suggests the contact may have been located near sensory neurons or connective tissue.
The results of the optimization were stable over time, with the amount of overlap between the included pairs of contacts remaining largely constant. This suggests that the nervecuff electrode provides a stable interface with the nerve, and that only occasional retuning of stimulation parameters would be required to account for changes in muscle strength during clinical implementation of selective stimulation.
The results of this study demonstrate that it is possible to efficiently determine optimal selective stimulation parameters for multicontact electrodes. Based on these results, it may be possible to use the selectivity achieved here to produce finer control of motor output by recruiting motor units from a small set of contacts for lowforce tasks and recruiting from additional contacts for highforce tasks. Additionally, by alternating stimulation between multiple independent populations, it may be possible to reduce the duty cycle of stimulation while maintaining a constant joint moment, and thereby delay the onset of fatigue.
Conclusions
This study presents a new method for optimizing stimulation parameters for multicontact peripheral stimulating electrodes. By collecting twitch responses to stimulation and fitting mathematical models to recruitment and overlap data, the method reduces the data requirements for characterizing and optimizing selective stimulation. The use of a cost function that includes terms representing both recruitment and pairwise overlap for all contacts within the electrode allows for maximization of the moments generated by all contacts while simultaneously minimizing the overlaps between all pairs of contacts. This method allows for an objective and automated means of selecting stimulation parameters for electrodes with high densities of contacts, where manual selection of stimulation parameters would be prohibitively time intensive.
The results of this study also suggest that it is possible to generate strong contractions with little or no overlap between contacts within a fourcontact spiral nervecuff electrode. Overlap in stimulation was also found to be stable over months after implantation, with little variation in overlap in response to the same stimulus parameters applied at multiple time points after implantation.
Abbreviations
 AICc:

corrected Akaike information criterion
 CWRU:

Case Western Reserve University
 FNS:

functional neuromuscular stimulation
 GOF:

goodnessoffit
 SCI:

spinal cord injury.
Declarations
Acknowledgements
We thank Drs. Musa Audu, Gerald Saidel, and Vira Chankong for their guidance and expertise in optimization theory. This work was supported in part by the National Institutes of Health under Grant NIH 5R01EB001889, Grant GCRC M01 RR000080, and Grant T32EB0431408.
Authors’ Affiliations
References
 Yoshida K, Horch K: Selective stimulation of peripheral nerve fibers using dual intrafascicular electrodes. IEEE Trans Biomed Eng. 1993, 40: 492494. 10.1109/10.243412.View ArticlePubMed
 Fisher LE, Miller ME, Bailey SN, Davis JA, Anderson JS, Rhode L, Tyler DJ, Triolo RJ: Standing after spinal cord injury with fourcontact nervecuff electrodes for quadriceps stimulation. IEEE Trans Neural Syst Rehabil Eng. 2008, 16: 473478.PubMed CentralView ArticlePubMed
 McDonnall D, Clark GA, Normann RA: Selective motor unit recruitment via intrafascicular multielectrode stimulation. Can J Physiol Pharmacol. 2004, 82: 599609. 10.1139/y04047.View ArticlePubMed
 Fisher L, Tyler D, Anderson J, Triolo R: Chronic stability and selectivity of fourcontact spiral nervecuff electrodes in stimulating the human femoral nerve. J Neural Eng. 2009, 6: 04601010.1088/17412560/6/4/046010.PubMed CentralView ArticlePubMed
 Schiefer MA, Polasek KH, Triolo RJ, Pinault GC, Tyler DJ: Selective stimulation of the human femoral nerve with a flat interface nerve electrode. J Neural Eng. 2010, 7: 2600610.1088/17412560/7/2/026006.PubMed CentralView ArticlePubMed
 McDonnall D, Clark GA, Normann RA: Interleaved, multisite electrical stimulation of cat sciatic nerve produces fatigueresistant, ripplefree motor responses. IEEE Trans Neural Syst Rehabil Eng. 2004, 12: 208215. 10.1109/TNSRE.2004.828425.View ArticlePubMed
 Gustafson KJ, Pinault GC, Neville JJ, Syed I, Davis JA, JeanClaude J, Triolo RJ: Fascicular anatomy of human femoral nerve: implications for neural prostheses using nerve cuff electrodes. J Rehabil Res Dev. 2009, 46: 973984. 10.1682/JRRD.2008.08.0097.PubMed CentralView ArticlePubMed
 Sunderland S: The intraneural topography of the radial, median and ulnar nerves. Brain. 1945, 68: 243299. 10.1093/brain/68.4.243.View ArticlePubMed
 Mahnam A, Hashemi SM, Grill WM: Measurement of the currentdistance relationship using a novel refractory interaction technique. J Neural Eng. 2009, 6: 03600510.1088/17412560/6/3/036005.PubMed CentralView ArticlePubMed
 Schiefer M, Polasek K, Triolo R, Pinault G, Tyler D: 31st Annual Internation Conference of the IEEE EMBS. Intraoperatice Demonstration of Selective Stimulation of the Common Human Femoral Nerve with a FINE. 2009, Minneapolis, 610613.
 Branner A, Stein RB, Normann RA: Selective stimulation of cat sciatic nerve using an array of varyinglength microelectrodes. J Neurophysiol. 2001, 85: 15851594.PubMed
 Yoshida K, Horch K: Reduced fatigue in electrically stimulated muscle using dual channel intrafascicular electrodes with interleaved stimulation. Ann Biomed Eng. 1993, 21: 709714. 10.1007/BF02368649.View ArticlePubMed
 Rainoldi A, Durfee W: 7th Vienna International Workshop on Functional Electrical Stimulation. Twitch summation with double stimulation. 2001, Vienna
 Rutten WL, van Wier HJ, Put JH: Sensitivity and selectivity of intraneural stimulation using a silicon electrode array. IEEE Trans Biomed Eng. 1991, 38: 192198.View ArticlePubMed
 Dowden BR, Frankel MA, Normann RA, Clark GA: Noninvasive method for selection of electrodes and stimulus parameters for FES applications with intrafascicular arrays. J Neural Eng. 2011, 9: 016006View ArticlePubMed
 Grill WM, Mortimer JT: Stability of the input–output properties of chronically implanted multiple contact nerve cuff stimulating electrodes. IEEE Trans Rehabil Eng. 1998, 6: 364373. 10.1109/86.736150.View ArticlePubMed
 Grill WM, Mortimer JT: Quantification of recruitment properties of multiple contact cuff electrodes. IEEE Trans Rehabil Eng. 1996, 4: 4962. 10.1109/86.506402.View ArticlePubMed
 Naples GG, Mortimer JT, Scheiner A, Sweeney JD: A spiral nerve cuff electrode for peripheral nerve stimulation. IEEE Trans Biomed Eng. 1988, 35: 905916. 10.1109/10.8670.View ArticlePubMed
 Smith B, Tang Z, Johnson MW, Pourmehdi S, Gazdik MM, Buckett JR, Peckham PH: An externally powered, multichannel, implantable stimulatortelemeter for control of paralyzed muscle. IEEE Trans Biomed Eng. 1998, 45: 463475. 10.1109/10.664202.View ArticlePubMed
 Bhadra N, Kilgore KL, Peckham PH: Implanted stimulators for restoration of function in spinal cord injury. Med Eng Phys. 2001, 23: 1928. 10.1016/S13504533(01)000121.View ArticlePubMed
 Tyler DJ, Durand DM: Chronic response of the rat sciatic nerve to the flat interface nerve electrode. Ann Biomed Eng. 2003, 31: 633642.View ArticlePubMed
 Tyler DJ, Durand DM: Functionally selective peripheral nerve stimulation with a flat interface nerve electrode. IEEE Trans Neural Syst Rehabil Eng. 2002, 10: 294303. 10.1109/TNSRE.2002.806840.View ArticlePubMed
 Durfee WK, MacLean KE: Methods for estimating isometric recruitment curves of electrically stimulated muscle. IEEE Trans Biomed Eng. 1989, 36: 654667. 10.1109/10.32097.View ArticlePubMed
 Akaike H: A new look at the statistical model identification. IEEE Trans Automat Control. 1974, 19: 716723. 10.1109/TAC.1974.1100705.View Article
 Kagaya H, Sharma M, Kobetic R, Marsolais EB: Ankle, knee, and hip moments during standing with and without joint contractures: simulation study for functional electrical stimulation. Am J Phys Med Rehabil. 1998, 77: 4954. 10.1097/0000206019980100000009. quiz 65–46View ArticlePubMed
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.