1. 摘要

神经影像数据分析和解释需要结合多学科的共同努力,不仅依赖于统计方法,而且越来越多地依赖于与其他脑源性特征相关的关联,如基因表达、组织学数据、功能和认知结构。在这里,我们介绍了BrainStat,它是一个工具箱,包括(i)在体素空间和皮层空间的神经影像数据集中的单变量和多变量线性模型,以及(ii)死后基因表达和组织学的空间图谱,基于任务的功能磁共振成像元分析,以及几个常见静息态功能磁共振成像大脑皮层模板在内的多模态特征关联。统计和特征关联结合成一个关键的工具箱简化了分析过程并加速了跨模态研究。工具箱用Python和MATLAB实现,这两种编程语言在神经影像和神经信息学领域中广泛使用的。BrainStat是公开提供的,并包括一个可扩展的文件。

2. 引言

神经影像可以在全脑范围内测量个体的大脑形态、微结构、功能和功能连接。通过越来越多的强大的图像处理技术,这些数据可以纳入标准化的参考空间,包括体素空间,如常用的MNI152空间、皮层表面基础空间,如fsaverage,MNI152表面,或grayordinates以及一些脑区分割模板。将神经影像数据注册到一个公共空间,我们能够允许应用统计分析在每个测量单位进行统计检验,包括的单变量广义线性模型和混合效应模型。通常,这种分析需要使用多种工具和程序来进行,减少了分析流程的可重复性并增加了人为错误的风险。在当前的论文中,我们提出了BrainStat,一个统一的工具箱,可以在一个统一的、透明的和开源的框架中实现这些分析。

神经影像研究的高级分析工作流程越来越依赖于以前获得的多种影像模式和非影像的数据集的可用性。当映射到与神经影像测量相同的参考框架时,这些数据集可以用于神经影像结果的背景化研究,并帮助解释和验证结果。例如,结果可能在已经建立的大脑功能架构中进行解释,如基于静息态功能磁共振成像内在功能社区或功能梯度。另一种常见的背景化方法是使用Neurosynth、NiMARE或BrainMap进行自动元分析。这些工具提供了对可能的数千项先前发表的功能磁共振成像研究进行荟萃分析的能力。将一个统计图谱与一个与认知术语相关的大脑激活地图的数据库联系起来,即所谓的元分析解读。因此,提供了一种定量的方法来推断与空间统计模式相关的可信认知过程。最后,将转录组学和组织学的死后数据集映射到一个共同的神经影像空间,使神经影像发现与基因表达和微观结构模式联系起来。这些发现可以提供关于大脑中分子和细胞特性的信息,这些特性在空间上与观察到的统计地图共同变化。通过结合这些特征关联技术,可以推断出神经影像学检查结果的功能、组织学和遗传学相关性。BrainStat提供了一个集成的解码引擎来执行这些多模态特性关联。

我们的工具箱有一个用Python和MATLAB并行实现的工具箱,这是神经影像研究领域中两种常见的编程语言。因此,BrainStat的一个关键设计选择是为了最大限度地同质化实现,从而提高工具的可访问性,并帮助用户在没有先验编程专业知识的情况下学习一种或两种编程语言。BrainStat依赖于一个简单的面向对象框架来简化分析工作流程。这个工具箱可以在https: //github.com/MICA-MNI/BrainStat上公开获得,还有一些文档可以在https://brainstat.readthedocs.io/上获得。我们将工具箱划分为两个主要模块:统计模块和背景化模块(图1)。在本文的其余部分中,我们描述了如何执行图1中所示的分析。

图1. BrainStat工作流程。工作流程被分为一个统计模块(深灰色)用于固定和混合效应线性模型,一个背景化模块(浅灰色)用于与外部数据集上下文化结果。要使用统计模块,用户必须提供体素、顶点或包裹水平的数据,必须指定固定或混合效果模型,并且必须指定模型对比度。一旦指定了这些,BrainStat计算t值并修正模型的p值。t值图或任何其他大脑图可以用于上下文模块,将统计结果与任务功能磁共振元分析的标记嵌入,并建立功能层次、遗传表达和组织学标记。

3. 统计模块

统计模块建立在SurfStat的基础上,这是一个经典但不维护的MATLAB包,用于实现固定和混合效应线性模型。为了在BrainStat实现中创建和拟合这样的模型,用户提供了一个主题×区域×响应变量的矩阵以及一个使用直观的模型公式框架创建的预测模型。这种方法允许将固定/随机效应直接定义为主要感兴趣的变量或控制协变量,促进了横向分析和纵向分析。对于混合效应模型建模,BrainStat使用了g侧规范,这样它就可以容纳多个随机效应作为独立效应,目前的拟合是通过严格的最大似然估计进行的。

为了比较感兴趣的变量(如健康/疾病、年龄)的影响,必须指定一个对比变量。BrainStat可以处理单变量或多变量响应变量数据,并提供两种广泛用于多重比较校正的分析选项,即错误发现率和随机场模型。错误发现率控制了数据中的特征变量(即顶点、体素、包裹)的假阳性比例,而随机场理论修正了曾经报告假阳性发现的概率(包括在峰值点水平和在团簇水平)。

为了说明统计模块,我们下载了微结构信息连接组学(MICAMICs)数据集中的70名被试的皮层厚度和人口统计学数据,其中12人是扫描两次的(图2A)。我们创建了一个线性模型,以年龄和性别以及它们的交互效应作为固定效应,被试作为随机效应(图2B)。这些都是使用固定效应和混合效应类(命名为如它可能包含随机和固定效应)来设置的。接下来,我们将对比数据定义为年龄,即正t值表示皮层厚度随着年龄的增长而减少。该模型采用单尾检验拟合皮层厚度数据。图2C绘制了由随机场理论得到的t值、聚类p值和峰值p值,以及经过错误发现率校正得到的顶点p值。我们发现,在MICAMICs数据集中,年龄对基于团簇水平的皮层厚度有影响。但在随机场理论中,这些团簇中没有显著的顶点级峰值,并且在顶点水平上没有边际显著性。这表明年龄对皮层厚度的影响覆盖了很宽的区域,而不是局部病灶。虽然我们在这个例子中使用了一个自由的团簇水平阈值(p<0.01),但我们通常建议使用一个更严格的阈值(p<0.001),特别是当使用很少的空间平滑的数据时。

图2.拟合年龄对皮层厚度影响的固定效应一般线性模型的Python代码。(A)BrainStat包含的MICA-MICS数据集包含皮层厚度和人口统计学数据。人口统计学数据包含散列的被试ID(SUB_ID)、扫描次数、年龄z分数(AGE_AT_SCAN)和性别(SEX)。(B)我们创建了一个以

BrainStat:一个用于全脑统计和多模态特征关联的工具箱相关推荐

  1. 一个用于Allen脑图谱基因数据的工具箱|abagen详细使用教程-获取基于脑区的基因表达矩阵(脑区*gene)

    艾伦人类脑图谱(Allen Human Brain Atlas) 艾伦人类脑图谱是一个由艾伦脑科学研究所(Allen Institute for Brain Science)开发的在线基因表达图谱数据 ...

  2. MICCAI 2022 | CLFC:基于对比学习的多模态脑肿瘤分割与单模态正常脑图像的特征比较

    MICCAI 2022 | CLFC基于对比学习的多模态脑肿瘤分割与单模态正常脑图像的特征比较 Multimodal Brain Tumor Segmentation Using Contrastiv ...

  3. nextvlad:一个用于大规模视频分类的聚集帧级特征的有效神经网络 NeXtVLAD: An Efficient Neural Network toAggregate Frame-level Fea

    nextvlad:一个用于大规模视频分类的聚集帧级特征的有效神经网络 NeXtVLAD: An Efficient Neural Network toAggregate Frame-level Fea ...

  4. EEGNet:一个小型的卷积神经网络,用于基于脑电的脑机接口

    脑机接口(BCI)利用神经活动作为控制信号,可以与计算机直接通信.这种神经信号通常从各种研究充分的脑电图(EEG)信号中选择.对于给定的脑机接口(BCI)范式,特征提取器和分类器是针对其所期望的脑电图 ...

  5. 从局部脑到全脑:时空域分层神经网络脑电情绪识别模型

    作者丨张玮玮 学校丨东北大学硕士生 研究方向丨情绪识别 引言 论文动机 近年来,脑电图 (EEG) 在处理情感识别方面得到了广泛的应用.虽然近年来人们提出了各种各样的脑电信号情感识别方法,但为了进一步 ...

  6. 实录丨戴琼海:深度学习遭遇瓶颈,全脑观测启发下一代AI算法

    2020-09-03 20:03:00 目前我们还无法精细到神经元级别的观测,只能从功能层面理解大脑,但这些成果也启发了很多经典的人工智能算法,例如卷积神经网络启发自猫脑视觉感受野研究,胶囊网络启发自 ...

  7. NPP: 重复经颅磁刺激产生抗抑郁效果的基础:全脑功能连接与局部兴奋度变化

    重复经颅磁刺激 (rTMS) 是重度抑郁症 (MDD) 常用的治疗方式,但是我们对经颅磁刺激产生抗抑郁效果的机制了解甚少.此外,我们也缺乏能够用以预测和追踪临床效果的大脑信号,而这些信号能够帮助进行分 ...

  8. Maplab:一个用于视觉惯性建图和定位研究的开源框架

    摘要 鲁棒且精确的视觉惯性估计是当今机器人领域的重要挑战.能够用先验地图(prior map)进行定位(localize)并获得准确且无漂移的姿态估计,可以推动该系统的适应性.然而,目前大多数可用的解 ...

  9. Cerebral Cortex:全脑功能动力学追踪抑郁症状的严重程度

    <本文同步发布于"脑之说"微信公众号,欢迎搜索关注~~> 20%的人口报告抑郁症状,与大规模脑网络功能完整性改变有关.瞬时的大脑功能和抑郁症状之间的联系,以及这些关系对 ...

最新文章

  1. Case study:数据库网页构建原理和实践
  2. Android Vector笔记
  3. keepalived配置文件详解   ​
  4. 【原创】Mac 工作协同之文件共享问题
  5. 让kbmmw 4.8 支持ios 64
  6. 【微波】【1】色散关系式与 TEM导波、TE导波、TM导波、混合波区别和特性
  7. 知网查重提交论文显示服务器错误,职称论文在进行知网查重时,经常出现的错误有哪些?...
  8. 设计师工资低?10大网站助你快速涨1万身价
  9. VASP+Phono3py计算声子linewidth
  10. 12. JLINK在线调试+软件调试方法与技巧
  11. 科技圈的“三座大山”正在崩塌中。。。
  12. 为什么选择嵌入式方向
  13. 【Matlab瑕疵检测】阙值瓶盖瑕疵检测【含源码 730期】
  14. idea出现decompiled .class file 解决方案
  15. 两相编码器的FPGA驱动
  16. 国产设计软件都有哪些,国产设计工具推荐
  17. 顺丰下单空运实际发陆运
  18. matlab画图高级,[转载]Matlab绘图系列之高级绘图
  19. wp主题模版xsbrand修改
  20. 风光储互补发电系统直流微网

热门文章

  1. 基于5G承载数通的PTN网络
  2. NXR75030H充电模块电源,直流充电桩整理模块
  3. 嵌入式Linux驱动开发笔记(未完待续。。。)
  4. java支付宝无密退款
  5. utools工具使用
  6. PD16 for Mac(支持M1芯片)v16.3.2(50531)中文版介绍
  7. 旧版本U232-P9 USB-Serial USB到串口转换器在Windows 10/11使用的问题
  8. 什么是委派模式(Delegate)?应用场景是什么?
  9. Web实时预览 界面组件Telerik——提高开发者工作效率的完美组合
  10. 当CComboUI控件处于不可见,SelectItem无效果的解决方法