SPM的统计分析基础——GLM模型

SPM的统计分析,操作相对简单,但其背后的原理却不容易。

image.png
                                    voxel-wise就是基于voxel水平的分析;ROI-wise是基于ROI水平分析

   在一系列的图像预处理之后,进入统计模型的建立及估计环节。 理论上,经过预处理的数据可以直接进行统汁分折,但是这样做会有很大的局限胜。我们需要考虑到实际信号(探测到的神经生理响应)包含有多方面的因素。如果我们可以把它分成感兴趣的部分、不感兴趣的部分以及误差项,这样以来,一方面可以用于复杂问题的研究,另一方面使有用信息更加突出。为了达到这个目的,SPM采用一般线性模型(GLM)对采集到的信号进行统计分析。

image.png

image.png

通俗的说,\color{red} {GLM基于这样一种假设}

  每个像素上的实验数据(用Y来表示,同一任务的时间序列或不同任务序列的),是一些未知参数(用x表示)的线性组合(用β表示未知参数x的系数)。这些未知参数与任务或时间有关、但与具体脑区(像素)无关的已知参数,它组成的矩阵通常又称为设计矩阵(\color{red}{Design matrix})。

  在经过GLM描述之后,分析的对象就发生了\color{red}{转变}。原本是对Y的统计分析,现在转向了拟合得到的β。SPM得到的脑功能激图实际上就是根据对参数β的统计推断而得到的。

\color{#4285f4}{对参数β的统计推断究竟是怎么一回事? }

  • 首先,是检查设计矩阵的正确性;
  • 其次,使用最小二乘法法拟合出参数β
  • 再次,对β系数做t检验或者F检验
  • 最后,根据t检验或者F检验统计量和阈值,对β系数做统计推断。

根据t或F值以及相应的阈值,可以得到该阈值对应置信度的脑功能激活图。由于它是对模型的参数做统计分析而得到的,所以被称为\color{red}{统计参数图}。与阈值对应的t或F值记为SPM{t},SPM{F}。

image.png

SPM还提供了与t或F值及它们相应的自由度所对应置信度的{0,1}正态分布的参数值z,这样也得到关于z图像,记为SPM{z}。此外,根据这些激活图,SPM还可以输出脑血流等生理信号的响应曲线以及不同状态的响应特性等。

t检验、F检验与β系数的关系 Null hypothesis:β=0.假定如果某一体素的检验结果达到显著性水平,我们有理由拒绝Null hypothesis。这就意味着这个点的β系数在组间存在异常差异,因此就表明该体素的实验数据Y值存在显著性差异。这里在逐个体素逐体素地进行t检验、方差分析等统计之后,SPM会进行一个多元校正的分析,常用的方法就是FDR、FWE校正。

\color{red}{后话:这一段我也看不懂,抄上来一起看看}
概率论里面有两个主要的研究内容,参数估计和假设检验。这里用到的是假设检验的内容。 这里补充一点假设检验的基础知识  假设检验的基本任务是根据样本所提供的信息,对未知总体分布某些方面的假设做出合理解释。假设检验的程序是,先根据实际问题的要求提出一个论断,称为零假设(null hypothesis)或原假设,记为H0(一般并列的有一个备择假设(alternative hypothesis),记为H1 )然后根据样本的有关信息,对H0的真伪进行判断,做出拒绝H0或不能拒绝H0的决策。  假设检验的基本思想是概率性质的反证法。概率性质的反证法的根据是小概率事件原理。该原理认为“小概率事件在一次实验中几乎是不可能发生的”。在原假设H0下构造一个事件(即检验统计量),这个事件在“原假设H0是正确的”的条件下是一个小概率事件,如果该事件发生了,说明“原假设H0是正确的”是错误的,因为不应该出现的小概率事件出现了,应该拒绝原假设H0 。

假设检验有两种方法:置信区间检验法(confidence interval approach)和显著性检验法(test of significance approach)。 显著性检验法中最常用的是t检验和F检验,前者是对单个变量系数的显著性检验,后者是对多个变量系数的联合显著性检验。实际上t检验法与置信区间法提供的结果是完全一样的。

假设检验的错误: 如果有一个零假设在5%的显著性水平下被拒绝了,有可能这个拒绝是不正确的,这种错误被称为第一类错误,它发生的概率为5%。另外一种情况是,我们得到95%的一个置信区间,落在这个区间的零假设我们都不能拒绝,当我们接受一个零假设的时候也可能犯错误,因为回归系数的真实值可能是该区间内的另外一个值,这一错误被称为第二类错误。 在选择显著性水平时人们面临抉择:降低犯第一类错误的概率就会增加犯第二类错误的概率。  P值度量的是犯第一类错误的概率,即拒绝正确的零假设的概率。P值越大,错误地拒绝零假设的可能性就越大;p值越小,拒绝零假设时就越放心。

voxel-wise就是基于voxel水平的分析;ROI-wise是基于ROI水平分析

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 210,914评论 6 490
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 89,935评论 2 383
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 156,531评论 0 345
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,309评论 1 282
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,381评论 5 384
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,730评论 1 289
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,882评论 3 404
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,643评论 0 266
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,095评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,448评论 2 325
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,566评论 1 339
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,253评论 4 328
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,829评论 3 312
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,715评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,945评论 1 264
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,248评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,440评论 2 348

推荐阅读更多精彩内容