A model using mutual influence of firing rates of corticomotoneurons for learning a precision grip task
© Boussaton et al; licensee BioMed Central Ltd. 2011
Published: 18 July 2011
As a part of a Brain-Machine Interface, we are currently defining a model for learning and forecasting muscular activity, given sparse brain activity in the form of action potential signals (spike trains). We have been working on the flexion of a finger during a trained precision grip performed by a monkey (macaca nemestrina), as she clasps a metal gauge with her finger and thumb. Experimentally, the activity of about a hundred neurons in the motor cortex can be recorded simultaneously with the help of a multielectrode array, see Figure A and  for more details about retrieving and filtering the data. Our method is based on a system of equations involving the firing rate of each recorded neurons, a set of thresholds, and Euclidian distances between averaged and current state at each time step. The firing rates are computed according to given time-windows, between 25 and 100 milliseconds. The thresholds used depend on these firing rates. The learning is done on a subset of the experiments and then evaluated on what remains. A raw estimation is done in order to be used as a reference for estimating the efficiency of each part of the learning formula brings to the final result. The complete improvement formula is divided into three stages and can be expressed as follows: p(t+1) = p(t) + A(t) + B(t) + C(t) where p(t) is the force exerted on the gauge at time t. The A part is the learning reference base of the method in which a straight matching is done between each neuronal code and the derivative of the observed force in the finger at each timestep. What we call a neuronal code is the vector of all the values of the firing rate functions at any given time and during the training stage, to any neuronal code is associated the average value of all the recorded derivatives of the force. The B part is an actuation made on the distance between the current activity of each neuron and its average activity over a former time window of the same length as the time-window used to compute the spike trains. Finally, the C part is a system of equations in which we suppose that every neuron is correlated to each other in a weighted way we optimize during the learning process. The purpose of this study is multifold, we want to estimate (i) the influence of the neurons on each other qualitatively, (ii) the efficiency of various easy-to-tract improvement techniques and (iii) the importance of thresholding the firing rates.
We are thankful to Pr. Lemon and his team for providing us with these and letting us use them.
- Maier MA, Bennett KMB, Hepp-Raymond MC, Lemon RN: Contribution of the monkey corticomotoneuronal system to the control of force in precision grip. Journal of neurophysiology. 1993, 18 (3): 772-785.Google Scholar
- Velliste M, Perel S, Chance S, Spalding A, Schwartz AB: Cortical control of a prosthetic arm for self feeding. Nature. 2008, 453: 1098-1101. 10.1038/nature06996.View ArticlePubMedGoogle Scholar
- Taira M, Georgopoulos AP: Cortical Cell types from spike trains. Neuroscience research. 1993, 17: 37-45.Google Scholar
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.