fMRI model specifictaion

GLM based

包括以下步骤:①明确GLM设计矩阵;②用经典或贝叶斯方法估计GLM参数;③利用对比向量检查结果,生成统计参数图(SPMs)以及后验概率图(PPMs)。

Design matrix:每行对应一个scan,每列对应一个影响变量或解释变量(回归量或解释变量)。关于GLM(一般线性模型)网络上有大量相关文章,

对于多被试分析,SPM通常分为两个阶段(level):

① 1st-level analysis:个体分析(within-subject)。通常有多少被试就要进行多少1st-level分析。使用 Specify fi rst level 和 Estimate 选项。这些选项的结果可以被用于做 case study。

② 2nd-level analysis:RFX analysis(Random effect analysis)。从被试推断被试所属主体的情况。

1st level analysis:

①Conditions:

name(1*n cell):condition name;

onsets(1*n cell):每个cell中为一保存condition开始时间的时间向量,每列对应一个condition,每行对应其出现时间(补充,心心水滴csddzh:任务刺激启动的扫描数(如1:14:70 ,代表任务从第 1个TR 开始,每14个TR 为一个周期,共70个TR));

durations:事件时长。

一般这些信息都保存在condition.mat文件中,使用的时候直接load进来,使用Multiple Conditions.

②Regressors

设计矩阵中用来和血液动力响应做卷积的一列。之后的leave-one-out实验要用。

暂时还没搞清楚是不是要用Multiple Regressors。

③Factorial design

设定该选项则SPM 可以直接生成主要效应和交互作用之间必要的test:F-test(被试个体(1st level)与被试间(2nd level)的contrast image)和T-test。

根据蜡笔的视频:被试根据三种商品颜色选择买或者不买,那么就有两个factor:color,3;choice,2。

SPM 默认第一个因素变化最慢,第二个因素变化第二慢...以此类推。

Factorial design一般在有复杂情况的时候才用,简单的时候可以不用设置。

其他因素没怎么看,因为在蜡笔的视频里都是默认。

=============================================

fMRI model estimation

模型参数可以用经典(ReML - Restricted Maximum Likelihood)或贝叶斯算法估计。

蜡笔视频中这一步直接使用了上一步生成(?)的SPM.mat之后用默认运行。

生成一SPM.mat,包含设计和估计的模型参数。

看结果:主面板Result。

=============================================

Factorial design specification(2nd level analysis)

配置设计矩阵,描述统计分析所需要的GLM,数据及必要的参数(保存在GLM中)。用于为2nd level analysis做准备。

Specify 2nd level analysis。

SPM12manual,统计部分(8-10)笔记的更多相关文章

  1. solr统计只返回10或者100个数据的解决办法

    因为我所在的公司为政府做的项目[风险管理系统],其中涉及大量的统计展示,多数以整个市的区划,行业部门等方式返回,在昨天,我发现听过填报单位的方式返回时,始终只有100个数据.通过对比发现,在前辈的代码 ...

  2. ubuntu 12.10 笔记

    笔记 more ec_unitouch.log |grep Thread-4 筛选日志 打开命令行终端 ctrl + alt + t     查看版本号 : sudo lsb_release -a t ...

  3. 斯坦福iOS7公开课10笔记及演示Demo

    这一节主要介绍了多线程中的串行队列以及滚动视图UIScrollView. 1 .多线程 这一节只是简单介绍了多线程的串行队列,即把任务加入线程队列后按顺序逐步执行. (1)目前iOS多线程提供的方法主 ...

  4. R语言与医学统计图形-【10】ggplot2图形映射

    ggplot2绘图系统--图形映射 颜色的映射. #aes中映射变量 ggplot()+geom_point(aes(x=carat,y=price,color='blue'),#color视为单一变 ...

  5. SQL 笔记 By 华仔

    -------------------------------------读书笔记------------------------------- 笔记1-徐 最常用的几种备份方法 笔记2-徐 收缩数据 ...

  6. 分享一下我研究SQLSERVER以来收集的笔记

    分享一下我研究SQLSERVER以来收集的笔记 前言 为什麽分享??因为像现在网上很多人攻城师那样,转行去卖水果,卖早餐,总有一日我也会离开这个行业的 由于本人不是在大公司上班工资很低,我希望有一天存 ...

  7. SQL统计

    --按周统计SELECT TOP 10DATENAME(year,AddDate) 年,DATENAME(week,AddDate) 周,COUNT(1) 单量,SUM(total) 总金额,AVG( ...

  8. C#+JQuery+.Ashx+百度Echarts实现全国省市地图和饼状图动态数据图形报表的统计

    在目前的一个项目中,需要用到报表表现数据,这些数据有多个维度,需要同时表现出来,同时可能会有大量数据呈现的需求,经过几轮挑选,最终选择了百度的echarts作为报表基础类库.echarts功能强大,界 ...

  9. HotApp小程序统计之自定义事件统计

    什么是自定义事件统计     官网:https://weixin.hotapp.cn/document 自定事件,就是自定统计任意事件的执行,灵活度最高. 用上图的云笔记说明想知道如下信息 (1)多少 ...

  10. 高性能Linux服务器 第10章 基于Linux服务器的性能分析与优化

    高性能Linux服务器 第10章    基于Linux服务器的性能分析与优化 作为一名Linux系统管理员,最主要的工作是优化系统配置,使应用在系统上以最优的状态运行.但硬件问题.软件问题.网络环境等 ...

随机推荐

  1. 转:查看sql语句执行时间/测试sql语句性能

    原文出处:http://www.cnblogs.com/qanholas/archive/2011/05/06/2038543.html 写程序的人,往往需要分析所写的SQL语句是否已经优化过了,服务 ...

  2. Hibernate 缓存介绍

    Hibernate中提供了两级缓存,一级缓存是Session级别的缓存,它属于事务范围的缓存,该级缓存由hibernate管理,应用程序无需干预:二级缓存是SessionFactory级别的缓存,该级 ...

  3. vmware克隆虚拟机eth0网卡无法启动

    概述: 通过vmware克隆安装好的虚拟机之后,出现了网卡未启动的问题. vmware安装虚拟机请看:<vmware快速安装linux虚拟机>. 定位过程: 1.通过ifocnfig命令只 ...

  4. 一条诡异的insert语句

    问题背景 有同事反馈在mysql上面执行一条普通的insert语句,结果报错, execute failed due to >>> Incorrect string value: ' ...

  5. XML 在SQLServer中的使用

    SQL Server对于XML支持的核心在于XML数据的格式,这种数据类型可以将XML的数据存储于数据库的对象中,比如variables, columns, and parameters.当你用XML ...

  6. 优化SqlServer--数据压缩

    数据压缩是对存储和性能优势的加强.减少数据库占用的磁盘空间量将减少整体数据文件存储空间,在一下几个方面增加吞吐量: 1.更好的I/O利用率,每个页面可以读写更多的数据. 2.更好的内存利用率,缓冲区可 ...

  7. Java代码获取NTP服务器时间

    apache的commons-net包下面有ntp相关的实现类,主要类有: 1  org.apache.commons.net.ntp.NTPUDPClient ? 1  org.apache.com ...

  8. 多元线性回归 ——模型、估计、检验与预测

    一.模型假设 传统多元线性回归模型 最重要的假设的原理为: 1. 自变量和因变量之间存在多元线性关系,因变量y能够被x1,x2-.x{k}完全地线性解释:2.不能被解释的部分则为纯粹的无法观测到的误差 ...

  9. python中mysqldb的用法

    1.引入MySQLdb库 import MySQLdb 2.和数据库建立连接 conn=MySQLdb.connect(host="localhost",user="ro ...

  10. iOS上传App Store报错:this action cannot be completed -22421 解决方案

    最近swift项目升了xcode8,提交版本时,遇到这个: this action cannot be completed -22421 瞬间懵逼,连具体报错原因都没有,只有一个代码 22421,找了 ...