人类进化历程中的近期自然选择:如何从分子遗传学中找证据-VG图的模式识别
VG图的模式识别
freesurfer是一个处理大脑3D结构像数据,进行自动皮层和皮下核团分割的工具,用起来非常方便。freesurfer wiki上的教程也非常详细,但是有一点,freesurfer的命令很复杂,很难准确地记住每个参数该怎么设置。本人比较懒,不愿记,也记不住,每次都需要打开wiki进行对照。由于wiki非常详尽,每次都是在一大篇英文中搜索命令。在这里弄一个简洁版,只把分析流程所用到的命令贴在这里,以便查阅。
刚从同学那里了解到spss.18之前的版本不能自动进行单因素多水平重复测量的Friedman检验,需要用spss syntax代码实现。而spss.18后就可以自己算了。在网上找到一段代码,并且也没有介绍怎么运行,可能对spss syntax代码不了解的同学用起来还是有些不便。自己琢磨了一下,所以在这里给大家一起分享了。
前面提到从scan中下来的数据需要用dcm2nii或者MNIConvert把原始的dcm数据转化为NifTI的hdr/img或者nii格式。再用dcm2nii转化的时候有一个选项orient,可以自动对3D数据进行方向的转换,转到与标准空间最相似的方向。在脑成像数据处理,图像的方向是非常重要的,同时也是很容易弄错的。数据处理时,一定要对数据的方向有很清楚的认识。目前常提到的方向有Neurological和Radiological两种方向。
在fsl和freesurfer的脑成像数据处理中,常常会用到sed命令,包括字符替换功能,提取文件内容等。比如说处理功能像数据的时候,经常在Feat中用一个被试的数据生成一个design.fsf,然后替换被试名字,生成所有被试的design_sub***.fsf文件,然后用for循环提交所有命令。有时候,又需要一行一行的读取某一个文件。所以sed非常实用,功能强大,但是也非常复杂,要完全掌握很难,一些常用的实用功能有必要了解,对于数据处理有很大的帮助。
在用 fsl和freesurfer处理数据的时候,除了需要看log文件是否有报错外,还常常需要check一下output文件数目是否对,复制数据和文件夹,文件和文件夹重命名等,并把这些结果输出到文本文件中。