CN114246551A - 一种诊断精神***症的***及其应用 - Google Patents
一种诊断精神***症的***及其应用 Download PDFInfo
- Publication number
- CN114246551A CN114246551A CN202111281024.5A CN202111281024A CN114246551A CN 114246551 A CN114246551 A CN 114246551A CN 202111281024 A CN202111281024 A CN 202111281024A CN 114246551 A CN114246551 A CN 114246551A
- Authority
- CN
- China
- Prior art keywords
- ppi
- alpha
- feature
- dfa
- frequency range
- 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.)
- Pending
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/165—Evaluating the state of mind, e.g. depression, anxiety
-
- 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]
- A61B5/377—Electroencephalography [EEG] using evoked responses
-
- 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]
- A61B5/377—Electroencephalography [EEG] using evoked responses
- A61B5/38—Acoustic or auditory stimuli
-
- 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/389—Electromyography [EMG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Psychiatry (AREA)
- Psychology (AREA)
- Neurology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Acoustics & Sound (AREA)
- Neurosurgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Child & Adolescent Psychology (AREA)
- Developmental Disabilities (AREA)
- Educational Technology (AREA)
- Hospice & Palliative Care (AREA)
- Social Psychology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明属于生物医学领域,具体涉及一种诊断精神***症的***及其应用。具体地,所述***中包括根据神经认知特征和/或电生理特征诊断精神***症的计算装置。所述神经认知特征是通过RBANS或Stroop测试方法测量得到的,所述电生理特征是通过脑电图检测后分析得到的。
Description
技术领域
本发明属于生物医学领域,具体涉及一种诊断精神***症的***及其应用。
背景技术
精神***症是最严重的精神疾病之一,全世界有2000万人受到影响。精神***症多在青壮年缓慢或亚急性起病,临床上往往表现为症状各异的综合征,涉及感知觉、思维、情感和行为等多方面的障碍以及精神活动的不协调。患者一般意识清楚,智能基本正常,但部分患者在疾病过程中会出现认知功能的损害。研究表明,认知缺陷是与精神***症相关的显著神经功能障碍的核心特征之一,并且通常与预后不良相关。此外,脉冲前抑制(PPI)被认为是反映精神***症患者信息处理缺陷的指标,脉冲前抑制是基于电生理测量得到的,PPI呈现中等效应大小。此外,脑电图(EEG)是一种非侵入性电生理测量方法,广泛应用于评估大脑对外部刺激的神经反应。EEG功率谱描述了每个频带的功率分布,通常用于精神***症研究。近年来出现了更先进的脑电分析方法,如DFA方法(detrended fluctuationanalysis)、LRTC(long-range temporal correlations)方法。
目前的精神***症的诊断金标准是在国际疾病分类第11版(ICD-11)或《精神疾病诊断和统计手册》第5版(DSM-5)的基础上建立的。这些诊断方法依赖于描述性精神病理学,在一定程度上受精神病医生的主观判断的影响。因此,客观地测量并诊断精神***症是临床医生的迫切需要。然而,由于病因和临床变异的异质性,我们仍然缺乏诊断精神***症的标志物。
发明内容
为提供更适于临床使用的精神***症诊断***,本发明基于Logistics算法、随机森林算法、XGBoost算法开发了包括神经认知特征和/或电生理特征的诊断***。
特征组合
一方面,本发明提供了一种诊断精神***症的特征组合,所述特征组合包括属于神经认知特征或电生理特征中的至少两种特征;
所述神经认知特征选自IMM(即刻记忆力,immediate memory)、LAN(语言能力,language functioning)、ATT(注意力,attention)、DEM(delayed memory,延迟记忆力)、INT-C(颜色干扰下时长,color interference time)、INT-W(字义干扰下时长,wordinterference time)。
所述电生理特征选自PPI,Abs-T(θ频段绝对功率,absolute theta power)、Abs-A(α频段相对功率)、Abs-AFp/AO(额枕区α频段绝对功比值,absolute power AFp/AOratio)、Abs-(D+T)/(A+B)(绝对功率δ+θ/α+β)、Rel-D(δ相对功率,relative deltapower)、Rel-T(θ相对功率,relative theta power)、Rel-A/B(相对功率α/β)、DFA-A(α频段DFA)、DFA-B(β频段DFA)。
优选地,所述神经认知特征是通过RBANS或Stroop测试方法测量得到的。
优选地,所述电生理特征是通过脑电图检测后分析得到的。
优选地,所述PPI包括PSC-PPI(知觉空间融合PPI,perceived spatial co-location prepulse inhibition)、PSS-PPI(知觉空间分离PPI,perceived spatialseparation prepulse inhibition)。
优选地,所述特征组合选自以下任意一组:
1)IMM、LAN、ATT、DEM、INT-C、INT-W;
2)PSC-PPI、PSS-PPI、Abs-T、Abs-A、Abs-AFp/AO、Abs-(D+T)/(A+B)、Rel-D、Rel-T、Rel-A/B、DFA-A、DFA-B;
3)1)和2)的合集。
优选地,所述IMM、LAN、ATT、DEM是通过RBANS测量的到的。
优选地,所述INT-C、INT-W是通过Stroop测试测量得到的。
优选地,所述PSC-PPI、PSS-PPI是通过震惊反射测试得到的。
优选地,所述Abs-T、Abs-A、Abs-AFp/AO、Abs-(D+T)/(A+B)、Rel-D、Rel-T、Rel-A/B、DFA-A、DFA-B是通过脑电图测试得到的。
***
另一方面,本发明提供了一种诊断精神***症的***,所述***包括根据任意一种诊断特征或任意一种诊断特征组合的检测结果诊断受试者是否患有精神***症的计算装置;
所述诊断特征包括IMM、LAN、ATT、DEM、INT-W、PPI、DFA-A、DFA-B;
所述诊断特征组合选自NSF、ESF、ASF;
所述NSF包括IMM、LAN、ATT、DEM、INT-C、INT-W;
所述ESF包括PSC-PPI、PSS-PPI、Abs-T、Abs-A、Abs-AFp/AO、Abs-(D+T)/(A+B)、Rel-D、Rel-T、Rel-A/B、DFA-A、DFA-B;
所述ASF是NSF和ESF的合集。
优选地,所述PPI包括PSC-PPI、PSS-PPI。
以上诊断特征的Cohen's d值皆大约0.5,AUC皆大于0.7,单独作为诊断特征时就可以有较好的诊断效果。
本发明所述“NSF”是neurocognitive selected features的缩写,也称为神经认知特征子集。
本发明所述“ESF”是electrophysiological selected features的缩写,也称为电生理特征子集。
本发明所述“ASF”是All selected features的缩写,也就是NSF和ESF的合集,也称为全特征集。
优选地,所述***还包括神经特征检测装置和/或电生理检测装置:
所述神经特征检测装置包括:
1)用于运行RBANS的第一神经特征检测装置;
2)用于运行Stroop测试的第二神经特征检测装置;
所述电生理检测装置包括:
1)用于运行人类震惊反射***检测的第一电生理检测装置;
2)用于运行脑电图***的第二电生理检测装置;
优选地,所述***还包括用于输入受试者信息的信息收集装置;
优选地,所述受试者信息包括任意一种诊断特征或任意一种诊断特征组合的检测结果、人口统计学信息、临床特征信息。
优选地,所述诊断结果包括患有精神***症、未患有精神***症。
优选地,所述***还包括使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法建模的建模装置。
所述Logistics算法、随机森林算法、XGBoost算法是指将前述神经认知特征子集、电生理特征子集或全特征集分别通过相关算法构建成模型。
优选地,所述***还包括用于发送诊断结果的结果发送装置,所述结果发送装置可以将诊断结果发送到患者或医护人员可以查阅的信息通信终端装置。
本发明所述“算法”包括Logistics算法、随机森林算法、XGBoost算法。
设备
另一方面,本发明提供了一种诊断精神***症的设备,所述设备包括:存储器和处理器;
所述存储器用于存储程序指令;
所述处理器用于调用程序指令,当程序指令被执行时,用于执行建立模型和/或诊断的操作:
优选地,所述建立模型的操作步骤是使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法处理诊断特征组合;
优选地,所述诊断的操作步骤是:
1)获得来着于受试者的任意一种诊断特征或任意一种诊断特征组合的检测结果;
2)根据任意一种诊断特征或任意一种诊断特征组合计算受试者是否患有精神***症;
3)显示诊断结果。
所述诊断特征包括IMM、LAN、ATT、DEM、INT-W。
所述诊断特征组合包括NSF、ESF、ASF;
所述NSF包括IMM、LAN、ATT、DEM、INT-C、INT-W;
所述ESF包括PSC-PPI、PSS-PPI、Abs-T、Abs-A、Abs-AFp/AO、Abs-(D+T)/(A+B)、Rel-D、Rel-T、Rel-A/B、DFA-A、DFA-B;
所述ASF是NSF和ESF的合集。
计算机可读存储介质
另一方面,本发明提供了一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现建立模型的方法和/或诊断方法:
所述建立模型的方法包括使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法处理诊断特征组合;
所述诊断方法包括以下步骤:
1)获得来着于受试者的任意一种诊断特征或任意一种诊断特征组合的检测结果;
2)根据任意一种诊断特征或任意一种诊断特征组合的检测结果计算受试者是否患有精神***症;
3)显示诊断结果。
应用
另一方面,本发明提供了IMM、LAN、ATT、DEM、INT-W、PSC-PPI、PSS-PPI、DFA-A、DFA-B在制备诊断精神***症的产品中的应用。
另一方面,本发明提供了前述特征组合、***、设备、计算机可读存储介质在制备诊断精神***症的产品中的应用。
方法
另一方面,本发明提供了一种建立诊断诊断精神***症的模型的方法,所述方法包括使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法处理诊断特征组合;
所述诊断特征组合选自神经认知特征子集、电生理特征子集、全特征集。
同时,本发明提供了该方法建立的模型在制备诊断精神***症的产品中的应用。
另一方面,本发明提供了一种诊断精神***症的方法,所述方法包括:
1)获得来着于受试者的前述任意一种诊断特征或任意一种诊断特征组合的检测结果;
2)根据任意一种诊断特征或任意一种诊断特征组合的检测结果计算受试者是否患有精神***症。
附图说明
图1是NSF、ESF、ASF使用Logistics算法构建的模型的诊断ROC曲线,纵坐标为敏感性,横坐标为特异性。
图2是NSF、ESF、ASF使用随机森林算法构建的模型的诊断ROC曲线,纵坐标为敏感性,横坐标为特异性。
图3是NSF、ESF、ASF使用XGBoost算法构建的模型的诊断ROC曲线,纵坐标为敏感性,横坐标为特异性。
图4是使用Logistics算法构建的模型的诊断散点图,NSF拟合值作为水平坐标,ESF拟合值作为纵向坐标。
图5是使用随机森林算法构建的模型的诊断散点图,NSF拟合值作为水平坐标,ESF拟合值作为纵向坐标。
图6是使用XGBoost算法构建的模型的诊断散点图,NSF拟合值作为水平坐标,ESF拟合值作为纵向坐标。
具体实施方式
下面结合实施例对本发明做进一步的说明,以下所述,仅是对本发明的较佳实施例而已,并非对本发明做其他形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更为同等变化的等效实施例。凡是未脱离本发明方案内容,依据本发明的技术实质对以下实施例所做的任何简单修改或等同变化,均落在本发明的保护范围内。
实施例1、影响因素筛选及验证
1、实验对象
这项研究中共69名精神***症患者和50名健康对照者参与。研究人员通过美国精神障碍诊断统计手册第四版(DSM-IV)结构化临床访谈(SCID)工具确认诊断。所有受试者均为右利手,听力评估(纯音测听,1000Hz)正常。精神***症患者登记的纳入标准如下:
1)所有临床稳定受试者均无神经***疾病或头部外伤史,
2)在过去六个月内无电休克治疗史,
3)无酒精/药物依赖或滥用史(烟草除外)。
由于病情不稳定或智商低于70,患者被排除在外。
在本研究期间,所有患者照常接受抗精神病药物治疗。健康对照组(CON)与精神***症组(SCZ)性别、年龄、受教育年限和吸烟史匹配。健康对照组的排除标准包括药物滥用、***风险、重大头部创伤和神经精神障碍。
在签署知情同意书之前,每位受试者都收到了参与研究的目的和程序的详细说明。北京安定医院独立伦理委员会批准了本研究。采用阳性和阴性症状量表(PANSS)对患者进行精神病理学评估。详细人口统计学和临床特征如表1。
表1.健康对照组和精神***症组的人口学和临床特征
注:a表示卡方检验的P值;b表示独立样本t检验的P值。
2、神经认知评估
我们使用可重复的成套神经心理状态测量(RBANS)对受试者的神经认知功能进行评估。RBANS评估了五个认知领域:IMM、DEM、VC、ATT、LAN。除RBANS外,还进行了Stroop测试。每个受试者被要求完成两个干扰任务,并记录颜色干扰下时长(INT-C)和字义干扰下时长(INT-W)。
在这部分实验中,共提取了5个RBANS特征(IMM、DEM、VC、ATT、LAN)和2个Stroop特征(INT-C、INT-W)。
3、电生理评估
3.1震惊反射测量(前脉冲抑制测量)
受试者舒适地坐在躺椅上,手臂以自然姿势完全放松。利用人类震惊反射定制***(北京天鸣宏远公司生产)通过记录右眼轮匝肌记录肌电信号测量听觉震惊反射。电极阻抗保持在<5kΩ。此外,肌电图经带通滤波至100-1000Hz,并放大10,000倍。听觉惊吓刺激通过两个耳机双耳呈现。利用声压校正仪(Larson Davis,AUDit and System 824)进行声压校正。
实验环境:听觉隔音屏蔽室中进行,保持一定光亮和温度。
声音素材:所需声音素材皆由使用MATLAB函数库中的“randn()”生成,抽样率为48kHz。生成后的高斯噪声通过512阶、截止频率为10kHz的低通数字滤波器得到宽带噪声,作为背景噪声。另外生成长度为150ms宽带噪声的前刺激声音,及长度为40ms宽带噪声的惊刺激声音。为避免能量飞溅现象,每个声音刺激都是在前脉冲刺激声音及惊跳刺激声音引入后再滤波。声音信号利用声卡(Audio CODEC'97)输入到Sennheiser HD600耳机呈现给被试。利用声压校正仪(Larson Davis,AUDit and System 824)进行声压校正。以下是各声音具体参数,a)背景声:白噪声,分为左声道或右声道领先3ms,持续时间15s,声压级60dB;b)前刺激:白噪声,分为左声道或右声道领先3ms,持续时间150ms,声压级65dB;c)惊刺激:白噪声,持续时间40ms,100dB。
测试流程:嘱被试坐在检查椅上,于右眼瞳孔下、外侧1.5cm贴2个Ag/AgCl电极记录眼轮匝肌电活动,于右乳突贴电极作为地线。眼睛距眼动仪显示器60cm。全身肌肉放松,嘱实验过程中尽量保持清醒,双眼注视屏幕,尽量保持头部固定不动。首先是对对各声音素材进行学习,每个声音素材学习3遍,对所有被试给予相同的指导语,背景声和前刺激需被试分辨左声道或右声道领先,对被试进行单眼3点刻度眼动校准,嘱被试注视前方十字,开始测试。
测试范式:前刺激持续时间为150ms,分为左领先及右领先。前刺激与背景声(左领先和右领先)呈现出2种知觉空间分离状态(前刺激左领先,背景声右领先;前刺激右领先,背景声左领先),每种状态重复5次。前刺激和惊刺激间隔120ms。仅有惊刺激的试次为7次。具体可见图2-1。
统计指标:刺激前500ms记为基线S0;仅含惊刺激记为S1,同时含前刺激记为S2;计算肌电变化百分比1-(S2-S0)/(S1-S0);肌电最大值潜伏期;肌电最大反应速度;持续注意水平(如持续盯住十字的时间)。
在这部分实验中,共提取了2个特征知觉空间分离PPI和知觉空间融合PPI(PSC-PPI,PSS-PPI)。
3.2脑电记录与处理
3.2.1 EEG数据预处理
受试者舒适地坐在躺椅上,闭上眼睛,保持放松和安静5分钟。使用脑电图***(EGI,Electrical Geodesics,Inc.,America.)以1000Hz采样率采集脑电数据,以Cz作为参考。在MATLAB(MATLAB Release 2017b,MathWorks,Inc.)中使用EEGLAB(v2019.1)和FieldTrip工具箱对EEG数据进行离线预处理。首先将EEG原始数据降采样至500Hz,并将带通滤波至0.5-45Hz.。对于每个受试者的脑电数据,使用EEGLAB通过独立成分分析(ICA,算法:runica)进行伪迹去除。ICA成分使用EEGLAB插件ICLabel工具进行分类。眼动、眨眼、心跳、肌肉活动或其他伪迹通过ICA算法分离,并插值坏导。然后,手动检查EEG数据以验证伪迹消除情况,剔除坏段。
最后,将所有电极重参考至全脑平均。
3.2.2功率谱特征
使用快速傅里叶变换(FFT、Welch方法、2秒滑动窗口、50%重叠、0.5Hz频率步长)评估每个电极的功率谱密度(PSD),得出范围为0.5-45赫兹。频带选择如下:δ(1.0-4.0Hz)、θ(4.0-8.0 9/19Hz)、α(8.0-14.0Hz)、β(14.0-30.0Hz)。D、T、A、B分别表示δ、θ、α和β频带。通过左半球(Fp1、F3、C3、P3、O1、F7、T3、T5)和右半球(Fp2、F4、C4、P4、O2、F8、T4、T6)α带的平均功率计算AL和AR。(D+T)L和(D+T)R的平均值为左右半球δ和θ的总和。AFp和AO通过平均Fp通道(Fp1,Fp2)和O通道(O1,O2)中的α带功率计算。对于每个电极,计算每个频带的绝对功率(Abs)和相对功率(Rel)。
在PSD中,共提取了20个特征:Abs-D、Abs-T、Abs-A、Abs-B、Abs-A/T、Abs-A/B、Abs-(D+T)/(A+B)、Abs-(D+T)L/(D+T)R、Abs-AL/AR、Abs-AFp/AO、Rel-D、Rel-T、Rel-A、Rel-B、Rel-A/T、Rel-A/B、Rel-(D+T)/(A+B)、Rel-(D+T)L/(D+T)R、Rel-AL/AR、Rel-AFp/AO,具体参见表2。
3.2.3去趋势波动(DFA)特征
DFA是一种基于无标度理论的分析方法,用于估计幂律形式的长程时间相关性(LRTC)。也就是说,如果时间序列数据具有缓慢衰减的自相关的非随机时间结构,则DFA可以量化这些相关性衰减的速度,如DFA幂律指数所示。一些证据表明DFA反映了大脑的成熟,可能被证明是神经发育障碍病理生理学的潜在生物标志物。
DFA的计算使用神经生理学生物标记物工具箱(NBT:http://www.nbtwiki.net/)。首先,所有电极分别以δ、θ、α和β振荡进行过滤。然后从每个频带生成振幅包络。
最后,每个参与者估计每个电极的DFA值,并分别存储每个频带的DFA值。DFA-D、DFA-T、DFA-A、DFA-B通过在δ、θ、α和β频带中平均所有电极来计算。
3.2.4分形维数(FD)特征
大脑复杂性可以描述为在纯随机性(如白噪声)和无变异性(恒定性或纯周期性)之间的中间情况下,在EEG信号中观察到的高度结构化的时间结构。我们使用了EEGLAB插件myFractal(https://github.com/ramicodes/myFractal)计算每个电极的FD。
最后,通过平均所有电极的FD值来提取FD特征。
4、统计分析
使用R软件(版本3.6.3)在RStudio(版本1.2.5033,RStudio,Inc.)中进行统计。人口统计学和基本临床数据,包括性别、年龄、受教育年限、吸烟史、病程、发病年龄、氯丙嗪等效剂量,以及PANSS总分、PANSS阳性分、PANSS阴性分和PANSS一般精神病理学分。除性别和吸烟史外,所有人口统计学和临床变量均以平均值±SDs表示。采用独立t检验和卡方检验评估CON和SCZ在人口统计学、临床变量、神经认知和电生理学方面的潜在差异。根据Benjamini-Hochberg方法,计算错误发现率(FDR)以调整多次测试的P值。P<0.05被认为具有统计学意义。
5、分类
所有分析均使用R3.6.3进行。为了选择用于SCZ和CON分类的最佳特征,我们首先评估了每个特征的分类能力。使用自编译函数计算每个特征的Cohen’s d值。每个特征的受试者操作特征曲线(ROC)分析是使用R中的pROC软件包创建。每个特征的ROC值包括准确性、敏感性、特异性和ROC曲线下面积(AUC)。
将所有特征分为两类:神经认知特征和电生理特征。使用rfe函数(R caret软件包)用于通过多元递归特征消除方法选择特征。rfe首先使用bagged树算法使模型适合所有特征。根据每个特征对模型的重要性对其进行排序。在特征选择的每次迭代中,保留已排序的特征,重新装配模型并评估性能。每组特征的最终选择基于10倍交叉验证。
然后得到了两个最佳子集(神经认知特征集NSF和电生理特征集ESF)和一个组合特征集ASF(包含前述神经认知特征集和电生理特征集)。
在确定以上三个特征集合后,使用Logistic算法(R统计数据包)、随机森林算法(R随机森林包)和XGBoost算法(R XGBoost包)从三个特征集中估计分类模型。
最后,使用10倍交叉验证方法来验证这些模型的性能。然后对验证结果进行平均。通过准确性、敏感性和特异性评估分类性能。此外,还使用ROC曲线评估了各模型的性能。
实验结果1:SCZ和CON之间所有特征的统计比较
表2.从神经认知和电生理评估中提取的所有特征的统计比较结果
注:*P<0.05;**P<0.01;***P<0.001。
所有特征的统计分析结果如表2所示。总的来说,SCZ和CON的所有神经认知特征在统计学上存在差异。还可以观察到,两组之间的所有PPI特征显著不同。在脑电功率谱特征中,只有Abs-T、Abs-AFp/AO、Rel-D、Rel-T有统计学意义的差异。此外,CON的θ波段DFA(DFA-T)低于SCZ,而在α和β波段CON显著高于SCZ(DFA-A,DFA-B)。SCZ组和CON组的FD无显著性差异。在所有功能中,IMM、LAN、ATT、DEM、PSS-PPI和DFA-B的P值均小于0.001。
实验结果2:单特征的诊断性能
为了选择区分SCZ和CON的最佳特征,我们首先评估每个特征的Cohen's d和ROC值。ROC值包括准确性(%)、敏感性(%)、特异性(%)和AUC(%)。
如表3所示,共有14个特征的Cohen's d值大于0.5达到了中等效应以上的水平,包括IMM、VC、LAN、ATT、DEM、INT-C、INT-W、PSC-PPI、PSS-PPI、Abs-T、Abs AFp/AO、Rel-T、DFA-A、DFA-B。其中,LAN、TT、DEM和PSS-PPI的Cohen’s d大于0.8,表示大效应水平;IMM的Cohen’s d达到1.42显示出非常大的效应水平。在ROC值方面,神经认知特征表现出比电生理特征更好的表现。IMM是最具神经认知特征的分类方法,准确度为84.03%,其AUC达到91.87%。
表3.Cohen's d与单特征分类性能
注:*Cohen's d值一列,灰色标注的表格代表Cohen's d值>0.5,达到中等效应以上;
*AUC值一列,灰色标注的表格代表AUC>70%,达到准确诊断的标准。
有5个来自神经认知特征(IMM、LAN、ATT、DEM、INT-W)的AUC大于70%,显示这些特征具有可用于开发辅助诊断精神***症工具的潜力。PSS-PPI是最佳的潜在电生理特征,准确度为80.67%,其AUC为84.32%。
有4个电生理特征(PSC-PPI、PSS-PPI、DFA-A、DFA-B)的AUC大于70%,显示这些特征具有可用于开发辅助诊断精神***症工具的潜力。
实验结果3:特征组合的分类性能
使用多元递归特征消除方法选择神经认知特征集和电生理特征集,以求能最大程度的区分CON和SCZ。此外,使用10倍交叉验证进行模型验证,以防止过度拟合。
最终得到的神经认知特征子集(NSF子集)包含IMM、LAN、TT、DEM、INT-C、INT-W。电生理特征子集(ESF子集)包含PSC-PPI、PSS-PPI、Abs-T、Abs-A、Abs-AFp/AO、Abs-(D+T)/(A+B)、Rel-D、Rel-T、Rel-A/B、DFA-A、DFA-B。ASF集是如上两个特征子集的合集。
然后,使用Logistics算法、随机森林算法、XGBoost算法,从NSF子集、ESF子集和ASF集建立分类模型。在此基础上,采用10倍交叉验证法对这些模型进行评估。最后,精神***症患者概率的拟合值介于0到1之间。通过将截止值设置为0.5,计算准确度、敏感性、特异性,然后根据拟合值评估ROC曲线并计算AUC。
这些模型评估指标(准确性、敏感性、特异性和AUC)如表4所示。图1-3显示了每个模型的ROC曲线。从表4和图1-3中可以看出,无论采用何种算法,通过包含所有选定特征集开发的模型都显示出较好的分类性能。NSF子集模型和ESF子集模型显示出大致相当的分类准确度,但在AUC值范围内,NSF子集模型显示出比ESF子集更好的性能。
此外,在所有比较的算法中,XGBoost算法在精度和AUC方面都表现得更加稳定和准确。采用XGBoost算法的ASF集合模型达到了最优的诊断效果。
为了更好地展示这三种算法中分类能力的差异,我们使用NSF子集拟合值作为水平坐标,ESF子集拟合值作为纵向坐标创建散点图(图4-6)。与Logistic算法相比,randomforest和XGBoost算法的拟合值在该图中分布更为密集。这也表明random forest和XGBoost算法在区分精神***症患者和健康人方面提供了更好的性能。
表4.特征组合的分类性能
Claims (10)
1.一种诊断精神***症的特征组合,所述特征组合包括属于神经认知特征或电生理特征中的至少两种特征;
所述神经认知特征选自即刻记忆力、语言能力、注意力、延迟记忆力、颜色干扰下时长、字义干扰下时长;
所述电生理特征选自PPI、θ频段绝对功率、α频段相对功率、额枕区α频段绝对功比值、绝对功率δ+θ/α+β、δ相对功率、θ相对功率、相对功率α/β、α频段DFA、β频段DFA;
优选地,所述PPI包括知觉空间融合PPI、知觉空间分离PPI。
2.如权利要求1所述的特征组合,其特征在于,所述特征组合选自以下任意一组:
1)即刻记忆力、语言能力、注意力、延迟记忆力、颜色干扰下时长、字义干扰下时长;
2)PPI、θ频段绝对功率、α频段相对功率、额枕区α频段绝对功比值、绝对功率δ+θ/α+β、δ相对功率、θ相对功率、相对功率α/β、α频段DFA、β频段DFA;
3)1)和2)的合集;
优选地,所述PPI包括知觉空间融合PPI、知觉空间分离PPI。
3.如权利要求1或2所述的特征组合,所述即刻记忆力、语言能力、注意力、延迟记忆力是通过RBANS测量得到的,
所述颜色干扰下时长、字义干扰下时长是通过Stroop测试方法测量得到的,
所述知觉空间融合PPI、知觉空间分离PPI是通过震惊反射***测试得到的,
所述θ频段绝对功率、α频段相对功率、额枕区α频段绝对功比值、绝对功率δ+θ/α+β、δ相对功率、θ相对功率、相对功率α/β、α频段DFA、β频段DFA是通过脑电图测试得到的。
4.一种诊断精神***症的***,所述***包括根据任意一种诊断特征或任意一种诊断特征组合的检测结果诊断受试者是否患有精神***症的计算装置;
所述诊断特征包括即刻记忆力、语言能力、注意力、延迟记忆力、字义干扰下时长、PPI、α频段DFA、β频段DFA;
所述诊断特征组合选自神经认知特征子集、电生理特征子集或全特征集;
所述神经认知特征子集包括即刻记忆力、语言能力、注意力、延迟记忆力、颜色干扰下时长、字义干扰下时长;
所述电生理特征子集包括知觉空间融合PPI、知觉空间分离PPI、θ频段绝对功率、α频段相对功率、额枕区α频段绝对功比值、绝对功率δ+θ/α+β、δ相对功率、θ相对功率、相对功率α/β、α频段DFA、β频段DFA;
所述全特征集是神经认知特征子集和电生理特征子集的合集;
优选地,所述PPI包括知觉空间融合PPI、知觉空间分离PPI。
5.如权利要求4所述的***,其特征在于,所述***还包括神经特征检测装置和/或电生理检测装置:
所述神经特征检测装置包括以下至少一种:
1)用于运行RBANS的第一神经特征检测装置;
2)用于运行Stroop测试的第二神经特征检测装置;
所述电生理检测装置包括以下至少一种:
1)用于运行震惊反射***的第一电生理检测装置
2)用于运行脑电图***的第二电生理检测装置。
6.如权利要求4或5所述的***,其特征在于,所述***还包括使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法建模的建模装置;
优选地,所述***包括以下装置:
1)使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法建模的建模装置;
2)神经特征检测装置和/或电生理检测装置;
3)任意一种诊断特征或任意一种诊断特征组合的检测结果诊断受试者是否患有精神***症的计算装置;
4)用于显示诊断结果的显示装置。
7.一种建立诊断诊断精神***症的模型的方法,所述方法包括使用算法中任意一种方法处理诊断特征组合;
所述诊断特征组合选自神经认知特征子集、电生理特征子集或全特征集;
所述神经认知特征子集包括即刻记忆力、语言能力、注意力、延迟记忆力、颜色干扰下时长、字义干扰下时长;
所述电生理特征子集包括知觉空间融合PPI、知觉空间分离PPI、θ频段绝对功率、α频段相对功率、额枕区α频段绝对功比值、绝对功率δ+θ/α+β、δ相对功率、θ相对功率、相对功率α/β、α频段DFA、β频段DFA;
所述全特征集是神经认知特征子集和电生理特征子集的合集;
优选地,所述算法包括Logistics算法、随机森林算法、XGBoost算法;
优选地,所述算法是随机森林算法或XGBoost算法。
8.一种诊断精神***症的设备,所述设备包括:存储器和处理器;
所述存储器用于存储程序指令;
所述处理器用于调用程序指令,当程序指令被执行时,用于执行建立模型和/或诊断的操作:
所述建立模型的操作步骤是使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法处理诊断特征组合;
所述诊断的操作步骤是:
1)获得来着于受试者的任意一种诊断特征或任意一种诊断特征组合的检测结果;
2)根据所述检测结果计算受试者是否患有精神***症;
3)显示诊断结果;
所述诊断特征包括即刻记忆力、语言能力、注意力、延迟记忆力、字义干扰下时长、知觉空间融合PPI、知觉空间分离PPI、Α频段DFA、Β频段DFA;
所述诊断特征组合选自神经认知特征子集、电生理特征子集或全特征集;
所述神经认知特征子集包括即刻记忆力、语言能力、注意力、延迟记忆力、颜色干扰下时长、字义干扰下时长;
所述电生理特征子集包括知觉空间融合PPI、知觉空间分离PPI、θ频段绝对功率、α频段相对功率、额枕区α频段绝对功比值、绝对功率δ+θ/α+β、δ相对功率、θ相对功率、相对功率α/β、α频段DFA、β频段DFA;
所述全特征集是神经认知特征子集和电生理特征子集的合集。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现建立模型的方法和/或诊断方法:
所述建立模型的方法包括使用Logistics算法、随机森林算法、XGBoost算法中任意一种方法处理诊断特征组合;
所述诊断方法包括以下步骤:
1)获得来着于受试者的任意一种诊断特征或任意一种诊断特征组合的检测结果;
2)根据所述检测结果计算受试者是否患有精神***症;
3)显示诊断结果;
所述诊断特征组合选自神经认知特征子集、电生理特征子集或全特征集;
所述神经认知特征子集包括即刻记忆力、语言能力、注意力、延迟记忆力、颜色干扰下时长、字义干扰下时长;
所述电生理特征子集包括知觉空间融合PPI、知觉空间分离PPI、θ频段绝对功率、α频段相对功率、额枕区α频段绝对功比值、绝对功率δ+θ/α+β、δ相对功率、θ相对功率、相对功率α/β、α频段DFA、β频段DFA;
所述全特征集是神经认知特征子集和电生理特征子集的合集。
10.即刻记忆力、语言能力、注意力、延迟记忆力、字义干扰下时长、PPI、Α频段DFA、Β频段DFA、权利要求1所述的特征组合、权利要求4所述的***、权利要求7所述的方法构建的模型、权利要求8所述的设备、权利要求9所述的计算机可读存储介质在制备诊断精神***症的产品中的应用;
优选地,所述PPI包括知觉空间融合PPI、知觉空间分离PPI。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111281024.5A CN114246551A (zh) | 2021-11-01 | 2021-11-01 | 一种诊断精神***症的***及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111281024.5A CN114246551A (zh) | 2021-11-01 | 2021-11-01 | 一种诊断精神***症的***及其应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114246551A true CN114246551A (zh) | 2022-03-29 |
Family
ID=80792254
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111281024.5A Pending CN114246551A (zh) | 2021-11-01 | 2021-11-01 | 一种诊断精神***症的***及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114246551A (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6434419B1 (en) * | 2000-06-26 | 2002-08-13 | Sam Technology, Inc. | Neurocognitive ability EEG measurement method and system |
CN102715903A (zh) * | 2012-07-09 | 2012-10-10 | 天津市人民医院 | 基于定量脑电图的脑电特征提取方法 |
US20130267866A1 (en) * | 2012-04-05 | 2013-10-10 | Sony Corporation | Electroencephalogram analysis apparatus, electroencephalogram analysis program, and electroencephalogram analysis method |
CN109480866A (zh) * | 2018-10-26 | 2019-03-19 | 首都医科大学附属北京安定医院 | 一种利用瞳孔测量ppi的方法 |
CN109671500A (zh) * | 2019-02-26 | 2019-04-23 | 上海交通大学 | 基于脑电图时域数据的精神***症辅助诊断分类方法 |
CN110063732A (zh) * | 2019-04-15 | 2019-07-30 | 北京航空航天大学 | 用于精神***症早期检测和风险预测*** |
CN110097930A (zh) * | 2019-04-04 | 2019-08-06 | 华南理工大学 | 精神***症医学数据处理方法、装置、***、服务器及介质 |
CN112259237A (zh) * | 2020-10-13 | 2021-01-22 | 阿呆科技(北京)有限公司 | 基于多情感刺激和多级分类模型的抑郁评估*** |
-
2021
- 2021-11-01 CN CN202111281024.5A patent/CN114246551A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6434419B1 (en) * | 2000-06-26 | 2002-08-13 | Sam Technology, Inc. | Neurocognitive ability EEG measurement method and system |
US20130267866A1 (en) * | 2012-04-05 | 2013-10-10 | Sony Corporation | Electroencephalogram analysis apparatus, electroencephalogram analysis program, and electroencephalogram analysis method |
CN102715903A (zh) * | 2012-07-09 | 2012-10-10 | 天津市人民医院 | 基于定量脑电图的脑电特征提取方法 |
CN109480866A (zh) * | 2018-10-26 | 2019-03-19 | 首都医科大学附属北京安定医院 | 一种利用瞳孔测量ppi的方法 |
CN109671500A (zh) * | 2019-02-26 | 2019-04-23 | 上海交通大学 | 基于脑电图时域数据的精神***症辅助诊断分类方法 |
CN110097930A (zh) * | 2019-04-04 | 2019-08-06 | 华南理工大学 | 精神***症医学数据处理方法、装置、***、服务器及介质 |
CN110063732A (zh) * | 2019-04-15 | 2019-07-30 | 北京航空航天大学 | 用于精神***症早期检测和风险预测*** |
CN112259237A (zh) * | 2020-10-13 | 2021-01-22 | 阿呆科技(北京)有限公司 | 基于多情感刺激和多级分类模型的抑郁评估*** |
Non-Patent Citations (3)
Title |
---|
LAL HUSSAIN: "Spatial Wavelet-Based Coherence and Coupling in EEG Signals With Eye Open and Closed During Resting State", 《 IEEE ACCESS》 * |
廖媛媛: "精神***症的脑电网络及其识别研究", 《中国优秀硕士学位论文全文数据库(医药卫生科技辑)》 * |
朱秀强: "缺陷型与非缺陷型精神***症患者的精神症状及认知功能损害的研究", 《中国处方药》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rogers et al. | Test-retest reliability of a single-channel, wireless EEG system | |
Russo et al. | Effects of background noise on cortical encoding of speech in autism spectrum disorders | |
Kujala et al. | The mismatch negativity as an index of temporal processing in audition | |
Hubl et al. | Competition for neuronal resources: how hallucinations make themselves heard | |
Hofmann et al. | The worried mind: autonomic and prefrontal activation during worrying. | |
Blumenthal et al. | Committee report: Guidelines for human startle eyeblink electromyographic studies | |
Hampton et al. | Non-linguistic auditory processing in stuttering: evidence from behavior and event-related brain potentials | |
JP4987338B2 (ja) | 認知能減退を予測するためのシステムおよび方法 | |
Clinard et al. | Aging alters the perception and physiological representation of frequency: evidence from human frequency-following response recordings | |
Bramon et al. | Mismatch negativity in schizophrenia: a family study | |
Shepherd et al. | Electrophysiological approaches to noise sensitivity | |
Rundle et al. | Contagious yawning and psychopathy | |
Scanlon et al. | Your brain on bikes: P3, MMN/N2b, and baseline noise while pedaling a stationary bike | |
Doan et al. | Predicting dementia with prefrontal electroencephalography and event-related potential | |
KR20130050817A (ko) | 뉴로-퍼지 네트워크 기반 심박 변이도를 이용한 우울증 진단 방법 | |
Berbano et al. | Classification of stress into emotional, mental, physical and no stress using electroencephalogram signal analysis | |
Teplan et al. | Spectral EEG features of a short psycho-physiological relaxation | |
Adochiei et al. | Complex Embedded System for Stress Quantification | |
RU2314028C1 (ru) | Способ диагностики и коррекции психоэмоционального состояния "нейроинфография" | |
Kotani et al. | The effect of stimulus discriminability on stimulus-preceding negativities prior to instructive and feedback stimuli | |
Buján et al. | Cortical auditory evoked potentials in mild cognitive impairment: Evidence from a temporal‐spatial principal component analysis | |
Sassi et al. | Reprint of: stuttering treatment control using P300 event-related potentials | |
Miller et al. | Auditory sensory gating predicts acceptable noise level | |
Guðmundsdóttir | Improving players' control over the NeuroSky brain-computer interface | |
Fan et al. | Acoustic change complex evoked by horizontal sound location change in young adults with normal hearing |
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 |