1、脑成像的元分析研究(ynji)：基于坐标点的元分析方法,报告人：封春亮单位(dnwi)：北京师范大学，脑与认知科学研究院,第一页，共四十四页。,脑成像元分析的背景(bijng)简介,第二页，共四十四页。,脑成像数据(shj)元分析的发展,元分析(fnx)(Glass,1976),“将现有研究整合并对其结果进行综合分析的统计(tngj)方法”,第三页，共四十四页。,脑成像数据(shj)元分析的发展,(胡传鹏等,2015),脑成像研究(ynji)的数量,第四页，共四十四页。,脑成像数据(shj)元分析的发展,元分析(fnx)对脑成像研究的意义(Wager et al.,2007;Eickhoff
2、et al.,2009),单个脑成像研究的被试数量普遍较少，结果不够稳定；单个脑成像的结果很可能受到特定的实验操作（如，扫描(somio)参数）的影响；单个脑成像研究对某个脑区功能的解释往往局限于使用的单个或几个任务。(e.g.,vmPFC=self-reflection?),-单个脑成像研究的局限性,(Wager et al.,2007),第五页，共四十四页。,脑成像数据元分析(fnx)的发展,脑成像技术：where?(位置(wi zhi)信息),坐标信息(x,y,z),最大激活点的坐标(x,y,z),第六页，共四十四页。,脑成像数据元分析(fnx)的发展,脑成像的元分析(fnx)：对不同研
3、究中位置信息的整合,早期脑成像数据的“元分析”：对不同研究(ynji)中报告的位置信息的主观描绘,(Peyron,2000),(Joseph,2000),第七页，共四十四页。,Activation likelihood estimation(激活(j hu)似然性估计，ALE,Turkeltaub et al.,2002),目标(mbio)：最大化地量化研究间位置信息的一致性；尽可能减小分析方法的主观性。,脑成像数据(shj)元分析的发展,第八页，共四十四页。,ALE的计算(j sun)(Turkeltaub et al.,2002),概念：把每个激活点拟合成一个概率分布 也即这个(zh ge
4、)激活点落在大脑每个voxel上的概率大小：与激活点越近的voxel，概率越大。,脑成像数据(shj)元分析的发展,文献中报告的激活点,第九页，共四十四页。,ALE的计算(j sun)(Turkeltaub et al.,2002),脑成像数据元分析(fnx)的发展,ALE：至少有一个(y)激活点落在某个voxel内的概率。,相关研究中报告的所有激活点,第十页，共四十四页。,ALE的统计(tngj)检验(Turkeltaub et al.,2002),Permutation test(置换(zhhun)检验),零假设：激活点没有重合，所有重合都是随机因素导致的假设元分析采集到了100个坐标（激
5、活点）；则产生100个随机坐标，然后根据这些随机坐标计算出大脑中每个voxel在零假设下的ALE值（如，重复1000次，然后求均值）产生零假设分布；通过(tnggu)比较实际ALE值与随机坐标产生的ALE值，来确定显著性。,脑成像数据元分析的发展,第十一页，共四十四页。,ALE的统计(tngj)检验(Turkeltaub et al.,2002),Permutation test(置换(zhhun)检验),在H0分布下，ALE为5.7210-3的概率为0.0001。因此，如果设置=0.0001为显著水平，那ALE值大于5.7210-3的voxel都视为显著（也即至少有一个激活点落在这个voxe
6、l内）。,脑成像数据(shj)元分析的发展,第十二页，共四十四页。,ALE的优点(yudin)(Turkeltaub et al.,2002),自动化的分析(fnx)除了位置信息外，也有对一致性水平的量化(ALE值)具备显著性检验，得到的结论有统计上的支持,脑成像数据(shj)元分析的发展,第十三页，共四十四页。,早期(zoq)ALE算法的缺点：,以激活点为单位，在每个voxel上，求出所有激活点对应概率(gil)的联合概率（即ALE值）。忽略了研究间的变异，属于固定效应的统计推断(fixed-effects inference):元分析的结论不能推广到该元分析未包含的研究(胡传鹏等,2015
7、;Eickhoff et al.,2009)。报告了较多激活点的研究对元分析的影响更大，甚至主导元分析的结果(Wager et al.,2007)。,脑成像数据(shj)元分析的发展,第十四页，共四十四页。,ALE的改进(gijn)算法(Eickhoff et al.,2009),脑成像数据元分析(fnx)的发展,第一步：在每个实验的内部(nib)先做一次整合,类比：把每个实验看做一个“被试”；每个voxel看做一个实验条件；每个实验报告出来的激活点看做一个采集到的数据点。,第一步的整合“”对每个被试各个条件下的所有数据点求和。,第一步整合的结果：每个实验在每个voxel上，都有只一个MA/A
8、LE值。,第十五页，共四十四页。,脑成像数据元分析(fnx)的发展,第二步：以实验(“被试”)为单位(dnwi)，计算每个voxel上的总ALE值，并做统计检验。,结果(ji gu)：每个voxel的总ALE值,置换检验：零分布：从每个实验的MA图中各自随机抽一个坐标点和对应的MA值，根据这些MA值算出一个总ALE值。重复1011次，每次能得到一个总ALE值。构建出零分布。根据不同ALE值在零分布中出现的概率，确定显著性水平。,ALE的改进算法(Eickhoff et al.,2009),第十六页，共四十四页。,ALE算法(sun f)进一步的优化：,脑成像数据元分析(fnx)的发展,Eick
9、hoff et al.,2012(NeuroImage):以MA值为单位代替以voxel为单位，减少(jinsho)了迭代次数，加快了运算速度；Turkeltaub et al.,2012(Human Brain Mapping):取最大值代替求和的方法，进一步减少单个实验对结果的影响。,第十七页，共四十四页。,脑成像元分析(fnx)的步骤,脑成像元分析(fnx)的步骤,第十八页，共四十四页。,脑成像元分析(fnx)的步骤,第十九页，共四十四页。,1.GingerALE:https:/www.brainmap.org/ale/,(邓沁丽,2015),脑成像元分析(fnx)的步骤,使用(shyn
10、g)到的软件：,第二十页，共四十四页。,脑成像元分析(fnx)的步骤,使用(shyng)到的软件：,2.Mango:rii.uthscsa.edu/mango,用于结果(ji gu)呈现,第二十一页，共四十四页。,第一步：选好主题，查找和筛选(shixun)文献,1.1.选好元分析(fnx)的主题,工作记忆(Nee et al.,2013);Go/No-go 任务(2008);奖赏加工(ji gng)(Bartra et al.,2013);社会决策(Feng et al.,2015).,共情与自身疼痛感受是否 基于共同的神经机制(Lamm et al.,2011);自传体记忆、展望未来、心理
11、理论以及默认网络是否基于共同的神经机制(Spreng et al.,2009);初级和次级奖赏的加工是否基于共同的神经机制(Sescousse et al.,2013)自己接受奖赏与观看别人接受奖赏是否基于共同的神经机制(Morelli et al.,2015);,“共性”,脑成像元分析的步骤,单个任务对应的神经机制：,多个任务共同的神经机制：,第二十二页，共四十四页。,第一步：选好主题(zht)，查找和筛选文献,1.2.查找和筛选(shixun)文献,查找途径：google学术(xush)(https:/http:/scholar.glgoo.org/)PubMed(http:/www.nc
12、bi.nlm.nih.gov/pubmed/)ISI Web of Science(http:/已查找文献的引用列表 相关领域综述的引用列表,文献筛选：被试群体：如，正常成人被试；是否使用fMRI技术；使用的是全脑(whole-brain)分析，而不是感兴趣区(ROI)分析；报告了感兴趣的统计结果；激活坐标报告在标准空间里(如MNI，Talairach).,脑成像元分析的步骤,查找内容：与主题相关的几个关键词 fMRI/functional fairness/fair/ultimatum game,第二十三页，共四十四页。,第二步：采集(cij)坐标点,需要提供的信息（规定(gudng)的格式
13、）：标准空间的类型（第一行）：MNI或者Talairach;文献信息：/作者名，年份等；/Subjects=N(被试量的信息)坐标点。每个点一行 不同文献之间需要有一个空行；注意MNI与Talairach空间之间的转换。,脑成像元分析(fnx)的步骤,第二十四页，共四十四页。,第三步：使用GingerALE进行(jnxng)元分析,3.1.计算(j sun)ALE(activation-likelihood-estimation),脑成像元分析(fnx)的步骤,Eickhoff et al.,2012(NeuroImage):Turkeltaub et al.,2012(Human Brain
14、 Mapping):,内置在GingerALE软件里的算法,第二十五页，共四十四页。,第三步：使用GingerALE进行(jnxng)元分析,3.2.矫正(jiozhng)结果：Settings,脑成像元分析(fnx)的步骤,第二十六页，共四十四页。,第三步：使用GingerALE进行(jnxng)元分析,3.2.矫正(jiozhng)结果,关于FDR的bug(在2.3.3以后的版本中更正)FDR矫正的原理：假设有N个比较，对应有N个p值。i).所有p值从小到大排序：p(1),p(2),p(3).p(N)；ii).若想控制FDR不能超过q，则只需找到最大的正整数i，使得 p(i)=(i*q)/
15、N.然后，挑选对应p(1),p(2),.,p(i)的比较做为有显著差异(chy)的比较，这样就能从统计学上保证FDR不超过q。,脑成像元分析的步骤,第二十七页，共四十四页。,第三步：使用(shyng)GingerALE进行元分析,3.3.选择输出结果(ji gu)的路径,Tools Preferences,脑成像元分析(fnx)的步骤,第二十八页，共四十四页。,第三步：使用(shyng)GingerALE进行元分析,3.3.点击(din j)“Compute”，进行运算。,脑成像元分析(fnx)的步骤,第二十九页，共四十四页。,第三步：使用GingerALE进行(jnxng)元分析,3.4.查
16、看(chkn)输出文件,ALE图像(t xin):*_ALE.niiP值图像:*_P.nii矫正后的ALE图像:*_ALE_pID05.niicluster图像:*_clust.nii,3.4.1.图像文件,脑成像元分析的步骤,第三十页，共四十四页。,第三步：使用GingerALE进行(jnxng)元分析,3.4.查看输出(shch)文件,cluster表格(*_cluster.xls)参数和统计信息(*_cluster.txt)1.使用的参数；2.每个cluster内，更详尽(xingjn)的label信息 3.每个文献对各个cluster的贡献值,3.4.2.文本文件,脑成像元分析的步骤,第三十一页，共四十四页。,第四步：使用Mango软件呈现(chngxin)元分析结果(rii.uthscsa.edu/mango),4.1.打开模板(mbn)图像：,open open Image选择(xunz)推荐的模板：Colin27_T1_seg_MNI.nii Colin27_T1_seg_MNI_2x2x2.nii,脑成像元分析的步骤,第三十二页，共四十四页。,第四步：使用Mango软件