Next Article in Journal
Genetic Pathways to Insomnia
Next Article in Special Issue
Impact of Spinal Manipulation on Cortical Drive to Upper and Lower Limb Muscles
Previous Article in Journal
Preconditioning of Spatial and Auditory Cues: Roles of the Hippocampus, Frontal Cortex, and Cue-Directed Attention
Previous Article in Special Issue
Body-Machine Interfaces after Spinal Cord Injury: Rehabilitation and Brain Plasticity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Teaching Signals for Motor Control in the Cerebellum during Real-World Robot Application

by
Ruben Dario Pinzon Morales
1,† and
Yutaka Hirata
1,2,*,†
1
Neural cybernetics laboratory, Department of Computer Science, Graduate School of Engineering, Chubu University, Kasugai 487-8501, Japan
2
Department Robotic Science and Technology, Graduate School of Engineering, Chubu University, Kasugai 487-8501, Japan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Brain Sci. 2016, 6(4), 62; https://doi.org/10.3390/brainsci6040062
Submission received: 6 July 2016 / Revised: 12 December 2016 / Accepted: 14 December 2016 / Published: 20 December 2016
(This article belongs to the Special Issue Motor Control and Brain Plasticity)

Abstract

:
Motor learning in the cerebellum is believed to entail plastic changes at synapses between parallel fibers and Purkinje cells, induced by the teaching signal conveyed in the climbing fiber (CF) input. Despite the abundant research on the cerebellum, the nature of this signal is still a matter of debate. Two types of movement error information have been proposed to be plausible teaching signals: sensory error (SE) and motor command error (ME); however, their plausibility has not been tested in the real world. Here, we conducted a comparison of different types of CF teaching signals in real-world engineering applications by using a realistic neuronal network model of the cerebellum. We employed a direct current motor (simple task) and a two-wheeled balancing robot (difficult task). We demonstrate that SE, ME or a linear combination of the two is sufficient to yield comparable performance in a simple task. When the task is more difficult, although SE slightly outperformed ME, these types of error information are all able to adequately control the robot. We categorize granular cells according to their inputs and the error signal revealing that different granule cells are preferably engaged for SE, ME or their combination. Thus, unlike previous theoretical and simulation studies that support either SE or ME, it is demonstrated for the first time in a real-world engineering application that both SE and ME are adequate as the CF teaching signal in a realistic computational cerebellar model, even when the control task is as difficult as stabilizing a two-wheeled balancing robot.

1. Introduction

Cerebellar research in the last few decades has been accompanied by general consensus about its role in motor control, motor learning and motor coordination, yet abundant questions regarding the fundamental operations handled by this neural structure remain elusive [1]. One of the points of avid debate is the mechanisms of plasticity in the cerebellum [2]. In the original work by Marr-Albus [3,4], it was suggested that the cerebellum behaves as an adaptive filter in which mossy fiber inputs to the cerebellum were analyzed into sparse components at the granular cell layer and then conveyed by parallel fibers to the Purkinje cell. In this model, the location of plasticity corresponds to synapses between the numerous parallel fibers and the Purkinje cell, which change their strength according to a teaching signal conveyed by the climbing fiber (CF) input. Although a large repertoire of cerebellar models has been proposed since the early days of Marr, parallel fiber synapses remained the most commonly suggested location for cerebellar plasticity [5,6,7,8,9,10,11,12]. The debate pivots around the nature of the teaching signal in the CF.
Two types of movement error have been recognized as candidates for teaching signals: one from the sensory origin encoding error in movement kinematics (sensory error (SE)), the other from premotor systems encoding error in the motor command (motor error (ME)). For example, in the case of a tracking eye movement called smooth pursuit that is initiated when we try to follow a smoothly-moving visual target with the eyes, the movement error is the deviation of the eye velocity from the target velocity. This movement error causes a slippage of the target image on the retina that is detected by the retinal neuronal network. Thus, the output of the retina (ganglion cell activities) encoding the image slip kinematics is SE. Once the image slip is detected, the smooth pursuit system works to reduce it by modifying motor commands sent to the extra-ocular muscles. This modified motor command is ME. There are physiological evidence and computational studies supporting each teaching signal. In support of the notion that CF carry SE, CF activities highly correlated with retinal slip during optokinetic visual stimuli have been reported in rabbits [13] and monkeys [14]. In contrast, support for CF carrying ME has been provided by single-unit recording experiments using rabbits [15] and monkeys [16].
Theories have been built around these two types of teaching signals. On the one hand, the feedback-error-learning algorithm proposed by Kawato and colleagues [17,18] posits that plasticity in the cerebellum is driven by ME, which is made available at CF by reference structures outside the cerebellum that transform SE into ME. To produce an adequate teaching signal, these reference structures must implement a partial inverse model of the plant under control, which is sometimes troublesome when the control plant includes several degrees of freedom. Feedback-error-learning seeks to solve this problem by using a simple approximation to the inverse model, such as a gain term to transform SE into ME, which has been shown to work properly in simple control tasks and robotic applications [6]. For biological plants, it has been argued that with multiple-input multiple-output redundant systems, such as the vestibulo-ocular reflex in three dimensions, the connectivity needed by the reference structure becomes unfeasibly complex. On the other hand, Porrill and colleagues proposed the recurrent architecture for an adaptive filter in which the learning is driven by SE [5,19]. The structure of this algorithm partially follows the connectivity of particular cerebellar microzones, such as the flocculus; however, the role of granular and molecular layer interneurons has been unaccounted for. This recurrent algorithm allows stable adaptive learning using only observable SE signals and has been shown to work properly in complex robotic and simulation scenarios.
To solve the apparent dichotomy about the teaching signal in the CF, computational models of the cerebellum widely used in engineering applications may provide a potential framework to evaluate different CF contents and evaluated them in appropriate movement control tasks. However, there has not been any attempt so far. Here, we conduct a real-world testing of our cerebellar neuronal network model (CNN) in robot control and compare the performance of the robot when the different teaching signals are used. We employ a revised version of our CNN that incorporates a realistic cerebellar network architecture and a learning algorithm whose validity has been proven in both simulation and real-world experiments [9,11]. For this CNN, we configure adequate control tasks in which SE and ME are clearly dissociated. We demonstrate that SE and ME yield equivalent performance when the control task is simple, while they yield comparable performance even when the control task is more complex. Possible mechanisms underlying these results in the cerebellar neural network are discussed.

2. Materials and Methods

2.1. Neuronal Network Model of the Cerebellum

Inspired by the neuronal circuit of the cerebellar cortex, in particular that of flocculus, known to be involved in the horizontal vestibulo-ocular reflex (Figure 1A), we developed a neural network model of the cerebellum [9,11] (Figure 1B). As the cerebellar cortex has a homogeneous structure all over its volume, its basic learning algorithm is sometimes called a generic algorithm of the cerebellum [20]. We focused on the flocculus and the vestibulo-ocular reflex because they have been widely studied in anatomy and physiology [20]. The model comprises cell types found so far for which their physiological and anatomical properties are well understood [20]. Those are granule cells, Golgi cells, basket/stellate cells and Purkinje cells. We do not differentiate basket and stellate cells in our model. The model receives two types of input carried by mossy fibers and CF as in the real cerebellum. Mossy fiber inputs are postulated to carry desired motion signals, the efference copy of motor commands and sensory error signals (desired trajectory minus actual trajectory) [21,22,23]. These inputs innervate granule cells and Golgi cells. Golgi cells receive excitatory input from granule cells while projecting back to granule cells via inhibitory synapses, forming a negative feedback loop to regulate the activity of granule cells. Axons of granule cells, called parallel fibers, bifurcate and innervate basket cells and Purkinje cells via excitatory synapses. Basket cells project to Purkinje cells via inhibitory synapses, forming a negative feedforward pathway from granule cells to Purkinje cells. Purkinje cells are the sole output cell in the cerebellar cortex sending inhibitory projections down to the vestibular nucleus. This model is available for download at [24] or by accessing the model database of the cerebellar platform of the Japanese Neuroinformatics Node [25].
A proportional and derivative (PD) controller that is a feedback controller widely used in industry and other applications [26] is employed in tandem with the CNN model (Figure 1B), representing the non-cerebellar pathway in the vestibulo-ocular reflex neuronal circuit in Figure 1A as in previous simulation studies [27,28]. It provides the non-adaptive motor command to the control plant. The PD controller also generates the ME signal in this context since it is a simple gain term that approximates the inverse model of the control plant transforming SE into ME as proposed in Kawato’s feedback error learning scheme [6]. The parameters of the PD controller (proportional gain k p and derivative gain k d ) were designed following optimal settings for automatic controllers [29], so that the PD controller alone can stably operate the control plant during a simple desired motion ( π sin ( 2 π 0.1 t ) ).
The structure of our CNN model and its connectivity with the control plant, the non-cerebellar pathway and the vestibular nucleus remain the same disregarding the content of the error signal in our experiments. However, to be inline with the adaptive filter architecture proposed by Porrill et al. [19], our model would require the modification of the cerebellar structure so that the output of the Purkinje cell is projected to the input of the non-cerebellar pathway. We do not modify our cerebellar structure, maintain the architecture of the model as reported in the flocculus and only change the CF content.

2.2. Equations of the Model

The CNN model includes 755 granule cells, 5 Golgi cells, 15 basket cells and 1 Purkinje cell. These numbers preserve the ratio of the cell population and the convergence/divergence ratios of each cell type close to the cerebellar neuronal circuit [28] (shown in Table 1) and a sufficient number of granule cells for adequate control of the plants considered here. For a detailed study on the effects of scaling the number of neurons in the CNN and the output performance, refer to [30]. Each neuron model is described as follows:
x GC = y MF ( M MF - GC · W MF - GC ) + y GO ( M GO - GC · W GO - GC )
y PF = 1 1 + e σ ( x GC μ )
x GO = y MF ( M MF - GO · W MF - GO ) + y PF ( M PF - GO · W PF - GO )
y GO = 1 1 + e σ ( x GO μ )
x BA = y PF ( M PF - BA · W PF - BA )
y BA = 1 1 + e σ ( x BA μ )
x PC = y PF W PF - PC + y BA W BA - PC
y PC = 1 1 + e σ ( x PC μ ) 0.5
where σ = 8 , μ = 1 / 2 , x GC R + 1 × N GC is the activity vector of all granule cells before being processed by the sigmoidal activation function y PF (Equation (2)), W MF - GC R + N MF × N GC , W GO - GC R N GO × N GC are the matrix of synaptic weights between mossy fibers-granule cells and Golgi-granule cells, M MF - GC N N MF × N GC , M GO - GC N N GO × N GC are masking matrices of zeros and ones that embed the convergence/divergence ratio for granule cells referring to Table 1. x GO R + 1 × N GO and x BA R + 1 × N BA , W MF - GO R + N MF × N GO , W PF - GO R + N GC × N GO , W PF - BA R + N GC × N BA , W PF - PC R + N GC × N PC , W BA - PC R N BA × N PC follow the same notation for Golgi cells, basket cells and the Purkinje cell, and N MF , N GC , N GO , N BA , N PC are the number of mossy fibers, Golgi, granule, basket and Purkinje cells, respectively; where R + = { x R : 0 < x 1 } and R = { x R : 1 < x 0 } for excitatory and inhibitory synapses, respectively. Firing rate vectors are all R + , but y C is R = { x : 0.5 < x 0.5 } , so that it can be subtracted from the motor command produced by the PD and sent to the plant.
In order to standardize the contribution of each mossy fiber input to the cerebellar model, pre-cerebellar gains were set as normalizing scaling constants for each of them. These values were taken from the maximum values in the description of the robot hardware. Table 2 shows the scaling gains for each control plant. Then, the scaled mossy fibers were passed through an activation function ( y MF = 1 / 1 + e σ ( x MF μ ) ) to obtain firing rates compatible with the model (i.e., R + ), where σ and μ are as defined previously.
The learning rule regulated by the activity in the CF input that modifies the weights of parallel fiber-Purkinje cell synapses in the model is computed every sampling time according to the following equations:
Δ W PF - PC = γ · ( y PF ( t ) 0.5 ) · C F ( t )
W PF - PC ( t + 1 ) = W PF - PC ( t ) + Δ W PF - PC
where Δ W PF - PC is the matrix of changes in the synaptic weights of parallel fiber-Purkinje cell synapses, W PF - PC ( t ) R + is the synaptic weight matrix in these synapses at sampling time t and γ is the learning rate. In this algorithm, when the signs of ( y PF ( t ) 0.5 ) and C F ( t ) are the same, the synaptic weights are decreased (long-term depression), whereas when their signs are different, the weights are increased (long-term potentiation). Note that the constant 0.5 operating y PF is required to convert the y PF firing rate from R + R . Learning rate γ was set to be 0.008 , which was chosen by considering a trade off between learning speed and convergence [9,11].

2.3. Control Plants

Two control plants, a brushed direct current motor and a two-wheeled balancing robot, are employed. The 2-W brushed motor (Figure 1C, RC-280SA, Mabuchi Motor, Chiba prefecture, Japan) generates a torque directly from the current supplied. The motor’s shaft is interfaced with an encoder circuit (ZMP Inc., Tokyo, Japan) for providing angular position information and a microcontroller board (e-nuvo CPU board, ZMP Inc., Tokyo, Japan) in charge of communication with the implementation computer via a serial interface. The mossy fiber inputs to the CNN model for this control are shown in Table 2. The PD controller for this plant is a position controller with k p = 0.8 and k d = 0.01 as the proportional and derivative constants, respectively. A virtual dynamical model simulation for this motor has been included in the repository of the CNN model as an example [24].
The two-wheeled balancing robot (Figure 1D, e-nuvo wheel, ZMP Inc., Tokyo, Japan) is an inverted pendulum system that is highly unstable and widely used in control engineering for testing control strategies [26]. It is equipped with a set of sensors, including a motor encoder and a gyroscope, which provide wheel angle (Figure 1D; ϕ ( t ) ) and body tilt angle (Figure 1D; θ ( t ) ), respectively. The robot is also equipped with a serial interface to allow communication with the computer (see its specification below) on which the model was implemented. The motion of the robot is driven by a single motor connected to the wheels of the robot that share the same shaft (single-input multiple-output system). The mossy fiber inputs for this control plant carry the signals described in Table 2. The PD controller in this control object is a parallel configuration of two controllers (body position controller: k p = 5 and k d = 0.5 ; and wheel position controller: k p * = 0.2 and k d * = 0.05 ) designed by following optimal settings for automatic controllers [26,29], so that the addition of both outputs (i.e., PD controller output) alone can stably operate the robot during a simple task ( ϕ d e s . ( t ) = π sin ( 2 π 0.1 t ) , where ϕ d e s . ( t ) is the desired wheel angular position).
The experiments were carried out using LabVIEW 2010 (National Instrument, Austin, TX, USA) running on a Windows computer (4 × 3.33 GHz Intel Core-i7 processor, memory: 16 GB) and communicating with the control plants via serial protocol at 57.6 kbps. The CPU board in both the motor and the two-wheeled balancing robot has a sampling time of t s = 10 ms that is the same time interval as used in the present configuration of the CNN model (for different setups and sampling times, see [31]).

2.4. Two Types of Error Information in CF

Sources of SE and ME, shown in Figure 1B, are generated from the non-cerebellar pathway. That is, the SE signal is computed from the difference between the desired and actual movement [20,32]. It carries position and velocity error components in kinematic coordinates (expressed in angle units: rad), whereas ME is computed as the output of the PD feedback controller, which represents error in motor command dynamic coordinates (expressed in electric current units: A) as proposed in the theory of the feedback-error-learning algorithm [6,17]. Additionally, the linear combination of the two types of signals (SE + ME) is also tested, as it has been suggested to be a plausible case in the real cerebellum [33]. The equations of SE and ME for each of the control plants are as follows:
c c SE DC motor = a 1 ϕ ( t ) e + a 2 ϕ ˙ e ( t )
SE robot = a 1 ϕ ( t ) e + a 2 ϕ ˙ e ( t ) + b 1 θ e ( t ) + b 2 θ ˙ e ( t )
ME DC motor = c 1 P D ( t )
where subscript “e” denotes error and corresponds to the difference between the desired and yielded motion (e.g., ϕ e ( t ) = ϕ d e s . ( t ) ϕ ( t ) ), constants a 1 = 0.5 rad 1 , a 2 = 0.02 rad/ s 1 , b 1 = 5 rad 1 , b 2 = 0.5 rad/ s 1 and c 1 = 0.4 A 1 are scaling values intended to equalize the contribution of each error component to the SE and ME signals so that all have the same importance. These values were calculated with the CNN model disabled (i.e., plant controlled only by the PD controller and desired motion ϕ d e s . ( t ) = π sin ( 2 π 0.1 t ) ). ME in Equation (13) remained the same for both control plants. The combination of SE and ME is computed directly from Equations (11)–(13).

2.5. Performance Assessment Methodology

The desired motion ϕ d e s . ( t ) for the control plants used to assess the performance of the cerebellar model is a sinusoidal motion. The sinusoidal wave was generated at frequencies ranging from 0.2 Hz– 0.4 Hz and an amplitude of π (maximum angular velocity: 7.89 rad/s). The frequencies of the stimulus were chosen below the maximum controllable velocity of the robot (9.82 rad/s). The desired body tilt angle θ d e s . ( t ) was set to zero degrees (90 degrees with the ground), so that the robot is commanded to remain vertical while following the desired wheel trajectory. Each stimulus was repeated up to 100 cycles. Considering that the random initialization of synaptic weights in the cerebellar neuronal network can be a source of variability in the performance of the model, five different initial sets of random synaptic weights were created for each set of weighted connections ( W PF - PC , W MF - GC , W MF - GO , W GO - GC , W PF - GO , W PF - BA , and W BA - PC ) sampled from a normal distribution in R + and R with a standard deviation of 1 and a mean value of 0.5 and −0.5 for excitatory and inhibitory synapses, respectively. Special care is required to avoid over inhibition of the Purkinje cell when W BA - PC is large or instability in the Purkinje cell output when the feedback loop formed by W PF - GO and W GO - GC is large. Excluding W PF - PC , all of the synaptic weights remained fixed during the experiments because they are considered the main synaptic plasticity locations of cerebellar motor learning [20,34]. The performance of the cerebellar neuronal network controller was measured as the root mean square error (RSE) of each control variable. In the case of the motor, there was one control variable (the shaft angle ϕ ( t ) ; Figure 1C), whereas there were two control variables for the two-wheeled balancing robot (wheel angle ϕ ( t ) and body tilt angle θ ( t ) ; Figure 1D).

3. Results

3.1. Simple Control Scenario

First, we tested our CNN model with different CF information contents in the easiest control scenario in our setup, that is the control of the angular position ϕ ( t ) of a metal shaft directly connected to a motor with a sinusoidal desired motion ( ϕ d e s . ( t ) = π sin ( 2 π 0.5 t ) ) (the experimental protocol is described in Materials and Methods). Three types of CF to the cerebellar model were tested separately; namely, CF carrying SE (Equation (11)), ME (Equation (13)) and the combination of the two, i.e., SE + ME. Figure 2 shows the behavioral consequences in control performance under each CF condition (SE, ME, SE + ME) in terms of the RSE of ϕ ( t ) (Figure 2A) and yielded shaft motions (Figure 2B). The performance obtained when the CNN model was disabled, i.e., the motor is controlled only by the PD, is included as “PD” in gray lines (Figure 2A,B). The RSE of ϕ ( t ) (Figure 2A) evidences the superior performance with the CNN model disregarding the CF content compared with the PD. The improvement was on average 0.015 rad or 30% of the initial value of 0.05 rad. Thus, using a CF with SE, ME or the combination of both produces similar performance in this simple control scenario. Not surprisingly, the temporal profiles of these CFs also look alike (Appendix Figure A1). Differences between the motion caused in the motor by the CNN model and that by the PD controller alone can be seen in the XY planes shown in Figure 2B. These XY planes are constructed by positioning the desired (X-axis) and yielded (Y-axis) motions in an XY plane rotated −45 degrees. In the ideal case, the desired and yielded motions are equal with a trajectory that lies in an horizontal line at y = 0 (Figure 2B, dashed lines). Erroneous motions in the clockwise (CW) and counterclockwise (CCW) rotations of the motor shaft are mapped in the planes y > 0 and y < 0 , respectively. To show the adaptation in the CNN model, the motions produced at Cycles #2, #50 and #95 are shown. Cycle #2 shows that the CNN model caused the motor to rotate in excess in the CW direction. This excess produced by the initialization conditions and the number of granule cells in the CNN model was corrected by Cycle #5. We have shown that increasing the number of granule cells brings robustness against initial conditions in a bi-hemispheric CNN [35]. XY planes corresponding to Cycles #50 and #95 show that the PD caused the motor to rotate in excess when transitioning from CW to CCW and from CCW to CW. In contrast, the CNN model caused the opposite effect in the shaft motion by reducing the rotation speed just before the transitions. As exemplified here, CF inputs carrying SE, ME or their combination are adequate error signals to drive plasticity at parallel fibers-Purkinje cell synapses in the CNN model and make little difference in the control performance if the control scenario is very simple (a shaft with 1 deg of freedom driven by a motor to follow a 0.5-Hz sinusoidal).

3.2. Difficult Control Scenario

After testing the easiest control scenario in our setup, the CF signals seemed to have similar temporal profiles, and thus, the relationship between the error signal used and the control performance attained is yet to be clarified. In this section, we tested our setup in a more difficult control scenario to elucidate this relationship.
The control scenario employs the two-wheeled balancing robot with sinusoidal desired motions for the wheel angle. The frequencies of the sinusoidal motion ϕ d e s . ( t ) range from 0.2 Hz– 0.4 Hz (the experimental protocol is described in Materials and Methods). Figure 3, in the same format as Figure 2, summarizes the control performance achieved when the frequency of the desired motion is 0.2 Hz. Similarly to the previous case, controlling the two-wheeled balancing robot with the CNN outperforms the PD controller alone. The wheel positions generated at Cycles #2, #50 and #95 are shown in Figure 3B. These trajectories evidence that the CNN progressively improved the yielded wheel motion so that the robot motion approached the desired motion ( ϕ d e s . ( t ) ). In contrast, the wheel motion generated by the PD shows considerably larger hysteresis, meaning that the yielded motion differs from the desired motion (Figure 3B, gray lines). Subtle differences can be seen in the control performance obtained by using SE, ME or their combination. Increasing the frequency of the desired motion from 0.2–0.3 and 0.4 Hz strengthened the differences between SE and ME. Figure 4 shows the temporal profile of the different CFs, the firing rate of the Purkinje cell and the average control performance at the three frequencies of ϕ d e s . ( t ) . The temporal profiles of CF carrying SE (Figure 4A, blue lines), ME (red lines) and SE + ME (green lines) show that at 0.2 Hz, all of the CFs look alike. However, increasing the frequency of ϕ d e s . ( t ) reveals that ME is delayed with respect to the other types of CFs (Figure 4A, black arrows).
The firing rates of the Purkinje cell and the vestibular nucleus (Figure 4B) also evidence clear differences when using SE or ME in the CF. At 0.4 Hz, the firing rate of the Purkinje cell caused by ME and SE is out of phase and presents distinctive shapes. The outputs of the vestibular nucleus due to these contributions by the Purkinje cell also differ depending on the CF types (Figure 4B, black arrows). Figure 4C summarizes the control performance in terms of the average RSE of ϕ ( t ) at each one of the frequencies tested at the beginning (Cycle #5) and the end (Cycle #95) of the experiments. This figure shows that the performance with ME deteriorates (Figure 4C, red bars) as the frequency increases, to such an extent that at 0.3 and 0.4 Hz, there is no improvement over the initial error (Figure 4C, asterisks). On the contrary, the performance with SE always reduces the initial error and is better than using the PD or the combination SE + ME. Thus, SE is the best error signal in our setup during the control of a highly unstable two-wheeled balancing robot at high frequencies (>0.3 Hz).

3.3. Neural Consequences of the Different CF Types

In this section, we analyze the neural consequences of using SE, ME or SE + ME in the CNN model during the control scenario with the two-wheeled balancing robot. We show the adaptation produced at parallel fibers-Purkinje cell synapses, which are the only plastic locus in our model, and the correlation of mossy fibers, granule cells and the CF used. The correlation values are shown to investigate the information in the granule cells that, in combination with the CF input, resulted in the learned parallel fibers-Purkinje cell synaptic weights.
Figure 5 shows the 755 parallel fibers-Purkinje cell synaptic weights ( W PF - PC ) during one of the experiments with CF carrying SE. W PF - PC with the other CF contents showed similar trends. The gray area in the figure shows the weights decreased by long-term depression from the initial value (black line). Those weights increased by long-term potentiation are shown above the black line. Some representative W PF - PC are presented in red lines. The number of potentiated and depressed weights varied depending on the types of CF. With CF carrying SE, ME and SE + ME, 65%, 71% and 68% of the W PF - PC weights were depressed, respectively. Across all of the experiments, the same trend was observed, namely the number of depressed synapses was larger by using ME than SE. The insets in Figure 5 show the normalized firing rate of the granule cells whose synaptic weights were the most potentiated (labeled as “a”) and most depressed (labeled as “b”) for each type of CF. The firing rate of potentiated granule cells are negatively correlated with the CF, whereas those depressed are positively correlated with the CF. Figure 6A shows the coefficient of correlation of the CF and the granule cell’s activity. The moving average (N = 13) is shown in bold lines. Granule cells have been sorted in the x-axis from the most depressed (x = 0) to the most potentiated (x = 755). The normalized synaptic weight of each sorted granule cell is shown (dashed line). Shadowed area shows those cells whose activity was depressed. The coefficient of correlation when using SE, ME or SE + ME shows that those granule cells potentiated are negatively correlated with the CF, meaning that increasing their activity reduces the error signal, whereas those GCs depressed are positively correlated with the error signal.
The correlation of the granule activities with the mossy fiber inputs reveals intrinsic differences when using SE (blue lines) and ME (red lines). Figure 6B shows the coefficient of correlation of the mossy fiber carrying wheel position error ( ϕ e ( t ) ) and the activities of granule cells. In this correlation, those granule cells positively correlated were preferably potentiated when using SE, contrary to when using ME. Figure 6C shows the coefficient of correlation of the mossy fibers carrying wheel desired position ( ϕ d e s . ( t ) ) and the activities of granule cells. In this correlation, those granule cells negatively correlated were preferably potentiated when using SE, contrary to when using ME. These two correlations evidence intrinsic differences when using SE and ME in our setup. Granule cells with firing activity in-phase with error components (i.e., ϕ e ( t ) and ϕ ˙ e ( t ) ) and those out of phase with desired motions (i.e., ϕ d e s . ( t ) and ϕ ˙ d e s . ( t ) ) were preferably engaged to produce the Purkinje cell activity when using SE. When using ME, the opposite relationship was more prominent among the granule cells potentiated. Other mossy fibers, such as the efference copy, did not show differences between SE and ME. Employing the combination of SE and ME seems to soften the correlations observed with SE and ME to an intermediate point (Figure 6, green lines); therefore, it is not a surprise that using SE + ME as CF produces intermediate control performance (Figure 4, green lines).

4. Discussion

Debate about the error information encoded in the CF input to the cerebellum has been going on for several years [2,5,16,36]. Two types of error have been greatly defended: SE and ME. Behavioral and neurophysiological support has been presented for SE in rabbits [13] and monkeys [14] and similarly for ME in rabbits [15] and monkeys [16]. Consensus has been difficult to reach because these behavioral and neurophysiological experiments employ simple movement tasks where the temporal patterns of SE and ME are alike. To the best knowledge of the authors, this study is the first direct comparison of the types of error information encoded in the CF input to a CNN and with the motor performance attained during a real-world engineering application. Our setup allowed us to configure the difficulty of the control task, change the type of error encoded by the CF and evaluated the yielded control performance, thus enabling us to effectively disassociate SE and ME.

4.1. Type of Error Signal and the Yielded Robot Control Performance

Our experimental results showed that both SE and ME, despite producing unique behavioral and neural changes in the CNN model, especially at the Purkinje cell output (Figure 6B), are adequate error signals to govern successfully the plants in a diversity of control tasks (simple and complex task) (Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6). The practical implications of these results are that under certain experimental setups, investigating the consequences in motor learning when using motor or error information in the CF might yield positive, but inconclusive results. Thus, the proper selection of control plants and stimuli that allow the dissociation of motor and control error contents is imperative. This is the first direct demonstration suggesting that both SE and ME can be an adequate error signal to teach a realistic cerebellar computational model in real-world engineering applications. These results are in agreement with biological evidence from the horizontal vestibulo-ocular reflex (VOR) and optokinetic reflex (OKR) systems, a feed-forward and a feedback system sharing the same controller (i.e., cerebellar microcomplex) and the same control object (i.e., eye plant) with CF carrying SE and ME, respectively [2]. However, further evaluation using the two-wheeled balancing robot with sinusoidal desired motion at different frequencies revealed quantitative differences in control performance caused by SE and ME (Figure 6). Namely, SE yielded better control performance than ME (Figure 6C). This might be trivial because we evaluated the control performance in terms of the RSE of SE (i.e., RSE of ϕ ( t ) ), and the CNN model was trained to reduce SE. Since the ultimate goal of a control task is to reduce the error between the desired motion and the yielded motion (i.e., reduce SE) in biological systems and engineering applications, we employed the RSE of SE as a measure of the goodness of the control performance. In addition, SE signals are observable. In the case of ME, the desired motor commands are not observable, and therefore, the performance cannot be computed directly in terms of ME.
The correlation of the firing rate activity of granule cells, CF and mossy fibers also revealed intrinsic differences caused by using SE, ME or their combination (Figure 6). In particular, granule cell activity preferably potentiated by LTP at parallel fibers-Purkinje cell synapses with SE and ME differed in the correlation with the mossy fiber inputs. This means that to produce the Purkinje cell activities shown in Figure 4 with SE and ME, granule cells with different mossy fiber inputs were engaged, thus the differences in control performance. In particular, granule cells carrying error information were preferentially potentiated by sensory error, and granule cells carrying desired motions were preferentially potentiated by motor error. In our setup, the relation between granule cells activity and the CF content holds true for the two control plants used; however, extrapolation to other systems is not evident since changing the number of mossy fiber inputs and the size of the population of granule cells would result in a different mapping of inputs to parallel fibers, which in turn also affects the information that reaches the Purkinje cell. Clarifying this relation requires further study of the mossy fiber inputs to each granule cell and their correlation with the CF error content. For a study in this direction with a bi-hemispherical version of the CNN, refer to [31].
We have employed a PD controller as the non-cerebellar pathway that produces ME not only because it is one of the most flexible, effective and popular feedback controllers [26,29], but because it has been proposed as a simple approximation to the inverse model of the controlled plant in the feedback-error-learning theory [18,37,38]. The PD may be actually physiologically plausible [19] since its computation requires sensory information and its first derivative (e.g., eye velocity and eye acceleration), which are available at brain stem [39] and the cerebellum [40]. However, as the complexity of the plant increases, this approximation is deficient, as we have shown in the experiments with the two-wheel balancing robot (Figure 6). We observed that increasing the frequency of the desired motion (i.e., increasing the complexity of the control task) produced a delayed ME signal, and performance weakened (Figure 6A). This is produced by the frequency characteristics of the PD controller and the fact that it is not a close approximation to the inverse model of the two-wheeled balancing robot. These results confirm a major drawback of using ME as the teaching signal in a realistic cerebellar model as proposed in the feedback-error-learning theory since the complexity of the reference structures transforming SE into ME must increase with the complexity of the control plant, which could become troublesome in highly elaborated control scenarios.
Our results also showed that a linear combination of SE and ME encoded in the CF is an adequate error signal in our CNN model to control the plants successfully. The performance attained was always better than when the CNN model was disabled (Figure 3, Figure 4 and Figure 5, magenta lines). These results are in-line with physiological evidence in monkeys during arm reaching experiments that has suggested that CF encodes both SE and ME at different times of the arm movement [33], and thus, CF carries the error signal required for real-time control and learning of movements.

4.2. Generalization and Limitation of the Current Results

The control performance with the CNN model was stable and consistent while considering a basic control plant (DC motor) and one of the most challenging plants in control engineering (the two-wheeled balancing robot) [26], combined with simple and complex stimulus (low and high frequency sinusoidal desired motion) within the operation boundaries of our setup (real-time operation at a sampling frequency of 10 ms). The characteristics of the control scenarios and our setup made it difficult to reproduce the delay of 100 ms (10× the sampling frequency) shown in the visual inputs in the real cerebellum [34], which correspond to ϕ ( t ) and θ ( t ) in our model. However, under different circumstances, we believe our CNN model would be able to learn and produce adequate motor commands. These control plants have one and two degrees of freedom, respectively. Therefore, we may not be able to simply extend the current results to those with higher numbers of degrees of freedom. There is a series of evidence in cerebellar control of eye movement showing that different zones in cerebellar flocculus control different single axes of rotation of the eye (i.e., yaw and pitch rotations) [41]. Thus, our current setup to control objects is inline with this view. To apply the model to controlling other objects with higher numbers of degrees of freedom, we may employ the same number of cerebellar models for controlling each axis of motion [35]. Evaluating the nature of error signals suitable for this situation is in our scope of future studies.
In the current configuration of our CNN model, we have only considered plasticity at parallel fibers-Purkinje cell synapses, which are regarded as the loci of fast adaptation in the cerebellum [42], as a minimum requirement for the motor learning algorithm [2,20]. Multiple sites of plasticity in the cerebellum have been reported, and their involvement in motor learning has been argued [10,43,44]. Therefore, evaluating the effects in motor performance during robot control when the CNN model includes multiple sites of plasticity remains a future study of our CNN model.

4.3. Comparison with Other Cerebellar Models

Computational models of the cerebellum and their successful interpolation into engineering applications have extensively been reported. In general, these models include a form of synaptic plasticity at parallel fibers-Purkinje cell synapses driven by the CF input as presented in our CNN model, and that corresponds to the basis for cerebellar motor learning. Examples of CF encoding SE include models used for simulation of the eye plant [45], control of pneumatic muscles [46] and control of robotic arms [10,47]. Similarly, examples of cerebellar models with ME in their CF input include robotic arms [18] and inverted pendulum systems [11,38]. We have demonstrated that the CF input in our CNN model can be configured to include SE or ME information for driving plasticity at parallel fibers-Purkinje cell synapses and control the robot plant successfully. To the best knowledge of the authors, this study is the first direct comparison of the type of teaching signal encoded in the CF input to a cerebellar model and with the motor performance attained during a real-world engineering application. In contrast to spiking models of the cerebellum [7,12,48,49], due to the level of abstraction in our CNN model (i.e., firing rate neuron models), spike patterns and temporal or spatial effects were impossible to evaluate. This would require the construction of a cerebellar network with spiking neuronal models that could endanger the real-time real-world application in control engineering and a network with realistic physical properties [28]. Nonetheless, large-scale cerebellum-like spiking model has also been demonstrated to run in real time [12], and its application into robotics has been suggested. We have also extended our CNN model to a bi-hemispherical configuration (biCNN) that presents additional benefits and reproduces a form of cerebellar asymmetrical motor adaptation using SE information in the CF input [30].

5. Conclusions

The current study demonstrated for the first time in a real-world engineering application that both sensory error (SE) and motor error (ME) are adequate as the climbing fiber teaching signal in a realistic computational cerebellar model not only when the control task is simple, but also when the control task is as difficult as stabilizing a two-wheeled balancing robot.

Acknowledgments

This study is supported by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan grant-in-aid (YH). RDPM is supported by the MEXT scholarship for postgraduate studies 2011.

Author Contributions

R.D.P.M. and Y.H. conceived of and designed the experiments. R.D.P.M. performed the experiments. R.D.P.M. and Y.H. analyzed the data. R.D.P.M. and Y.H. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Climbing fiber activity during the control of the DC motor. Five cycles superimposed of the different CF activities during the control of the DC motor with the CNN model.
Figure A1. Climbing fiber activity during the control of the DC motor. Five cycles superimposed of the different CF activities during the control of the DC motor with the CNN model.
Brainsci 06 00062 g007

References

  1. Manto, M.; Taib, N.O.B. The contributions of the cerebellum in sensorimotor control: What are the prevailing opinions which will guide forthcoming studies? Cerebellum 2013, 12, 313–315. [Google Scholar] [CrossRef] [PubMed]
  2. Ito, M. Error detection and representation in the olivo-cerebellar system. Front. Neural Circuits 2013, 7, 1. [Google Scholar] [CrossRef] [PubMed]
  3. Marr, D. A theory of cerebellar cortex. J. Physiol. 1969, 202, 437–470. [Google Scholar] [CrossRef] [PubMed]
  4. Albus, J.S. A theory of cerebellar function. Math. Biosci. 1971, 10, 25–61. [Google Scholar] [CrossRef]
  5. Dean, P.; Porrill, J.; Ekerot, C.F.; Jorntell, H. The cerebellar microcircuit as an adaptive filter: Experimental and computational evidence. Nat. Rev. Neurosci. 2010, 11, 30–43. [Google Scholar] [CrossRef] [PubMed]
  6. Kawato, M. Internal models for motor control and trajectory planning. Curr. Opin. Neurobiol. 1999, 9, 718–727. [Google Scholar] [CrossRef]
  7. Hofstotter, C.; Mintz, M.; Verschure, P.F.M.J. The cerebellum in action: A simulation and robotics study. Eur. J. Neurosci. 2002, 16, 1361–1376. [Google Scholar] [CrossRef] [PubMed]
  8. Carrillo, R.R.; Ros, E.; Boucheny, C.; Coenen, O.J.M. A real-time spiking cerebellum model for learning robot control. Biosystems 2008, 94, 18–27. [Google Scholar] [CrossRef] [PubMed]
  9. Tanaka, Y.; Ohata, Y.; Kawamoto, T.; Hirata, Y. Adaptive control of 2-wheeled balancing robot by cerebellar neuronal network model. In Proceedings of the 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Buenos Aires, Argentina, 31 August–4 September 2010; pp. 1589–1592.
  10. Garrido Alcazar, J.A.; Luque, N.R.; DAngelo, E.; Ros, E. Distributed cerebellar plasticity implements adaptable gain control in a manipulation task: A closed-loop robotic simulation. Front. Neural Circuits 2013, 7, 159. [Google Scholar]
  11. Pinzon-Morales, R.D.; Hirata, Y. Cerebellar-inspired bi-hemispheric neural network for adaptive control of an unstable robot. In Proceedings of the Biosignals and Biorobotics Conference (2013 BRC), Rio de Janeiro, Brazil, 18–20 February 2013; pp. 1–4.
  12. Yamazaki, T.; Igarashi, J. Realtime cerebellum: A large-scale spiking network model of the cerebellum that runs in realtime using a graphics processing unit. Neural Netw. 2013, 47, 103–111. [Google Scholar] [CrossRef] [PubMed]
  13. Graf, W.; Simpson, J.I.; Leonard, C.S. Spatial organization of visual messages of the rabbit’s cerebellar flocculus. II. Complex and simple spike responses of Purkinje cells. J. Neurophysiol. 1988, 60, 2091–2121. [Google Scholar] [PubMed]
  14. Stone, L.S.; Lisberger, S.G. Visual responses of Purkinje cells in the cerebellar flocculus during smooth-pursuit eye movements in monkeys. II. Complex spikes. J. Neurophysiol. 1990, 63, 1262–1275. [Google Scholar] [PubMed]
  15. Frens, M.A.; Mathoera, A.L.; van der Steen, J. Floccular complex spike response to transparent retinal slip. Neuron 2001, 30, 795–801. [Google Scholar] [CrossRef]
  16. Kobayashi, Y.; Kawano, K.; Takemura, A.; Inoue, Y.; Kitama, T.; Gomi, H.; Kawato, M. Temporal firing patterns of purkinje cells in the cerebellar ventral paraflocculus during ocular following responses in monkeys II. Complex spikes. J. Neurophysiol. 1998, 80, 832–848. [Google Scholar] [PubMed]
  17. Kawato, M.; Gomi, H. The cerebellum and VOR/OKR learning models. Trends Neurosci. 1992, 15, 445–453. [Google Scholar] [CrossRef]
  18. Kawato, M.; Gomi, H. A computational model of four regions of the cerebellum based on feedback-error learning. Biol. Cybern. 1992, 68, 95–103. [Google Scholar] [CrossRef] [PubMed]
  19. Porrill, J.; Dean, P.; Anderson, S.R. Adaptive filters and internal models: Multilevel description of cerebellar function. Neural Netw. 2013, 47, 134–149. [Google Scholar] [CrossRef] [PubMed]
  20. Ito, M. The Cerebellum: Brain for an Implicit Self; FT Press Science: Upper Saddle River, NJ, USA, 2011. [Google Scholar]
  21. Blazquez, P.M.; Hirata, Y.; Heiney, S.A.; Green, A.M.; Highstein, S.M. Cerebellar signatures of vestibulo-ocular reflex motor learning. J. Neurosci. 2003, 23, 9742–9751. [Google Scholar] [PubMed]
  22. Hirata, Y.; Highstein, S.M. Acute adaptation of the vestibuloocular reflex: Signal processing by floccular and ventral parafloccular purkinje cells. J. Neurophysiol. 2001, 85, 2267–2288. [Google Scholar] [PubMed]
  23. Huang, C.C.; Sugino, K.; Shima, Y.; Guo, C.; Bai, S.; Mensh, B.D.; Nelson, S.B.; Hantman, A.W. Convergence of pontine and proprioceptive streams onto multimodal cerebellar granule cells. eLife 2013, 2, e00400. [Google Scholar] [CrossRef] [PubMed]
  24. Atlassian Bitbucket. Available online: https://bitbucket.org/rdpinzonm/the-bicnn-model/downloads (accessed on 18 December 2016).
  25. Cerebellar Platform. Available online: https://cerebellum.neuroinf.jp/ (accessed on 18 December 2016).
  26. Li, Z.; Yang, C.; Fan, L. Advanced Control of Wheeled Inverted Pendulum Systems; Springer: London, UK, 2013. [Google Scholar]
  27. Ito, M. Cerebellar learning in the vestibulo-ocular reflex. Trends Cogn. Sci. 1998, 2, 313–321. [Google Scholar] [CrossRef]
  28. Solinas, S.; Nieus, T.; DAngelo, E. A realistic large-scale model of the cerebellum granular layer predicts circuit spatio-temporal filtering properties. Front. Cell. Neurosci. 2010, 4, 12. [Google Scholar] [CrossRef] [PubMed]
  29. Ziegler, J.G.; Nichols, N.B. Optimum settings for automatic controllers. Trans. ASME 1942, 64, 759–768. [Google Scholar] [CrossRef]
  30. Pinzon-Morales, R.; Hirata, Y. A bi-hemispheric neuronal network model of the cerebellum with spontaneous climbing fiber firing produces asymmetrical motor learning during robot control. Front. Neural Circuits 2014, 8, 131. [Google Scholar] [CrossRef] [PubMed]
  31. Pinzon-Morales, R.; Hirata, Y. An stand-alone and portable bi-hemispherical neuronal network model of the cerebellum for engineering applications. In Proceedings of the 2014 IEEE International Conference on Robotics and Biomimetics (ROBIO), Bali, Indonesia, 5–10 December 2014; pp. 1–4.
  32. Hirata, Y.; Lockard, J.M.; Highstein, S.M. Capacity of vertical VOR adaptation in squirrel monkey. J. Neurophysiol. 2002, 88, 3194–3207. [Google Scholar] [CrossRef] [PubMed]
  33. Kitazawa, S.; Kimura, T.; Yin, P.B. Cerebellar complex spikes encode both destinations and errors in arm movements. Nature 1998, 392, 494–497. [Google Scholar] [CrossRef] [PubMed]
  34. Raymond, J.L.; Lisberger, S.G. Neural learning rules for the vestibulo-ocular reflex. J. Neurosci. 1998, 18, 9112–9129. [Google Scholar] [PubMed]
  35. Pinzon-Morales, R.D.; Hirata, Y. A realistic bi-hemispheric model of the cerebellum uncovers the purpose of the abundant granule cells during motor control. Front. Neural Circuits 2015, 9, 18. [Google Scholar] [CrossRef] [PubMed]
  36. Thach, W.T. On the specific role of the cerebellum in motor learning and cognition: Clues from PET activation and lesion studies in man. Behav. Brain Sci. 1996, 19, 411–433. [Google Scholar] [CrossRef]
  37. Wolpert, D.M.; Miall, R.; Kawato, M. Internal models in the cerebellum. Trends Cogn. Sci. 1998, 2, 338–347. [Google Scholar] [CrossRef]
  38. Ruan, X.; Chen, J. On-line NNAC for a balancing two-wheeled robot using feedback-error-learning on the neurophysiological mechanism. J. Comput. 2011, 6, 489–496. [Google Scholar] [CrossRef]
  39. Blazquez, P.M.; Hirata, Y.; Highstein, S.M. Chronic changes in inputs to dorsal Y neurons accompany VOR motor learning. J. Neurophysiol. 2006, 95, 1812–1825. [Google Scholar] [CrossRef] [PubMed]
  40. Takemura, A.; Inoue, Y.; Gomi, H.; Kawato, M.; Kawano, K. Change in neuronal firing patterns in the process of motor command generation for the ocular following response. J. Neurophysiol. 2001, 86, 1750–1763. [Google Scholar] [PubMed]
  41. Lisberger, S.G.; Pavelko, T.A.; Broussard, D.M. Neural basis for motor learning in the vestibuloocular reflex of primates. I. Changes in the responses of brain stem neurons. J. Neurophysiol. 1994, 72, 928–953. [Google Scholar] [PubMed]
  42. Kassardjian, C.D.; Tan, Y.F.; Chung, J.Y.J.; Heskin, R.; Peterson, M.J.; Broussard, D.M. The site of a motor memory shifts with consolidation. J. Neurosci. 2005, 25, 7979–7985. [Google Scholar] [CrossRef] [PubMed]
  43. Gao, Z.; Beugen, B.J.V.; De Zeeuw, C.I. Distributed synergistic plasticity and cerebellar learning. Nat. Rev. Neurosci. 2012, 13, 619–635. [Google Scholar] [CrossRef] [PubMed]
  44. McElvain, L.E.; Bagnall, M.W.; Sakatos, A.; du Lac, S. Bidirectional Plasticity gated by hyperpolarization controls the gain of postsynaptic firing responses at central vestibular nerve synapses. Neuron 2010, 68, 763–775. [Google Scholar] [CrossRef] [PubMed]
  45. Kettner, R.E.; Mahamud, S.; Leung, H.C.; Sitkoff, N.; Houk, J.C.; Peterson, B.W.; Barto, A.G. Prediction of complex two-dimensional trajectories by a cerebellar model of smooth pursuit eye movement. J. Neurophysiol. 1997, 77, 2115–2130. [Google Scholar] [PubMed]
  46. Lenz, A.; Anderson, S.; Pipe, A.; Melhuish, C.; Dean, P.; Porrill, J. Cerebellar-inspired adaptive control of a robot eye actuated by pneumatic artificial muscles. IEEE Trans. Syst. Man Cybern. B Cybern. 2009, 39, 1420–1433. [Google Scholar] [CrossRef] [PubMed]
  47. Eskiizmirliler, S.; Forestier, N.; Tondu, B.; Darlot, C. A model of the cerebellar pathways applied to the control of a single-joint robot arm actuated by McKibben artificial muscles. Biol. Cybern. 2002, 86, 379–394. [Google Scholar] [CrossRef] [PubMed]
  48. Inagaki, K.; Hirata, Y.; Blazquez, P.; Highstein, S. Computer simulation of vestibuloocular reflex motor learning using a realistic cerebellar cortical neuronal network model. In Neural Information Processing; Ishikawa, M., Doya, K., Miyamoto, H., Yamakawa, T., Eds.; Springer: Berlin/Heidelberg, Germany, 2008; Volume 4984, pp. 902–912. [Google Scholar]
  49. Luque, N.R.; Garrido, J.A.; Carrillo, R.R.; Tolu, S.; Ros, E. Adaptive cerebellar spiking model embedded in the control loop: Context switching and robustness against noise. Int. J. Neural Syst. 2011, 21, 385–401. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Neural circuit of the horizontal vestibulo-ocular reflex (VOR) and cerebellar neuronal network (CNN) model. (A) The neural circuit of the horizontal VOR. CF, climbing fiber; GO, Golgi cell; GC, granule cell; HC, horizontal canal; IO, inferior olive; LR, lateral rectus muscle of the eye; MF, mossy fiber; MN, motor neuron; BA, basket/stellate cell; PC, parallel fibers; VN, vestibular excitatory/inhibitory neurons. Blue and red lines indicate excitatory and inhibitory action, respectively. The dashed line shows the CF. (B) The structure of the CNN model with a proportional and derivative (PD) feedback controller that represents the non-cerebellar pathway depicted in (A). Labeled dashed lines show the origin of the sensory error (SE) (blue dashed line) and motor error (ME) (red dashed line) sources of CF, respectively. (C) A 2-W direct current motor with 1 degree of freedom (DOF). The control variable is shaft angular position ϕ ( t ) . (D) Two-wheeled balancing robot with 2 DOFs. The control variables are wheel angular position ϕ ( t ) and body tilt angle θ ( t ) .
Figure 1. Neural circuit of the horizontal vestibulo-ocular reflex (VOR) and cerebellar neuronal network (CNN) model. (A) The neural circuit of the horizontal VOR. CF, climbing fiber; GO, Golgi cell; GC, granule cell; HC, horizontal canal; IO, inferior olive; LR, lateral rectus muscle of the eye; MF, mossy fiber; MN, motor neuron; BA, basket/stellate cell; PC, parallel fibers; VN, vestibular excitatory/inhibitory neurons. Blue and red lines indicate excitatory and inhibitory action, respectively. The dashed line shows the CF. (B) The structure of the CNN model with a proportional and derivative (PD) feedback controller that represents the non-cerebellar pathway depicted in (A). Labeled dashed lines show the origin of the sensory error (SE) (blue dashed line) and motor error (ME) (red dashed line) sources of CF, respectively. (C) A 2-W direct current motor with 1 degree of freedom (DOF). The control variable is shaft angular position ϕ ( t ) . (D) Two-wheeled balancing robot with 2 DOFs. The control variables are wheel angular position ϕ ( t ) and body tilt angle θ ( t ) .
Brainsci 06 00062 g001
Figure 2. Performance during the control of the direct current motor with sinusoidal desired motion. (A) Control performance in terms of the root mean squared error (RSE) of ϕ ( t ) for the different types of climbing fiber (CF) averaged across the five experiments with different sets of initial synaptic weights. The mean value (thick line) and raw performance are shown. For reference, the performance with the proportional and derivative (PD) controller alone is shown (gray lines). (B) Yielded ( ϕ ( t ) , Y-axis) versus desired ( ϕ d e s . ( t ) , X-axis) shaft angle position (rotated −45 degrees) when the different CFs are employed. The trajectories at Cycles #2, #50 and #95 are shown. Black scale lines are 0.05 rad.
Figure 2. Performance during the control of the direct current motor with sinusoidal desired motion. (A) Control performance in terms of the root mean squared error (RSE) of ϕ ( t ) for the different types of climbing fiber (CF) averaged across the five experiments with different sets of initial synaptic weights. The mean value (thick line) and raw performance are shown. For reference, the performance with the proportional and derivative (PD) controller alone is shown (gray lines). (B) Yielded ( ϕ ( t ) , Y-axis) versus desired ( ϕ d e s . ( t ) , X-axis) shaft angle position (rotated −45 degrees) when the different CFs are employed. The trajectories at Cycles #2, #50 and #95 are shown. Black scale lines are 0.05 rad.
Brainsci 06 00062 g002
Figure 3. Performance during the control of the two-wheeled balancing robot with sinusoidal desired motion. (A) Control performance in terms of the root mean squared error (RSE) of ϕ ( t ) for the different types of climbing fiber (CF) averaged across the five experiments with different sets of initial synaptic weights. The mean value (thick line) and raw performance are shown. For reference, the performance with the proportional and derivative (PD) controller alone is shown (gray lines). (B) Yielded ( ϕ ( t ) , Y-axis) versus desired ( ϕ d e s . ( t ) , X-axis) wheel angle position (rotated −45 degrees) when the different CFs are employed. The trajectories at Cycles #2, #50 and #95 are shown. Notation as in Figure 2. Scaled black lines mark 0.05 rad.
Figure 3. Performance during the control of the two-wheeled balancing robot with sinusoidal desired motion. (A) Control performance in terms of the root mean squared error (RSE) of ϕ ( t ) for the different types of climbing fiber (CF) averaged across the five experiments with different sets of initial synaptic weights. The mean value (thick line) and raw performance are shown. For reference, the performance with the proportional and derivative (PD) controller alone is shown (gray lines). (B) Yielded ( ϕ ( t ) , Y-axis) versus desired ( ϕ d e s . ( t ) , X-axis) wheel angle position (rotated −45 degrees) when the different CFs are employed. The trajectories at Cycles #2, #50 and #95 are shown. Notation as in Figure 2. Scaled black lines mark 0.05 rad.
Brainsci 06 00062 g003
Figure 4. Firing rates of climbing fiber (CF) and the Purkinje cell (PC) during the control of the two-wheeled balancing robot. Three frequencies of ϕ d e s ( t ) are shown, 0.2, 0.3 and 0.4 Hz. (A) CF activity when carrying sensory error (SE) (blue), motor error (ME) (red) and SE + ME (green). Arrows show that CF carrying ME is delayed with respect to the other types of CF. (B) Firing rate of the Purkinje cell and vestibular nucleus (VN, black) at the different frequencies tested. Arrows show that the vestibular nucleus activity produced with CF carrying SE and ME is clearly different at 0.4 Hz. (C) Averaged control performance in terms of the root mean squared error (RSE) of ϕ ( t ) at Cycles #5 (light colors) and #95 (dark colors) at the different frequencies of the desired motion.
Figure 4. Firing rates of climbing fiber (CF) and the Purkinje cell (PC) during the control of the two-wheeled balancing robot. Three frequencies of ϕ d e s ( t ) are shown, 0.2, 0.3 and 0.4 Hz. (A) CF activity when carrying sensory error (SE) (blue), motor error (ME) (red) and SE + ME (green). Arrows show that CF carrying ME is delayed with respect to the other types of CF. (B) Firing rate of the Purkinje cell and vestibular nucleus (VN, black) at the different frequencies tested. Arrows show that the vestibular nucleus activity produced with CF carrying SE and ME is clearly different at 0.4 Hz. (C) Averaged control performance in terms of the root mean squared error (RSE) of ϕ ( t ) at Cycles #5 (light colors) and #95 (dark colors) at the different frequencies of the desired motion.
Brainsci 06 00062 g004
Figure 5. Synaptic weights between granule cells and the Purkinje cell ( W PF - PC ) during the control of the two-wheeled balancing robot. Red lines show some examples of W PF - PC . The horizontal line separating the shadowed areas of long-term depression (LTD) and long-term potentiation (LTP) is the initialization synaptic weight. Insets show the normalized firing rate (NFR) of granule cells whose W PF - PC synapses were predominately potentiated “a” or depressed “b” when climbing fiber (CF) carried sensory error (SE) (blue lines), motor error (ME) (red lines) or the combination of both (green lines).
Figure 5. Synaptic weights between granule cells and the Purkinje cell ( W PF - PC ) during the control of the two-wheeled balancing robot. Red lines show some examples of W PF - PC . The horizontal line separating the shadowed areas of long-term depression (LTD) and long-term potentiation (LTP) is the initialization synaptic weight. Insets show the normalized firing rate (NFR) of granule cells whose W PF - PC synapses were predominately potentiated “a” or depressed “b” when climbing fiber (CF) carried sensory error (SE) (blue lines), motor error (ME) (red lines) or the combination of both (green lines).
Brainsci 06 00062 g005
Figure 6. Correlation of granule cell (GCs) firing rate, climbing fiber (CF) and mossy fibers (MFs) in the CNN model. The granule cells in the x-axis have been sorted according to their W PF - PC value, so that at the origin is the most depressed and at x = 755 is the cell most potentiated during the control of the two-wheeled balancing robot. Shadowed areas show the granule cells depressed. The moving average (N = 13) is shown in bold colors. The synaptic weight is shown with a dashed lined. (A) Coefficient of correlation between the CF and granule cells; (B) coefficient of correlation between the mossy fiber carrying wheel position error ( ϕ e ( t ) ) and the granule cells; (C) coefficient of correlation between the mossy fiber carrying desired wheel position ( ϕ d e s ( t ) ) and the granule cell.
Figure 6. Correlation of granule cell (GCs) firing rate, climbing fiber (CF) and mossy fibers (MFs) in the CNN model. The granule cells in the x-axis have been sorted according to their W PF - PC value, so that at the origin is the most depressed and at x = 755 is the cell most potentiated during the control of the two-wheeled balancing robot. Shadowed areas show the granule cells depressed. The moving average (N = 13) is shown in bold colors. The synaptic weight is shown with a dashed lined. (A) Coefficient of correlation between the CF and granule cells; (B) coefficient of correlation between the mossy fiber carrying wheel position error ( ϕ e ( t ) ) and the granule cells; (C) coefficient of correlation between the mossy fiber carrying desired wheel position ( ϕ d e s ( t ) ) and the granule cell.
Brainsci 06 00062 g006
Table 1. Convergence and divergence synaptic ratio in the CNN (cerebellar neuronal network) model.
Table 1. Convergence and divergence synaptic ratio in the CNN (cerebellar neuronal network) model.
No. of CellsDivergenceConvergence
Mossy fibers (MF)7/5 *
Golgi (GO)5
Granule (GC)755
Basket/stellate (BA)15
Purkinje (PC)1
Parallel fibers (PF)755
MF → GC 1:1504:1
MF → GO 1:550:1
GC (parallel fibers)→ GO 1:4150:1
GO → GC 1:6003:1
GC (PF)→ BA 1:3050:1
GC (PF)→ PC 1:1755:1
BA → PC 1:115:1
* Number of mossy fiber inputs to the cerebellar model changes with the control plant; seven for the two-wheeled balancing robot, five for the motor in the present study.
Table 2. Characteristics of the control plants.
Table 2. Characteristics of the control plants.
ObjectOutputs (Sensors)Mossy FibersScaling Gain
Motor ϕ ( t ) : shaft ang. pos. (rad)1-des. shaft ang. pos.0.1 rad 1
ϕ ˙ ( t ) : shaft ang. vel. (rad/s)2- des. shaft ang. vel. 0.19 s/rad
3- shaft ang. pos. error 0.5 rad 1
4- shaft ang. vel. error 0.07 s/rad
5- efference copy1 A 1
Robot ϕ ( t ) wheel angle (rad)1- des. wheel ang. pos.0.03 rad 1
ϕ ˙ ( t ) wheel ang. vel. (rad/s)2- des. wheel ang. vel.0.04 s/rad
θ ( t ) body tilt ang. pos. (rad)3- body tilt ang. pos. error1 rad 1
θ ˙ ( t ) body tilt ang. vel. (rad/s)4- body tilt ang. vel. error0.5 s/rad
5- wheel ang. pos. error0.1 rad 1
6- wheel ang. vel. error0.2 s/rad
7- efference copy0.5 A 1
Angular (ang.); desired (des.); velocity (vel.); position (pos.).

Share and Cite

MDPI and ACS Style

Pinzon Morales, R.D.; Hirata, Y. Evaluation of Teaching Signals for Motor Control in the Cerebellum during Real-World Robot Application. Brain Sci. 2016, 6, 62. https://doi.org/10.3390/brainsci6040062

AMA Style

Pinzon Morales RD, Hirata Y. Evaluation of Teaching Signals for Motor Control in the Cerebellum during Real-World Robot Application. Brain Sciences. 2016; 6(4):62. https://doi.org/10.3390/brainsci6040062

Chicago/Turabian Style

Pinzon Morales, Ruben Dario, and Yutaka Hirata. 2016. "Evaluation of Teaching Signals for Motor Control in the Cerebellum during Real-World Robot Application" Brain Sciences 6, no. 4: 62. https://doi.org/10.3390/brainsci6040062

APA Style

Pinzon Morales, R. D., & Hirata, Y. (2016). Evaluation of Teaching Signals for Motor Control in the Cerebellum during Real-World Robot Application. Brain Sciences, 6(4), 62. https://doi.org/10.3390/brainsci6040062

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop