信息速览
手术台上的麻醉医生盯着一块监护屏,屏幕上真正反映大脑的只有一个数字:双频指数、状态熵,或者某个厂商自己起的名字。这个数字背后是一整条被压缩掉的脑电曲线,而那条曲线里藏着 GABA 能麻醉药在什么时候把大脑从一种无意识推进另一种无意识的答案。更要命的是,这两种无意识状态各自都被研究得很透,唯独中间那段路没人画出来:慢波与 α 振荡怎么一步步让位给爆发抑制,剂量降下来之后又怎么走回去。PNAS 上那篇论文找到了两个能跟踪这段路的指标,分析了四十个人,公开到 PhysioNet 上的只有四个人、四段会话。
少的三十六个人不是失误。这份数据只放出论文图 3 到图 6 对应的那四条记录,每条配一张全程频谱图、两条二阶 Butterworth 滤波后的波形、一段给药记录,外加一个写明会话总时长与分析窗起止时刻的小文件。它不是用来训练模型的,它是让你把论文里的两条调制曲线在原始维度上重新算一遍的。四例、两种药、四条时间轴,数量少得让人第一眼就皱眉,稀缺性却也正在这里:一个健康志愿者在计算机控制的丙泊酚输注下按五个台阶先升后降地走完全程,这种给药形状的数据,公开场合几乎找不到第二份。先看清它能给你什么、给不了什么,再往下读。一屏速览见下表。
| 维度 | 内容 |
|---|---|
| 官方全名 | Electroencephalogram dynamics during unconsciousness mediated by GABAergic-anesthetics |
| 发布方 | PhysioNet(由 MIT 计算生理学实验室维护托管) |
| 贡献者 | Elie Adam、Ohyoon Kwon、Karla Alejandra Montejo(MIT);Emery Brown(MIT Picower Institute for Learning and Memory) |
| 首次发布 | 2023 年 7 月 21 日 |
| 当前版本 | 1.0.0,版本日期 2023 年 7 月 21 日;官方 Release Notes 写明这是该数据集的首发版本 |
| 版本 DOI | 10.13026/dx44-kw30(v1.0.0);最新版本 DOI 10.13026/7ygq-3f91 |
| 配套论文 | Adam、Kwon、Montejo、Brown(2023)PNAS 120(30) e2300058120,PMID 37467269 |
| 公开受试者数 | 4 名(1 名健康志愿者、3 名手术患者),分别对应论文图 3、图 4、图 5、图 6 |
| 论文分析队列 | 4 组各 10 例,合计 40 例;公开的是每组各 1 例 |
| 目录组织 | 4 个文件夹:Propofol_Volunteer、Propofol_OR_strongAlpha、Propofol_OR_weakAlpha、Sevoflurane_OR |
| 每文件夹文件 | 6 个 CSV(三个频谱图文件、两条波形、一个给药文件)加 duration_start_end 时窗文件 |
| 分析通道 | 志愿者 Fp1(64 通道采集,分析只取一路);手术患者 Fp2(SedLine 额部监护) |
| 采集参数 | 志愿者 5000 Hz、带宽 0.016–1000 Hz、0.5 µV 每最低有效位;手术丙泊酚 178 Hz;手术七氟烷 250 Hz、前置放大器带宽 0.5–92 Hz、16 位 29 nV |
| 共享波形采样率 | 官方 Data Description 段写 α 波与慢波为 250 Hz 或 178 Hz,与 Methods 段的 5000 Hz 存在张力(见 §2.3) |
| 频带定义 | 慢 δ 0.3–4 Hz;α 8–14 Hz;均为二阶 Butterworth 带通结果 |
| 文件格式 | 全部 CSV |
| 许可 | PhysioNet Contributor Review Health Data License 1.5.0,配套 Health Data Use Agreement 1.5.0 |
| 访问门槛 | 需 PhysioNet 认证账号、签署 DUA、完成 CITI 的 Data or Specimens Only Research 培训,且个别研究需贡献者审阅 |
| 典型用途 | 复现 α 调制指数与慢调制指数、研究慢波-α 与爆发抑制之间的过渡动力学、开发麻醉深度监测新指标 |
| 官方明写的局限 | 手术会话中合并用药会改变动力学;手术操作会造成脑电瞬时改变;两条波形被截断在调制指标计算区间内;频谱图覆盖清醒、诱导、苏醒等未被分析的时段 |
§1 它要解决的具体麻烦
1.1 一个被压缩成标量的器官
麻醉监护里存在一个长期的不对称:心脏给一条曲线,大脑给一个数字。心电、血压、血氧都以波形或趋势上屏,唯独脑电先被厂商的算法吃进去,再吐出一个 0 到 100 的标量。压缩的代价在教科书里被反复念叨,最典型的是同一个指数值可以对应完全不同的频谱形态,尤其在老年患者身上,指数正常而脑电已经接近平线的情形并不罕见。
为什么不能直接用频谱去看?可以,但要看懂它,得先知道这两个状态之间发生了什么。GABA 能麻醉药维持的无意识状态有两个主角:0.3–4 Hz 的慢 δ 振荡,与 8–14 Hz 的 α 振荡,两者合称慢波-α 复合体,英文缩写 SDA。剂量再往上抬,两者会解体为爆发抑制,也就是等电位平线与爆发段交替出现的深度失活。临床上的通行建议是把剂量压在不出现爆发抑制的位置,因为后者与术后认知障碍的关联在老年人群里被反复报道。
问题就卡在这里:慢波-α 是一回事,爆发抑制是另一回事,两者之间的转换过程没有可用的刻画。官方页面的 Background 段把这一点写得很直白:这篇工作报告一种刻画两个无意识状态之间过渡的生理调制,而这种调制同时出现在 α 频带与慢 δ 频带上。这份 GABA 麻醉脑电调制数据集的全部价值都建立在回答这一句上。
1.2 两条被分开研究的曲线
把两者分开看是有历史原因的。慢波-α 主要被当作神经生理过程:α 是丘脑与前额皮层之间的高度同步节律,慢 δ 则是皮层大面积进入上下状态(up-down state)的外在表现。爆发抑制则被建模为神经代谢过程,ATP 供给跟不上,神经元维持不住膜电位,于是平线插进了原本连续的振荡之间。一个讲环路,一个讲能量,研究群体与方法论天然分成两拨。
分开研究的后果很具体:医生在床旁拿到的是一个「要么还没到、要么已经过了」的二元判断,中间那段连续变化没有任何可操作的坐标轴。论文给出的答案是两个调制过程,α 波幅度调制与慢波频率调制,并由此构造 α 调制指数与慢调制指数,用它们连续地追踪从慢波-α 到爆发抑制再回到慢波-α 的整条路径。
这份数据放出的正是构造这两个指数所需的原料:一条被滤波出来的 α 波形、一条慢波波形,以及一整段会话的多锥谱图。想验证或者推翻这两个指数,不需要找作者要原始文件,也不需要那四十个人。四个人,够画图,也够复现方法。
1.3 三条常见的错误期待
第一条,以为它能拿来训练分类器。 四名受试者、四个会话、每条会话一个连贯的时间序列,这不是一个有样本量的建模场景。它连最基本的训练集与测试集划分都无从下手,因为同一个人不能同时在两侧,而这里根本没有多余的人去分组(§4.1)。
第二条,以为它是原始脑电。 它不是。官方分享的是三条派生结果:频谱图的三个坐标文件、两条滤波后的窄带波形,以及给药记录。原始 64 通道信号、原始 SedLine 波形都不在里面。想换一种参考电极、想试别的谱估计方法、想看 Fp1 以外的第二个额头通道,做不到。
第三条,以为拿文件不需要门槛。 它是 Contributor Review 级别:注册、签署数据使用协议、完成 CITI 培训是三道路径上的固定关卡,官方 Access 段还额外写明使用者的个别研究须经贡献者审阅。这不是下载一个开放心电库那种随手可得的东西,从决定用到真正拿到第一份 CSV,中间隔着人(§5.1)。
1.4 它能给到的与给不到的
| 你的需求 | 它给的 | 够不够 |
|---|---|---|
| 复现论文的两个调制指数 | α 波与慢波波形,滤波口径在论文中写明 | 够,这正是它被发布的目的 |
| 研究剂量变化与调制深度的关系 | 志愿者的效应室浓度曲线;手术患者的输注速率与推注量 | 志愿者的形状干净,手术患者的粗糙但真实 |
| 比较丙泊酚与七氟烷的差别 | 三条丙泊酚记录、一条七氟烷记录 | 只够定性对照,不够统计推断 |
| 训练一个端到端的麻醉深度模型 | 无 | 不够,只有四个人 |
| 换一种谱估计方法重算频谱图 | 无原始波形 | 不够 |
| 做多通道空间分析(相干、源定位) | 每例只给一路额部通道的派生结果 | 不够 |
| 拟合一例药代动力学模型的参数 | 志愿者给的是效应室浓度而非血浆浓度 | 部分够,取决于模型输入要求 |
| 研究 α 波为何在部分患者身上缺失 | 专门收了一名 70 岁的弱 α 手术患者 | 够做单例剖析,不够做归因 |
最后一行是这份数据最容易被忽略的一格。弱 α 被单独做成一个文件夹,等于承认「不是每个人都按教科书长那样」这件事本身值得研究,而绝大多数公开麻醉脑电数据在发布时就把这类患者剔掉了。
1.5 你的处境与它能帮到的事
| 你的处境 | 它能帮你什么 | 下一步该做什么 |
|---|---|---|
| 你在审一篇用调制指数做麻醉深度监测的论文,想自己复算 | 原料与论文口径一致,幅度调制与频率调制都能重取 | 先读 §2.3 的单条解剖,再按 §5.4 把几列读到内存里算一遍 |
| 你在做麻药闭环给药控制,需要一个可参考的目标轨迹 | 志愿者会话给了先升五个台阶再关泵的完整剂量形状 | 把这条曲线当工况模板,注意只有一例 |
| 你在写综述,需要一个「过渡过程也有生理学」的实例 | 四张图覆盖四条典型路径:志愿者、强 α、弱 α、七氟烷 | 引用论文图 3 至图 6 与这份 DOI,不必再找作者 |
| 你想比较计算机控制输注与人工滴定的差别 | 同一套频带口径下有一条靶控记录与两条人工滴定记录 | 只看形状不看绝对强度,因为给药记录形式不同 |
| 你在教麻醉监测,需要一个能操作的脑电示例 | CSV 直接可读,无需专用软件 | 用 §9.2 的流程先跑通志愿者那一份 |
| 你想在爆发抑制检测上做基准测试 | 频谱图上能看出高剂量段,但没有逐段标签 | 不够,这不是标注数据集(§4.4) |
1.6 一句话判断与不适合谁
它不是通常意义上的数据集,而是一篇论文的可运行附录:只有四例,但每一例都带着完整的剂量轴、频谱轴与两条调制波形,够你判真伪,不够你做统计。
四条边界,写清楚之后就该掉头:
不适合当作机器学习素材。 四例的数量级只够做单例分析与方法验证,任何交叉验证都无从设计。
不适合当作临床人群的代表样本。 志愿者是 18–35 岁的健康年轻人,三名手术患者分别 44、70、54 岁,四例一共只覆盖有限几档剂量与两种给药方式。
不适合当作多药相互作用的证据来源。 官方 Usage Notes 明确提醒,手术会话中合并使用了镇痛药等其他药物,这些药会改变主要麻醉药的效果与动力学。
不适合当作干净 α 节律的数据源。 与之相反,发布者刻意保留了一个 α 波消散的 70 岁患者;把这一例剔掉再谈规律,才是对这份数据的误用。
§2 数据构成与规格
2.1 四个文件夹与七类文件
官方 Data Description 段把结构交代得很规整:每名受试者一个文件夹,四个文件夹分别是 Propofol_Volunteer、Propofol_OR_strongAlpha、Propofol_OR_weakAlpha、Sevoflurane_OR。每个文件夹的内部构成如下表。
| 文件名模式 | 维度 | 物理量 | 单位 | 说明 |
|---|---|---|---|---|
| t_spec_*.csv | 一维数组 | 频谱图时间轴 | s | 频谱图的横轴 |
| f_spec_*.csv | 一维数组 | 频谱图频率轴 | Hz | 频谱图的纵轴 |
| val_spec_*.csv | 二维数组 | 频谱图取值 | dB | 行列归属需与两个轴对形状确认 |
| alphaWave_*.csv | 一维数组 | α 波段波形 | µV | 8–14 Hz 带通结果,已截断到分析窗 |
| slowWave_*.csv | 一维数组 | 慢波段波形 | µV | 0.3–4 Hz 带通结果,已截断到分析窗 |
| propofolConcentration_*.csv | 两列 | 丙泊酚效应室浓度 | s / µg 每毫升 | 仅 Propofol_Volunteer 有 |
| propofolInfusionRate_*.csv | 两列或三列 | 丙泊酚输注方案 | s / mg 每分钟 / mg | 三个手术文件夹都有,列数不同 |
| SevofluraneConcentration_*.csv | 两列 | 七氟烷呼气末浓度 | s / % | 仅 Sevoflurane_OR 有 |
| duration_start_end_*.csv | 三项 | 时长与起止时刻 | s | 依次为会话总时长、α 波起止、慢波起止 |
把这张表和 §1.3 的第二条连起来读,会出现一个关键判断:文件总数在二十几个,比大多数人以为的「四份原始脑电」要少,但也比大多数人以为的「只有波形」要多。频谱图提供全程视野,波形提供可算指数的窄带序列,给药文件提供那条自变量轴。
2.2 字段字典
CSV 没有统一的列名约定,官方是用「第几列是什么」的方式说明的。下面逐条列出,这是唯一能对齐语义的地方。
| 文件 | 第 1 列 | 第 2 列 | 第 3 列 |
|---|---|---|---|
| t_spec_*.csv | 频谱图时间点,单位 s | 无 | 无 |
| f_spec_*.csv | 频谱图频率点,单位 Hz | 无 | 无 |
| val_spec_*.csv | 频谱图 dB 值矩阵 | 无 | 无 |
| alphaWave_*.csv | α 滤波波形,单位 µV | 无 | 无 |
| slowWave_*.csv | 慢波滤波波形,单位 µV | 无 | 无 |
| propofolConcentration_*.csv | 时间,s | 效应室浓度,µg 每毫升 | 无 |
| propofolInfusionRate_*(手术丙泊酚) | 时间,s | 持续输注速率,mg 每分钟 | 单次推注量,mg(以 1 分钟内的 mg 每分钟表示) |
| propofolInfusionRate_*(七氟烷那例) | 时间,s | 单次推注量,mg(以 1 分钟内的 mg 每分钟表示) | 无 |
| SevofluraneConcentration_*.csv | 时间,s | 呼气末浓度,% | 无 |
| duration_start_end_*.csv | 会话总时长,s | α 波相对频谱图的起止时刻 | 慢波相对频谱图的起止时刻 |
注意第三行与倒数第四行的差别:同样是 propofolInfusionRate 这一族文件,在七氟烷那个文件夹里只有两列,且第 2 列语义变成了推注量,因为那一例以七氟烷维持为主,丙泊酚只用于推注。文件名相同而列语义不同,这是最容易写出静默错误的地方(§5.6)。
2.3 一条记录的完整解剖
以 Propofol_Volunteer 为例走一遍。这个文件夹对应论文图 3,是一名 18–35 岁的健康志愿者,在计算机控制的靶控输注下先接受递增的五个目标效应室浓度,随后关泵让其自然下降。官方 Methods 段还记下了一个与数据表强相关的细节:受试者被要求全程闭眼,以避免眨眼伪迹进入脑电。
采集端是全头 64 通道 BrainVision MRI Plus 系统,采样率 5000 Hz,带宽 0.016–1000 Hz,分辨率 0.5 µV 每最低有效位。但官方 Data Description 段却写着,α 波与慢波的采样频率为 250 Hz 或 178 Hz,如 Methods 段所述。这两句话之间存在一处官方未解释的张力:共享出来的波形,显然不可能是 5000 Hz 的原始分辨率。
| 解读方式 | 支持它的官方文字 | 反对方的理由 |
|---|---|---|
| 志愿者波形被降采样到 250 Hz | Data Description 段明确写 250 或 178 Hz | 官方没有写明降采样这一步 |
| 250 与 178 分别对应七氟烷例与手术丙泊酚例,志愿者另有处理 | Methods 段只出现 5000、250、178 三个数 | 那么 Data Description 的「或」就没有涵盖志愿者 |
| Methods 段的 5000 Hz 只描述采集端,与派生文件无关 | 两句都成立 | 缺少让两者自洽的中间说明 |
三种解读在官方文字里都站得住,也都缺一环(§10 列为未核验项)。落地时的防御写法很简单:从 duration_start_end 的第 1 项读出会话总时长,用它和波形长度反推实际采样率,再据此建立时间轴,而不是照抄任何一个写在网页上的数字。
2.4 映射:从文件夹到论文章节
| 文件夹 | 论文出处 | 受试者 | 年龄 | 药物与给药方式 | 分析通道 | 给药物件 |
|---|---|---|---|---|---|---|
| Propofol_Volunteer | 图 3 | 健康志愿者 | 18–35 岁 | 丙泊酚,计算机控制输注,先升后降 | Fp1 | 效应室浓度 |
| Propofol_OR_strongAlpha | 图 4 | 手术患者 | 44 岁 | 丙泊酚,人工滴定 | Fp2 | 持续输注速率加推注量 |
| Propofol_OR_weakAlpha | 图 5 | 手术患者 | 70 岁 | 丙泊酚,人工滴定,α 波消散 | Fp2 | 持续输注速率加推注量 |
| Sevoflurane_OR | 图 6 | 手术患者 | 54 岁 | 七氟烷,人工滴定 | Fp2 | 呼气末浓度加丙泊酚推注量 |
这张表同时暴露三个边界条件。第一,年龄只有四个点,且志愿者与三名手术患者的年龄段不重合。第二,通道不统一,Fp1 与 Fp2 之间不可直接比较空间信息。第三,给药记录的语义不统一,人工滴定只有速率与推注,计算机控制输注才有效应室浓度。
2.5 与同类资源的字段差异
把这份数据和库内另外几份麻醉相关的记录放在一起看,差别不在谁更好,而在谁保留了哪一维。
| 对比项 | 本数据集 | eeg-power-anesthesia(rid 1317) | propofol-anesthesia-dynamics(rid 1316) | jsmf-access(rid 1006) |
|---|---|---|---|---|
| 主体形态 | 全程频谱图加两条窄带波形 | 两秒窗的多锥功率谱 | 逐拍心率变异加皮电事件 | 术中额部脑电的衍生指数特征 |
| 是否含原始波形 | 否 | 否 | 否 | 否 |
| 样本规模 | 4 例 | 54 例 | 9 人 10 阶段 | 646 例 |
| 标签体系 | 无逐段标签 | 逐两秒意识标签加信号质量标记 | 意识消失与恢复时刻 | 术后谵妄等结局变量 |
| 剂量信息 | 有,逐例给出浓度或输注曲线 | 本页未见同等粒度 | 有,靶控给药分阶段 | 有,药物剂量与时点 |
| 分析通道 | 单通道,Fp1 或 Fp2 | 单通道,近似 Fp1 或 Fp2 | 以非脑电通道为主 | 额部四电极 |
| 典型任务 | 方法复现与动力学建模 | 意识状态分类 | 自主神经动态分析 | 结局预测 |
读这张表的正确顺序是从第二行开始:四份全都不给原始波形。这不是巧合,而是麻醉脑电公开数据的常态,能连续采集脑电的场合本来就不是以研究为目的的门诊病房,愿意公开原始波形的机构极少。差别只在于把原始信号压成什么形态:功率谱、带通波形,还是监护仪的衍生指数。
2.6 通道的物理含义
| 通道 | 位置属性 | 在麻醉脑电里的角色 | 本数据集里的归属 |
|---|---|---|---|
| Fp1 | 左额极,国际 10-20 系统 | α 振荡在额区表达强,闭眼状态下尤其明显 | Propofol_Volunteer |
| Fp2 | 右额极,国际 10-20 系统 | SedLine 等额部监护仪的标准安放位点之一 | 三个手术患者文件夹 |
选额区不是随意的。维持麻醉无意识时出现的那条 α 节律是丘脑与前额皮层环路高度同步的产物,它在额区的表达强于枕区,这就是所谓前额 α 与清醒闭眼时枕区 α 的区别。慢 δ 振荡则是皮层大规模上下状态的体现,额区同样看得清楚。
但这里有一个必须写出来的物理限制:一个单极额部通道无法分离源与参考。Fp1 与 Fp2 之间、额区与参考电极之间的每一次体位变动都会反映在同一条曲线上;志愿者被要求闭眼,手术患者没有这个条件。官方 Usage Notes 直接点明,术中操作的瞬间改变会体现在脑电里。所以跨文件夹比较绝对幅度时要格外小心,同一文件夹内部的时序比较才是可靠的部分。
§3 它是谁的数据
3.1 贡献者与机构
官方页面署名四位作者:Elie Adam、Ohyoon Kwon、Karla Alejandra Montejo 三人隶属 MIT,Emery Brown 的挂靠单位明确写为 Picower Institute for Learning and Memory, MIT。论文的署名单位则展开为三条:MIT 的 Picower 学习与记忆研究所与脑认知科学系、麻省总医院麻醉重症与疼痛医学科、哈佛医学院麻醉学系。三条并列的署名解释了这份数据的双面身份,前半是计算神经科学,后半是临床麻醉。
这条双通道在数据形态上留下了痕迹。志愿者那条记录有着严格设计的剂量方案,五个递增台阶之后关泵,这是典型的实验设计口吻;三名手术患者的记录则是临床常规监测的副产物,给药记录是输注泵实际输出的东西,没有受控的台阶。同一份数据里塞了两种采集文化,用之前要先认出自己拿的是哪一半。
资助方的信息同样写在明处:JPB 基金会、Picower 研究所、数位个人捐赠者与麻醉计划基金,以及 NIH 的两项资助 P01 GM118269 与 R01 NS123120。利益冲突栏只写了一句话,专利待审。论文页面的披露更细:通讯作者持有麻醉状态监测与控制相关的专利,持有一家生理监测初创公司的创始人权益,并从经麻省总医院许可给 Masimo 的知识产权收取版税,相关利益由麻省总医院与 Mass General Brigham 按政策管理。这条披露对使用者是有分量的,Masimo 恰恰是三名手术患者所用 SedLine 监护仪的制造商。
3.2 时间线
| 时间 | 事件 | 出处与意义 |
|---|---|---|
| 2023-01-06 | 论文收稿 | PNAS 页面标注 received January 6, 2023 |
| 2023-04-26 | 论文被接收 | PNAS 页面标注 accepted April 26, 2023 |
| 2023-07-19 | 论文在线发表,PNAS 120 卷 30 期 e2300058120 | 论文 DOI 10.1073/pnas.2300058120 |
| 2023-07-21 | 数据集在 PhysioNet 上线,版本 1.0.0 | 官方页发布日期与 DOI 10.13026/dx44-kw30 |
| 2023-07-25 | 印刷卷期出刊 | PubMed 记录日期 2023 年 7 月 25 日 |
| 至今 | 无更新版本 | 官方 Versions 区块仅列出 1.0.0;Release Notes 写这是首发版本 |
论文上线两天后数据上线,间隔短得像刻意为之,说明发布动作从属于论文的可复现承诺,而不是独立的资源建设。同时也说明一件事:三年过去仍只有一个版本,指望它扩容到四十例,大概率会落空。
3.3 采集场景:两个完全不同的现场
现场一:研究病房。 健康志愿者接上一台 64 通道 BrainVision MRI Plus 放大器。剂量由计算机按目标效应室浓度台阶式推进,五个浓度层逐级抬高,然后关泵。期间有人每隔一段时间给出语言提示或声音咔哒,用正确反应概率判定意识消失与意识恢复。受试者全程闭眼。
现场二:手术室。 三名患者额头上贴着 Masimo SedLine 脑功能监测仪的传感器,采样率 178 或 250 Hz,前置放大器带宽 0.5–92 Hz。麻醉由人工滴定完成,给药记录是输注泵实际输出的速率与推注事件。除了主要麻醉药,官方提醒手术中还合并使用了镇痛药等其他药物。
| 维度 | 研究病房(志愿者) | 手术室(三名患者) |
|---|---|---|
| 分析通道 | Fp1,取自 64 通道全头系统 | Fp2,取自额部监护仪 |
| 采样率 | 采集端 5000 Hz | 178 Hz(丙泊酚)或 250 Hz(七氟烷) |
| 给药控制 | 计算机靶控,五个限定目标浓度 | 人工滴定,速率与推注随时调整 |
| 合并用药 | 无 | 有,镇痛药等 |
| 干扰源 | 限于眨眼,已被闭眼要求消除 | 手术操作、电刀、体位变动 |
| 伦理口径 | 受试者签署知情同意 | 无专门同意,属常规诊疗采集且已去标识 |
| 会话边界 | 一个完整的升到降的周期 | 一次真实手术全程,含清醒、诱导、苏醒 |
最后一行是理解两组记录差别的关键:频谱图覆盖整次会话,包括未被分析的时段。志愿者记录也多于分析窗,但它的剂量形状是受控的;手术记录连起止边界都不受控。
3.4 数据来源的三个层次
| 层次 | 内容 | 是否公开 | 影响 |
|---|---|---|---|
| 第一层:原始脑电 | 64 通道 5000 Hz 原始信号;SedLine 原始波形 | 否 | 无法换参考、换谱估计算法 |
| 第二层:派生表示 | 频谱图三件套、两条窄带波形 | 是 | 可复算两个调制指数,可重画论文图 |
| 第三层:给药与元数据 | 浓度、速率、推注、呼气末浓度、时窗 | 是 | 提供自变量轴,但语义逐文件夹不同 |
第二层与第三层之间没有任何机器可读的联结键,duration_start_end 里那三个数字是唯一的桥梁,而它们写的是相对频谱图的起止时刻。想把波形的局部对齐到频谱图的某个时间点,只能照官方说明手工换算(§5.5)。
§4 它最独特的东西
4.1 稀缺性:四十个人的分析,四个人能拿到
这是这份数据最鲜明的特点,也最容易被当成缺点。论文 Results 段的原话是分析了四组不同患者的脑电数据,四组各十例:10 名 18 至 35 岁志愿者、10 名 24 至 82 岁丙泊酚手术患者、10 名 50 至 90 岁且 α 波减弱的丙泊酚手术患者、10 名 53 至 72 岁七氟烷手术患者,合计四十例。全部结果在补充材料图 S1 至 S40,统计汇总在图 S41。数据集只放出每组一例,恰好对应论文图 3 至图 6。
把这两个数字并列放在一起,意图很清楚:它的定位是可复现的实例,不是可抽样的总体。这一判断决定了几乎所有用法。说它能复现论文的图是一种读法,说它四例没法做统计是另一种读法,两句都对,只取决于你站在哪一边。
为什么不多放?合理的推测是伦理与贡献者审阅机制的双重限制:这是 Contributor Review 级别的数据,访问申请要经贡献者逐项审阅,发布越多受试者,审阅与去标识的成本越高。这个推测没有官方文字支持,但与 §3.2 的时间线一致。
4.2 差异化证据:一条按控制程度排成的队列
四个文件夹的差异不是简单的三种病人,而是三种可以互相对照的实验条件。
| 条件轴 | Propofol_Volunteer | Propofol_OR_strongAlpha | Propofol_OR_weakAlpha | Sevoflurane_OR |
|---|---|---|---|---|
| 给药控制 | 计算机靶控 | 人工滴定 | 人工滴定 | 人工滴定 |
| 药物 | 丙泊酚 | 丙泊酚 | 丙泊酚 | 七氟烷,合用丙泊酚推注 |
| α 振荡表达 | 有 | 强 | 消散 | 有 |
| 年龄 | 18–35 岁 | 44 岁 | 70 岁 | 54 岁 |
| 设备 | BrainVision 64 通道 | SedLine | SedLine | SedLine |
| 采样率(采集端) | 5000 Hz | 178 Hz | 178 Hz | 250 Hz |
中间两列的对照是这张表里最有价值的一行:同一药物、同一种滴定方式、同一台监护仪,只是一个 44 岁且 α 强、一个 70 岁且 α 消散。官方专门把这两例拆成两个文件夹而不是合并成一个手术组,等于承认有无 α 表达是一个独立变量。真正想研究为什么有些人没有额区 α 的人,公开场合能拿到的对照组基本只有这一对。
4.3 具体案例:弱 α 文件夹为什么值得单独看
七十岁、丙泊酚人工滴定、脑电显示 α 波正在消散,这三个词组放在一起,指向一个临床上扎手的难题。麻醉深度监护依赖的就是额区脑电的频谱特征,厂商指数也好,α 功率也好,当患者的 α 表达本身就很弱甚至缺失时,这些间接指标的适用边界都会随之改变。
这一例给的东西是:一条慢波波形、一条被标注为弱的 α 波形、一段人工滴定的输注速率与推注记录,以及一张覆盖全局的频谱图。拿到手之后能做的事有三件。第一,在同一时间轴上叠两条波形,看调制是否在慢波一侧保留。第二,用 duration_start_end 的第 2、3 项确认两条波形各自对应的时间区间,检查它们是否等长。第三,把它与 44 岁那例在同一 dB 尺度上放在相邻两张图里,肉眼判断差别有多大。
这三件事都不需要那四十个人,也不需要联网,所有字段都在文件里。这就是四例也够用这句话的具体含义:够回答问题,不够给答案打包票。
4.4 标注粒度:这里没有标签,只有曲线
必须先说清楚,这份数据不是标注数据集,没有一个状态标签列。它的标注只有两种。
| 层级 | 形式 | 位置 | 能做什么 |
|---|---|---|---|
| 时间窗标注 | 会话总时长、α 波起止、慢波起止三个数值 | duration_start_end_*.csv | 定位分析的起点与终点 |
| 派生量 | 频谱图与两条窄带波形本身 | 各 spec 与 Wave 文件 | 自行计算两个调制指数,或自绘时频图 |
两句话足以把这个事实讲完。想要清醒或无意识的逐段标签的人,应该去看库内的 eeg-power-anesthesia(rid 1317),那里给的是逐两秒的意识标签;想要意识消失与恢复时刻的人,propofol-anesthesia-dynamics(rid 1316)给的是这两个时刻。三份数据的关系是分工而不是重复(§6.1)。
把没有标签这件事说清楚格外重要,因为 PhysioNet 上绝大多数公开脑电数据,尤其大额的长程癫痫库,都以标注为其主要卖点。带着标注在哪里的预设打开这四个文件夹,会直接得出这份数据什么都没有的错误结论。它有的是一条从量变到质变都不中断的连续时间轴,而在慢波-α 到爆发抑制的过渡研究里,这比稀疏标签有用得多。
4.5 边界条件
| 条件 | 是否适用 | 判断依据 |
|---|---|---|
| 复现论文的两个调制指数 | 适用,且是首选数据源 | 官方 Usage Notes 明确指向这两个指数 |
| 研究剂量台阶式的受控过渡 | 仅志愿者例适用 | 手术患者为人工滴定,剂量不连续 |
| 比较静脉与吸入麻醉药的差异 | 勉强 | 七氟烷只有一例,且该例合用了丙泊酚推注 |
| 建立慢 δ 与 α 频段的相位幅度耦合基准 | 适用 | 两条波形同源同窗,正是这类分析所需 |
| 训练麻醉深度回归模型 | 不适用 | 四人四会话 |
| 研究手术室中的多药相互作用 | 不适用,官方明说这是干扰而非研究对象 | Usage Notes 提到合并用药改变动力学 |
| 儿童或青少年麻醉脑电 | 不适用 | 最年轻的是 18 至 35 岁的志愿者队列,下限 18 岁 |
| 多通道空间分析 | 不适用 | 每例只给出一路额部通道的派生结果 |
§5 怎么拿到
5.1 目录结构与访问门槛
| 路径 | 内容 | 备注 |
|---|---|---|
| https://physionet.org/content/eeg-gaba-anesthesia/1.0.0/ | 项目首页,含摘要、方法、使用说明、许可 | 公开可读 |
| https://physionet.org/files/eeg-gaba-anesthesia/1.0.0/ | 文件目录 | 需认证;未通过审查时返回 403 |
| 该目录下 Propofol_Volunteer 子目录 | 志愿者会话的六份 CSV 加时窗文件 | 需通过访问审查 |
| 该目录下 Propofol_OR_strongAlpha 子目录 | 44 岁手术患者 | 同上 |
| 该目录下 Propofol_OR_weakAlpha 子目录 | 70 岁手术患者 | 同上 |
| 该目录下 Sevoflurane_OR 子目录 | 54 岁七氟烷手术患者 | 同上 |
三道门槛必须按顺序走完:注册 PhysioNet 账号并完成认证,在 CITI 完成 Data or Specimens Only Research 课程,在项目页签署 Health Data Use Agreement 1.5.0。官方 Access 段还有一条附加条款,使用者的个别研究需经贡献者审阅。也就是说,即使账号层面已经解锁,某项具体用途仍可能被拦下。审稿与解锁由人工执行,等待时间不可预测,申请时把研究计划写清楚能省掉一轮往返。
5.2 批量下载
#-- 前提:已注册认证、完成 CITI 培训并签署 DUA 1.5.0
#-- 未通过审查时直接访问文件目录会拿到 403
wget -r -N -c -np --user='你的账号' --ask-password \
https://physionet.org/files/eeg-gaba-anesthesia/1.0.0/
#-- 也可以把凭据写进 ~/.netrc,避免在命令行历史里残留密码
cat > ~/.netrc <<'EOF'
#-- machine physionet.org login 你的账号 password 你的密码
EOF
chmod 600 ~/.netrc
wget -r -N -c -np https://physionet.org/files/eeg-gaba-anesthesia/1.0.0/
第一段里密码不会出现在 shell 历史中。第二段里 netrc 的权限位必须是 600,否则 wget 会拒绝读取。递归时请不要去掉 -np,否则会把相邻项目的目录一起拉下来。
5.3 读取方式
import numpy as np
#-- base 指向解开后的 1.0.0 目录,sub 指定受试者文件夹
base = 'eeg-gaba-anesthesia/1.0.0'
sub = 'Propofol_Volunteer'
#-- 频谱图三件套:时间轴、频率轴、dB 矩阵
t_spec = np.loadtxt(f'{base}/{sub}/t_spec_{sub}.csv', delimiter=',')
f_spec = np.loadtxt(f'{base}/{sub}/f_spec_{sub}.csv', delimiter=',')
val_spec = np.loadtxt(f'{base}/{sub}/val_spec_{sub}.csv', delimiter=',')
#-- 两条窄带波形,单位 µV
alpha = np.loadtxt(f'{base}/{sub}/alphaWave_{sub}.csv', delimiter=',')
slow = np.loadtxt(f'{base}/{sub}/slowWave_{sub}.csv', delimiter=',')
#-- 会话总时长与两条波形的起止时刻
dur = np.loadtxt(f'{base}/{sub}/duration_start_end_{sub}.csv', delimiter=',')
print(t_spec.shape, f_spec.shape, val_spec.shape)
print('会话总时长 s =', float(dur[0]))
val_spec 的行列归属不能假设。正确做法是用形状去匹配它到底是频率数乘时间数还是转置,用一次布尔判断解决,比猜要省事。此外 loadtxt 读到一维还是二维取决于文件是否含换行分隔的行矩阵,读取后先看 shape 再决定后续操作,不要默认。
5.4 转表
import numpy as np
import pandas as pd
#-- 把志愿者会话的 α 波与给药曲线对齐到同一张表
base, sub = 'eeg-gaba-anesthesia/1.0.0', 'Propofol_Volunteer'
alpha = np.loadtxt(f'{base}/{sub}/alphaWave_{sub}.csv', delimiter=',')
conc = np.loadtxt(f'{base}/{sub}/propofolConcentration_{sub}.csv', delimiter=',')
#-- 用会话时长反推采样率,不照抄网页上写死的数值
total = float(np.loadtxt(f'{base}/{sub}/duration_start_end_{sub}.csv',
delimiter=',')[0])
fs_alpha = len(alpha) / total
t_axis = np.arange(len(alpha)) / fs_alpha
df = pd.DataFrame({'t_s': t_axis, 'alpha_uV': alpha})
#-- 给药文件第 1 列时间、第 2 列效应室浓度,按最近邻对齐
idx = np.abs(conc[:, 0][None, :] - t_axis[:, None]).argmin(axis=1)
df['ce_ug_per_mL'] = conc[idx, 1]
print(df.head())
这段里有两处是有意写笨的。一是采样率从文件长度反推(参见 §2.3 的张力),二是对齐采用最近邻而不是插值。最近邻在这里足够,因为要看的是趋势对照;如果下游是药代动力学拟合,那就要换成插值并显式记录插值口径。
5.5 时窗映射:把三个数字用起来
import numpy as np
#-- duration_start_end 的三项:会话总时长、α 波起止、慢波起止
root = 'eeg-gaba-anesthesia/1.0.0/Sevoflurane_OR'
d = np.loadtxt(f'{root}/duration_start_end_Sevoflurane_OR.csv',
delimiter=',')
total_s = float(d[0]) #-- 全程时长,秒
a_start, a_end = d[1], d[2] #-- α 波窗,相对频谱图时间轴
s_start, s_end = d[3], d[4] #-- 慢波窗
#-- 在频谱图时间轴上切出被分析的那一段
t = np.loadtxt(f'{root}/t_spec_Sevoflurane_OR.csv', delimiter=',')
mask = (t >= a_start) & (t <= a_end)
print('窗内时间点', int(mask.sum()), '/ 全程时间点', len(t))
官方说频谱图覆盖全会话,包含清醒、诱导、苏醒这些未被分析的时段,而分析窗起止写在 duration_start_end 里。想只画被分析的那一段,也就是复现论文图的样子,必须用这两个边界去切。注意 α 波窗与慢波窗是分开给的两组值,两者的起止不一定相同,切两条波形时要各用各自的那一组。
5.6 五个最容易出错的地方
| 位置 | 错误写法 | 后果 | 正确做法 |
|---|---|---|---|
| 采样率 | 一律按 250 Hz 处理 | 志愿者那条时间轴全错 | 由 duration_start_end 第 1 项与波形长度反推 |
| val_spec 方向 | 默认行是时间 | 频谱图上下颠倒 | 用形状对 len(f_spec) 与 len(t_spec) 判明 |
| propofolInfusionRate | 统一按三列读 | 七氟烷那一例报维度错 | 先读维度再分支处理 |
| 推注量单位 | 当成瞬时速率 | 剂量形状被放大 | 官方写明 mg 以 1 分钟的 mg 每分钟表示 |
| 分析窗 | 用整段波形当分析窗 | 混入清醒与苏醒段 | 用 duration_start_end 的第 2、3 项切窗 |
第五行的后果最隐蔽。直接拿整条 α 波去算调制指数,而前面那段是诱导前、后面那段是苏醒后,得到的指数会被拉偏,且偏的方向不可预测。
§6 与库内条目的关系
6.1 关系图谱
下面列出的是与这份数据在主题上最接近的一批在架记录。判定是否为同一份发布的依据统一为三条:官方 URL 路径不同、DOI 不同、样本构成不同。三条同时成立即可排除撞库。
| 在架记录 | rid 与链接 | 主体数据类型 | 是否同一份发布 | 判定依据 |
|---|---|---|---|---|
| eeg-power-anesthesia | rid 1317 | 54 例、0 至 50 Hz 一百个频点的两秒窗功率谱加逐两秒意识标签 | 否 | 不同的 PhysioNet 路径与 slug;样本 54 例对 4 例;成品是两秒窗的谱,这里是整段会话的谱加滤波波形 |
| propofol-anesthesia-dynamics | rid 1316 | 9 人 10 阶段的逐拍心率变异与皮电事件加意识消失恢复时刻 | 否 | 路径不同;数据主体不是频谱与滤波波形;且该记录为公开署名许可,这份为受限许可 |
| multimodal-surgery-anesthesia | rid 1313 | 101 台手术的逐五秒生理特征加伤害性刺激标注 | 否 | 样本、数据类型、标注体系均不同,官方 URL 各自独立 |
| jsmf-access | rid 1006 | 646 例术中额部脑电的衍生指数特征加术后谵妄结局 | 否 | 四个医学中心组成的独立联盟,规模与用途不同 |
| chb-mit-scalp-eeg | rid 72 | 686 个文件的癫痫头皮脑电与发作标注 | 否 | 任务域是癫痫而非麻醉,人群为儿科,与这份无任何共享记录 |
| siena-scalp-eeg | rid 1070 | 14 例癫痫患者 47 次发作的视频脑电 | 否 | 任务域与人群不同,场景是癫痫外科评估 |
| auditory-eeg | rid 1429 | 听觉注意解码的头皮脑电 | 否 | 清醒被试的听觉任务范式,与药理学无关 |
| eda-rest-sedation | rid 1487 | 镇静前后的皮肤电活动 | 否 | 信号模态为皮电,场景为门诊镇静而非全身麻醉 |
| calcium-imaging-sleep-state | rid 1465 | 小鼠宽场钙成像的睡眠状态 | 否 | 物种为小鼠、模态为钙成像;与这份混淆只因两者都研究状态切换 |
最后一行值得单说。查重时按惯例要扫一遍小鼠与啮齿类相关的 KEY,而这份数据与动物实验毫无关系。官方全名写的是 GABA 能麻醉药介导的无意识期间的脑电动态,Methods 段点明受试者是人类志愿者与手术患者,伦理段落点明由 Mass General Brigham 人类研究委员会审查,确认这一点只需读三段官方文字。
6.2 检索路径
| 你要找的东西 | 检索词 | 会落到哪里 |
|---|---|---|
| 同一课题组、同一篇论文的数据 | PNAS e2300058120 | 这份数据,DOI 10.13026/dx44-kw30 |
| 有逐段标签的麻醉脑电 | eeg-power-anesthesia | rid 1317 |
| 有意识消失与恢复时刻的丙泊酚实验 | propofol-anesthesia-dynamics | rid 1316 |
| 术中剂量与结局成对的数据 | jsmf-access | rid 1006 |
| 与这份数据常被混为一谈的动物数据 | 钙成像睡眠状态 | rid 1465 |
6.3 分工判据
| 任务 | 该用谁 | 理由 |
|---|---|---|
| 复现慢波-α 到爆发抑制过渡的调制指数 | 这份数据 | 唯一同时提供成对 α 与慢波波形和连续剂量轴的组合 |
| 训练脑电意识状态分类器 | rid 1317 | 有逐两秒标签与 54 例规模 |
| 研究丙泊酚对自主神经的影响 | rid 1316 | 有逐拍心率变异与明确的阶段划分 |
| 预测术后谵妄等临床结局 | rid 1006 | 有 646 例与结局变量 |
| 研究术中镇痛深度 | rid 1313 | 有伤害性刺激标注与镇痛监护仪对照 |
| 做法医学式的麻醉深度鉴定 | 都不合适 | 四份临床数据各有其伦理访问限制 |
§7 质量评估与局限性
7.0 阅读本节前的一个前置校正
把只有四例当成质量问题,是评价这份数据时最容易犯的错误。它的设计目标从一开始就不是提供一个可抽样的样本:论文自己分析了四十例并给出了统计汇总,发布的部分只承担让别人能重算那两条曲线这一件事。用数据量的尺子去量它,等于用体重秤量身高。
真正该问的质量问题另有两条:派生过程逆向可用的程度有多大,以及四个受试者点在多大程度上能暴露方法的脆弱性。下面的评分围绕这两条展开。
7.1 五维评分
按常用维度打分,每项 1 至 5 分,分数针对作为麻醉脑电动力学研究素材这一用途。
| 维度 | 分数 | 理由 |
|---|---|---|
| 数据规模 | 1 | 4 例、4 会话,任何统计推断都无从下手 |
| 信号保真度 | 3 | 给的是滤波后的窄带波形与频谱图,原始信号不可得;但滤波器口径写明为二阶 Butterworth,频带为 0.3 至 4 Hz 与 8 至 14 Hz |
| 任务适配 | 5 | 成对 α 与慢波加连续剂量轴的组合,在公开数据里近乎唯一 |
| 文档完整 | 4 | Methods 段交代设备、采样率、分辨率与通道选择;细节上留有 §2.3 的未说明断层 |
| 可获取性 | 2 | 三道门槛加贡献者逐一审阅,明显区别于 PhysioNet 上多数开放访问的数据 |
合计 15 分,满分 25。这个分数结构的含义很直接:拿它当主数据源会很痛,拿它当对照与复现工具几乎无损。
7.2 信息的两层保真
第一层是派生过程的可逆性。频谱图由多锥谱分析得到,两条窄带波形由二阶 Butterworth 带通得到,频带边界明确给出。这意味着重做一遍这件事在滤波这一层是可行的:拿到波形后在同一套频带上跑同样的滤波器,就能得到同类的结果。可逆性到此为止,再往回推那套 64 通道原始信号,做不到。
第二层是时窗的保真。duration_start_end 用三个数字交代了全程多长、α 波从哪到哪、慢波从哪到哪。这是官方提供的唯一位置信息,粒度粗到只有起止,但足以把分析窗从频谱图上切出来。缺少的是逐点质量标记。这条波形在哪个时间点上被电刀污染过,文件里没有任何提示,而官方 Usage Notes 又明确说手术操作会造成脑电瞬时改变。这一层的缺失必须自己填补,补的办法只能是目测频谱图舍弃可疑时段,并在方法节写清舍弃规则。
7.3 可枚举的偏差
| 偏差 | 表现 | 影响 |
|---|---|---|
| 样本量偏差 | 4 例,而论文中是 40 例 | 任何统计结论都不可由此单独得出 |
| 选择偏差 | 公开的是每组典型的一例 | 不能据此推断各组内部的变异范围 |
| 设备混淆 | 志愿者与手术患者用不同放大器与采样率 | 跨文件夹比较绝对值时设备差异无法剥离 |
| 通道混淆 | Fp1 与 Fp2 混用 | 空间信息不可比 |
| 年龄偏差 | 最年轻 18 岁,三名手术患者分别 44、70、54 岁 | 覆盖不到儿科与高龄两极 |
| 合并用药 | 手术会话中使用镇痛药等 | 动力学改变无法归因于单一药物 |
| 手术操作 | 脑电的瞬时改变可能来自操作而非药物 | 短时段分析易受污染且无质量标记 |
| 分析窗截断 | 两条波形被截断在指标计算区间内 | 看不到这一区间之外的波形,不能外推 |
| 利益冲突 | 官方写明专利待审;论文披露通讯作者持相关专利、持初创公司权益、从许可给 Masimo 的知识产权收取版税 | 解释任何与设备有关的差异时应作为背景信息计入 |
最后一行不是质量问题而是诚信披露的一部分:发布方把利益冲突写在明处,而 Masimo 恰恰是这批手术患者所用监护仪的制造商。这一事实本身不需要评价,但任何用这批数据做设备比较的研究都应该在方法节里提到。
7.4 外部验证矩阵
| 验证角色 | 做法 | 讲得通的地方 | 风险 |
|---|---|---|---|
| 方法复现 | 按论文口径重建两个调制指数 | 原料与论文一致 | 没有第三方公开参考实现,只能比图形不比数值 |
| 单例剖析 | 拿弱 α 那例研究 α 缺失 | 公开场合少见的对照 | 一例无法区分年龄效应与个体效应 |
| 教学演示 | 用志愿者会话演示剂量与频谱的关系 | 曲线干净、可直接读取 | 容易把一次实验当成普遍规律 |
| 设备比较 | 志愿者 5000 Hz 对手术 178 或 250 Hz | 恰好覆盖两端 | 同时混淆了设备、人群、给药方式三个变量 |
第三行的风险要特别点出。正因为志愿者那条曲线太干净,它特别适合演示,也特别容易被误读成规律就是如此。计算机靶控的健康年轻志愿者不等于真实手术场景,这句话应当出现在每一张用它画出的图的图注里。
7.5 它为何还能用下去
它没有社区、没有更新、没有大规模传播,官方页面 Discovery 区块显示的项目浏览次数为个位数。它能继续存在的理由只有一条:研究慢波-α 与爆发抑制之间的过渡动力学时,公开数据里再没有第二处同时给了对称的两条窄带波形与连续的剂量轴。只要有研究组想复核那两个调制指数,或者想在真实手术数据上试试自己的滤波器与谱估计,这四个文件夹就是入口。
§8 生态位与历史
8.1 沿革:论文两天之后的产物
它的诞生方式在公开数据里并不常见:先有论文,再有数据发布,中间只隔两天。多数 PhysioNet 数据要么先建库再写文,要么由一个项目同时给出两者;这里是论文在线发表两天后,作者把支撑图 3 至图 6 的原料挂上了平台。这个时间差说明发布动作从属于论文的可复现承诺,而不是独立的资源建设。
后果写在 §3.2 的版本表里:三年过去仍只有一个版本,Release Notes 那句首发说明至今仍是唯一一条。它不是那种会因为用户反馈而扩容的资源,也没有相应的维护承诺。
8.2 与挑战赛的关系:没有关系
这一点必须明确否定。PhysioNet 上多数脑电数据会通过挑战赛扩散,比如 CinC 系列、I-CARE,癫痫发作检测也是有明确赛事的。这份数据没有出现在任何 PhysioNet 挑战赛的来源清单里,官方页面也没有挑战赛链接,它不是经由比赛扩散的那类资源,传播路径只有论文引用一条。
这个否定有实际用处:有些人会从 PhysioNet 的脑电数据联想到是不是某个赛事的基线数据,进而期待带有官方划分的训练集与测试集、榜单与评测代码。这里没有划分、没有榜单、没有基线,只有四个文件夹。
8.3 生态位
| 角色 | 训练素材 | 复现基点 | 临床结局库 |
|---|---|---|---|
| 规模要求 | 数十到数百例 | 一两例即可 | 数百例加结局变量 |
| 标签要求 | 逐段或逐点标签 | 可以没有标签 | 结局必须完整 |
| 代表情形 | rid 1317(54 例带标签) | 这份数据(4 例,重方法不重样本) | rid 1006(646 例带结局) |
| 失败模式 | 标签口径不统一 | 单点故障,一例不可用即全盘失效 | 混杂因素过多 |
这份数据坐在中间那一列,而且是这一列里比较少见的依附于单篇论文的成员。中间列的成员通常稀少,因为它们对宿主论文的依赖太重;好处是方法节可以写得很具体,坏处是论文一旦不再被引,数据也随之沉默。
8.4 谁在用,用来做什么
从可核实的线索看,使用者集中在三类群体。第一类是想复核那两个调制指数的人,任何人质疑这种调制是否真实存在,都能用这四个文件夹算一遍,这是它被发布时唯一明确的用途,官方 Usage Notes 的第二句写的就是这个。第二类是做麻醉深度监测新方法的人,开发新指标时需要一条已知的剂量轴作为参照,志愿者会话的五个浓度台阶恰好提供了这样一条参考曲线。第三类是教学的人,四条典型路径正好构成一套完整的教学样例库,且全部是 CSV,不需要任何专用软件。
PNAS 页面显示的论文使用度,包括约 7920 次浏览、25 次被引、被 18 家新闻机构报道,是论文侧的传播数据,不等于数据集的使用量。两者之间隔着那道 Contributor Review 的访问门槛,这也解释了为什么这种依附型数据的引用数通常远低于论文。
§9 使用指引
9.1 三种角色的路径
| 角色 | 第一天 | 第一周 | 注意事项 |
|---|---|---|---|
| 复核者 | 读论文 Methods,跑通 §5.3 | 重建两个调制指数,与论文图对比 | 先判明 val_spec 的行列方向 |
| 方法开发者 | 取志愿者会话,建立剂量与调制的对照 | 在三个手术会话上测试鲁棒性 | 三个手术文件夹的给药语义不同 |
| 教学者 | 下载并画出四张频谱图 | 用强 α 与弱 α 两例做对照板书 | 年龄与设备差异要在图上标注 |
9.2 一条最短的可用流程
#-- 1. 完成 CITI 培训与 DUA 签署后,确认文件目录可访问
curl -s -o /dev/null -w '%{http_code}\n' \
https://physionet.org/files/eeg-gaba-anesthesia/1.0.0/
#-- 2. 从目录索引里确认四个受试者文件夹都在
curl -s https://physionet.org/files/eeg-gaba-anesthesia/1.0.0/ \
| grep -o 'Propofol_[A-Za-z_]*\|Sevoflurane_[A-Za-z_]*' | sort -u
import numpy as np
#-- 3. 读志愿者会话,建立时长与三条主键的关系
sub = 'Propofol_Volunteer'
root = f'eeg-gaba-anesthesia/1.0.0/{sub}'
d = np.loadtxt(f'{root}/duration_start_end_{sub}.csv', delimiter=',')
al = np.loadtxt(f'{root}/alphaWave_{sub}.csv', delimiter=',')
sl = np.loadtxt(f'{root}/slowWave_{sub}.csv', delimiter=',')
print('会话总时长 s =', float(d[0]))
print('α 波样本数 =', len(al), '慢波样本数 =', len(sl))
print('由时长反推的采样率 =', len(al) / float(d[0]))
第二步的 grep 只是确认目录结构,不是解析正式清单,官方没有提供汇总式RECORDS文件,文件列表只能从目录索引上读。第三段最后一行输出是整个流程里最关键的一个数:它告诉你到底该用哪个采样率,而这一点在官方页面上是含糊的(§2.3)。
9.3 引用方式
| 引用对象 | 内容 |
|---|---|
| 数据集 | Adam, E., Kwon, O., Montejo, K. A., & Brown, E. (2023). Electroencephalogram dynamics during unconsciousness mediated by GABAergic-anesthetics (version 1.0.0). PhysioNet. RRID:SCR_007345. https://doi.org/10.13026/dx44-kw30 |
| 配套论文 | Adam, E., Kwon, O., Montejo, K. A. & Brown, E. N. (2023) Modulatory dynamics mark the transition between anesthetic states of unconsciousness. PNAS 120(30), e2300058120 |
| 平台 | PhysioNet 平台综述(2026,Nature Health),按项目页 Cite 按钮生成 |
| 建议额外注明 | 使用了哪几个文件夹;志愿者与手术记录是否分别处理;分析窗是否按 duration_start_end 的第 2、3 项切出 |
最后一行的三项不是格式要求,而是让结果可比的必要条件。同一个分析在切窗与不切窗之间可能给出不同的调制深度,评审无从通过一句「用了某版本的数据」反推出你做了什么。
§10 AI 使用声明
| 项目 | 声明 |
|---|---|
| 文本生成 | 本文由 AI 辅助撰写,结构与事实由编辑把关 |
| 数字来源 | 规模、结构、边界三类数字取自 PhysioNet 官方页面(版本 1.0.0)与 PNAS 论文页面明写的内容 |
| 文件核验 | 受限访问:未认证且未签署 DUA 时文件目录返回 403,未取得任何原始文件;文件名与列语义均引自官方 Data Description 段原文 |
| 不一致说明 | 共享波形的采样率在 Methods 段(5000、250、178 Hz)与 Data Description 段(250 或 178 Hz)之间存在官方未解释的张力,已在 §2.3 与 §5.4 说明,正文未自行取值 |
| 论文信息 | PNAS 页面标题、卷期、DOI、PMID、收稿与接收日期、并列作者标注、使用度指标均已核对;论文 Materials and Methods 中的完整公式未读全文 |
| 未核验内容 | 四条会话的具体时长;每个文件夹对应的实际采样率;val_spec 的行列方向;七氟烷那例的输注文件是否确为两列;duration_start_end 三个项目组在扁平数组中的确切下标;四十例队列是否另有公开计划 |
| 图像 | 无配图;封面为生成式提示词,不含真实脑电与计算公式 |
| 免责 | 本文不构成临床建议;数据使用须遵守 Health Data Use Agreement 1.5.0,并经过贡献者审阅 |
常见误读(避坑清单)
| 说法 | 实际 | 证据 |
|---|---|---|
| 它是四十个人的数据 | 论文分析了四十例,数据集只公开四例,每组一例 | 论文 Results 段列出四组各十例;数据页写明共享四名受试者的数据 |
| 它是动物实验的脑电 | 全部为人类受试者,由 Mass General Brigham 人类研究委员会审查 | Methods 段写明志愿者与手术患者;Ethics 段写明伦理委员会 |
| 它是原始脑电,可以换参考电极或换标度 | 没有原始波形,只有频谱图三件套与两条滤波后的窄带波形 | Data Description 段只列出这三类文件 |
| 四个文件夹的分析通道一致 | 志愿者取 Fp1,三名手术患者取 Fp2 | Methods 段逐个写明所选通道 |
| 它有清醒或无意识的逐段标签 | 没有标签,只有三个数字说明会话时长与两条波形的起止 | Data Description 段;需要标签应看向 rid 1317 |
| 采样率写死是 250 Hz | Data Description 段写 250 或 178,Methods 段写志愿者采集 5000 Hz | 两段原文并存且官方未解释;应由文件长度反推 |
| 四个文件夹的给药文件列数一样 | 志愿者给效应室浓度,手术丙泊酚给速率加推注,七氟烷那例只有推注 | Data Description 段逐文件夹写明列语义 |
| 频谱图只覆盖被分析的那一段 | 频谱图覆盖整次会话,含清醒、诱导与苏醒;分析窗另写在 duration_start_end 里 | Data Description 段末段原话 |
| 它是公开库,注册就能下 | Contributor Review 级别:认证、CITI、DUA 三道关,外加逐研究审阅 | Access 段与许可区块 |
| 志愿者会话能代表真实手术场景 | 志愿者的剂量形状受控、闭眼、无合并用药与手术干扰 | Methods 段与 Usage Notes 段的对照 |
FAQ
Q1:为什么论文分析了四十个人,公开的数据只有四个人?
官方并没有解释这个取舍。可以确认的是两件事:论文四组各十例的结果放在补充材料图 S1 到 S40,统计汇总在图 S41;数据集覆盖的是论文图 3 到图 6,也就是每组各一例。把它理解成可复现的样例而不是抽样总体,是唯一不会踩坑的读法。
Q2:它的许可到底是什么,能不能直接下载?
不能直接下载。它是 PhysioNet Contributor Review Health Data License 1.5.0,配套 Health Data Use Agreement 1.5.0,还要求完成 CITI 的 Data or Specimens Only Research 培训。Access 段另有一条,使用者的个别研究需经贡献者审阅。这几项缺一不可,缺任一项时文件目录直接返回 403。
Q3:共享的 α 波和慢波是什么采样率?
官方文字在这里不一致:Data Description 段写 α 波与慢波的采样率是 250 Hz 或 178 Hz,Methods 段写志愿者采集为 5000 Hz、手术丙泊酚为 178 Hz、七氟烷为 250 Hz。稳妥做法是从 duration_start_end 的第 1 项读出会话总时长,用波形长度除以它反推,不要照抄任一处的数值。
Q4:它有让我可以直接做监督学习的标签吗?
没有。这份数据不包含任何逐段或逐点的状态标签,唯一的标注是 duration_start_end 里交代会话时长与两条波形起止的三个数字。需要标签的话,库内的 eeg-power-anesthesia(rid 1317)给的是逐两秒的意识标签加信号质量标记,propofol-anesthesia-dynamics(rid 1316)给的是意识消失与恢复时刻。
Q5:志愿者会话和手术会话能放在一起比较吗?
能比较形状,不能比较绝对值。志愿者用 64 通道系统、分析通道 Fp1、计算机靶控给药、全程闭眼;手术患者用 SedLine、分析通道 Fp2、人工滴定、且有合并用药与术中操作干扰。设备、人群、给药方式三个变量同时变化,任何跨组的定量差异都无法归因。
Q6:我只想画一张像论文那样的图,最少要读几个文件?
三个:时间轴文件、频率轴文件、取值矩阵文件,外加 duration_start_end 里的第 2、3 项用来把分析窗切出来。如果还想叠加两条波形,再加两个 Wave 文件,但要先把它们的采样率按 Q3 的办法反推出来。
术语表
| 术语 | 解释 |
|---|---|
| GABA 能麻醉药 | 通过增强 γ-氨基丁酸受体抑制发挥作用的全身麻醉药,这里涉及丙泊酚与七氟烷 |
| SDA | slow-delta–alpha 的缩写,指慢 δ 振荡与 α 振荡并存的无意识状态 |
| 慢 δ 振荡 | 0.3–4 Hz 的皮层节律,对应上下状态下的神经元活动抑制 |
| α 振荡 | 8–14 Hz 的节律;麻醉无意识时的额区版本是丘脑与前额环路的同步产物 |
| 爆发抑制 | 等电位平线与爆发段交替出现的深度脑失活状态,常由剂量过高导致 |
| 上下状态 | 皮层神经元在去极化与超极化之间慢速切换的状态,对应慢 δ 振荡 |
| 多锥谱分析 | 用多个正交窗函数估计功率谱的谱估计方法,这里的频谱图由此得到 |
| Butterworth 滤波器 | 通带内最平坦的一类滤波器,这里两个频带均用二阶版本 |
| 调制指数 | 论文构造的 α 调制指数与慢调制指数,用于连续跟踪两个状态之间的过渡 |
| 效应室浓度 | 药物在作用部位而非血浆的推定浓度,是靶控输注的控制目标 |
| 靶控输注 | 按药代模型自动调节输注速率以逼近目标浓度的给药方式 |
| 呼气末浓度 | 呼气末期呼出气中的药物浓度,临床上用于估算吸入麻醉药的分压 |
| SedLine | Masimo 公司的额部脑功能监测仪,是这里三名手术患者的数据来源设备 |
| Fp1 与 Fp2 | 国际 10-20 系统中的左、右额极电极位点 |
| 意识消失与恢复 | 由行为正确反应概率判定的两个时刻,英文缩写为 LoC 与 RoC |
| Contributor Review | PhysioNet 的一档受限访问级别,需要 DUA、CITI 且逐研究经贡献者审阅 |
| DUA | Data Use Agreement,数据使用协议 |
| CITI | 提供研究伦理培训的在线课程平台 |
| RRID | 研究资源标识符,这里 PhysioNet 自身的标识为 SCR_007345 |
附录 A:核心数字核对表
| 数字 | 出处 |
|---|---|
| 版本 1.0.0,版本日期 2023 年 7 月 21 日 | 项目页 Versions 区块与标题区 |
| 版本 DOI 10.13026/dx44-kw30 | 项目页 Discovery 区块 |
| 最新版本 DOI 10.13026/7ygq-3f91 | 项目页 Discovery 区块 |
| 公开受试者 4 名 | 项目页 Methods 段与 Data Description 段 |
| 论文分析队列 4 组共 40 例 | PNAS 论文 Results 段 |
| 志愿者年龄 18–35 岁 | 项目页 Methods 段 |
| 三名手术患者年龄 44、70、54 岁 | 项目页 Methods 段 |
| 论文各队列年龄区间 24–82、50–90、53–72 岁 | PNAS 论文 Results 段 |
| 志愿者采集 64 通道、5000 Hz、带宽 0.016–1000 Hz、0.5 µV 每最低有效位 | 项目页 Methods 段 |
| 手术丙泊酚 178 Hz | 项目页 Methods 段 |
| 七氟烷例 250 Hz、带宽 0.5–92 Hz、16 位、29 nV | 项目页 Methods 段 |
| 共享 α 与慢波采样率 250 或 178 Hz | 项目页 Data Description 段,与 Methods 段存在张力 |
| 频带边界 0.3–4 Hz 与 8–14 Hz,二阶 Butterworth | 项目页 Methods 段 |
| 四个文件夹名称 | 项目页 Data Description 段 |
| 每文件夹六份 CSV 加时窗文件 | 项目页 Data Description 段 |
| 论文 PNAS 120 卷 30 期 e2300058120,PMID 37467269 | 项目页 References 段与 PubMed |
| 论文收稿 2023-01-06、接收 2023-04-26、在线 2023-07-19 | PNAS 论文页面 |
| 许可 Health Data License 1.5.0 与 Use Agreement 1.5.0 | 项目页 License 与 Data Use Agreement 区块 |
| 所需培训 CITI Data or Specimens Only Research | 项目页 Required training 区块 |
| 伦理审查方 Mass General Brigham Human Research Committee | 项目页 Ethics 段 |
| 资助 NIH P01 GM118269 与 R01 NS123120 | 项目页 Acknowledgements 段 |
| 利益冲突标注「专利待审」 | 项目页 Conflicts of Interest 段 |
| 项目浏览次数为个位数 | 项目页 Discovery 区块 |
相关数据集导航
以下为站内 AI-Ready 数据集百科中与本词条共享多个主题标签的相关数据集,按相关度降序排列:
- eeg-power-anesthesia—共享标签:生理信号 / 脑电信号 / 麻醉与镇静
- propofol-anesthesia-dynamics—共享标签:生理信号 / 麻醉与镇静 / 丙泊酚
- multimodal-surgery-anesthesia—共享标签:生理信号 / 麻醉与镇静 / 术中监测
- jsmf-access—共享标签:生理信号 / 脑电信号 / 麻醉与镇静
- eda-rest-sedation—共享标签:生理信号 / 麻醉与镇静
- calcium-imaging-sleep-state—共享标签:神经科学 / 状态切换
- auditory-eeg—共享标签:生理信号 / 脑电信号
- chb-mit-scalp-eeg—共享标签:生理信号 / 脑电信号 / 时序信号
- siena-scalp-eeg—共享标签:生理信号 / 脑电信号 / 时序信号
- tuh-eeg-corpus—共享标签:生理信号 / 脑电信号 / 神经科学
导航说明:本章节由全站统一标签体系自动计算生成(标签重合度算法),双向可达;点击链接可跳转至对应数据集词条。

