当前位置：首页 > news >正文

9. 从个体到群体：Cartool静息态EEG微状态全流程实战解析

news 2026/6/19 13:17:16

1. 静息态EEG微状态分析入门指南

第一次接触静息态EEG微状态分析时，我和大多数初学者一样感到困惑。这种通过脑电图信号识别大脑功能状态的技术，听起来既神秘又复杂。但经过多次实战后，我发现只要掌握几个关键概念，就能轻松入门。

微状态分析的核心思想很简单：大脑在静息状态下会自发产生特定的电活动模式，这些模式会持续几十到几百毫秒，然后切换到另一种模式。就像观看一部电影，大脑也在不断切换不同的"画面"。通过识别这些典型模式（即微状态），我们可以了解大脑的功能组织方式。

Cartool软件是这个领域的利器。它提供了完整的分析流程，从个体水平聚类到组水平分析，再到最终的模板匹配。我特别喜欢它的可视化界面，即使没有编程基础也能快速上手。不过要注意，软件对数据格式有严格要求，必须使用Brain Vision（BP）格式，这是很多新手容易踩的第一个坑。

2. 数据预处理的关键步骤

2.1 格式转换与基础处理

在开始微状态分析前，数据预处理是重中之重。我习惯在EEGLAB中完成基础预处理（去噪、去伪迹等），然后导出为BP格式。这里有个实用技巧：导出时务必检查电极名称和顺序是否一致。曾经有一次，我因为不同被试的电极顺序不一致，导致后续分析全部出错，白白浪费了两天时间。

Cartool对数据量比较敏感。如果电极数过多或采样率太高，软件很容易崩溃。我的经验是：将电极数控制在64个以下，采样率降到500Hz左右。虽然会损失一些信息，但稳定性大幅提升。对于静息态分析来说，这样的设置已经足够。

2.2 静息态专用预处理

静息态微状态分析需要一些特殊处理。我通常会做2-20Hz的带通滤波，这能保留微状态分析最相关的频段。全脑平均参考也是必须的，它能消除共同参考的影响。有个细节需要注意：如果数据已经分段，要确保所有被试的分段长度一致，否则后续聚类会出现问题。

3. 个体水平聚类实战

3.1 参数设置技巧

进入Cartool后，在Tools菜单选择"Segmentation of EEG files"。第一次操作时，我被众多参数选项搞得晕头转向。经过多次尝试，我总结出几个关键设置：

在File Presets中选择"First Stage"，这是个体水平聚类的标志。加载数据时，建议一次不要超过20个被试，否则容易导致软件无响应。Output选项中，只需勾选ep文件，它包含了各类地形图的关键信息。

3.2 聚类算法选择

Computation Presets部分需要特别注意。对于静息态数据，我强烈建议选择"Using only GFP Peaks Data"。GFP（全局场功率）峰值处的信噪比最高，能获得更稳定的聚类结果。算法方面，k-means和T-AAHC各有优劣：k-means速度更快，但可能陷入局部最优；T-AAHC更稳定，但计算量较大。我通常先试用k-means，如果结果不理想再换T-AAHC。

聚类数范围保持默认的1-12即可。记得勾选"Common Best Clustering Directory"，它会自动保存每个被试的最优聚类结果。运算时间取决于数据量，30个被试大概需要3-5小时，建议在晚上开始计算。

4. 组水平聚类深度解析

4.1 从个体到群体的过渡

完成个体水平聚类后，接下来是组水平分析。在File Presets中选择"Second Stage"，然后加载所有被试的ep文件。这里有个常见误区：有人会误加载原始dat文件，这会导致分析完全错误。正确的做法是加载RSWhole Sub.BestClustering文件夹中的ep文件。

Data Preprocessing部分要选择"Using Whole Data"，因为现在是对模板图进行聚类，而不是原始EEG数据。算法选择最好与个体水平一致，这样结果更具可比性。聚类数范围可以扩大到1-15，给软件更多选择空间。

4.2 结果解读与验证

运算完成后，Cartool会弹窗显示最优聚类数。有趣的是，结果常常显示5类最优，而非文献中常见的4类。这不是错误，可能反映了被试群体的特殊性。我建议将得到的模板图与经典图谱（A/B/C/D类）对比，确定对应关系。

绘制模板图时，可以使用EEGLAB的topoplot函数。记得检查电极位置文件是否正确，这是影响可视化效果的关键因素。如果发现某些模板图难以归类，可能是数据质量或参数设置问题，需要重新检查。

5. 模板匹配与指标提取

5.1 精细调整匹配参数

最后一步是将组水平模板匹配到个体数据。在"Fitting templates to EEG files"中加载模板ep文件和所有被试的dat文件。Labeling Parameters部分需要特别注意："Labeling at low correlations"选项会影响结果解读。如果勾选，当某个时刻的地形图与所有模板的相关性都低于阈值（默认0.5）时，该时刻会被标记为无类别。这意味着各类别的总时间占比可能不到100%。

Temporal postprocessing中的平滑操作要谨慎使用。早期研究普遍推荐平滑，但最新版本的Cartool反而建议不做平滑。我的经验是：对于高质量数据可以不做平滑，但噪声较大的数据，适度的平滑（不超过30ms）能提高结果稳定性。

5.2 关键指标计算

在Variables To Extract中，最常用的三个指标是：Mean Duration（平均持续时间）、Time Coverage（时间占比）和Segment Count Density（出现频率）。注意Mean Duration的单位是采样点数，需要根据采样率转换为毫秒。例如，对于500Hz的数据，采样点乘以2就是毫秒数。

马尔可夫链分析能揭示微状态间的转换规律，但计算量很大。我通常先计算基础指标，确认数据质量后再单独进行马尔可夫分析。生成的转换概率矩阵很有价值，可以反映大脑状态切换的模式特征。