原创,转载请指明出处!
首先找到一篇微博,理解了下BOLD signal change http://blog.sina.com.cn/s/blog_92f37ae20102vsqp.html
翻译一下就是:
- 找个每个ROI 的平均 timecourse
- 找到baseline 的平均timecourse
- 用ROI 的timecourse / baseline timecourse
- 每个条件下 onsets 的时间点的 平均这个条件的开始的时间点的 平均timecourse
- 遍历时间 差不多30秒的样子 其实20秒差不多了
- 例子
然后翻开论文看看人家是怎么画的图
人家的图是这样的
还有这样的
那么我们就要学习一下人家怎么 求 怎么画图
这里还要让 PyMVPA教程还要 上场 http://www.pymvpa.org/ 我的垃圾翻译版本 https://blog.csdn.net/cocosion/article/details/95978031
教你把你的nii 4D文件倒进来(三维空间中的点 加一维的时间 )时间可以自己写个文件导入具体volume下对应的条件是啥,
看了人家的教程 抽个时间点算(建议先算4秒的,依据HRF函数) ,把 baseline 和 具体你要的条件下的 在4秒下的事件全部抽出来,先做平均,再利用下面的公式算 signal change,算出来就是第一幅图人家那个样子;再遍历时间点即可画出第二张图人家那个样子(2秒算一次)
公式:
其实可以不用 乘100-100,直接 做除法之后-1 即可
参考文献即链接再附一遍:
[1] BOLD signal change http://blog.sina.com.cn/s/blog_92f37ae20102vsqp.html
[2] 图片1来源:Astrid M. Rosenthal-von der Pütten, Krmer, N. C. , Maderwald, S. , Brand, M. , & Grabenhorst, F. . (2019). Neural mechanisms for accepting and rejecting artificial social partners in the uncanny valley. The Journal of Neuroence, 39(33).
[3] 图片二以及公式来源:Albers, A. M. , Meindertsma, T. , Toni, I. , & De Lange, F. P. . (2017). Decoupling of bold amplitude and pattern classification of orientation-selective activity in human visual cortex. Neuroimage, S1053811917307942.
[4] PyMVPA教程来源:http://www.pymvpa.org/