CN108742517B - 一种基于Stacking单导联脑电睡眠自动分期方法 - Google Patents
一种基于Stacking单导联脑电睡眠自动分期方法 Download PDFInfo
- Publication number
- CN108742517B CN108742517B CN201810259825.3A CN201810259825A CN108742517B CN 108742517 B CN108742517 B CN 108742517B CN 201810259825 A CN201810259825 A CN 201810259825A CN 108742517 B CN108742517 B CN 108742517B
- Authority
- CN
- China
- Prior art keywords
- sleep
- algorithm
- classification
- feature
- entropy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000007958 sleep Effects 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 95
- 230000006870 function Effects 0.000 claims abstract description 27
- 238000010801 machine learning Methods 0.000 claims abstract description 21
- 230000008667 sleep stage Effects 0.000 claims abstract description 21
- 238000012216 screening Methods 0.000 claims abstract description 20
- 210000004556 brain Anatomy 0.000 claims abstract description 15
- 230000009467 reduction Effects 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- 238000007635 classification algorithm Methods 0.000 claims description 15
- 238000007637 random forest analysis Methods 0.000 claims description 15
- 238000003066 decision tree Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000000354 decomposition reaction Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 8
- 230000033764 rhythmic process Effects 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012549 training Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000007477 logistic regression Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 5
- 238000011160 research Methods 0.000 description 9
- 238000000605 extraction Methods 0.000 description 8
- 230000001020 rhythmical effect Effects 0.000 description 6
- 238000012706 support-vector machine Methods 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 230000003860 sleep quality Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000282412 Homo Species 0.000 description 1
- 206010049816 Muscle tightness Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005291 chaos (dynamical) Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 230000004630 mental health Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 230000008452 non REM sleep Effects 0.000 description 1
- 230000037053 non-rapid eye movement Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000004461 rapid eye movement Effects 0.000 description 1
- 230000036385 rapid eye movement (rem) sleep Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 208000019116 sleep disease Diseases 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4806—Sleep evaluation
- A61B5/4812—Detecting sleep stages or cycles
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Psychiatry (AREA)
- Psychology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种基于Stacking单导联脑电睡眠自动分期方法,属于机器学习算法领域。该方法包括步骤:S1:睡眠脑电信号预处理;S2:睡眠脑电多特征提取和筛选;S3:机器学习分类;S4:睡眠自动分期。本发明能够获取自适应阈值的小波函数和IIR滤波函数结合的滤波方法对脑电信号进行降噪处理,有效地提高信号的信噪比;能够对特征算法进行优化筛选获取一种新的特征参数组,作为睡眠分期的特征;能够获取一种新的准确率高的多特征与集成学习算法组合作为睡眠分期方法。
Description
技术领域
本发明属于机器学习算法领域,涉及一种基于Stacking单导联脑电睡眠自动分期方法。
背景技术
睡眠研究对于人的身心健康和日常工作具有重要意义。睡眠分期是了解睡眠状况和睡眠质量评价的关键,基于睡眠分期的睡眠评分同时也是精神病学和神经病学最重要的诊断方法。
Rechtschaffen和Kates根据人体睡眠过程中脑电图(Eletroencephalogram,EEG)、眼电图(Electrooculogram,EOG)、肌电图(Electromyography,EMG)中的眼球运动情况和肌肉张大的变化等因素,提出R&K分期标准,把睡眠分为觉醒期、非快速眼动睡眠期(NREM)和快速眼动睡眠期(REM),其中NREM又可分为睡眠Ⅰ期、睡眠Ⅱ期、睡眠且Ⅲ期和睡眠Ⅳ期。2007年AASM协会,推出了AASM睡眠标准,AASM标准将睡眠分为五期,分别是清醒期、非快速眼动期(NREM,又分为1期,2期,3期)、快速眼动期(REM)。
目前对于睡眠自动分期研究主要在于脑电信号。脑电信号幅度非常微弱。频率范围0.5-50hz。脑电信号幅值一般在50μV左右,高于100μV一般可以看做是噪音。脑电信号是一种随机性很强的非平稳信号。影响脑电的因素很多,表现的很多规律和意义至今仍没有认识。不得不利用统计处理技术如方差、均值等来检测和估计它的特征。脑电信号同时也是非线性信号。传统意义上的信号处理方法都是建立在线性的理论上的,因此,基于传统信号处理,会不可避免的丢失很多原始信号所带的信息。所以,非线性动力学、信息论、混沌理论发展起来迅速应用到脑电的分析研究中。非线性的研究越来越在脑电信号中。所以基于脑电信号特征,研究者在对睡眠脑电信号进行特征提取中,一般从时域特征、频域特征、时频域特征和非线性特征四个方面来提取单个特征或者多个特征。
睡眠脑电特征通过特征算法提取出来之后,需要利用机器学习算法或模型进行进一步的分析和处理。机器学习按照输入情况大致可以分为两类:监督式学习和非监督式学习。睡眠脑电分期大多采用监督式学习,即输入部分(训练集和测试集)包括特征数据列表和特征数据标签值列表。特征数据标签是睡眠分期状态,一般由人工专家提供。机器学习按照输出情况可分为两类:分类算法和回归算法。睡眠分期状态是离散的,属于分类问题。目前常用于睡眠分期的方法有支持向量机(Support vector machine,SVM)、人工神经网络(Artificial Neural Network,ANN)、决策树(decision tree,DT)、随机森林(Randomforest,RF)。
Baha等利用C3-A2的单导联脑电信号,再利用小波系数、排列熵和赫斯特指数作为特征算法,使用随机森林分类算法对睡眠进行六期分类,算法的准确率为97.03%。Arthur等利用C3和C4的脑电数据,再利用复杂度和相关系数,使用隐马尔科夫(HMM)模型对睡眠进行三期分类,算法的准确率为80%。Luay等利用Pz-Oz脑电信号,再利用小波系数,使用决策树分类算法对睡眠进行六期分类,算法的准确率为75%。Salih等利用C4-A1脑电信号,再利用韦尔奇分析作为特征算法,使用决策树分类算法对睡眠进行六期分类,算法的准确率82.15%。Thiago等利用Pz-Oz脑电信号,再利用方差、偏度、峰度,使用随机森林分类算法对睡眠进行六期分类,算法的准确率为90.5%。Farideh等利用Pz-Oz脑电信号,再利用小波包系数,使用人工神经网络分类算法对睡眠进行五期分类,算法的准确率为93%。Zhu等利用Pz-Oz脑电信号,再利用不同可视图(Difference Visibility Graph,DVG)作为特征算法,使用支持向量机(Support Vector Machine,SVM)分类算法对睡眠进行六期分类,算法的准确率为87.5%。kaveh等利用Fpz-Cz脑电信号,再利用RDSTFT作为特征算法,使用随机森林进行4期睡眠分类,算法准确度为92.5%。刘志勇等利用Pz-Oz脑电信号,再利用符号熵指数、去趋势波动指数、δ频率能量作为特征,使用SVM对睡眠进行五期分类,算法准确度为92.87%。林秀晶等利用Pz-Oz脑电信号,再利用均值、方差、频率能量、频谱熵,使用SVM对睡眠进行五期分类,算法准确度为80.1%。
现有技术的不足包括以下几点:
(1)睡眠脑电信号的预处理
脑电是一种随机性,非平稳性很高的生理信号,不同的节律、不同心理改变都会影响脑电波的变化。因此,EEG具有很高的时变敏感性,其信号极易被无关噪音污染。从而形成各种EEG噪音和伪迹。常见的伪迹有包括来及仪器的干扰和来自人体的如眼电伪迹(ocularartifact)、肌电伪迹(muscle artifact)。如何对常见的噪音进行有效的消除和尽可能避免在信号中产生眼电伪迹是预处理的关键。
(2)睡眠分期算法准确度
可用于睡眠脑电信号的特征很多,如时域上的均值、方差、时频域上的不同节律波的小波系数、样本熵、Reyi系数、排列熵等。基于时域和时频域上的特征参数很容易选择。但是基于非线性的算法的如样本熵、关联维数、最大Lyapunov熵等进行算法提取时如嵌入维数、延迟时间等特征参数的选择十分重要。目前没有统一的原则针对于信号是睡眠脑电信号时特征参数选择。特征选择之后,需要对特征进行优化和筛选。不同的选择原则有不同特征选择方案。如针对于分期算法高效性原则,特征算法中不应有关联维数、最大Lyapunov熵等运算周期长的特征算法。如针对于分期算法的准确度原则,需要对特征进行特征筛选工作。其目的在于找出和原信号相关系数高的特征算法。对特征算法进行如何的筛选方法决定了分期算法的准确度。
可应用于睡眠分期的机器学习算法很多,如支持向量机、随机森林、神经网络等。如何选择合适的机器学习分类器来对脑电信号进行分类。且在每个机器学习算法内部参数很多,如何利用一定准则对参数进行调优,从而提高算法准确度,是目前研究的关键。
发明内容
有鉴于此,本发明的目的在于提供一种基于Stacking单导联脑电睡眠自动分期方法,提高睡眠分期方法的准确率和分期效率,再利用睡眠分期结果设计睡眠质量评价模型,提供一种科学的参考依据,使个体和医师能够利用本发明得到睡眠分期的结果了解睡眠质量状况和某些睡眠疾病发病可能性。
为达到上述目的,本发明提供如下技术方案:
一种基于Stacking单导联脑电睡眠自动分期方法,包括以下步骤:
S1:睡眠脑电信号预处理;
S2:睡眠脑电多特征提取和筛选;
S3:机器学习分类;
S4:睡眠自动分期。
进一步,所述步骤S1具体为:使用自适应阈值的小波函数和IIR滤波函数结合的滤波方法对脑电信号进行降噪处理,提高信号的信噪比,小波函数使用自适应阈值和软阈值函数,小波降噪的参数为小波基选择‘db4’,小波分解层数为5层;IIR滤波函数使用20阶巴特沃斯滤波函数,滤波频率为5-30HZ。
进一步,所述步骤S2具体为:采用特征算法,得出特征数据组;特征算法来自于时域、时频域和非线性域三个领域;时域采用的特征算法包括:基于Hjorth三个参数活动性,移动性,复杂性,峰度、偏度、基于α、β、θ、δ方差、基于α、β、θ、δ最大幅值;时频域采用的特征算法包括:基于α、β、θ、δ的小波能量,基于δ/α的小波能量比;非线性采用的特征算法包括:模糊熵、样本熵、分形维度,基于α、β、θ、δ的样本熵,复杂度、赫斯特指数、Tsallis熵,最大李雅科夫指数、排列熵;
特征中α、β、θ、δ节律波由小波分解进行五层分解得到;模糊熵和样本熵参数为嵌入维度m和相似容限r;利用相关系数方法找出模糊熵、样本熵和排列熵最优的参数;其中模糊熵的m和r分别为1和0.15,样本熵的m和r分别为2和0.15,排列熵参数中嵌入维数为2,延迟时间为10;利用C-C办法计算最大李雅科夫指数中时间延迟为17、平均周期256为和嵌入维数为1.5;
特征筛选算法使用Fisher Score特征筛选算法,利用Fisher Score选择特征的基本流程,即先求出每个特征的Fisher Score分数,然后设定一个阈值,若某个特征的FisherScore分数大于阈值,则选择这个特征;否则,舍弃该特征,阈值设定为0.3;
从特征数据,选择模糊熵、Hjorth的复杂性(complexity)、基于β和δ的节律波样本熵、分形维度、排列熵、基于β、θ和δ节律波的小波能量、基于β节律波的方差。
进一步,所述步骤S3具体为:采用的分类算法是基于Stacking集成学习算法,通过元分类器或元回归聚合多个分类或回归模型;基础层次模型(level model)基于完整的训练集进行训练,然后元模型基于基础层次模型的输出进行训练;基础层次由不同的学习算法组成,stacking集成是异构的;
基础层采用的元学习算法包括K最近邻(K-NerestNeighbor,KNN)分类算法、随机森林(RandomForest,RF)、极端随机树(Extremely Randomized Trees,ERT)、梯度提升决策树(Gradient Boosting Decision Tree,GBDT)和多层感知机(Multi-layer Perceptron,MLP);
第二层采用通用的线性逻辑回归分类方法,即特征信号先通过基础层五个异构的元学习器,各自得出分类结果;在第二层中输入基于第一层产生的结果,通过线性分类算法得出最后的分类结果。
进一步,所述步骤S4具体为:利用机器学习算法对睡眠状态进行分期后,对分期的结果进行评估,判定分期方法的优劣;评估参数有混淆矩阵、分类准确度和Kapper系数;
混淆矩阵基本形式是基于二分类,同时适用多分类情况;真阳性(True Positive,TP)是正类预测为正类数,假阴性(False Negative,FN)是正类预测为负类数,假阳性(False Positive,FP)是将负类预测为正类,真阴性(True Negative,TN)是将负类预测为负类;基于R&K睡眠分期标准,将睡眠分为六期:Awake、N-REM stage 1、N-REM stage 2、N-REM stage 3、N-REM stage 4和REM;
分类准确度(Accuracy,ACC)定义准确分类的样本数与样本总数的比值,从混淆矩阵得出:
TP为真阳性样本数,TN为真阴性样本数,FN为假阴性样本数,FP为假阳性样本数;
利用kappa系数用衡量专家分期和本方法分期结果的吻合程度,计算公式为:
ACC为分类准确度,机会水平ACC0=1/NY,NY为分类数目;Kappa系数为0时,表示分类准确度为机会水平程度,Kappa系数为1表示最好的分类效果。
本发明的有益效果在于:
(1)获取自适应阈值的小波函数和IIR滤波函数结合的滤波方法对脑电信号进行降噪处理,有效地提高信号的信噪比;
(2)对特征算法进行优化筛选获取一种新的特征参数组,作为睡眠分期的特征;
(3)获取一种新的准确率高的多特征与集成学习算法组合作为睡眠分期算法。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明主要流程图;
图2为原始信号与降噪后信号对比图;
图3为本发明中采用的Stacking模型。
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述。
图1为本发明主要流程图,解决的技术问题:
(1)脑电信号的预处理
脑电信号的特点是高度随机性、高度的非平衡性、高度的非线性性,脑电信号的降噪处理是脑电信号处理关键。需要对传统的巴特沃斯滤波器和目前流行的小波降噪算法进行对比,同时通过比较不同的小波基函数、小波的阈值函数、小波层数、小波阈值大小。研究出新的适合睡眠脑电降噪算法。
获取自适应阈值的小波函数和IIR滤波函数结合的滤波方法对脑电信号进行降噪处理,有效的提高信号的信噪比。
(2)睡眠分期算法设计与评估
①如何选择特征参数和如何对特征参数进行筛选
对于得出的多特征值,需要对特征值进行筛选,理由如下:特征值并不是都能体现出脑电信号的真实信息,需要找出主要的、关联性好的特征参数。特征算法提取的难度、时间效率有所区别,基于时域的特征算法提取速度很快,但是特征不能很好的表征睡眠信息,而关联维数、分型维数等可以很好的表征睡眠信息但是提取速度较慢。需要对其进行平衡处理。本发明使用Fisher score特征筛选算法对特征进行筛选,选择的标准在于是否提高了分期算法的准确度。
②如何选取机器学习算法进行分期
通过特征向量,利用机器学习算法来对其进行分类,目前的机器学习算法很多,如神经网络、支持向量机、决策树、随机森林、梯度提升决策树(Gradient Boosting DecisionTree,GBDT)等。本发明使用Stacking集成学习算法进行分类,利用机器学习常见的评估方法进行比较,得出最适合本方法的机器学习方法。常用的评估参数有混淆矩阵、分类准确度、和Kapper系数。
本发明主要研究内容包括以下几个方面:
(1)睡眠脑电信号预处理
(2)睡眠脑电多特征提取和特征筛选
(3)机器学习算法的选择和性能评估
(4)睡眠自动分期
1、数据选择
本方法验证的数据来自于PhysioNet website公共数据库。数据名为Sleep-EDF。利用61个多导睡眠仪(polysomnograms,PSGs)在1987-1991和1994分别两次记录睡眠信号。第一次采集睡眠信号时,采集了20个无睡眠方面病史的正常成年人20个小时睡眠信号,其中年龄在23到25之间,男女各占一半。信号采样频率为100Hz。使用脑电导联Fpz-Cz和Pz-Oz代替传统的C4-A1和C3-A2。数据中带有专家的人工分期标准,每30s判断一次睡眠分期状态。本研究采用Pz-Oz单脑电导联,把睡眠按照R&K准则分为六期,分别为Awake、N1、N2、N3、N4、REM。从10个人提取到本次使用的原始数据如下表1所示:
表1本研究采用的睡眠数据分布图
2、信号预处理
使用自适应阈值的小波函数和IIR滤波函数结合的滤波方法对脑电信号进行降噪处理,有效地提高信号的信噪比,小波降噪的主要参数:小波基选择‘db4’,小波分解层数为5层、使用自适应阈值和软阈值函数。IIR滤波函数使用20阶巴特沃斯滤波函数,滤波频率为5-30HZ。
图2为原始信号与降噪后信号的对比图。
PER值为降噪信号的方差/原始信号的模。ERR为降噪信号和原始信号差的模。本方法中PER值为0.93.ERR值为142.7.从数值上看,降噪算法达到了很好的滤波效果。
从图上可以看出原始信号上,局部上混杂较为严重,通常这是噪音的标志。在降噪后的信号中,有效的降低了噪音。与原信号相比,细节更为具体。降噪函数同时有效的消除了原信号具有的较多的孤立尖峰点。由此可以看出,降噪函数可以明显的提高信号的纯净度。
3、特征提取和特征筛选
本研究使用15个特征算法,得出30个特征数据组。特征算法来自于时域、时频域和非线性域三个领域。表2为30特征算法的分布。本文在这里不详细对每个特征算法原理和提取细节进行阐述,只对重要的影响特征算法优劣的参数进行叙说。特征中α、β、θ、δ节律波由小波分解进行五层分解得到。模糊熵和样本熵参数为嵌入维度m和相似容限r。本研究利用相关系数方法找出模糊熵、样本熵和排列熵最优的参数。其中模糊熵的m和r分别为1和0.15,样本熵的m和r分别为2和0.15,排列熵参数中嵌入维数为2,延迟时间为10。利用C-C办法计算最大李雅科夫指数中时间延迟为17、平均周期256为和嵌入维数为1.5。
表2本研究采用的特征算法
特征筛选算法使用常见的Fisher Score特征筛选算法,利用Fisher Score选择特征的基本流程是,先求出每个特征的Fisher Score分数,然后设定一个阈值,若某个特征的Fisher Score分数大于阈值,则选择这个特征;否则,舍弃该特征。Fisher Score特征筛选计算方法简单,准确性高、可操作性强、节省运算时间。故本研究采用该特征筛选办法,阈值设定为0.3。
从30个特征数据,选择模糊熵、Hjorth的复杂性(complexity)、基于β和δ的节律波样本熵、分形维度、排列熵、基于β、θ和δ节律波的小波能量、基于β节律波的方差。
4、机器学习分类
本研究采用的分类算法是基于Stacking集成学习算法。Stacking是一种集成学习技术,通过元分类器或元回归聚合多个分类或回归模型。基础层次模型(level model)基于完整的训练集进行训练,然后元模型基于基础层次模型的输出进行训练。基础层次通常由不同的学习算法组成,因此stacking集成通常是异构的。本研究采用的Stacking模型如图3所示。
基础层采用五个异构的元学习算法包括K最近邻(K-NerestNeighbor,KNN)分类算法、随机森林(RandomForest,RF)、极端随机树(Extremely Randomized Trees,ERT)、梯度提升决策树(Gradient Boosting Decision Tree,GBDT)和多层感知机(Multi-layerPerceptron,MLP)。
第二层采用通用的线性逻辑回归分类方法。特征信号先通过基础层五个异构的元学习器,各自得出分类结果。在第二层中输入基于第一层产生的结果,通过线性分类算法得出最后的分类结果。
5、分类评估
利用机器学习算法对睡眠状态进行分期后,需要对分期的结果进行评估,来判定分期算法的优劣。常用的评估参数有混淆矩阵、分类准确度和Kapper系数。
混淆矩阵基本形式是基于二分类,同时适用多分类情况。真阳性(True Positive,TP)是正类预测为正类数(准确判断为正类)。假阴性(False Negative,FN)是正类预测为负类数。假阳性(False Positive,FP)是将负类预测为正类。真阴性(True Negative,TN)。是将负类预测为负类(准确判断为负类)。基于R&K睡眠分期标准,将睡眠分为六期:Awake,N-REM 1,N-REM 2,N-REM 3,N-REM 4,REM。建立基于睡眠六期的混淆矩阵如表3。
表3睡眠六期的混淆矩阵
从对睡眠的混淆矩阵中可以看出,算法对于各期的睡眠准确率都大于90%,对于清醒期和睡眠期主要的分期,算法达到了98.27%。可以认为算法进行高效准确的分期。同时,在比较难分期的N-REM stage 2和REM期,算法在这两期主要的错分类数据也集中在于此。但是算法在这个期间的分期准确度是可观的,超过了90%。大于一般的睡眠分期算法。
分类准确度(Accuracy,ACC)定义了准确分类的样本数与样本总数的比值,可以从混淆矩阵得出
本分类算法分期准确度为96.67%,高于最近的研究现状中睡眠六期的准确度。本文利用kappa系数用来衡量专家分期和本研究算法分期结果的吻合程度。计算公式为:
ACC为分类准确度,机会水平ACC0=1/NY,NY为分类数目。Kappa系数为0时,表示分类准确度为机会水平程度,Kappa系数为1意味着最好的分类效果。本方法的kappa系数为0.96。说明本算法进行的睡眠分期的效果具有很高的一致性,睡眠分期准确度具有可靠的参考价值。
本方法提出的特征筛选方法同样适用于其他医学生理信号,有助于发现高效的特征算法。本研究算法可以为其准备一个可视化界面应用,增加适用性。本方法的成功表明了特征参数的多样性会有效地提高脑电信号及其他生理信号的识别精度,具有广泛的应用前景。
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。
Claims (1)
1.一种基于Stacking单导联脑电睡眠自动分期方法,其特征在于:该方法包括以下步骤:
S1:睡眠脑电信号预处理;
S2:睡眠脑电多特征提取和筛选;
S3:机器学习分类;
S4:睡眠自动分期;
所述步骤S1具体为:使用自适应阈值的小波函数和IIR滤波函数结合的滤波方法对脑电信号进行降噪处理,提高信号的信噪比,小波函数使用自适应阈值和软阈值函数,小波降噪的参数为小波基选择‘db4’,小波分解层数为5层;IIR滤波函数使用20阶巴特沃斯滤波函数,滤波频率为5-30HZ;
所述步骤S2具体为:采用特征算法,得出特征数据组;特征来自于时域、时频域和非线性域三个领域;时域采用的特征包括:基于Hjorth三个参数活动性,移动性,复杂性,峰度、偏度、基于α、β、θ、δ方差、基于α、β、θ、δ最大幅值;时频域采用的特征包括:基于α、β、θ、δ的小波能量,基于δ/α的小波能量比;非线性采用的特征包括:模糊熵、样本熵、分形维度,基于α、β、θ、δ的样本熵,复杂度、赫斯特指数、Tsallis熵,最大李雅科夫指数、排列熵;
特征中α、β、θ、δ节律波由小波分解进行五层分解得到;模糊熵和样本熵参数为嵌入维度m和相似容限r;利用相关系数方法找出模糊熵、样本熵和排列熵最优的参数;其中模糊熵的m和r分别为1和0.15,样本熵的m和r分别为2和0.15,排列熵参数中嵌入维数为2,延迟时间为10;利用C-C办法计算最大李雅科夫指数中时间延迟为17、平均周期为256、嵌入维数为1.5;
特征筛选算法使用Fisher Score特征筛选算法,利用Fisher Score选择特征的基本流程,即先求出每个特征的Fisher Score分数,然后设定一个阈值,若某个特征的FisherScore分数大于阈值,则选择这个特征;否则,舍弃该特征,本算法阈值设定为0.3;
利用特征筛选算法从30个特征数据进行筛选,选择模糊熵、Hjorth的复杂性(complexity)、基于β和δ的节律波样本熵、分形维度、排列熵、基于β、θ和δ节律波的小波能量、基于β节律波的方差;
所述步骤S3具体为:采用的分类算法是基于Stacking集成学习算法,通过元分类器或元回归聚合多个分类或回归模型;基础层次模型(level model)基于完整的训练集进行训练,然后元模型基于基础层次模型的输出进行训练;基础层次由不同的学习算法组成,stacking集成是异构的;
基础层采用的元学习算法包括K最近邻(K-NerestNeighbor,KNN)分类算法、随机森林(RandomForest,RF)、极端随机树(Extremely Randomized Trees,ERT)、梯度提升决策树(Gradient Boosting Decision Tree,GBDT)和多层感知机(Multi-layer Perceptron,MLP);
第二层采用通用的线性逻辑回归分类方法,即特征信号先通过基础层五个异构的元学习器,各自得出分类结果;在第二层中输入基于第一层产生的结果,通过线性分类算法得出最后的分类结果;
所述步骤S4具体为:利用机器学习算法对睡眠状态进行分期后,对分期的结果进行评估,判定分期方法的优劣;评估参数有混淆矩阵、分类准确度和Kapper系数;
混淆矩阵基本形式是基于二分类,同时适用多分类情况;真阳性(True Positive,TP)是正类预测为正类数,假阴性(False Negative,FN)是正类预测为负类数,假阳性(FalsePositive,FP)是将负类预测为正类,真阴性(True Negative,TN)是将负类预测为负类;基于R&K睡眠分期标准,将睡眠分为六期:Awake、N-REM stage 1、N-REM stage 2、N-REMstage 3、N-REM stage 4和REM;
分类准确度(Accuracy,ACC)定义准确分类的样本数与样本总数的比值,从混淆矩阵得出:
TP为真阳性样本数,TN为真阴性样本数,FN为假阴性样本数,FP为假阳性样本数;
利用kappa系数用衡量专家分期和本方法分期结果的吻合程度,计算公式为:
ACC为分类准确度,机会水平ACC0=1/NY,NY为分类数目;Kappa系数为0时,表示分类准确度为机会水平程度,Kappa系数为1表示最好的分类效果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810259825.3A CN108742517B (zh) | 2018-03-27 | 2018-03-27 | 一种基于Stacking单导联脑电睡眠自动分期方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810259825.3A CN108742517B (zh) | 2018-03-27 | 2018-03-27 | 一种基于Stacking单导联脑电睡眠自动分期方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108742517A CN108742517A (zh) | 2018-11-06 |
CN108742517B true CN108742517B (zh) | 2023-12-29 |
Family
ID=63980421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810259825.3A Active CN108742517B (zh) | 2018-03-27 | 2018-03-27 | 一种基于Stacking单导联脑电睡眠自动分期方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108742517B (zh) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109567748B (zh) * | 2018-11-09 | 2021-09-28 | 速眠创新科技(深圳)有限公司 | 睡眠分期的确定方法、装置、计算机设备和存储介质 |
CN109363670A (zh) * | 2018-11-13 | 2019-02-22 | 杭州电子科技大学 | 一种基于睡眠监测的抑郁症智能检测方法 |
CN109700434A (zh) * | 2018-12-28 | 2019-05-03 | 武汉中旗生物医疗电子有限公司 | 基于心电图诊断的元学习模型训练方法、系统和设备 |
CN109674468A (zh) * | 2019-01-30 | 2019-04-26 | 大连理工大学 | 一种单导脑电自动睡眠分期方法 |
CN109770892A (zh) * | 2019-02-01 | 2019-05-21 | 中国科学院电子学研究所 | 一种基于心电信号的睡眠分期方法 |
CN110327040A (zh) * | 2019-04-24 | 2019-10-15 | 武汉理工大学 | 基于云平台的睡眠分期方法及系统 |
CN110584596B (zh) * | 2019-07-15 | 2022-05-27 | 天津大学 | 基于双输入卷积神经网络的睡眠阶段分类方法及应用 |
CN110558934B (zh) * | 2019-08-02 | 2022-06-28 | 杭州连帆科技有限公司 | 睡眠质量评估方法、装置、电子设备及介质 |
CN110432870A (zh) * | 2019-08-13 | 2019-11-12 | 重庆邮电大学 | 一种基于1d cnn-lstm的睡眠信号自动分期方法 |
CN110859593A (zh) * | 2019-08-22 | 2020-03-06 | 西安八水健康科技有限公司 | 一种基于宽度神经网络的睡眠阶段自动划分算法 |
CN111067513B (zh) * | 2019-12-11 | 2022-08-23 | 杭州电子科技大学 | 一种特征权重自学习的睡眠质量检测关键脑区判定方法 |
CN111053537A (zh) * | 2019-12-27 | 2020-04-24 | 电子科技大学 | 基于Heaviside核函数的生理信号非平衡性分析方法 |
CN111436904B (zh) * | 2020-03-17 | 2022-11-15 | 深圳赛安特技术服务有限公司 | 智能化睡眠分期方法、装置及计算机可读存储介质 |
CN111493822B (zh) * | 2020-03-23 | 2023-01-31 | 济南国科医工科技发展有限公司 | 一种基于睡眠脑电的快速眼动期睡眠行为障碍分类方法 |
CN111466910B (zh) * | 2020-04-30 | 2023-11-21 | 电子科技大学 | 一种睡眠监测方法、系统、存储介质、计算机程序、装置 |
CN111714118B (zh) * | 2020-06-08 | 2023-07-18 | 北京航天自动控制研究所 | 一种基于集成学习的脑认知模型融合方法 |
CN112690759B (zh) * | 2020-12-24 | 2022-04-19 | 江苏爱谛科技研究院有限公司 | 一种自动睡眠分期的建立方法及其应用 |
CN112842342B (zh) * | 2021-01-25 | 2022-03-29 | 北京航空航天大学 | 一种结合希尔伯特曲线和集成学习的心电磁信号分类方法 |
CN112827039A (zh) * | 2021-03-01 | 2021-05-25 | 中山大学 | 麻醉呼吸机的控制方法、装置及可调控麻醉呼吸机系统 |
CN113100756A (zh) * | 2021-04-15 | 2021-07-13 | 重庆邮电大学 | 一种基于Stacking的帕金森震颤检测方法 |
CN113116307A (zh) * | 2021-04-26 | 2021-07-16 | 西安领跑网络传媒科技股份有限公司 | 一种睡眠分期方法、计算机可读存储介质及程序产品 |
CN113807299B (zh) * | 2021-09-26 | 2023-08-08 | 河南工业大学 | 基于平行频域脑电信号的睡眠阶段分期方法及系统 |
CN113876339A (zh) * | 2021-10-11 | 2022-01-04 | 浙江工业大学 | 一种睡眠状态脑电特征信号特征集的构建方法 |
CN114366038B (zh) * | 2022-02-17 | 2024-01-23 | 重庆邮电大学 | 基于改进的深度学习算法模型的睡眠信号自动分期方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104970773A (zh) * | 2015-07-21 | 2015-10-14 | 西安交通大学 | 一种基于双重特征筛选的自动睡眠分期方法 |
WO2016154298A1 (en) * | 2015-03-23 | 2016-09-29 | Temple University-Of The Commonwealth System Of Higher Education | System and method for automatic interpretation of eeg signals using a deep learning statistical model |
CN106388778A (zh) * | 2016-09-21 | 2017-02-15 | 广州视源电子科技股份有限公司 | 睡眠状态分析中的脑电信号预处理方法和系统 |
CN106503799A (zh) * | 2016-10-11 | 2017-03-15 | 天津大学 | 基于多尺度网络的深度学习模型及在脑状态监测中的应用 |
CN106709469A (zh) * | 2017-01-03 | 2017-05-24 | 中国科学院苏州生物医学工程技术研究所 | 基于脑电和肌电多特征的自动睡眠分期方法 |
CN106778865A (zh) * | 2016-12-13 | 2017-05-31 | 重庆邮电大学 | 一种多域特征下的半监督脑电信号睡眠分期方法 |
CN107007278A (zh) * | 2017-04-25 | 2017-08-04 | 中国科学院苏州生物医学工程技术研究所 | 基于多参数特征融合的自动睡眠分期方法 |
CN107205652A (zh) * | 2014-12-05 | 2017-09-26 | 新加坡科技研究局 | 具有特征生成和自动映射的睡眠分析系统 |
CN107495962A (zh) * | 2017-09-18 | 2017-12-22 | 北京大学 | 一种单导联脑电的睡眠自动分期方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108697323A (zh) * | 2015-11-06 | 2018-10-23 | 生命Q全球有限公司 | 应激水平的非侵入式生理量化 |
-
2018
- 2018-03-27 CN CN201810259825.3A patent/CN108742517B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107205652A (zh) * | 2014-12-05 | 2017-09-26 | 新加坡科技研究局 | 具有特征生成和自动映射的睡眠分析系统 |
WO2016154298A1 (en) * | 2015-03-23 | 2016-09-29 | Temple University-Of The Commonwealth System Of Higher Education | System and method for automatic interpretation of eeg signals using a deep learning statistical model |
CN104970773A (zh) * | 2015-07-21 | 2015-10-14 | 西安交通大学 | 一种基于双重特征筛选的自动睡眠分期方法 |
CN106388778A (zh) * | 2016-09-21 | 2017-02-15 | 广州视源电子科技股份有限公司 | 睡眠状态分析中的脑电信号预处理方法和系统 |
CN106503799A (zh) * | 2016-10-11 | 2017-03-15 | 天津大学 | 基于多尺度网络的深度学习模型及在脑状态监测中的应用 |
CN106778865A (zh) * | 2016-12-13 | 2017-05-31 | 重庆邮电大学 | 一种多域特征下的半监督脑电信号睡眠分期方法 |
CN106709469A (zh) * | 2017-01-03 | 2017-05-24 | 中国科学院苏州生物医学工程技术研究所 | 基于脑电和肌电多特征的自动睡眠分期方法 |
CN107007278A (zh) * | 2017-04-25 | 2017-08-04 | 中国科学院苏州生物医学工程技术研究所 | 基于多参数特征融合的自动睡眠分期方法 |
CN107495962A (zh) * | 2017-09-18 | 2017-12-22 | 北京大学 | 一种单导联脑电的睡眠自动分期方法 |
Non-Patent Citations (4)
Title |
---|
Comparison of classification accuracy using Cohen’s Weighted Kappa;Arie Ben-David;《Expert Systems with Applications》;20080229;第34卷(第2期);825-832页 * |
Reranking for Stacking Ensemble Learning;Tang, Buzhou, 等;《NEURAL INFORMATION PROCESSING: THEORY AND ALGORITHMS, PT I》;20101231;575-584页 * |
唐军.Kappa评价.《中华医学统计百科全书 管理与健康统计分册[M]》.2013, * |
运动想象脑机接口的在线分类与识别算法研究;周伟伟;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180215(第2期);论文第三章 * |
Also Published As
Publication number | Publication date |
---|---|
CN108742517A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108742517B (zh) | 一种基于Stacking单导联脑电睡眠自动分期方法 | |
Hussein et al. | Robust detection of epileptic seizures based on L1-penalized robust regression of EEG signals | |
Acharya et al. | Automated diagnosis of epileptic EEG using entropies | |
CN110897639A (zh) | 一种基于深度卷积神经网络的脑电睡眠分期方法 | |
Ranjan et al. | A fuzzy neural network approach for automatic K-complex detection in sleep EEG signal | |
CN109674468A (zh) | 一种单导脑电自动睡眠分期方法 | |
Kumari et al. | Seizure detection in EEG using time frequency analysis and SVM | |
CN109758145B (zh) | 基于脑电信号因果关系的自动睡眠分期方法 | |
CN114010171B (zh) | 一种基于心跳数据的分类器设置方法 | |
Yücelbaş et al. | Automatic detection of sleep spindles with the use of STFT, EMD and DWT methods | |
CN113343805A (zh) | 一种基于rr间期心电数据和集成学习的房颤节律识别方法 | |
Held et al. | Extracting fuzzy rules from polysomnographic recordings for infant sleep classification | |
CN108567418A (zh) | 一种基于PCANet的脉搏信号亚健康检测方法及检测系统 | |
Rangappa et al. | Classification of cardiac arrhythmia stages using hybrid features extraction with k-nearest neighbour classifier of ecg signals | |
Sabor et al. | Detection of the interictal epileptic discharges based on wavelet bispectrum interaction and recurrent neural network | |
Liu et al. | Automatic sleep arousals detection from polysomnography using multi-convolution neural network and random forest | |
Raiesdana | Automated sleep staging of OSAs based on ICA preprocessing and consolidation of temporal correlations | |
Oleiwi et al. | Efficient ECG Beats Classification Techniques for The Cardiac Arrhythmia Detection Based on Wavelet Transformation. | |
Zhang et al. | Automatic Sleep Staging using Multi‐dimensional Feature Extraction and Multi‐kernel Fuzzy Support Vector Machine | |
Yuan et al. | Automatic seizure detection using logarithmic Euclidean-Gaussian mixture models (LE-GMMs) and improved deep forest learning | |
CN111603135B (zh) | 一种基于主从支持向量机的低功耗癫痫检测电路 | |
Alrawashdeh et al. | Enhanced epileptic seizure diagnosis using EEG signals with support vector machine and bagging classifiers | |
CN115399735A (zh) | 基于时频双流增强的多头注意机制睡眠分期方法 | |
Romaissa et al. | Epileptic seizure detection from imbalanced EEG signal | |
Abdel-Maboud et al. | Comparative study of machine learning techniques based on TQWT for EMG signal classification |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |