fMRI: how to do ROI analysis

42 sec read

ROI stands for region of interest. The region is predefined usually from glm contrast. For example, you find visual cortex is activated when a flash is presented through glm analysis, you may then do more analysis on this particular region.

ROI analysis usually means a plot of BOLD signal in this region against time. In some cases, you may also want to do correlation analysis. I will explain how to do these analysis below. I will also give an example.

Before you do any ROI analysis, you should have the following data in hand:
1. Preprocessed image files which contains the BOLD signals. The file names are usually of the format swaf*.img and swaf*.hdr.
2. The onset time of each event of interest.
3. The coordinates of the voxels in the region of interest, or a mask image.
4. If you want to find the correlation between two regions, you need another set of 3.
5. If you want to find the correlation between BOLD signal and the stimulus property (such as the intensity of flash), you need data of this property.

The following script is written in MatLab. It explains how to do basic ROI analysis. download roi.m

第三十期 fNIRS Journal Club 通知 2022/05/28,10am

三个臭皮匠顶个诸葛亮背后的机理是什么?浙江大学潘亚峰研究员将为大家讲解他们最近发布的一篇用三人超扫描研究合作学习的文章。热烈欢迎大家参与讨论。 时间: 北京时间2022年5月28日周六上午10点地点: https://zoom.com房间号: 869 5374 8016 密码: 383720 Pan, Cheng, Hu (2022) Three heads are better than one: cooperative learning brains wire together...
Xu Cui
6 sec read

第二十九期 fNIRS Journal Club 视频

光子在大脑中是如何散射的?如何用计算机模拟的方法了解fNIRS的测量基础?来自波士顿大学的程小峻博士为大家分享他们近期在 J Biomed Opt 发表的文章。 Youtube: Youku: 相关资源: 文献鸟(追踪科学文献)文献大分析(两分钟了解一个领域)
Xu Cui
5 sec read

第二十九期 fNIRS Journal Club 通知 2022/4/30,10:00am

光子在大脑中是如何散射的?如何用计算机模拟的方法了解fNIRS的测量基础?来自波士顿大学的程小峻博士将为大家分享他们近期在 J Biomed Opt 发表的文章。 热烈欢迎大家参与讨论。 时间: 北京时间2022年4月30日上午10:00地点: https://zoom.com房间号: 868 4070 8135密码: 655816 程小峻要讲解的文章如下: Cheng, Chen, Sie, Marsili, Boas (2022) Development of a Monte...
Xu Cui
9 sec read

8 Replies to “fMRI: how to do ROI analysis”

  1. I am attempting to run your ROI script, but I always receive errors for the line: coord = mask2coord(‘mymask.img’). Matlab does not undersand the mask2coord command or function. I have to ten .img files for the ROI. How can I find the coordinates so that they work with your script?

  2. Dear Xu,
    What is the difference of the ROI analysis between the method you introduced here and the way carried out marsbar?

  3. Hello Xu,

    I am doing a MSc internship in a neuroimaging lab and being asked to run over 100 ROI analysis for every patient we are scanning and also at group level. Our regions have all been defined a-priori by reading the literature in particular a meta-analysis of the disorder. Nonetheless these amount to over 50 activation region masks (sometimes multiple masks per region), and around 20 deactivations and 20 correlations with little overlap.

    My basic question is how many ROI analysis is it ok to run? Is there a guideline somewhere? Most studies seem to concentrate on no more than 10 regions. Is my supervisor being unreasonable and not rigorous enough? Is this frowned upon?

    I know this comment is not specifically on your program but I cannot email the SPM list because my supervisor will know who I am and I am trying to understand his reasoning behind ROI analysis (which he is not helping me understand). I would be really grateful if you could give me some advice or point me to someone who could.


  4. @Alex

    During data exploration phase, one can do a lot of analysis, including tons of ROIs. Of course, only a small fraction are included in the final manuscript. However, the more ROI you do, the more likely you will discover something false positive. You can use FDR correction though.

    Keep in mind the time spent on 1 ROI is not much different from the time on 100 ROI if you have a matlab script to do it automatically (Feel free to adopt my code to your own use). I would do what your advisor asked for, summarize the findings, and let your advisor decide which to include in the final manuscript.

    Good luck!

  5. Hi,
    I have 2 groups of data,controls & patients.
    I have fmri datas of them in auditory task.
    now I need extracting ROIs of them( related to activity),
    How can I do it?


Leave a Reply

Your email address will not be published. Required fields are marked *