Segundo J P, Sugihara G, Dixon P, Stiber M, Bersier L F
Department of Neurobiology, University of California, Los Angeles 90095-1763, USA.
Neuroscience. 1998 Dec;87(4):741-66. doi: 10.1016/s0306-4522(98)00086-4.
This communication describes the new information that may be obtained by applying nonlinear analytical techniques to neurobiological time-series. Specifically, we consider the sequence of interspike intervals Ti (the "timing") of trains recorded from synaptically inhibited crayfish pacemaker neurons. As reported earlier, different postsynaptic spike train forms (sets of timings with shared properties) are generated by varying the average rate and/or pattern (implying interval dispersions and sequences) of presynaptic spike trains. When the presynaptic train is Poisson (independent exponentially distributed intervals), the form is "Poisson-driven" (unperturbed and lengthened intervals succeed each other irregularly). When presynaptic trains are pacemaker (intervals practically equal), forms are either "p:q locked" (intervals repeat periodically), "intermittent" (mostly almost locked but disrupted irregularly), "phase walk throughs" (intermittencies with briefer regular portions), or "messy" (difficult to predict or describe succinctly). Messy trains are either "erratic" (some intervals natural and others lengthened irregularly) or "stammerings" (intervals are integral multiples of presynaptic intervals). The individual spike train forms were analysed using attractor reconstruction methods based on the lagged coordinates provided by successive intervals from the time-series Ti. Numerous models were evaluated in terms of their predictive performance by a trial-and-error procedure: the most successful model was taken as best reflecting the true nature of the system's attractor. Each form was characterized in terms of its dimensionality, nonlinearity and predictability. (1) The dimensionality of the underlying dynamical attractor was estimated by the minimum number of variables (coordinates Ti) required to model acceptably the system's dynamics, i.e. by the system's degrees of freedom. Each model tested was based on a different number of Ti; the smallest number whose predictions were judged successful provided the best integer approximation of the attractor's true dimension (not necessarily an integer). Dimensionalities from three to five provided acceptable fits. (2) The degree of nonlinearity was estimated by: (i) comparing the correlations between experimental results and data from linear and nonlinear models, and (ii) tuning model nonlinearity via a distance-weighting function and identifying the either local or global neighborhood size. Lockings were compatible with linear models and stammerings were marginal; nonlinear models were best for Poisson-driven, intermittent and erratic forms. (3) Finally, prediction accuracy was plotted against increasingly long sequences of intervals forecast: the accuracies for Poisson-driven, locked and stammering forms were invariant, revealing irregularities due to uncorrelated noise, but those of intermittent and messy erratic forms decayed rapidly, indicating an underlying deterministic process. The excellent reconstructions possible for messy erratic and for some intermittent forms are especially significant because of their relatively low dimensionality (around 4), high degree of nonlinearity and prediction decay with time. This is characteristic of chaotic systems, and provides evidence that nonlinear couplings between relatively few variables are the major source of the apparent complexity seen in these cases. This demonstration of different dimensions, degrees of nonlinearity and predictabilities provides rigorous support for the categorization of different synaptically driven discharge forms proposed earlier on the basis of more heuristic criteria. This has significant implications. (1) It demonstrates that heterogeneous postsynaptic forms can indeed be induced by manipulating a few presynaptic variables. (2) Each presynaptic timing induces a form with characteristic dimensionality, thus breaking up the preparation into subsystems such that the physical variables in each operate as one
本通讯描述了通过将非线性分析技术应用于神经生物学时间序列可能获得的新信息。具体而言,我们考虑从经突触抑制的小龙虾起搏神经元记录的串列脉冲间隔Ti序列(“时间安排”)。如先前报道,通过改变突触前脉冲序列的平均速率和/或模式(意味着间隔离散度和序列)可产生不同的突触后脉冲序列形式(具有共享特性的时间安排集)。当突触前序列为泊松分布(独立指数分布间隔)时,形式为“泊松驱动”(未受干扰和延长的间隔不规则地相继出现)。当突触前序列为起搏器型(间隔几乎相等)时,形式要么是“p:q锁定”(间隔周期性重复)、“间歇性”(大多几乎锁定但不规则中断)、“相位遍历”(具有更短规则部分的间歇性),要么是“杂乱的”(难以预测或简洁描述)。杂乱的序列要么是“无规律的”(一些间隔正常而其他间隔不规则延长),要么是“结巴式”(间隔是突触前间隔的整数倍)。使用基于时间序列Ti中连续间隔提供的滞后坐标的吸引子重建方法分析单个脉冲序列形式。通过试错程序根据众多模型的预测性能对其进行评估:最成功的模型被视为最能反映系统吸引子的真实性质。每种形式根据其维度、非线性和可预测性进行表征。(1) 通过为系统动力学建模所需的最少变量数(坐标Ti)来估计潜在动态吸引子的维度,即通过系统的自由度来估计。测试的每个模型基于不同数量的Ti;其预测被判定成功的最小数量提供了吸引子真实维度的最佳整数近似(不一定是整数)。三到五的维度提供了可接受的拟合。(2) 通过以下方式估计非线性程度:(i) 比较实验结果与线性和非线性模型数据之间的相关性,以及(ii) 通过距离加权函数调整模型非线性并确定局部或全局邻域大小。锁定与线性模型兼容,结巴式处于边缘状态;非线性模型最适合泊松驱动、间歇性和无规律的形式。(3) 最后,将预测准确性与越来越长的间隔预测序列作图:泊松驱动、锁定和结巴式形式的准确性不变,揭示了由不相关噪声引起的不规则性,但间歇性和杂乱无规律形式的准确性迅速下降,表明存在潜在的确定性过程。对于杂乱无规律和某些间歇性形式能够实现的出色重建尤其重要,因为它们的维度相对较低(约为4)、非线性程度高且预测准确性随时间衰减。这是混沌系统的特征,并提供了证据表明相对较少变量之间的非线性耦合是这些情况下明显复杂性的主要来源。这种对不同维度、非线性程度和可预测性的证明为基于更启发式标准先前提出的不同突触驱动放电形式的分类提供了严格支持。这具有重要意义。(1) 它表明通过操纵少数突触前变量确实可以诱导异质的突触后形式。(2) 每个突触前时间安排诱导出具有特征维度的形式,从而将制剂分解为子系统,使得每个子系统中的物理变量作为一个整体运作