基于DFA方法的健康人与癫痫病人EEG数据分析附代码

引言

       DFA分析方法是由C.-K提出的一种研究时间序列波动长时相关性的方法。主要用来区别复杂系统本身产生的波动和由外界及环境刺激作用在系统上产生的波动。外部刺激产生的变化假设引起了局部效应，而系统本身内部的动力学产生的变化假设展示了长时的相关性。DFA是尺度不变行为的另一测度方法，因为它估计了所有尺度的显示分形特性的趋势。DFA计算中包括强调由非平稳性引起的相关性减少的局部趋势，和量化代表系统本身性质的长时分形关联性特征。

       DFA提供给临床医生一种研究系统本身内在机制产生的生理信号长时相关性的方法，该研究不包括由外界产生的与系统本身不相关的刺激信号。该方法利用整个序列来进行计算而且尺度自由，所以能够提供区别生理信号的有用信息。理论上，尺度指数

       从0.5（随机序列）变化到1.5（随机步速），生理信号的尺度指数接近1。大于1的尺度指数表示长时尺度相关性的丢失和机体本身的病理改变。该技术最初应用来探测DNA序列的长时相关性，后来广泛的应用在生理时间序列的分析。

       DFA方法已经被利用来估计癫痫的发生与诊断，老年人，癫痫病人都显示了分形尺度的丢失。在心脏研究中通过对2小时ECG纪录进行分析发现，DFA能够提供传统的时域、频域分析所不能提供的信息。在与其他分析方法相比较分析24小时癫痫病发后病人的癫痫数据中，尺度指数的下降，明显的预测了脑部活动失常的猝发。在分析睡眠无呼吸病人的脑部活动变化中DFA方法明显优于谱分析方法。在癫痫病发后预测死亡率的研究中，短时尺度指数的下降是个很好的预测参数。总之，脑部波动信号中DFA尺度指数的变化能够对癫痫的发生和诊断提供辅助诊断和预测信息。

       我在本文中运用DFA方法分析了不同生理病理状态的脑电图的自相似特性，得出了一些比较结论，可以对临床诊断有重要的提示作用。

1  原理

       脑电图（electroencephalogram，简称EEG）作为脑部生理电活动的图形记录，反映脑部兴奋的产生、传导和恢复过程中的电变化，疾病对脑部神经细胞电兴奋传导的影响也会在脑电图中表现出来，如图一。因此采用脑电图对脑部生理电活动进行检测和分析一直是医学临床实践中心脏功能检测和诊断的最重要方法和手段。

图一 正常与癫痫病人的脑电图

       为了克服生理数据时间序列的高度非平稳性，提出了处理随机游走的一种改进的均方根分析方法（去趋势波动分析，简称DFA）以分析生物医学数据的自相似特性。它比传统方法（如频谱分析和Hurst分析）的优点是，它可以检测出一个似乎是非平稳时间序列的内在的自相似性，同时可避免杂散地检测可明显看出来的可能是一个人为的外在趋势自相似性。DFA算法更适用于某些慢变趋势非平稳时间序列。

具体算法如下：

       针对如图1所示的某个心电时间序列(总长度为N)，首先进行求和

         

其中，Bi是第i个数据，Bave而是所分析的心电时间序列的平均值。这种求和步骤可以映射原始时间序列到一个自相似的过程。

       接着，量度求和后的时间序列的垂直特征尺度。方法是：求和后的时间序列被分成等长为正在上传…重新上传取消的许多个小片段。画出每一个长度为的小片段的最小二乘拟和直线（它用来代表该片段中的趋势）（图2）。直线段的y坐标标记为ynk。

图二 局域的去趋势图示

（垂直线显示片段长度为10000，图中直线段代表着该求和后的时间序列的趋势由对每一个片段进行线性最小二乘拟合画出）

       其次，对求和后的时间序列进行去趋势，即在每一个片段中把yk减去局域趋势ynk。对一个给定的片段长度，这种求和并且去趋势后的时间序列的涨落的特征尺寸可以由下式计算：

 

       在所有时间尺度（片段大小）上重复上面的计算，可以得出和片段大小的关系曲线。对的斜率决定尺度指数（自相似参量）。

2  方法

       应用DFA方法来分析两种典型的生理时间序列。我们采集了5个健康人、5个癫痫病人的EEG数据.在每组数据里取10000个点，然后计算该时间序列并进行分析。

3  结果

       为了进一步探讨DFA与尺度的关系，在这次实验中，我选择了4、8、20、50、100、200、500、1000八个不同的尺度；并且，为了让实验具有普遍性，避免实验初患者的不适应带来的生理反应的影响，我们随机选择了正常人和癫痫患者EEG信号的中部信号。将这部分信号加载进python进行实验。为了研究正常人和癫痫病人的EEG信号，我们不妨选择了C3导联所采集的信号作为主要的研究对象。