Reference slice in slice timing in SPM

I saw at least two web pages saying that the reference slice can be chosen during slice timing correction. For example, http://www.fil.ion.ucl.ac.uk/spm/doc/manual/spatial.htm. I disagree.

Let’s take the first volume as example. Assume TR=2 and number of slices are 31. During data analysis (GLM or ROI analysis), we assume that the first volume is at time 0. In reality, the first acquired slice in this volume is at time 0, the 2nd is at time 2/31, the 3rd 4/31, and so on. The 31st slice is acquired at time 60/31. So except for the first acquired slice, the timing of all other slices should be corrected (or shifted) to time 0. That means, only the 1st acquired slice can be the “reference slice”. Reference slice is not a choice in data preprocessing step.

In SPM’s convention, slice #1 is the most bottom slice which may or may not be the first acquired slice depending on the slice order. Assume the slice order is descending (meaning data are acquired from top to bottom). Then the first acquired slice is slice #31; and the last acquired is slice #1. So in this case, the reference slice should be 31, and the slice order is 31, 30, 29, …, 1.

Choosing a wrong reference slice can affect the timing of your signal. Below I plotted the fMRI signal time series in motor cortex in a finger tapping experiment (tapping for 15s, rest 20s). The top plot is for reference slice = 1, the middle plot is for reference slice = 31. For the bottom plot, I did not do slice timing correction at all. I simply plot the time series in this region (which is at slice #23, or 0.5 seconds away from #31). The vertical line indicates the onset of tapping instruction. Click to zoom the figure.

You can see from the above plot that the rising timing of the signal can be off by 2s if difference reference slice is chosen. Does this matter? It depending on if timing is important to your experiment. For a block design, it’s probably fine; but for event design, or when the question you address is about timing, then this is important.

(The only preprocessing done on the images are reslicing to the first volume. No normalizing, smoothing etc ).

Don't want to miss new papers in your field? Check out Stork we developed:

fNIRS Journal Club 视频

546个被试的大型实验是怎么做的？近红外超扫描技术如何揭开群体冲突的神经机制？北京时间2020年5月29日周五上午11点，北京师范大学的马燚娜教授（Yina Ma）为大家讲解了她们组刚刚在Nature Neuroscience发表的文章。下面是报告视频。第一个是在Youtube，第二个是在youku。 Youtube：https://www.youtube.com/watch?v=4qZ7zP-BGz4Youku：https://v.youku.com/v_show/id_XNDY5MTc0MjAxMg==.html 她讲的文献如下： Within-group synchronization in the prefrontal cortex associates with intergroup conflict. Nature neuroscience https://www.storkapp.me/pubpaper/32341541 中文摘要（人工智能翻译，仅供参考）：陷入群体的个人有时会失去自己的个性，冒着通常会避免的风险，以无端的敌对态度与外界接触。在这项研究中，我们确定了右侧背外侧前额叶皮层（rDLPFC）和右侧颞顶交界处（rTPJ）的组内神经同步是组间敌对性的潜在机制。我们将546个人组织为91个三对三人小组间比赛，诱导了组内亲和，并使用功能性近红外光谱仪测量了神经活动和组内同步。在组内亲和之后，个人给组内成员的钱比给组外成员的钱多，并且捐出更多的钱来击败竞争对手。组内亲和减少了rDLPFC的活动，并增加了rDLPFC和rTPJ之间的功能连接。尤其是在组外攻击期间，组内亲和还增加了rDLPFC和rTPJ中的组内同步，并且组内rDLPFC同步与组间敌对性正相关。组内同步减少前额活动可能可以解释组内联结如何导致对外界的冲动和集体敌视。

Remote fNIRS

Everytime I met Dr Leanne Hirshfield, I am impressed by her energy and passion with the fNIRS technology. It was a pleasure listening to...

fNIRS Journal Club 通知 2020/5/29, 11am

546个被试的大型实验是怎么做的？近红外超扫描技术如何揭开群体冲突的神经机制？北京时间2020年5月29日周五上午11点，北京师范大学的马燚娜教授（Yina Ma）将为大家讲解她们组刚刚在Nature Neuroscience发表的文章。 欢迎大家参加并参与讨论。 时间: 北京时间2020年5月29日周五上午11点 地点: https://zoom.com/j/84320310196房间号: 843 2031 0196 密码: 600516 她要讲的文献如下： Within-group synchronization in the prefrontal cortex associates with intergroup conflict....

7 Replies to “Reference slice in slice timing in SPM”

1. Ting Xiang says:

Agree. I also found that using a different reference slice changed my GLM results (event-related analysis).

2. Jian says:

well, I actually disagree on this issue. First of all, the ref slice selection is a subject preference. If you select the 1st or last slice as ref slice, you also need to do corresponding event timing adjustment when you specify the design matrix in SPM. This latter procedure ensures that the ref slice maps with the event timing well.
Also, give the slice-timing correction tends to get more inaccurate with longer time delays, when TR = 2s, for example, people tended to choose the middle slice so that the maximum correction is 1 sec not 2 sec if you choose the 1st or last slice as ref slice.

3. Xu Cui says:

@Jian
Yes, if you do event timing adjustment then you can choose whatever ref slice. But it seems many people are not aware of this.

about 2nd point: I thought about this and think there is no “better” ref point (slice). If you choose the middle slice as ref, then you need to adjust the 1st and last slice by 1s; if you choose the 1st slice as ref, you need to adjust the middle slice by 1s, but you don’t need to adjust the last one by 2s, because the last slice is very close to the 1st slice in the next TR so you only need to adjust by TR/n second (n is number of slices). The periodicity of the scanning, I think, makes no slice a better ref slice in terms of accuracy of temporal interpolation.

4. Jian says:

@xu
I am glad we can communicate this way! :), again, congrats for the new daddy!
as to the 2nd point, I still think choose the middle slice as ref slice would save people a lot of hassle since every slices with that specific TR is aligned well, and you only need to remember one lag. Instead, if you choose the 1st slice as ref slice, and decide to align the last slice with the first slice (ref slice) of next TR, then you need to worry about that for this specific slice, you have a different ref point, and I don’t think SPM or any conventional Statistic package allow you to deal with slices that are “acquired” at different time and automatically adjust that.

5. Ting Xiang says:

for the event timing adjustment in the design matrix, do you mean the hrf with time derivative? But, I was suggested not to use this option…

6. evp says:

@Ting Xiang
Some people adjust the onsets of events in their designs, but the correct SPM way for doing this is adjusting microtime resolution and onset. Sorry for posting at such an outdated discussion, but someone may occasionly find it on google just like me 🙂

7. Gaha says:

How is the microtime computed in SPM? Is it an upsampling via interpolation?