Print

Print


Hello Raphael,
                           This may be due to a current issue with the -e option in randomise - we will hopefully be releasing FSL 4.1.8. with a patched version of randomise in the next few days, I would strongly recommend rerunning the analysis with this new version.

Many Regards

Matthew

> Dear FSL-experts,
> 
> I`m very surprised by the difference between two tbss analyses i ran. In fact, I forget to include the DTI images of one subject into a longitudinal analysis with 3 points in time and 2 groups. Thus I reran the analysis with altogether 13 subjects in each group and surprisingly, results differed extremly. Testing for a main effect of time and an interaction of groups with all 26 subjects included in the model all the skeleton shows up when looking at tfce-corrected fstat-images at p < .05 (> .95 in fsl-terminology). However, the same results look way more plausibel (two cluster at the expected location) when I only consider 25 subjects, that is to say with the subject that I initally forgot to incorporate in the analysis. I visually inspected the results of the preprocssing steps before subjecting data to tbss analysis as well as the results of the tbss analysis (e.g. the all_FA_skeletonized image, etc.), results look fine.
> 
> I attached pictures of the design matrices and contrasts as well as the exchangeability block (EB) or to be precise ".grp" file/s. As can be seen from the design files, I only added one subject and of course adjusted the EB file accorindgly. The first 26 columns of my design matrix code subject factors, the last four columns the experimental group and control group at T1 and T2. The columns are respectively labeled in the design matrix. The last four columns look somewhat scrabbled which is due to the odering of subjects in the all_FA file which makes the membership of subjects to groups most likely only obvious to me.
> 
> Randomise was used with the following commands:
> 
> randomise -i all_FA_skeletonized.nii.gz -m -o ... -m mean_FA_skeleton_mask.nii.gz -d ... .mat -t ... .con -f ... .fts -e ... .grp --fonly -n 5000 --T2 -V.1
> 
> I would vey much appreciate any help on this as I´m a little bit desperate on how to explain this differences and on how to proceed.
> 
> Thank you very much in advance.
> 
> Raphael
> 
> 
> <13vs12.zip><13vs13.zip>