Greetings fmrib,
I want to run false discovery rate thresholding on TBSS/Randomise output and
came up with the analysis steps mentioned below. When you have the time,
could someone double-check the following. Thank you yet again. Greg
After TBSS and Randomise using default parameters.
1.) Turn unthresholded tstat into logp image:
avwmaths tbss_tstat1 -mul 0 -add 1 dummyvar
ttologp dummyvar tbss_tstat1 1000
avwmaths logps -exp ptstat
2.) FDR:
fdr -i ptstat -m mean_FA_skeleton_mask -q 0.05
3.) Threshold p values according to FDR output
avwmaths ptstat -mul -1 -add 1 -thr {(1 - FDR OUTPUT)} resultname
|