Update: I used SPM to coregister my lesion maps to the MNI 152 template in the bcbtoolkit directory and everything worked like a charm!
Thanks again Mich for your help!
1 Like
My pleasure. Happy to hear that it is working. Always glad to receive feedback.
Regarding the templates, I’ve been in the same situation as yours, here is my summary of the MNI templates story:
While the coordinates system is supposed two be roughly the same between one MNI and another, some small differences typically crash analyses. We are using the MNI from FSL, which is an older version of the MNI from the Montreal Neurological Institute (http://www.bic.mni.mcgill.ca/ServicesAtlases/ICBM152NLin2009). While the ICBM template looks better, we prefer to use its older version because not everybody has a top-notch MRI in the stroke community, and an older template might be more representative of your T1 images.
SPM uses another MNI template cropped and smoothed for various reasons. The cropping is an issue as the images register in SPM MNI will not match the size of the images of tractotron or disconnectome templates. Finally, the Colin template from MRIcron is a single brain registered to the MNI (the brain of Colin Holmes!) but is missing one slice compared to the MNI we use (which will crash the analysis as well).
We could have modified the code so that any MNI works, but results might have subsequently been biased in some cases because of errors in the alignment.
I hope this helps!
Kind regards
michel
Hi!
I am trying to use the BCBtoolkit to create disconnectome maps, but the program is producing an error that says it cannot read my roi files. I am using a lesion in MNI 2x2x2mm space. How can I resolve this error?
Thanks,
William
Hi William
Thanks for your email.
BCBtoolkit is optimised for 1mm MNI
So you can either reslice your lesions to 1x1x1mm (the easiest, I can even do it for you if you don’t know how to) or I can provide you with the disconnectome 2x2x2 mm.
I usually use fsl and bash for this:
flirt -interp nearestneighbour -in M2006_lesion.nii -ref MNI152_T1_1mm_brain.nii.gz -out wM2006_lesion.nii -usesqform -applyxfm
If you use something else be especially careful with the interpolation rule you use. E.g. trilinear will give you continuous values. You want to use nearestneighbour in order to keep your lesions images binary.
Hope this makes sense.
We can zoom this out if you run in any trouble.
Hello,
I have a few questions related to the AnaCOM2 software in BCBtoolkit and the clever way you utilized it in your paper “Advanced lesion symptom mapping analyses and implementation as BCBtoolkit”.
In your paper, you use disconnectome maps instead of lesion masks as input when running AnaCOM2 to identify clusters associated with the performance in a category fluency test. I wish to use a similar approach for our sample of 328 MRI-scanned MS patients in relation to various clinical and cognitive measures.
- Is a comparison with a control group a requirement for this analysis even if I only wish to run comparisons between disconnected vs. spared?
- If inclusion of a normative control value is required (we do not have scores from healthy controls), how is this normative value included in a disconnected vs. spared comparison?
- Finally, a technical support question: When attempting to run AnaCOM2 it immediately stops and only produces an empty csv-file called copypatients.csv without showing any error messages.
- Lesion directory contains disconnectome maps in nii.gz-format (I have also tried running with lesions in nii-format).
- Result directory is specified as an empty location.
- Patient names and scores are included in a csv-file.
- I have tried different formats for the patient name including full path and with/without file extension.
- I run the software on macOS Catalina Version 10.15.3
- I have tried the suggested fix for Catalina users (http://toolkit.bcblab.com/): sudo spctl –master-disable
- I have successfully run Disconnectome Maps, FSL and R on the same computer earlier.
Best,
Henning
Hi Henning,
Thanks for your interest in the BCBtoolkit.
Regarding Anacom2, this has been designed for small samples (less than 50 participants) otherwise it is going to take an eternity for you to get your results.
For big samples like yours, I would recommend you to use fsl randomise instead as we did for Pacella et al. Elife (2019) https://elifesciences.org/articles/46075
Best of luck with your analyses. Kind regards. mich
Just wanted to post I got a script error and Michel kindly pointed out my files were MNI 2mm and not 1mm — easy fix 
thanks Michel!
1 Like
Hello BCBToolkit Stars,
I have an issue running AnaCOM2 on disconnectome data that I have created with the BCBToolkit.
The error provided by AnaCOM2 states that the behavioural data contains scores that are <= 0. I have checked the file and there are no scores in the file that are <=0.
To solve this problem I have tried;
- Using commas and colons to separate the score and filename columns and receive the error regardless.
- Specifying the filename with and without extensions (errors either way).
- To check the filenames are correct, I have also written out the filenames of the disconnectome maps in the folder [using dir(’*.nii.gz’)] to csv to verify that the image filenames match the filenames in the behavioural data.
- I’ve looked at the file in TextEdit and see no hidden characters or empty rows
- I have selected a smaller subset of the scans (just frontal patients) and still receive the error (n=88).
I’m running version 4.1.0 on a macbook pro (mid 2017) with catalina 10.15.6. I have run the suggested command line functions to give permission to bcbtoolbox to run.
I’ve attached a copy of the csv (as a txt) as an example in case I have missed something obvious.
PLSM_behaviouraldataf1.txt (694 Bytes)
Thank you in advance,
Bronson
Hi Bronson, Apologies for the delay, I am in the middle of the grant season for my students.
As far as I remember my colleague ran through the same issue here.
A few reasons that might help:
- The name of the file must be the same as the name of the lesion. I see that you used F01, F02, F03… please use the same name as you lesion file (maybe F01.nii.gz)
- No headers in the CSV file should be used.
- Another potential reason is that some MACOSX distributions of Excel create weird CSV with hidden characters (that even TextEdit can’t see)
Try with the attached file to see whether it works now. (check that I used the correct extension fo your lesions though! and change the .txt for .csvscores.txt (1.2 KB) )
Kind regards
mich
3 Likes
Hello Chris and Michel.
I just downloaded BCBtoolkit in a MacBook Pro Catalina. I am a newbie with MRI analysis and wanna use tractotron and disconnectome to analyze my MRI lesions. However, I am having trouble running tractotron itself. Even when I put the Lesion directory as the default one provided with the software, the program gets stuck about 1/5th of the way through.
Just no response at all (everything else is running ok at the same time). I have run the sudo spctl --master-disable as suggested in the BCBtoolkit website. Is there anything else I need to run to start with the software. Sorry for such a beginners question. Thanks in advance.
Hi Lushun,
Thanks for your interest in bcbtoolkit.
That is weird. Are you sure you did not use space characters in your directory?
When you run bcbtoolkit a log file is created in your folder. Would you have the kindness to share your log file with us? This will help us identify the problem.
Cheers
Mich
Hi Mich,
Thank you very much for getting back; and I apologise for the delayed response. Yes it was the space in the path and also some privacy issues with my Mac. Solved it.
However, I am having difficulty defining the lesion. I used MRICron to manually draw roi over my lesions. Do you have any suggestion on how to properly define the lesion as 1 and rest of the image as 0, as suggested in the bcbtoolkit webpage, with MRIcron? Do I also need to create a mask of the native MRI and overlap it with the ROI that I define or do I use an MNI template and use my ROI on that template itself?
Or can you suggest me the simplest way to adopt my image file for bcbtoolkit tractotron analysis using MRI image of the patients? I am experienced in MRI diagnosis itself so can manually draw the lesions correctly but a novice in computational analysis.
Lushun
Hi Lushun,
You can convert your ROI in nifti in MRIcron with 1 for the lesion an 0 for the non lesion.
Let me know if you don’t find the option I can send you some snapshots (currently in holiday with family :))
Kind regards,
Mich
Hi Chris and Michel,
I am trying to get started with the Tractotron part of the BCBToolkit and I am encountering an error:
++(70)/Users/selmalugtmeijer/Documents/scripts/BCBToolKitOSX/Tools/binaries/bin/fslstats /Users/selmalugtmeijer/Documents/scripts/BCBToolKitOSX/Tools/tmp/multresh/multresh_Anterior_Commissure 2.nii -R
Image Exception : #22 :: ERROR: Could not open image /Users/selmalugtmeijer/Documents/scripts/BCBToolKitOSX/Tools/tmp/multresh/multresh_Anterior_Commissure
+(70)max=
+(70)rm -rf /Users/selmalugtmeijer/Documents/scripts/BCBToolKitOSX/Tools/tmp/multresh
+(70)exit 1
+(71)printf ‘%s\t’
++(74)/Users/selmalugtmeijer/Documents/scripts/BCBToolKitOSX/Tools/binaries/bin/fslstats /Users/selmalugtmeijer/Documents/scripts/BCBToolKitOSX/Tools/tmp/multresh/tmpAnterior_Commissure 2.nii -V
++(74)awk ‘{print $1}’
Image Exception : #22 :: ERROR: Could not open image /Users/selmalugtmeijer/Documents/scripts/BCBToolKitOSX/Tools/tmp/multresh/tmpAnterior_Commissure
I get this same error when selecting my own lesion files and when selecting the lesion file that came with the toolbox.
Any thoughts on this?
Thanks in advance,
Selma
Hi Selma,
Thanks for your interest and sorry for your trouble. Would you have the kindness to share the log file that should have been created in the path you used for your analyses?
Thanks!
Hi,
can I send the log file to an e-mail address? I cannot upload files as new member.
Kind regards,
Selma
1 Like
Hi Chris and Michel.
Thank you for tour help. Premise: i am completely new to this field.
I recently installed your BCBtoolkit on my macbook pro running mac os big sur 11.1;
Here my problem: i’ve tried to use the Normalisation feature in order to use the diconnectome maps feature afterwards. The Normalisation process starts but it never ends, I waited more than an hour and at some point i decided to stop the process.
How long is this process usually? Is this normal? Have I done something wrong?
I can not attach the log file that was created during the process since I am a new user.
Thank you for your help!
Regards,
Aron
Hi Aron,
The process of normalisation is quite long actually in BCBtoolkit. Just let it run overnight 
Cheers
Mich
Hi Mich,
thank you for your quick answer. I guess it’s gonna be a lot of nights then 
Cheers.
Aron
You can run them all in a row in one night