VBM后的双样本t检验

article/2025/10/16 0:31:15

上一篇文章写到做完了VBM,做完后因为数据一般都是患者组和HC组,这两个组之间需要进行比较,那么我们就要进行双样本t检验。
这里介绍双样本t检验的做法。
依然使用的是SPM-fMRI。

1.第一步是选择Specify 2nd-level

在这里插入图片描述

打开以后我们可以看到这个界面,Directory就选择自己准备保存数据的文件夹。
注意:我们要做双样本t检验,打开这个界面我们会发现,Design下面默认的是单样本t检验,我们点击Design,选择Two-sample t-test
选择完以后,就会发现有以下这些需要选择。Group 1 scans和Group 2 scans就分别选择HC组和患者组的灰质图像即可。就是上一篇文章中最后得到的smwc1xxxx图像。
在这里插入图片描述
之前也提到了一般是将颅内体积作为协变量来进行t检验。那么这里的Multiple covariates就需要选择New:Covariates来添加一个协变量。这里的File(s)选择时可以是.mat或者是.txt文档。我这里使用的是.txt.

txt文件的创建:
这个txt文件里应包含所有被试的颅内体积,例如前面共选择了10例HC和10例患者,那么.txt文件里应该有20个数据,并且注意顺序,应该与之前的选择顺序一致。HC在前或者患者在前。
在这里插入图片描述
其余的不用动,运行以后会得到一个 SPM.mat文件

2.estimation

对上一步生成的.SPM文件进行估计,从首页选择estimate或者是直接在batch上方找到都可以。
在这里插入图片描述
在这里插入图片描述
这一步很简单,直接在select SPM.mat选择上一步得到SPM.mat文件就可以。
这一步可能会弹出是否要覆盖之前得到的SPM.mat,选择继续就可以。

3.对比

这个我在首页没有看到,在batch里可以直接选择。在这里插入图片描述
SPM.mat的选择就选上一步estimate后的结果。在Contrast Sessions里选择T-contrast。
Name可以随便命名,我建议是你比较的是哪两组就写哪两组的名字,例如:HC-xx患者这种。
重要的是Weghts vector这里,我是输入的[-1 1]。这里的意思实际上就是怎么比较,
如果输入【1 -1】,得到的结果如果是正值就是A-B>0,以灰质体积为例,就可以理解为B组的灰质体积相对A组减少。如果是负值,那么意味着A-B<0,理解为B组的灰质体积相对A组增加。个人认为这样不清楚,不如写【-1 1】,结果为正就是体积增加,结果为负就是体积减少。当然这也和A和B组的顺序相关,看个人怎么选择。
在这里插入图片描述
其余的就没有什么了,可以运行得到结果。
这个spmT图就是最终我们得到的双样本t检验结果。这是一个nii图像,可以使用xjview打开,然后对其进行FDR校正或者FWE校正,看自己需要。xjview也可以report激活的脑区,还是很方便的。对图像进行校正后,可以使用xjview把图保存下来,再使用brainnet viewer看看,也可以得到比较好看的图。至于brainnet viewer的使用我前面有写过,可以参考一下。
在这里插入图片描述


http://chatgpt.dhexx.cn/article/HDXnPbCO.shtml

相关文章

VBM后的配对t检验以及xjview使用

之前写了VBM后的双样本t检验&#xff0c;再记录一下配对t检验。 配对t检验和双样本t检验的过程基本一致。包括以下三个步骤。 第一步输入两组被试时&#xff0c;应该成对输入&#xff0c;共有几个被试就有几个pair。 但是这里我在做的过程中没有加协变量&#xff0c;不知道会不…

Visual Basic

目录 一&#xff0c;Visual Basic 二&#xff0c;控制台程序 三&#xff0c;可视化程序 1&#xff0c;IDE 2&#xff0c;实例——加法计算器 一&#xff0c;Visual Basic Visual Basic是可视化的Basic&#xff0c;简称VB VB是第一个可视化编程语言。 二&#xff0c;控制…

VBM法MRI图像处理——记第一次使用cat12

1.环境 MATLAB 2015b SPM12 CAT12 2.SPM部分 命令行输入 spm 出现 以及 点击Toolbox 出现 3.CAT部分 点击上图 设置请根据自己需求 多分割了一种surface皮层数据&#xff0c;当做皮层统计分析SBM时需要提取surface皮层指标时会用到。 我本意只是获得灰质、白质的体积…

VBM_DARTEL算法对灰质变化的计算

根据一些文献得知&#xff0c;VBM目前比较新的算法是DARTEL算法&#xff0c;这一算法被集成在SPM里&#xff0c;这里记录一下做法。 VBM是对T1像进行分割得到灰质等。所以要有结构T1加权像数据。 整个流程应该是这样&#xff1a; 1.手动调整前联合&#xff08;AC&#xff09; …

基于cat12搞定VBM的ROI分析——vertex水平和ROI水平的双样本T检验

前言 本来上周要更新此篇的&#xff0c;但由于本身有问题没有解决清楚&#xff0c;再加上导师给了数据处理的任务下来了&#xff0c;两下耽搁&#xff0c;就等到现在了。上回说到&#xff0c;做了VBM和SBM的指标提取及双样本T检验的统计分析&#xff0c;那接下来我们还可以做什…

基于cat12和SPM12进行VBMSBM数据分析笔记1——数据预处理

前言 今年是小编步入研究生生活的第一年&#xff0c;研究方向待定&#xff0c;但主要以磁共振成像为主&#xff0c;以后会不断地总结这方面的知识&#xff0c;涉及MRI&#xff0c;数据分析基础方法理论&#xff0c;软件操作教程&#xff0c;文献阅读笔记总结等等&#xff0c;用…

基于cat12和SPM12进行大脑VBM数据分析笔记2——统计分析

小白一顿操作猛如虎——拿下VBM的双样本t检验 前言 培训的deadline到了&#xff0c;小编需要交作业了&#xff0c;所以就在此先把笔记记下来。作业需要提交一个报告&#xff0c;but这个报告我应该如何写呢&#xff0c;培训总结&#xff1f;哎&#xff0c;现在还是不清楚&…

磁共振影像分析之: 基于FSL的VBM分析(1)

磁共振影像分析之: 基于FSL的VBM分析(1) 基本概念 在开始之前, 先给出几个基本概念: 磁共振成像 磁共振成像(Magnetic Resonance Imaging, MRI)是一种基于核磁共振原理和射频波探测组织结构和功能的医学影像技术. 值得一提的是, 磁共振成像技术没有辐射, 对人体没有影响. …

MRI图像处理:VBM原理和步骤

VBM是voxel-based morphometry的缩写&#xff0c; 是对被试之间灰质体素粒度统计分析。VBM可以得到人群中volume和gyrification的不同&#xff0c;对clinical score进行相关性分析&#xff0c;VBM不仅仅可以用于结构数据&#xff08;T1&#xff09;&#xff0c;还可以用在fmri和…

MRI脑影像分析从哲学到技术:一文搞懂VBM预处理基本原理(全网最详细解析)

基于体素的形态学方法&#xff08;voxel-based morphometry, VBM&#xff09;&#xff0c;是分析大脑解剖学&#xff08;结构&#xff09;差异最常用方法之一。 其通过给大脑volume逐体素打标签&#xff08;分类&#xff09;的方式来进行组织分割&#xff08;segmentation&…

LabVIEW程序测试

LabVIEW程序测试 工程师经常不太关注测试&#xff0c;将更多时间用于其他开发。通过一定程度的测试&#xff0c;可以保证节省时间。 开发人员必须清楚地了解对测试的期望程度。此外&#xff0c;还必须标准化测试方法并跟踪测试结果。在开发需求和设计规范时&#xff0c;还要制…

史上最详细的小程序测试用例模板,以后再也不怕测试用例不会写了

目录 1. 文档介绍 1.1 文档目的 1.2 文档范围 1.3 读者对象 1.4 术语与缩写解释 2. 功能测试用例 2.1 被测试对象的介绍 2.2 测试范围与目的 2.3 测试环境与测试辅助工具的描述 2.4 功能测试用例 附录&#xff1a;评审意见 1. 文档介绍 1.1 文档目的 设计测试用例…

软件测试人员一定要会的微信小程序测试点

微信小程序&#xff1a; 不需要下载安装即可使用的应用&#xff0c; 实现了应用“触手可及”的梦想&#xff0c;用户扫一扫或者搜一下即可打开应用。 体现了“用完即走”的理念&#xff0c;用户不用关心是否安装太多应用的问题。 应用将无处不在&#xff0c;随时可用&#…

小程序测试流程

流程分为两个方面&#xff0c;解包可以挖掘信息泄露问题、隐藏的接口&#xff0c;抓包可以测试一些逻辑漏洞、API安全问题。两者结合起来就可以边调试边进行测试&#xff0c;更方便于安全测试。 搜索目标小程序 目标搜索不能仅仅局限于主体单位&#xff0c;支撑单位、供应商、…

【软件测试】

系列文章目录 文章目录 系列文章目录前言第四章 单元测试4.1 软件测试过程概述4.2 什么是单元测试4.2.1 单元测试的定义4.2.2 单元测试的重要性4.2.3 单元测试原则 4.3 单元测试的目标和任务4.3.1 单元测试的目标&#xff1a;单元模块被正确编码4.3.2 单元测试的主要任务 4.4 单…

软件测试之测试程序开发

软件测试培训之测试程序开发 脚本维护&#xff0c;测试数据准备&#xff0c;测试监控准备; 测试程序开发 性能测试场景设计和性能测试脚本设计中需要注意以下几个问题: 场景真实性: 每个脚本的具体操作步骤,是和真实环境操作类似的,每个场景内的测试程序之间的逻辑关系,压力比重…

程序员的测试课

git项目地址&#xff1a;GitHub - dreamhead/geektime-todo: Geektime Todo is a demo todo project for Geektime column. 1、实现一个Todo应用 设计规范 1、对于输入参数的检测&#xff0c;由入口部分代码进行处理。如空字符串。 2、Repository 的问题以运行时异常 的形式抛…

软件测试重点知识

软件测试的定义&#xff1a; 利用人工或自动化的手段来运行或测定某个系统&#xff0c;其目的在于测它是否满足规定的需求或检测预期结果与实际结果之间的差异。 软件测试的目的&#xff1a; 发现软件中存在的问题或缺陷&#xff0c;检查软件系统是否满足要求。 软件测试的对象…

跑批利器--读取文件

在上篇文章中已经对示例的基本业务和流程进行了解,同时也知道了SpringBatch的三个接口ItemReader,ItemProcessor,ItemWriter.接下来我们主要是通过示例来具体处理一个文本文件,将里面的数据和内容通过SpringBatch的加工处理来存到新的数据库中. SpringBatch提供专门的类去读取文…

金融任务实例实时、离线跑批,Apache DolphinScheduler 在新网银行的三大应用场景与五大优化...

在新网银行&#xff0c;每天都有大量的任务实例产生&#xff0c;其中实时任务占据多数。为了更好地处理任务实例&#xff0c;新网银行在综合考虑之后&#xff0c;选择使用 Apache DolphinScheduler 来完成这项挑战。如今&#xff0c;新网银行多个项目已经完成了实时与准实时的跑…