大家好,我在做SBM的anova的报错如图1,basic model 和estimate 都没有问题,到roi analysi 就报错了(图1),做两样本t检验就没问题,spm和cat的工具包也重新添加路径了,依旧是如上报错;另外在提取surface的指标的过程中,发现部分病人运行过程中是红色的,请问这个是错误吗?但是最后也显示done,担心是由于原始数据的问题,请大家帮我看下,谢谢!感激不尽!!
我确认一下你的分析步骤,你是做皮层的统计分析,先设置模型和估计模型,有查看结果(Results模块)吗?然后你再做的ROI的统计分析,是这样吗?
提醒一下,上传图片尽量避免一些隐私信息(比如被试姓名),以防可能的风险。
你好,谢谢提醒,我会注意的。是的,做的皮层的统计分析,设置和估计模型都成功了,然后设置矩阵:f检验,输入:
1 -1 0
1 0 -1
之后选择p为0.05,图谱选择为DK40,选择指标thickness,之后就出现这样的报错,烦请您指导一下这个要怎么解决
虽然你这里用的F contrast是可以的,但是好像不是CAT12使用手册上使用的。所以可以试一试。
对了,我忘了问了,vertex-wise的统计结果还是正常显示的,只是说ROI-wise的分析会报错,对吧?
老师你好,vertes水平的统计是正常的,不管是用您刚刚建议的矩阵还是 1 -1 0这样的矩阵,都可以成功;
然后按照您的建议更改了矩阵,报错和之前一样,图1-3;之后尝试了重新添加cat12至工作路径,报错改变了,如图四
另外老师看到您这个回复自己不太理解,我理解的cat12做anova只能得到三组有差异的脑区(一个脑图),不能体现两组间的差异,需要把这些差异脑区的值提出来,再做两两比较(比如箱图之类的),不知道和您的理解是否一样?谢谢您!
是的,有两种做法,一种是先做F-test,然后把差异脑区提取出来两两比较。另一种是,在voxel-level去做两两比较,如果是这种,我建议就不做F-test了。
目前我的问题是vertex水平结果不能通过矫正,所以才做的roi分析,三组的话是需要做anova,但是这个anova怎么一直报错?如果分开做两两的t检验会被argue吧?
你这里应该不存在多个版本的CAT吧(即使旧版本文件夹改名了也不行)?
CAT的报错我其实也看不懂,我自己的经验是,如果CAT报错,可能换一个版本就好了,因为CAT更新特别频繁,小bug不断,特别是关于surface的分析。
另外,你也可以考虑不用CAT做统计,其实如果是ROI的数据,我感觉用SPSS或者R要更容易。
最佳的解决办法还是去SPM邮件列表里提问,CAT的作者Christian Gaser的回复一般都很及时。
我觉得不会,你只要两两分开做的时候控制一下多重比较就行了。比如,如果是三组,任意两组做T检验,就比较了三次,显著水平卡在0.05/3就行了。
当然,你不用完全相信我说的,我只是根据我的理解给你提供一些思路,你可以多查一查其他来源的资料。这些统计方面的基本问题应该有很多资料可以参考。(你如果查到有不同的看法,也欢迎和我讨论)
是的,我不存在多个版本的cat,下载的是最新版本的;
对于roi的数据spss是可以的,但是理论上要先筛出有意义的脑区对吧?另外我想得到脑图,不知道老师还有其他可以做统计的软件建议吗?
spm邮件列表里我也会去提问的,辛苦老师!
另外,你的spm12版本也可以尝试换一下最新版本的。
对于ROI的分析,因为涉及到多重比较校正,确实还是应该先做F检验,并对F检验的结果进行校正。不过这个肯定是可以在SPSS里做的。
我自己是习惯于用R做统计分析。画图的话,可以尝试ggseg(https://ggseg.github.io/ggseg/)或者enigma toolbox(https://enigma-toolbox.readthedocs.io/en/latest/pages/12.visualization/index.html),前者基于R,后者可以用MATLAB。这两个工具的可视化效果都是文献里常见的。
老师,我的spm12和cat12都是重新下载的是用的最新版本
eniogma toolbox是那种输入图像就可以做统计分析的工具包吗?已经被这个统计搞了好几天了
好的,那可能的因素都排除掉了。
enigma toolbox我只用来做过可视化,其他功能没用过,也不清楚还有哪些功能。如果是nii的数据,每个常用软件的统计模块几乎都可以做,但是对于CAT12的皮层分析,我不知道除了CAT本身还有什么软件可以做。
从前面的讨论下来,我可能倾向于认为你遇到的错误是CAT12的Bug,除了SPM的邮件列表,你也可以直接给CAT作者写邮件,回复会更快一些。