Neural control on multiple time scales: Insights from human stick balancing

The time-delayed feedback control mechanisms of the nervous system are continuously subjected to the effects of uncontrolled random perturbations (herein referred to as noise). In this setting the statistical properties of the fluctuations in the controlled variable(s) can provide non-invasive ins...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2006
Hauptverfasser: Cabrera, J.L., Luciani, C., Milton, J.
Format: Artikel
Sprache:English
Veröffentlicht: Інститут фізики конденсованих систем НАН України 2006
Schriftenreihe:Condensed Matter Physics
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/121324
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Neural control on multiple time scales: Insights from human stick balancing / J.L. Cabrera, C. Luciani, J. Milton // Condensed Matter Physics. — 2006. — Т. 9, № 2(46). — С. 373–383. — Бібліогр.: 42 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-121324
record_format dspace
spelling irk-123456789-1213242017-06-15T03:03:10Z Neural control on multiple time scales: Insights from human stick balancing Cabrera, J.L. Luciani, C. Milton, J. The time-delayed feedback control mechanisms of the nervous system are continuously subjected to the effects of uncontrolled random perturbations (herein referred to as noise). In this setting the statistical properties of the fluctuations in the controlled variable(s) can provide non-invasive insights into the nature of the underlying control mechanisms. We illustrate this concept through a study of stick balancing at the fingertip using high speed motion capture techniques. Experimental observations together with numerical studies of a stochastic delay differential equation demonstrate that on time scales short compared to the neural time delay (“fast control”), parametric noise provides a non-predictive mechanism that transiently stabilizes the upright position of the balanced stick. Moreover, numerical simulations of a delayed random walker with a repulsive origin indicate that even an unstable fixed point can be transiently stabilized by the interplay between noise and time delay. In contrast, on time scales comparable to the neural time delay (“slow control”), feedback and feedforward control mechanisms become more important. The relative contribution of the fast and slow control mechanisms to stick balancing is dynamic and, for example, depends on the context in which stick balancing is performed and the expertise of the balancer. 2006 Article Neural control on multiple time scales: Insights from human stick balancing / J.L. Cabrera, C. Luciani, J. Milton // Condensed Matter Physics. — 2006. — Т. 9, № 2(46). — С. 373–383. — Бібліогр.: 42 назв. — англ. 1607-324X PACS: 89.75.-k, 87.19.St, 02.30.Ks, 05.45.-a, 02.50.-r DOI:10.5488/CMP.9.2.373 http://dspace.nbuv.gov.ua/handle/123456789/121324 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description The time-delayed feedback control mechanisms of the nervous system are continuously subjected to the effects of uncontrolled random perturbations (herein referred to as noise). In this setting the statistical properties of the fluctuations in the controlled variable(s) can provide non-invasive insights into the nature of the underlying control mechanisms. We illustrate this concept through a study of stick balancing at the fingertip using high speed motion capture techniques. Experimental observations together with numerical studies of a stochastic delay differential equation demonstrate that on time scales short compared to the neural time delay (“fast control”), parametric noise provides a non-predictive mechanism that transiently stabilizes the upright position of the balanced stick. Moreover, numerical simulations of a delayed random walker with a repulsive origin indicate that even an unstable fixed point can be transiently stabilized by the interplay between noise and time delay. In contrast, on time scales comparable to the neural time delay (“slow control”), feedback and feedforward control mechanisms become more important. The relative contribution of the fast and slow control mechanisms to stick balancing is dynamic and, for example, depends on the context in which stick balancing is performed and the expertise of the balancer.
format Article
author Cabrera, J.L.
Luciani, C.
Milton, J.
spellingShingle Cabrera, J.L.
Luciani, C.
Milton, J.
Neural control on multiple time scales: Insights from human stick balancing
Condensed Matter Physics
author_facet Cabrera, J.L.
Luciani, C.
Milton, J.
author_sort Cabrera, J.L.
title Neural control on multiple time scales: Insights from human stick balancing
title_short Neural control on multiple time scales: Insights from human stick balancing
title_full Neural control on multiple time scales: Insights from human stick balancing
title_fullStr Neural control on multiple time scales: Insights from human stick balancing
title_full_unstemmed Neural control on multiple time scales: Insights from human stick balancing
title_sort neural control on multiple time scales: insights from human stick balancing
publisher Інститут фізики конденсованих систем НАН України
publishDate 2006
url http://dspace.nbuv.gov.ua/handle/123456789/121324
citation_txt Neural control on multiple time scales: Insights from human stick balancing / J.L. Cabrera, C. Luciani, J. Milton // Condensed Matter Physics. — 2006. — Т. 9, № 2(46). — С. 373–383. — Бібліогр.: 42 назв. — англ.
series Condensed Matter Physics
work_keys_str_mv AT cabrerajl neuralcontrolonmultipletimescalesinsightsfromhumanstickbalancing
AT lucianic neuralcontrolonmultipletimescalesinsightsfromhumanstickbalancing
AT miltonj neuralcontrolonmultipletimescalesinsightsfromhumanstickbalancing
first_indexed 2025-07-08T19:39:22Z
last_indexed 2025-07-08T19:39:22Z
_version_ 1837108892577824768
fulltext Condensed Matter Physics 2006, Vol. 9, No 2(46), pp. 373–383 Neural control on multiple time scales: Insights from human stick balancing J.L.Cabrera1∗, C.Luciani2†, J.Milton3‡ 1 Centro de Fı́sica, Instituto Venezolano de Investigaciones Cientı́ficas, Apartado 21827, Caracas 1020A, Venezuela 2 Departamento de Fı́sica, Universidad Simón Bolı́var, Caracas 1080–A, Venezuela 3 Joint Science Department,The Claremont Colleges, Claremont, CA 91711, USA Received April 10, 2006, in final form May 12, 2006 The time-delayed feedback control mechanisms of the nervous system are continuously subjected to the effects of uncontrolled random perturbations (herein referred to as noise). In this setting the statistical prop- erties of the fluctuations in the controlled variable(s) can provide non-invasive insights into the nature of the underlying control mechanisms. We illustrate this concept through a study of stick balancing at the fingertip using high speed motion capture techniques. Experimental observations together with numerical studies of a stochastic delay differential equation demonstrate that on time scales short compared to the neural time delay (“fast control”), parametric noise provides a non-predictive mechanism that transiently stabilizes the upright position of the balanced stick. Moreover, numerical simulations of a delayed random walker with a repulsive origin indicate that even an unstable fixed point can be transiently stabilized by the interplay between noise and time delay. In contrast, on time scales comparable to the neural time delay (“slow control”), feedback and feedforward control mechanisms become more important. The relative contribution of the fast and slow control mechanisms to stick balancing is dynamic and, for example, depends on the context in which stick balancing is performed and the expertise of the balancer. Key words: fluctuation phenomena, time delay, inverted pendulum, feedback control, noise PACS: 89.75.-k, 87.19.St, 02.30.Ks, 05.45.-a, 02.50.-r 1. Introduction It is becoming clear that biological systems tuned at or close to a stability boundary can exhibit qualitatively the same kinds of complex dynamics [1,2] that are observed for physical systems at or close to a critical point [3], including power laws, intermittency, and critical slowing down. A case in point arises in the neural control of stick balancing at the fingertip. In the control literature, the studies of the inverted planar pendulum figure prominently in discussions of the stabilization of an unstable fixed point [4]. Since the equations of motion take the form [5] θ̈(t) + γ m θ̇(t)− g ` sin θ(t) + F (θ(t)) = 0, (1) where ` is the length of the stick, g is the gravitational constant, γ is the damping coefficient, and θ is the vertical displacement angle (θ = 0 corresponds to the upright position, hence the “–” sign), the control problem is how to design the feedback, F , to stabilize an unstable saddle point. However, in real applications, such as those that arise in the context of neural control, it is necessary to take into account the time required to detect a deviation in θ, and then effect a corrective movement. Consequently the feedback is time-delayed and (1) becomes θ̈(t) + γ m θ̇(t)− g ` sin θ(t) + F (θ(t− τ)) = 0, (2) ∗Email: jlc@ivic.ve & juluisca@gmail.com †Email: cluciani@gmail.com ‡Email: jmilton@jsd.claremont.edu c© J.L.Cabrera, C.Luciani, J.Milton 373 J.L.Cabrera, C.Luciani, J.Milton where τ is the time delay, and θ(t), θ(t − τ)) refer, respectively, to the values of θ at times t and t − τ . Since the upright position for (2) can be stable provided that τ and ` are appropriately chosen [6–8], the goal is to figure out the nature of the control strategy used by the nervous system from an analysis of the observed fluctuations. This neuro-biological inverse problem has attracted great interest [9,10]. It is widely held that an understanding of how the human nervous system addresses the problem of stick balancing at the fingertip may provide useful insights into solving the problems ranging from minimizing the risk of falling in the elderly to the developing of stable two-legged robots [11]. The observation that longer sticks are easier to balance at the fingertip than shorter ones is a reflection of the fact that the neural latency is finite (∼ 100 − 200 msec): as the stick becomes longer the speed of its movements eventually becomes slow enough so that the nervous system has time to make a corrective movement. However, it is found that no matter how long the stick or how expert the individual may be, the stick always eventually falls. Thus the intriguing question for stick balancing is not why the stick can be balanced at the fingertip, but to understand why the balanced stick eventually falls! Here we review our studies on stick balancing at the fingertip [2,6,12–16]. We demonstrate that it is possible to gain insight into the nature of the control strategy used by the nervous system by carefully studying the statistical properties of the fluctuations made in the controlled variable, i.e. θ, and the controller (i.e. the movements of the fingertip). 2. Survival functions The simplest hypothesis to account for the observation that a stick balanced at the fingertip falls is that the dimensions of the basin of attraction associated with the stable upright position are small compared to the intensity of the random perturbations, herein referred to as noise. In this situation, for example, a sufficiently large perturbation occurring at the right time could move the trajectory outside the basin of attraction and hence produce a fall. If the hypothesis were correct, then the survival function, i.e. the fraction of trajectories remaining in the basin of attraction at time t, might be anticipated to be exponential. Indeed, exponential survival curves have been associated with the times to cross a threshold in chaotic dynamical systems [17,18] and in stochastic dynamical systems subjected to the effects of additive noise [19], i.e. the effects of noise are independent of the state variable. An example of a stochastic dynamical system with additive noise can be obtained by adding a noise term to the right-hand side of (2). 10 0 10 1 10 2 10 3 t (sec) 10 -2 10 -1 10 0 P (t es c> t) (a) (b) Figure 1. (a): Survival function, P (tesc > t), measured for three different stick lengths, 39 cm (•), 62 cm (4) and 120 cm (¦). The P (tesc > t) generated by (3) (solid line) is compared to that measured experimentally for the 62 cm stick length. Parameters: m = 35 g, k = 35, γ = 40 sec−1, r0 = 0.97 , σ = 0.15. Data is from [12]. (b): Distribution of escape times for a 55 cm stick. The solid line is for the Weibul distribution with β = 1.8105. 374 Neural control on multiple time scales: Insights from human stick balancing Figure 1 (a) shows the survival curve, i.e. the probability that stick is balanced at the fingertip at time t, P (tesc > t), as a function of t, for two different stick lengths. A longer stick having the same mass is easier to balance and hence the survival curve is shifted to the right. As can be seen the survival curve is clearly not an exponential, but has the form of a Weibul function (figure 1(b)), P (tesc > t) ∼ exp(−λt)β , where β > 1, λ are constants. Weibul-type survival functions have been observed previously in the context of the analysis of the failure times of ball bearings [20] and in survival statistics of patients with certain forms of cancer [17]. One mechanism that can generate survival functions having this form occurs when the effects of noise are parametric, i.e. the effects of noise enter through a parameter that is multiplied by the state variable [12]. Vertical displacements in the pivot point of an inverted pendulum enter the equations of motion through a parameter [5]; hence parametric noise is particularly relevant to the study of stick balancing at the fingertip [21]. If we expand F as a Taylor series F (θ(t− τ)) ≈ r0θ(t− τ) + · · · and assume that the gain, r0 varies noisily, then (2) becomes θ̈(t) + γ m θ̇(t)− g ` sin θ(t) + r0(t)θ(t− τ) = 0, (3) where r0(t) = R0 + ξ(t) and R0 is a constant and ξ(t) is Gaussian distributed white noise with zero mean and variance σ. The survival curves in figure 1 can be reproduced by discretizing (3) to obtain the second-order stochastic discrete dynamical system with state-dependent noise (compare triangles with the solid line for the right survival curve in figure 1) xt+1 = xt + yt , yt+1 = ( 1− γ m ) yt − k m xt + rn m xt−τ , (4) where rt = r0 + ξt and where r0 is a constant, k/m = g/`, and ξt is a random variable chosen from a Gaussian distribution of discrete variables having mean equal to zero and variance, σ. Thus, the observations in figure 1 point to a dynamical system in which parametric noise plays an important role in shaping the observed dynamics. Further insight into the mechanisms of the control of stick balancing can be obtained by examining the fluctuations made in the controlled variables. 3. Stick and fingertip motions The fluctuations in the vertical displacement angle, θ, and the controller, i.e. the fingertip, can be measured with high accuracy using motion capture technology [6,12,13]. Infrared light reflected from markers attached to each end of the stick are detected by two specialized motion capture cameras (Qualisys, Inc): the image projected onto the CCD of each camera determines two of the spatial coordinates and the third is determined using triangulation techniques involving both cameras. Measurements can be easily made with a temporal resolution of 1 msec and a spatial resolution of 4µm. Figure 2 (a) shows the fluctuations in θ. It should be noted that θ is determined solely from the position of the two markers in the vertical (z) plane, i.e. cos θ = ∆z/`. The fluctuations in θ are characterized by intervals with small amplitude fluctuations that alternate with intervals of larger amplitude fluctuations. Figure 2 (b) shows the movements of the fingertip 375 J.L.Cabrera, C.Luciani, J.Milton (a) (b) Figure 2. (a) Temporal series for ∆z/l for a 62 cm stick balanced at the fingertip. The horizon- tal dashed line depicts the threshold position. (b) Two-dimensional representation of the path traveled by the fingertip during stick balancing. The digitization timestep was 1 ms. Movements in the anterior-posterior direction are shown on the y axis; side to side movements on the x axis. Data is taken from [13]. (i.e. the marker closest to the fingertip) in the (x, y)-plane during stick balancing. These movements are also characterized by two components: 1) relatively slow translations of the pivot point for the inverted pendulum; and 2) small amplitude corrective movements. Taken together the observations in figure 2 suggest that it is convenient to consider the mecha- nisms for neural control on two time scales: 1) fast control, i.e. the time between successive corrective movements occur on time scales short compared to the neural time delay; and 2) slow control, i.e. the times between successive corrective movements occur on time scales that are comparable to the neural time delay. 4. Fast components for balance control The times between successive crossings of a threshold (dashed horizontal line in figure 2 (a)) in the upward, or corrective, direction were measured. These times correspond to the laminar phases of a dynamical system that exhibits intermittency. Figure 3 (a) shows a double-log plot of the probability that the laminar phases have length δt, P (δt), as a function of δt. There are two important observations [6]. First, > 98% of the laminar phases are shorter than the measured neural latency. Second, the fluctuations exhibit two statistical features of on-off intermittency [22], namely, the laminar phases exhibit a –3/2 power law (solid line in figure 3 (a) and the power spectrum of the fluctuations in θ exhibits a –1/2 power law (figure 3 (b)). The power laws that characterize on-off intermittency can be reproduced using (3) provided that the parameters are chosen to place the deterministic version of this equation close enough to the stability boundary so that the r0(t) can be stochastically forced back and forth across it [6]. These observations suggest that parametric noise plays an important role in the fast component of the control of stick balancing. It should be noted that the zero crossings for a random walk also scale as –3/2 [23]. Thus the beneficial effects of on-off intermittency for stick balancing likely arise because the fluctuations in θ resemble a random walk for which the mean value of θ is approximately zero, i.e. the balance position is statistically stabilized [24]. In summary, the fast corrective movements for stick balancing exhibit three characteristics: 1) 376 Neural control on multiple time scales: Insights from human stick balancing they are intermittent and ballistic; 2) the time between successive corrective movements is much shorter than the neural latency; and 3) they are open-loop, i.e. non-predictive (see below). It has been recently suggested that ballistic fluctuations that scale as –3/2 are associated with optimal control [25]. Fast components having these same properties have also been observed in other human balancing tasks, namely postural sway during quiet standing [26,27] and ankle movement controlled balancing of an inverted pendulum [28]. Although it has been suggested that the mechano-elastic properties of limbs and joints could provide a “zero delay” feedback control mechanism for balance [27,28], our observations on stick balancing suggest that the effects of parametric noise may be a more plausible explanation for these fast controlling movements. In addition, these observations suggest that methods used to characterize intermittent behavior in turbulent flow [29], such as the flatness factor, may be useful for the study of skill acquisition in humans. (a) (b) Figure 3. (a) Log-log plot of the normalized probability of having laminar phases of length δt, P (δt) for a balanced stick of length 62 cm. The dotted line has slope −3/2 and the physiological delay is depicted by the vertical line. (b) Power spectrum, S(f), of the fluctuations in θ showing the presence of two power laws. The slopes of the solid lines are, respectively, −0.5 and −2.5. Data is from [6]. 5. Stabilizing with noise and delay An important question concerns whether the control system is tuned slightly to the stable or to the unstable side of the stability boundary, the proximity to the stability boundary depending on the intensity of the state-dependent noise. Two experimental observations suggested that the system may actually be tuned slightly to the unstable side of the boundary: 1) power laws identical to those observed for stick balancing are observed in a virtual stick balancing task that does not possess a stable point [14] and 2) the measured survival functions can be best fit by (4) when the parameters are chosen so that the deterministic version possesses only an unstable fixed point (see figure 1). The effect of partial or even full stabilization of unstable states by noise is well known to occur in stochastic dynamical parametric systems [30]. Furthermore, recently it has been shown that additive noise can induce transient stability by delaying the escape time of unstable systems in periodically modulated systems [31] and in pure unstable damped systems [32]. In the case of time delayed dynamical systems the stabilizing effect of additive and parametric noise by postponement of bifurcations has been also reported for specific systems [33]. Thus it seems appropriate to ask ourselves whether there exist cooperative effects between random perturbations and time delay 377 J.L.Cabrera, C.Luciani, J.Milton feedback at a more fundamental level responsible for improved stabilization in simple unstable systems. To address this question we have considered the case of the repulsive delayed random walk [15] which, in turn, is a special case of the delayed random walk (DRW) [34]. In the DRW the walker takes a discrete step of unit length per unit time in a direction determined by a set of conditional probabilities that depend on the position of the walker at a given time in the past, τ . The conditional probabilities for the walker to take a step to the right are P (X(t + 1) = X(t) + 1|X(t− τ) > 0) = p, P (X(t + 1) = X(t) + 1|X(t− τ) = 0) = 1 2 , P (X(t + 1) = X(t) + 1|X(t− τ) < 0) = 1− p, (5) where X(t) is the position of the walker at time t and 0 < p < 1. This walk differs from the well known brownian walk not only because its dynamics depend on the past but also because its origin can behave attractively or repulsively depending on the value of p. In particular the origin, X = 0, is attractive if p < 1/2, unbiased for p = 1/2 (as in the case of the brownian random walk) and repulsive if p > 1/2. Models of delayed random walks with an attractive origin enable a number of important properties of the delayed Langevin equation to be obtained analytically, including the autocorrelation function [34]. However, analytical insights in the repulsive case are not immediate and thus we restrict our exposition to a numerical approach. Figure 4. Mean first passage time as a function of the delay time for X∗ = 120 steps and different values of the bias p. (a) From top to bottom: p = 0.7, 0.8, 0.9, 1.0. (b) From top to bottom: p = 0.425, 0.45, 0.475, 1/2. Data obtained with 100000 realizations. Computer simulations of a DRW with a repulsive origin for a given τ require the definition of an initial function, Φ(s), where s ∈ [−τ, 0]. We generated this function by using a brownian random walk with p = 1/2 and τ = 0. The origin of the delayed random walk corresponds to Φ(0). For t > 0 the position of the walker evolves while being influenced by its history at τ , X(t− τ). To study the stability properties in the repulsive regime we concentrate our attention on the escape behavior. Thus the walker is said to escape at a first passage time L if at such time the walker crosses an arbitrarily established threshold ±X∗. The question we ask is how long the walker can be kept around the repulsive origin without escaping? Figure 4 (a) depicts the mean first passage time (MFPT), 〈L〉, for values of 1/2 < p 6 1. The MFPT shows a non-monotonous dependence on the time delay showing a maximum at a particular value of the time delay, τ∗. When p increases from 1/2 this optimal delay time moves towards smaller values. The behavior observed in the repulsive regime should be contrasted to that obtained for a DRW in the attractive regime, i.e. for p < 1/2. In this situation the MFPT also depends on the time delay. In particular, in the zero delay case the walker remains trapped around the origin for 378 Neural control on multiple time scales: Insights from human stick balancing long times. As the time delay increases the MFPT rapidly and monotonously decreases indicating that the inclusion of delay in the dynamics destabilizes the walker. Moreover, the MFPT tends to decrease with increasing p and, as expected, it reaches a constant value when p = 1/2. The mechanism by which this transient stabilization of the unstable origin occurs is still under investigation. Preliminary observations suggest that it reflects the interplay between two opposing tendencies: 1) the stabilizing effects of the past history of the walker that delay escape, and 2) the greater the probability of escape the longer the walk is. 6. Slow control The above mechanism for fast control does not readily explain why, for example, a stick cannot be balanced at the fingertip so well with the eyes closed as with the eyes open. Obviously closed- loop feedback control mechanisms that depend, in part, on visual and other sensory inputs should also be involved in this task. We refer to these unidentified controlling mechanisms as the “slow components” for control since they cannot be expected to respond quicker than the related neural conduction times. The identification of the nature of these slow control mechanisms has been more difficult [10]. The observation that closed-loop feedback control of voluntary movements become of lesser importance as expertise develops [35], suggests that it might be possible to gain insights into their role by identifying the controlling movements that occur as expertise develops. -1 0 1 ∆V ( m/sec) 10 -5 10 -4 10 -3 10 -2 P (∆ V ) Figure 5. Comparison of the probability distribution of ∆V of hand movements during stick balancing for the same person at two different skill levels [13]: low skill (•) and high skill (4). The stick length was 62 cm. The reflector closest to the fingertip was used to track the hand movements. The vertical dashed lines show the truncation point: low skill at 0.26 and high skill at 0.41. Similar changes in the probability distributions with the development of stick balancing expertise have been obtained in four individuals (in preparation). Data is taken from [13]. Suppose we model the movements of the fingertip during stick balancing as a random walk. The step size of the random walker per unit time is equal to the speed, V . In order to focus on the fast components of these movements (figure 2 (b)) it is necessary to high pass filter the time series. This is equivalent to calculating the change in speed, ∆V [13]. Figure 5 compares the distribution in the changes of the distribution of ∆V , P (∆V ),that occur when a subject was a novice stick balancer and that obtained when they became more expert after 8 days of intensive practice. The increase in expertise coincides with a change in the shape of P (∆V ), i.e. the distribution develops fatter tails at the higher skill level. Surprisingly the Lévy exponent was found to be ∼ 0.9 [13]. Lévy flights with α ≈ 1 have the characteristics of super-diffusion and have been shown to be optimal for a variety of random search patterns [36]. This observation suggests that the fast components of the control of stick balancing resemble a random foraging strategy. A simple explanation for the changes shown in figure 5 is that they reflect differences in trun- 379 J.L.Cabrera, C.Luciani, J.Milton cation [37,38]. If we write P (∆V ) = { c1Lα(∆V )f(∆V ) if |∆V | > Π, c2Lα(∆V ) if |∆V | 6 Π, where Lα(∆V ) describes the Lévy distribution in ∆V ) with Lévy index α, and c1, c2 are normali- zation constants. The truncation function, f(∆V ), is typically a decreasing function of ∆V with compact support and the truncation threshold, Π, is the critical value of ∆V at which the distribu- tion begins to deviate from Lα(∆V ). Truncation accounts for the fact that ∆V is limited by both the biomechanical properties of the musculoskeletal system (which are unlikely to change over the training period of 1–2 weeks) and the effects of time-delayed (reactive) feedback mechanisms: the more active the feedback control, the more truncated [38]. No skilled-related changes in the Lévy index, α, or the Lévy scale factor were observed [13]. Thus the slow components for the control of stick balancing that change as the expertise develops, are likely to either involve a reactive feedback at the level of local sensorimotor loops (e.g. [7]) or possibly the mechanisms based on feedforward control involving higher neural centers (e.g. [10]). 7. Discussion The parametric noise mechanism we have described for the fast component of the control for stick balancing will obviously be most important when the movements of the stick are quick compared to the neural time delay. Although this observation would argue that this mechanism should be of lesser importance once the stick movements become sufficiently sluggish, e.g. long and heavier sticks, experimental observations indicate that the situation is far from being that simple. For example, novice stick balancers are capable of improving performance, i.e. shift the stick survival curve to the right, by concurrently moving a leg rhythmically [15]. The fact that a similar improvement in performance can be obtained by having the subject think about moving the leg rhythmically, but not actually moving it [16], demonstrates that this improvement is not simply due to mechanical vibrations of the fingertip induced by leg movement. The simplest explanation is that the increase in performance is the result of the diversion of the intention for voluntary corrective movements away from the stick balancing task. In other words the improvement in performance arises because a corrective movement guided by volition is not necessarily optimal and hence their role should be minimized since they may do more harm than good. At any instant in time the control of a voluntary movement by the nervous system reflects the interplay between multiple controlling effects. Closed-loop neural feedback loops in which informa- tion obtained from sensory receptors to adjust corrective movements are organized in a hierarchical manner [39]: local reflex control loops are organized into cooperative synergies under the control of neural centers located at progressively higher neural centers (e.g. spinal cord, brainstem, cortex). As the environment and the demands of the task become more predictable, it becomes possible to use feed-forward mechanisms to program the neural activations before the movements should be made. Moreover, it follows directly from the equations of motion for every motor task, that control cannot be completely due to the effects of neural activations alone, but should also take into account the biomechanical properties of the musculoskeletal system and the changing nature of the environment in which the movements occur [40]. Often neglected is the fact that the care- ful identification of noise sources and their effects on variability can provide us with important non-invasive methods to gain insights into the structure of the underlying control networks [2,41]. What are the rules that govern the interplay between all of these controlling effects? To what extent are these rules dependent on factors such as the nature of the task, the expertise of the performer, and the demands placed on critical, but limited resources, e.g. attention? Stick balancing at the fingertip provides an important and accessible experimental paradigm to investigate these issues using non-invasive techniques. The development of virtual stick balancing tasks that involve the interplay between a human and a computer [10,14] will enable these studies to be extended into the functional neuroimaging suite and possibly to even the study of non-human primates (e.g. [42]). 380 Neural control on multiple time scales: Insights from human stick balancing Thus obtaining an understanding of how the nervous system controls complex voluntary tasks may be no further away than the tip of our finger! Acknowledgements JLC thanks Dr. Bertrand Berche and Dr. Ricardo Paredes for their invitation to participate in the school. Part of the material presented here was done in collaboration with Dr. Toru Ohira, Dr. Christian Eurich, Dr. Tadaaki Hosaka and Dr. Ronald Bormann. JLC acknowledges research standard support from IVIC (Venezuela). JM acknowledges support from the William R. Kenan Charitable Trust. References 1. Kelso J.A.S. Dynamic Patterns: The Self-organization of Brain and Behavior. MIT Press, Cambridge, Mass. 2. Milton J.G., Foss J., Hunter J.D., Cabrera J.L. Quantitative Neurosciences: Models, Algorithms, Diag- nostics, and Therapeutic Applications, (Pardolos P.M., Sackellares J.C., Carney P.R., Iasemidis L.D., eds.). Kluwer Academic Publishers, Boston 117, 2004. 3. Hohenberg P.C., Halperin B.I., Rev. Mod. Physics, 1977, 49, 435. 4. Bradshaw A., Saho J.D., Robotica, 1996, 14, 397; Huang S.J., Huang C.L., IEEE Trans. Ind. Appl., 2000, 36, 452; Xu Y.C., Iwase M., Furuta K., J. Dyn. Sys. Mod. and Control, 2001, 123, 518. 5. Baker G.L., Blackburn J.A. The Pendulum: A Case Study in Physics. Oxford University Press, New York, 2005. 6. Cabrera J.L., Milton J.G., Phys. Rev. Lett., 2002, 89, 158702. 7. Stépán G. Retarded Dynamical Systems: Stability and Characteristic Functions. Pitman Research Notes in Mathematics Series, Wiley, New York, 1989; Stépán G., Kollár L., Math Comp. Modeling, 2000, 3, 199. 8. Landry M., Campbell S.A., Morris K., Aguilar C.O., SIAM J. Appl. Dyn., 2005, 4, 333. 9. de Guzman G.C. Coordination Dynamics: Issues and Trends, ( Jirsa V.K., Kelso J.A.S., eds.). Springer, New York, 2004; Foo P., Kelso J.A.S., de Guzman C.G., J. Exp. Psychol: Human Percept. Perform., 2000, 26, 1281; Treffner P.J., Kelso J.A.S., Ecol. Psychol., 1999, 11, 103. 10. Mehta B., Schaal S., J. Neurophys., 2002, 88, 942. 11. Moss F., Milton J., Nature, 2003, 425, 911. 12. Cabrera J.L., Milton J.G., Nonlinear Studies, 2004, 11, 305. 13. Cabrera J.L., Milton J.G., Chaos, 2004, 14, 691. 14. Bormann R., Cabrera J.L., Milton J.G., Eurich C.Sw., Neurocomputing, 2004, 58–60C, 517; Cabr- era J.L., Bormann R., Eurich C., Ohira T., Milton J., Fluctuation Noise Lett., 2004, 4, L107. 15. Hosaka T., Luciani C., Ohira T., Cabrera J.L., Milton J.G., Prog. Theor. Phys. in press; Luciani C., Hosaka T., Ohira T., Milton J.G., Cabrera J.L., in preparation. 16. Milton J., Cummins J., Gunnoe J., Tolefson M., Luciani C., Cabrera J.L., Hosaka T., Ohira T., Bull. Amer. Phys. Soc., in press. 17. Mackey M.C., Milton J.G., J. Math. Biol., 1990, 28, 22. 18. Bauer W., Bertsch G.F., Phys. Rev. Lett., 1990, 65, 2213. 19. Kramers H.A., Physica, 1940, 7, 284; Legrand O., Sornette D., Phys. Rev. Lett., 1991, 65, 2172; Foss J., Moss F., Milton J., Phys. Rev. E, 1997, 55, 4536. 20. Kalbfleisch J.D., Prentice R.L. The Statistical Analysis of Failure Time Data. Wiley, Toronto, 1980. 21. Bogdanoff J.L., J. Acoust. Soc. Amer., 1962, 34, 1055; Bogdanoff J.L., Citron S.J., J. Acoust. Soc. Amer., 1964, 38, 447. 22. Yamada T., Prog. Theor. Phys., 1985, 74, 918; Platt N., Spiegel E.A., Tresser C., Phys. Rev. Lett., 1993, 70, 279; Heagy J.F., Platt N., Hammel S.N., Phys. Rev. E, 1994, 49, 1140; Venkataramani S.C., Antonsen T.M., Ott E., Sommerer J.C., Physica D, 1996, 96, 66. 23. Ding M., Yang W., Phys. Rev. E, 1995, 52, 207. 24. Cabrera J.L., Physica A, 2005, 356, 25. 25. Eurich C.W., Pawelzik K., Lecture Notes Computer Sci., 2005, 3697, 365. 26. Collins J.J., De Luca C.J., Phys. Rev. Lett., 1994, 73, 764. 381 J.L.Cabrera, C.Luciani, J.Milton 27. Zarsiorsky V.M., Duarte M., Motor Control, 1999, 3, 28; Zarsiorsky V.M., Duarte M., Motor Control, 2000, 4, 185; Botarro A., Casadio M., Morasso P.G., Sanguineti V., Human Movement Sci., 2005, 25, 588. 28. Loram I.D., Lakie M., J. Physiol. (London), 2002, 540, 1111. 29. McComb W.D. The Physics of Fuild Turbulence. Oxford University Press, Oxford, 1999. 30. Horsthemke W., Lefever R. Noise Induced Transitions. Springer-Verlag, Berlin, 1984. 31. Dayan I., Gitterman M., Weiss G.H., Phys. Rev. A, 1992, 46, 757; Gitterman M., Weiss G.H., J. Stat. Phys., 1993, 70, 107; Casado J.M., Morillo M., Phys. Rev. E., 1994, 49, 1136; Mantegna R.N., Spag- nolo B., Phys. Rev. Lett., 1996, 76, 563. 32. Agudov N.V., Malakhov A.N., Int. J. Bif. and Chaos, 1995, 5, 531; Agudov N.V., Phys. Rev. E, 1998, 57, 2618; Agudov N.V., Malakhov A.N., Phys. Rev. E, 1999, 60, 6333; Malakhov A.N., Pankratov A.L., Physica C, 1996, 46, 269. 33. Longtin A., Milton J., Bos J.E., Mackey M.C., Phys. Rev. A, 1990, 41, 6992; Longtin A., Phys. Rev. A, 1991, 44, 4801. 34. Ohira T., Milton J.G., Phys. Rev. E, 1995, 52, 3277; Ohira T., Yamane T., Phys. Rev. E, 2000, 61, 1247. 35. Pew R.W., J. Exp. Psychology, 1966, 71, 764. 36. Viswanathan G.M., Afanazyev V., Buldyrev S.V., Murphy E.J., Stanley H.E., Nature, 1996, 381, 413; Viswanathan G.M., Buldyrev S.V., Havlin S., da Luz M.G.E., Raposo E.P., Stanley H.E., Nature, 1999, 401, 911; D. Brockman, Giesel T., Neurocomputing, 2000, 32–33, 643; Buldyrev S.V., Havlin S., Kazakov A.Ya., Lu M.G.E., Raposo E.P., Stanley H.E., Viswanathan G.M., Phys. Rev. E, 2001, 64, 041108; Brockman D., Giesel T., Phys. Rev. Lett., 2003, 90, 170601. 37. Koponen I., Phys. Rev. E, 1995, 52; Mantegna R.N., Stanley H.E., Nature, 1995, 46. 38. Gupta M., Campanha J.R., Physica A, 2000, 275, 531. 39. Bernstein N. The Coordination and Regulation of Movement. Pergamon, London, 1967. 40. Chiel H.J., Beer R.D., TINS, 1997, 20, 553; Milton J., Small S., Solodkin A., J. Clin. Neurophys., 2004, 21, 134. 41. Frank T.D., Beek P.J., Friedrich R., Phys. Lett. A, 2004, 328, 219; Frank T.D., Friedrich R., Beek P.J., Phys. Lett. A, 2005, 338, 74. 42. Wessberg J., Stambaugh C.R., Kralik J.D., Beck P.D., Laubach M., Chapin J.K., Kim J., Biggs S.J., Srinivasan M.A., Nicolelis M.A.L., Nature, 2000, 408, 361. Questions and answers Q (Dragi Karevski): About the distribution of big jumps: are they rare (extreme) events? A We need to consider that in the task of stick balancing there are two important distributions where large fluctuations (i.e., big jumps) can be observed: 1) In the case of the distribution of laminar phases (figure 3 (a)) and 2) In the distribution of fingertip speed changes (figure 5). However the first one corresponds to time intervals (corrective movements) and do not reflect necessarily big fluctuations in the stick. Meanwhile large fluctuations in case 2) have low probability as can be seen from (figure 5) indicating that they are quite rare extreme events. These big changes in speed are not strictly related with big changes in the stick vertical deviation which are a persistent signature of its dynamics thus not necessaryly rare because the sticks always fall. 382 Neural control on multiple time scales: Insights from human stick balancing Ìåõàíiçì êîíòðîëþ íåðâîâî¿ ñèñòåìè íà ðiçíèõ ÷àñîâèõ ìàñøòàáàõ: áàëàíñóâàííÿ ïàëèöåþ íà êií÷èêó ïàëüöÿ Õ.Ë.Êàáðåðà1, Ê.Ëó÷iàíi2, Äæ.Ìiëüòîí3 1 Öåíòð ôiçèêè, Âåíåñóåëüñüêèé Iíñòèòóò íàóêîâèõ äîñëiäæåíü, ïîøò. ñêðèíüêà 21827, Êàðàêàñ 1020–À, Âåíåñóåëà 2 Ôiçè÷íèé ôàêóëüòåò, Óíiâåðñèòåò Ñiìîíà Áîëiâàðà, Êàðàêàñ 1080–À, Âåíåñóåëà 3 Ôàêóëüòåò ïðèðîäíè÷èõ äèñöèïëií, Êëåðåìîíò êîëåäæ, Êëåðåìîíò, ÑÀ 91711 ÑØÀ Îòðèìàíî 10 êâiòíÿ 2006 ð., â îñòàòî÷íîìó âèãëÿäi – 12 òðàâíÿ 2006 ð. Çàïiçíåíèé â ÷àñi ìåõàíiçì êîíòðîëþ çâîðîòíüîãî çâ’ÿçêó íåðâîâî¿ ñèñòåìè ïîñòiéíî çàçíຠâïëèâó íåêîíòðîëüîâàíèõ âèïàäêîâèõ çáóäæåíü (ÿêi íàçèâàòèìåìî øóìîì). Ç öiº¿ òî÷êè çîðó, ñòàòèñòè÷íi âëàñòèâîñòi ôëóêòóàöié êîíòðîëüîâàíèõ çìiííèõ ìîæóòü ñëóæèòè äæåðåëîì iíôîðìàöi¿ ïðî ïðèðî- äó öüîãî ìåõàíiçìó êîíòðîëþ. Ìè iëþñòðóºìî öå òâåðäæåííÿ âèâ÷åííÿì áàëàíñóâàííÿ ïàëèöåþ íà êií÷èêó ïàëüöÿ, âèêîðèñòîâóþ÷è òåõíiêó øâèäêiñíîãî ôiêñóâàííÿ ðóõó. Åêñïåðèìåíòàëüíi ñïîñòåðå- æåííÿ ïîºäíàíi ç ÷èñëîâèì àíàëiçîì äèôåðåíöiàëüíîãî ðiâíÿííÿ iç ñòîõàñòè÷íèì çàïiçíåííÿì ïî- êàçóþòü, ùî íà ìàñøòàáàõ ÷àñó ìàëèõ ïîðiâíÿíî ç ÷àñîì çàïiçíåííÿ íåðâîâî¿ ñèñòåìè (“øâèäêèé êîíòðîëü”) ïàðàìåòðè÷íèé øóì ñïðè÷èíÿº íåïåðåäáà÷óâàíèé ìåõàíiçì, ÿêèé ñòàáiëiçóº áàëàíñîâà- íó ïàëèöþ ó âåðòèêàëüíîìó ïîëîæåííi. ×èñëîâi ñèìóëÿöi¿ âèïàäêîâîãî áëóêàííÿ iç çàïiçíåííÿì òà ç âiäøòîâõóâàëüíèì ïî÷àòêîì âiäëiêó ïîêàçàëè, ùî íàâiòü íåñòiéêà íåðóõîìà òî÷êà ìîæå ñòàáiëiçóâà- òèñÿ çàâäÿêè ñïiëüíîìó âïëèâó øóìó i ÷àñîâîãî çàïiçíåííÿ. I íàâïàêè, íà ìàñøòàáàõ ÷àñó, ïîðiâíÿíèõ ç ÷àñîì çàïiçíåííÿ íåðâîâî¿ ñèñòåìè (“ïîâiëüíèé êîíòðîëü”) ìåõàíiçìè êîíòðîëþ çâîðîòíüîãî çâ’ÿç- êó i ïåðåäáà÷åííÿ ñòàþòü âàæëèâiøèìè. Âiäíîñíèé âêëàä øâèäêîãî i ïîâiëüíîãî ìåõàíiçìó êîíòðîëþ â ïðîöåñ áàëàíñóâàííÿ ïàëèöåþ º äèíàìi÷íèì i çàëåæèòü, íàïðèêëàä, âiä óìîâ, â ÿêèõ çäiéñíþºòüñÿ öåé ïðîöåñ òà âiä äîñâiäó òîãî, õòî çäiéñíþº áàëàíñóâàííÿ. Êëþ÷îâi ñëîâà: ôëóêòóàöiéíi ÿâèùà, ÷àñîâå çàïiçíåííÿ, îáåðíåíèé ìàÿòíèê, êîíòðîëü çâîðîòíüîãî çâ’ÿçêó, øóì PACS: 89.75.-k, 87.19.St, 02.30.Ks, 05.45.-a, 02.50.-r 383 384