我有个问题请教,还请您春节里抽空给予解答:单样本t检验,参数估计都很好,但结果显示时,选中FDR校正之后问P=0.05点回车确定,没有&extent threshold选项就直接蹦出下面的结果.
并且MATLAB command window里提示:
========================================================================
"Warning: No voxels survive height threshold",意思是在给定的阈值下,没有得到差异.你的单样本t检验,检验的具体是什么指标(越具体越好)?样本量多大?
非常感谢臧老师的回答,我们进行EPI序列扫描,了解海洛因成瘾者在静息态时脑内的激活情况,预实验已进行了十多人,分析数据时,我按照严超赣博士的视频(Data_Process_of_Rest_fMRI.mp4),首先利用SPM进行了slice timingrealignmentnormalization及smoothing(ReHo是后做的smoothing),再利用REST工具做的detrend和filter;这些预处理工作完成后,先后做了ALFF,ReHo和FC,将mALFF,smReHo减1后联合zFC经单样本t检验,参数估计过程都很好,但结果显示时,选中FDR校正之后问P=0.05点回车确定,没有&extent threshold选项就直接蹦出没有任何激活的玻璃脑结果,按理说做完FDR校正后应该是有&extent threshold选项的.&extent threshold我的理解应该是显示体素的阈值,没有设置该选项所以没有差异,不知能否这样理解,如果是这样请问怎么实现结果输出(产生静息态下的激活图)和统计检验.
注:我是为了验证自己做法是否正确,只先分析了3个被试(n=3).
非常感谢您能在春节期间抽时间解答.
3个被试做统计检验,几乎是不可能出结果的.请把十几个被试全部分析.请注意,ALFF的one-sample t-test结果中,脑池,外侧裂t值一般也会非常高.
感谢臧老师的回答,另外我想补充一点,&extent
感谢臧老师的回答,另外我想补充一点,&extent threshold选项,SPM设置的是只有在检验之后有体素存活的情况下,才会出现.您做的实验中经过FDR校正之后,没有体素存活下来,所以SPM就跳过了&extent threshold选项这一步.如果你仅仅是想做个验证,可以提高FDR的P值,比如将0.05改成0.5,或者更大,可以见到&extent threshold选项.