大家好,请问做完TBSS 之后是用cluster 这个命令汇报结果吗?如果想知道是哪个纤维束或者脑区的白质具有统计学差异,应该怎么做呢?谢谢!
cluster用于生成显著的团块。如果想知道这些团块的属于什么结构或纤维束,需要使用atlasquery这个命令,详细用法参见FSL官网:https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Atlasquery
谢谢,我已经尝试过atlasquery 这个命令了,用的是cluster 得到的vox 坐标,但是atlasquery 总是说显著区域团块unclassified 这种情况该怎么办呢?
这说明FSL自带的分区模板里没有包含这个坐标位置的标签。一个解决方法是更换一个分区模板,FSL 6.0.2后有一个新的XTRACT atlas包含的白质区域更广泛一些。另外,你不用一个voxel坐标去查询,试试用整个cluster mask查询。
谢谢,我在FSL邮件列表也查到了一个相似的命令autoaq,并且采用您所说的XTRACT atlas之后的确有label出来了,然后我发现有的cluster比较大,3w+的体素,汇报结果的时候是不是如同xjview那般,除了峰值坐标所在的区域,该cluster中比较大的部分所在的白质区域也应该汇报呢?谢谢
是的。如果cluster比较大,只报告峰值会漏掉很多信息。
你在cluster 命令中加上 --mm选项,生成的坐标是毫米坐标,有些模板我用voxel坐标的时候就查不出来,用mm坐标就可以。
谢谢,我也是用mm坐标查出来了。我还有个疑问,不知是否可以解答。通常查询团块的坐标只会报告峰值坐标,但是我的团块很大,我也想汇报其他团块,请问如何查询非峰值区域的那些坐标呢?
我记得cluster命令可以输出质心坐标,就是加权平均后的中心坐标,你可以看一下这个。另外,如果你的cluster比较大的话,那应该比较少吧,可以用fsleyes视觉查看。导入标准控件图像和统计结果,在fsleyes里面打开atlas panel,就可以了。可以看一下我的博客哈哈 https://www.wetiqe.xyz/quirying-atlas-using-fsleyes-or-command-line-tool-f0d51934dae940729147d7f688604cd6#260e8ca6e12b49268eb14128d0cfe4c0
还有一点,fslwiki上面写的说,cluster很大的话只是代表这个cluster有区别,并不能说明cluster覆盖的地方都有区别。如果cluster很大的话,最好结合着原始数据就是tstat文件看一下
谢谢你的回复,中心坐标和峰值坐标我都试过了,我是觉得因为cluster太大,里面包含的信息很多,如果只是汇报峰值坐标或中心坐标所在的纤维束,那就会遗漏很多信息,所以我觉得里面包含比较多的体素所在的纤维束也应该汇报
您好,我的汇报结果是1个大的cluster,但是会在下面标注很多团块,后面的数值应该怎么理解呢?我看相关文献也有对年龄和性别进行协变量回归,但是没有找见对应的命令,向您请教。
没懂你的意思,能给个截图之类的吗?第二个问题能描述得详细一些吗?
你是说后面这些MAX什么的是什么意思吗?这些就是峰值和峰值坐标。如果cluster太大,可以提高一下显著性水平,让大cluster分裂成小cluster,这样报告结果会更清晰一些。
应该是属于这个脑区的概率,autoqc
调用atlasquery
,所以直接查看atlasquery
的文档即可。