脑电图动态中的拓扑递归：一种关于睡眠分期组织的从觉醒到快速眼动睡眠分级的补充性描述符

《Neuroscience Informatics》：Topological Recurrence in EEG Dynamics: A Complementary, Wake-to-REM-Graded Descriptor of Sleep-Stage Organisation

　　

背景：睡眠分期脑电图（EEG）分析传统上依赖频谱与复杂性测度。几何组织是否提供补充信息尚未被检验。本研究检验源自持续同调（persistent homology）的拓扑递归描述符是否标定一个在统计上区别于频谱与复杂性测度的睡眠分期EEG组织维度。新方法：K0为单通道额中线（Fpz–Cz）EEG时延嵌入（m=10，τ=2）持续同调的被试内z分数汇总，计算于Sleep-EDF Expanded公开多导睡眠库的100名受试者197夜数据。混合效应模型、二项广义线性模型（GLM）与留一被试分类将K0与频带功率、频谱熵、排列熵及Lempel–Ziv复杂性基准比对。结果：K0跨睡眠分期系统变化，与频谱内容统计独立（与对数频带功率Pearson r为?0.49至?0.25；快速眼动（REM）–觉醒方向相反），其构念效度由独立队列复现、生命周期两倍衰减及被试内苯二氮?不变性支持。REM–觉醒分离大（+1.01 SD；p=2.62×10?182）且跨嵌入保持。与现有方法比较：二项GLM确认K0携带统计独立于频带功率的信息（REM对主动眼动觉醒ΔAIC=?744.6，OR=5.04）。留一被试ROC下面积（AUC）为0.79（K0）与0.91（频带功率）；联合增益小（ΔAUC≈+0.004，配对bootstrap p<10?3）。结论：K0是睡眠分期EEG组织中拓扑有别、频谱非冗余的描述符。

研究背景与动机

睡眠分期脑电图（EEG）长期依赖频谱功率、振荡特征与统计描述符区分觉醒、非快速眼动（NREM）与快速眼动（REM）睡眠。深NREM以高幅慢振荡为特征，REM呈低幅宽频活动，与安静觉醒相似，导致常规单通道EEG难以区分REM与觉醒。神经生理学证据表明REM伴丘脑皮质与神经调制特异配置（胆碱能升高、单胺类下降），提示脑状态差异可能不在频谱内容而在底层动力学组织。时延嵌入（time-delay embedding）可重构标量序列的状态空间，持续同调（persistent homology，TDA子框架）以尺度不变方式量化其中的循环（loop）状递归结构（H₁同调特征）。该研究在《Neuroscience Informatics》发表，旨在以构念效度框架检验单通道EEG时延嵌入的H₁总持续性（total persistence）归一化指标K₀能否作为睡眠分期组织的补充描述符。

主要技术方法

研究人员采用公开Sleep-EDF Expanded多导睡眠库，纳入100名受试者197夜（Sleep-Cassette亚队153夜、Sleep-Telemetry亚队44夜含替马西泮交叉设计）。单通道Fpz–Cz经0.5–40 Hz零相位Butterworth滤波、降采样50 Hz、30秒分帧；每帧做m=10、τ=2样本时延嵌入，嵌入云标准化后二折子采样得25 Hz点云。以Ripser计算Vietoris–Rips filtration的H₁持续图，取有限死亡特征的生命期和（L₁≡TP₁），按夜内z分数得K₀。对照测度包括Welch频带功率（δ/θ/α/σ/β）、频谱熵、排列熵（m=5）、Lempel–Ziv复杂性。推断用线性混合效应模型（被试随机截距、觉醒参照）、二项GLM（REM vs 静息/主动眼动觉醒，含log频带功率+K₀）、留一被试逻辑回归/随机森林AUC、Cohen d_z与Holm校正；并做EOG回归校正、Pz–Oz多通道、独立队列复现、年龄分层、替马西泮药控、嵌入参数网格稳健性。

研究结果

3.1 原始与归一化拓扑递归（描述性）

原始L₁=TP₁跨期池化显示REM（93.53）与N1（91.68）高于觉醒（71.58）、N2/N3居中；因受增益尺度影响不可跨人比。夜内z分数得K₀后，各期均值：觉醒?0.67、N1 +0.28、N2 +0.03、N3 ?0.01、REM +0.37，REM与N1为高递归态，觉醒最低，呈现W→N1→REM抬升、N2/N3回落的形态。

3.2 统计对比与睡眠分期分化（推断）

混合效应整体分期效应χ²(4)=723.13（p=3.41×10^?155）。REM vs 觉醒配对d_z=1.92（每夜n=197时MDE d_z=0.20，观测达9.6倍），分离+1.01 SD（p=2.62×10^?182）；该正向分离在m∈{6,10,12}、τ∈{2,11}网格中保持。K₀与对数频带功率Pearson r ?0.49~?0.25，方向异于REM–觉醒频谱关系，证实非冗余。

2.11 全配对分期对比与单调性

十对分期对比Holm校正后均显著；按W→N1→N2→N3→REM赋秩的夜内Spearman ρ均值经单样本t检验>0（p极小），说明K₀沿觉醒—REM轴呈单调分级而非离散跳变。

2.13 增量二项GLM

REM vs 安静觉醒：加K₀较纯频带功率模型ΔAIC=?744.6，K₀系数OR=5.04（p≈0）；REM vs 主动眼动觉醒同方向稍弱。表明K₀在控制频谱后仍需统计显著，携独立信息。

2.14 留一被试分类

REM vs 觉醒AUC：K₀单独0.79，频带功率0.91，联合0.914（ΔAUC≈+0.004，配对bootstrap p<10^?3）；REM vs NREM、REM vs 其余期结论同向。K₀单用不及频谱但联合有微小显著提升。

2.10/2.16 独立队列复现与人口学分层

Sleep-Cassette与Sleep-Telemetry两亚队分别重拟合REM–觉醒对比，效应方向与量级一致。年龄分层（18–29至75+）显示K₀的REM–觉醒分离随龄衰减约两倍，与文献REM组织年龄相关改变吻合；性别分层无交互。Telemetry替马西泮vs安慰剂交叉中，REM–觉醒K₀对比在药内不变（β_REM,T不显著），说明K₀对被试内药理扰动不敏感、具状态纯粹性。

2.8/2.9 伪迹与通道控制

EOG水平通道自身算K₀负控制效应远小于EEG；Fpz–Cz经全记录OLS减EOG回归后REM–觉醒分离不消。Pz–Oz复算分离同向但略弱，表明递归结构以前中线为主但不只限局。

讨论与结论翻译

讨论指出，K₀不替代频谱分期而补其未捕维度：在REM与觉醒频谱近似时，以重构状态空间H₁总持续性标定“内部闭合递归”程度。其与被试内苯二氮?不变、随龄衰减、跨嵌入稳健、EOG不驱动，共同支持作为睡眠动力学组织合法描述符。结论原文意译：K₀是睡眠分期EEG组织中拓扑有别、频谱非冗余的递归描述符；它沿觉醒→REM轴分级变化，在REM相对觉醒呈约1 SD正向分离，该分离独立于频带功率、排列熵与Lempel–Ziv复杂性，并在独立队列、生命周期与药控下复现，可作为现有临床睡眠分期频谱/复杂性测度的补充构念。