基于微状态的轻度认知障碍脑电信号特征研究.pdf
《基于微状态的轻度认知障碍脑电信号特征研究.pdf》由会员分享,可在线阅读,更多相关《基于微状态的轻度认知障碍脑电信号特征研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 47 卷 第 4 期燕山大学学报Vol.47 No.42023 年 7 月Journal of Yanshan UniversityJuly 2023 文章编号:1007-791X(2023)04-0370-07基于微状态的轻度认知障碍脑电信号特征研究刘泽达1,2,崔 冬1,2,李小俚3,王 磊4,尹世敏4(1.燕山大学 信息科学与工程学院,河北 秦皇岛 066004;2.燕山大学 河北省信息传输与信号处理重点实验室,河北 秦皇岛 066004;3.北京师范大学 认知神经科学与学习国家重点实验室,北京 100875;4.中国人民解放军火箭军特色医学中心 神经内科,北京 100088)收稿日
2、期:2023-01-18 责任编辑:唐学庆基金项目:国家自然科学基金资助项目(62173291);河北省重点实验室项目(202250701010046)作者简介:刘泽达(1999-),女,河北保定人,硕士研究生,主要研究方向为脑电信号处理;通信作者:崔冬(1978-),女,吉林白城人,博士,教授,主要研究方向为神经信号建模与处理,Email:cuidong 。摘 要:静息态脑电微状态是一种探究认知活动时间动力学的有效工具。为了研究轻度认知障碍(MCI)患者中是否存在特定的微状态变化,本文将 15 名 MCI 患者和 10 名对照者的脑电图通过 k 均值聚类算法解析为四种典型的微状态,并分析其时
3、间参数在两组间的显著性差异以及与神经心理学量表得分间的相关性,又对微状态序列的复杂性、转换特点以及各状态间是否具有依赖性做了分析。研究发现,MCI 组在微状态 A 的持续时间以及 AD 和 BD 的转移概率上与对照组有着显著性差异。微状态 A 的持续时间与 AVLT2 分数成显著性正相关。MCI 组微状态序列的复杂度要高于对照组,且在微状态的转换过程中尚未发现某种特殊结构,未来时刻的状态与过去无关。关键词:轻度认知障碍;脑电信号处理;微状态;聚类中图分类号:TP391 文献标识码:A DOI:10.3969/j.issn.1007-791X.2023.04.0100 引言轻度认知障碍(Mild
4、 Cognitive Impairment,MCI)是指正常衰老和早期痴呆认知变化之间的过渡状态1-2。MCI 患者的显著特征为认知以及执行功能存在一定程度下降,例如不符合年龄规律的记忆力衰退,但依旧拥有进行日常基本生活的能力,65 岁及以上老人约有 20%患有 MCI,与同年龄和性别相匹配的老年人群相比,MCI 患者发展成阿尔茨海默症(Alzheimer Disease,AD)的转换率要高几十倍,确诊后一年内约 10%30%转化为痴呆,其中绝大多数转化为 AD3-4。因此,研究 MCI的发病机制,为 MCI 的临床诊断提供有效的客观指标,具有重要科学价值和临床意义。Lehmann 等人最早提
5、出了多通道静息态 EEG信号的 频段(813 Hz)可以解析为有限数量的不同准稳态,称为“微状态”5。通过聚类方法,微状态分析将脑电地形图解析为几类最主要的微状态原型,进一步表达脑电时间序列的时空特征6。近几年,微状态分析已经开始应用于分析 MCI7、AD8、精神分裂症9以及自闭症10等多种脑功能疾病,在这些患者中 EEG 微状态均有显著性的变化。由于空间地形可以定义于任意时间点,独立于之前或后续的地形,具有毫秒级分辨率,因此微状态更适合检测大规模神经认知网络中的快速动态活动,在一定程度上弥补了 EEG 在空间分辨率和分析大规模脑网络异常方面的不足11。本文利用脑电微状态分析,研究与对照组相比
6、,轻度认知障碍患者中是否存在特定的微状态变化。将 EEG 数据每个时刻处的脑电地形图根据空间相似度指标,聚类得到 4 个最具代表性的微状态原型地图,通过此微状态原型重新表达脑电时间序列特征,统计分析了 MCI 组和健康对照组在平均持续时间、平均发生频率、覆盖率以及转第 4 期刘泽达 等 基于微状态的轻度认知障碍脑电信号特征研究371 移概率 4 个时间参数上的差异,利用皮尔森相关分析探究了这些时间参数与神经心理学量表得分之间的关联,利用样本熵算法计算了两组微状态序列的复杂性,最后对两组微状态序列做了随机化测试及马氏性检验,探究微状态间的转换特点。1 方法1.1 受试者和诊断标准本文被试来自于中
7、国人民解放军火箭军特色医疗中心神经内科。所有的受试人员具有一定视觉和听觉分辨力,可以独立完成神经心理学检测内容,包括简易智力状态检查量表(Mini-Mental State Examination,MMSE)、蒙特利尔认知评估(Montreal Cognitive Assessment,MoCA)、听觉言语学习测验(Auditory Verbal Learning Test,AVLT)的即刻回忆 AVLT1、30 min 延迟回忆 AVLT2、长时延迟再认 AVLT3、Boston 命名测试、语言流畅性测试、轨迹测试 A 和 B 以及数字广度测验和日常生活量表(Frequently Asked
8、 Questions,FAQ)。并排除了:1)具有额叶痴呆,锥体外系病变,脑外伤,脑肿瘤,癫痫,神经性衰退性疾病等;2)严重酗酒和抑郁症患者;3)存在临床大血管并发症(如心绞痛、心肌梗死病史、有临床症状的脑梗死及外周血管栓塞);4)入组前 1 个月内使用过具有镇静、催眠效果的药物及影响中枢神经系统的药物。将所有受试者按照划分标准进行分类12,分为 MCI 组15 人(男性 6 名,女性 9 名)和认知功能正常的对照组 10 人(男性 6 名,女性 4 名)。使用 SPSS 26.0 软件对 MCI 组和对照组神经心理学量表测试得分与人口特征数据进行独立样本 T 检验分析,以均值标准差的方式呈现
9、结果,如表 1 所示,表示 p0.05,表示 p0.01,表示 p0.001,p0.05 认为存在统计学意义的显著差异。结果显示人口统计学数据在两组之间没有显著差异。所有受试者在年龄、性别和教育水平方面均匹配。神经心理学测试显示 MCI 组和对照组的 MoCA 分数和 AVLT2 分数具有显著差异性(p0.05),表明 MCI 组在认知功能尤其是记忆功能方面有 显 著 的 损 伤。而 在 MMSE 分 数、AVLT1、AVLT3、Boston 命名测试、语言流畅性测试、轨迹测试、数字广度测验、FAQ 等方面两组之间没有出现明显的差别。表 1 神经心理学量表得分Tab.1 Table of ne
10、uropsychological scale score因素MCI 组对照组p 值年龄69.877.0472.507.840.216受教育程度14.002.7013.302.500.520MMSE28.471.8128.900.740.978MoCA22.801.6626.700.820.000AVLT16.232.106.941.670.378AVLT24.734.038.403.170.024AVLT310.933.8613.001.490.177Boston19.271.0319.700.480.461语言18.273.7718.004.060.868轨迹测试 A63.9221.3160
11、.9018.100.723轨迹测试 B109.6934.77102.3040.360.642数字广度测验11.803.3813.702.710.062FAQ2.143.660.701.380.5851.2 数据采集和预处理脑电信号采用 EGI 的 GES300 型 128 通道脑电记录仪记录,试者处于闭眼静息状态进行静息EEG 数据采集,以双侧乳突为参考电极,阻抗小于10 k,采样频率 1 000 Hz,时长 5 min。使用EEGLAB 对数据做预处理,包括通过凹陷滤波去除 50 Hz 工频干扰;进行 130 Hz 的带通滤波消除低频漂移和高频噪声;平均重参考电极。1.3 微状态分析首先计算
12、全局场功率(Global Field Power,GFP)来表示大脑电场的瞬时强度5,公式如下:g=(Ki=1(Vi(t)-Vmean(t)2)/K,(1)其中,Vi(t)为第 i 路电极信号强度,K 为电极数目,Vmean为多路信号强度的平均值。脑电地形图在 GFP 峰值附近信噪比最高且最稳定13,因此选取 GFP 峰值即局部极大值处所对应的脑电数据作为时间样本点。设置最小峰距为10 ms,以便更好地区分相邻峰值;每个受试者进入分割的 GFP 峰值数量为 1 000,保证聚类结果能够收敛的前提下,使用较少的样本减少计算量;剔除超过 GFP 标准差 1 倍的峰值,防止极端 GFP 处的高振幅非
13、神经源伪影带来的不利影响。通过 k372 燕山大学学报2023均值聚类算法,预设 4 种微观状态,随机初始化的次数为 50,最大迭代次数为 1 000,收敛阈值设为10-6。以空间相似性为指标,全部样本地图分成 4个类群,对群内所有样本地图取平均,得到 4 种典型的微状态地图,即 A、B、C、D 四个状态。把每个被试的所有脑电地形图,依旧根据空间相似性指标替换成与其最相似的微状态原型地图,得到 A、B、C、D 四个原型交替出现的微状态时间序列。为了避免噪声短段以及极性变化所带来的虚假影响,去除长度小于 30 ms 的微状态片段。1.4 微状态时间参数统计分析提取 MCI 组和对照组微状态的四类
14、时间参数:1)平均持续时间:各微状态出现并保持短暂稳定的平均时间长度;2)平均发生频率:各微状态在记录期间内每秒发生的平均次数;3)覆盖率:各个微状态存在于总的时间中所占的比例;4)转移率:各个微状态之间的转换概率。利用独立样本 T 检验统计分析 MCI 组和对照组的微状态在平均持续时间、平均发生频率、覆盖率和转移率之间是否具有显著性差异。利用皮尔森相关分析探究神经心理学量表得分与微状态时间参数之间的相关性差异。1.5 微状态序列复杂性分析利用样本熵算法对两组被试微状态时间序列的复杂性进行分析。研究两组微状态序列的熵值大小随嵌入维数 m 的变化趋势,其中,嵌入维数 m设为 1 到 11,阈值
15、r 设为 0.2 倍的标准差。1.6 微状态转换结构与依赖性分析为了探究微状态序列中是否存在某特定转移结构以及各状态间是否具有依赖性,分别对两组微状态的转移概率进行了随机化测试和马氏性检验。随机化测试先假设从一个微状态到下一个微状态的转换是随机发生的,观测转移概率将与微状态的相对发生率成比例,可从微状态序列提取了观测转移概率,又基于各微状态的相对发生率计算出预期转移概率,并使用二者的卡方距离 D0来评估差异,将观测转移概率和预期转移概率随机进行 5 000 次重新分配,得到 5 000 个卡方距离的经验分布,若大于 D0的有 q 个,则估计的 h 值为 q/5 000,当 h2(s-1)2),
16、可认为微状态序列具备马氏性,反之,则不具备马氏性。2 结果MCI 组和对照组的 A、B、C、D 四个微状态原型地图如图 1 所示。不同颜色及颜色深浅代表采集电极所记录的电压正负及其大小。在地形结构上,A 类微状态地形图的特点为右前-左后;B 类微状态地形图为左前-右后;C 类微状态地形图为前-后;D 类微状态地形图为前-中末。四种微状态原型图彼此之间有很大区别,在以往的研究中具有高度的可重复性。图 1 MCI 组和对照组微状态原型地图Fig.1 Map of microstate prototype of MCI group and control group2.1 时间参数分析结果微状态时间
17、参数的统计分析结果如表 2 所示,表中包含了 MCI 组与对照组两组微状态的平均持续时间、发生频率和覆盖率,并以均值标准差的方式呈现,可以看出,MCI 组微状态 A、B 和 C的平均持续时间小于对照组,其中微状态 A 具有显著差异;而 MCI 组微状态 D 的平均持续时间高于对照组,但没有显著差异。MCI 组微状态 A 和D 的发生频率增高,B 和 C 的发生频率降低,均未有显著差异。MCI 组微状态 A、B 和 C 的覆盖率降第 4 期刘泽达 等 基于微状态的轻度认知障碍脑电信号特征研究373 低,微状态 D 的覆盖率增高,也均未有显著差异。表 2 持续时间、发生频率以及覆盖率统计分析Tab
18、.2 Table of statistical analysis of duration,occurrence and coverage微状态MCI 组对照组p 值持续时间/msA69.585.2375.368.410.044B72.375.2075.625.350.144C80.5911.2883.5810.390.367D88.0325.0674.6715.350.071发生频率/(次/s)A2.890.512.800.930.723B3.140.373.240.360.338C3.310.683.820.510.057D3.240.482.870.510.080覆盖率A0.200.050
19、.210.090.709B0.230.040.250.030.233C0.270.080.320.080.129D0.300.110.220.070.063 表 3 为 MCI 组与对照组的微状态转移概率统计分析结果,可以看出,MCI 组微状态 A、B 和 C 转移到 D 的概率、C 和 D 转移到 A 的概率以及 D 到B 的转移概率增高,其中 AD 和 BD 的转移概率具有显著差异。MCI 组微状态 A、B 和 D 转移到 C 的概率、A 和 C 转移到 B 的概率以及 B 到 A的转移概率降低,但尚未具有显著性。表 3 转移概率统计分析Tab.3 Table of statistical
20、 analysis of transfer probability转移方向转移概率MCI对照组p 值AB0.3150.0410.3420.0560.184AC0.3470.0930.4090.0620.077AD0.3380.0930.2520.0500.007BA0.2840.0510.3050.1200.599BC0.3410.0840.4280.0990.080BD0.3760.1070.2670.0820.012CA0.3010.0710.2960.1080.904CB0.3250.0390.3630.0560.056CD0.3740.0910.3410.1090.412DA0.293
21、0.0390.2470.0930.238DB0.3360.0710.3080.0490.428DC0.3710.0720.4440.1210.1062.2 相关性分析结果图 2 是微状态时间参数与神经心理学量表得分之间的显著性相关图。结果显示,微状态 A 的持续时间与 AVLT2 得分成显著性正相关。该结果表明微状态的持续时间与神经心理学存在着一定关联,认知障碍程度越大,微状态 A 的持续时间则越短。图 2 神经心理学量表得分与微状态时间参数相关性Fig.2 Correlation between neuropsychological scale score andmicrostate tim
22、e parameters2.3 复杂性分析结果图 3 是两组微状态序列的样本熵值随嵌入维数 m 的变化趋势图,实线表示 MCI 组样本熵的均值,虚线表示对照组样本熵的均值,图中的阴影部分表示标准差。MCI 组在不同嵌入维数下的熵值都高于对照组,说明 MCI 组微状态序列的复杂度要高于对照组,MCI 患者在处理认知相关的任务时可能需要调用更多的认知资源。对照组中,样本熵值随着嵌入维数 m 的增加而降低,表明微状态的转换序列在多个转换中是自相似的;而在MCI 组中,样本熵值没有随着嵌入维数的增加而降低,微状态序列缺乏自相似性。样本熵的计算使我们能够量化微状态转变的可预测性。随着嵌入维数 m 的增加
23、,微状态序列的样本熵在对照组中减少,但在 MCI 患者中没有,这表明来自对照组被试的 EEG 数据比来自患者的数据更加可预测。2.4 转换结构与依赖性分析结果随机化测试得出对照组的 h 值为 0.998 6,MCI 组的 h 值为 0.985 0。即接受原假设,微状态序列中不同状态的转换是较为随机的,尚未发现某种特定转移结构存在。马氏性检验中,经查表得,20.05(9)=16.919,374 燕山大学学报2023在 MCI 组中2=579.99820.05(9),在对照组中2=599.71220.05(9),具备马氏性,说明在微状态序列中,未来的状态与过去的状态互相独立,不具有依赖性。图 3
24、样本熵随嵌入维数变化图Fig.3 Change of sample entropy with the increase of embedding dimension3 讨论本文对 MCI 患者与对照组相比,是否存在特定的微状态变化进行了分析,证实了 MCI 患者和对照组在微状态 A 的平均持续时间以及 AD、BD 的转移概率上具有显著性差异;并报告了神经心理学量表得分与微状态时间参数之间的相关性,得知微状态 A 的持续时间与 AVLT2 得分成显著性正相关;探究了两组微状态序列的复杂度差异,发现 MCI 患者的微状态序列复杂度更高,序列缺乏自相似性;最后检验出在状态转换过程中不存在特定转移结构
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 状态 轻度 认知 障碍 电信号 特征 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。