Dear Sara, When you say you run melodic, does that mean you run melodic with the same mask, same resolution and smoothing you use in the toolbox? Or do you mean an unmasked ICA? In the latter case it would not be surprising to get completely different results because of different intrinsic dimensionalities of the different brain regions. Unfortunately, we had to remove that part of the results from our toolbox paper prior to publication. To understand this problem, I need more information. Are you running split-half or test-retest reproducibility analysis? If the first, how many repetitions did you choose? if it was only one, this may explain the crash since there is nothing to calculate the correlation coefficients from. Try running more repetitions then. The ICA should not have convergence problems for all samplings. The no_convergence_error.png file should then be in the top folder of the analysis (together with an output "corr_mean.png" that is essentially the figure 2A of our paper). If it is not there, that means that the analysis did not finish properly.

If everything else went smoothly, you can use
python '/Path_to_your_toolbox_installation/py/ic_corr /Path_to_folder_containing_the_samples_folder/ number_of_samplings dimensionality_range 1'
and re-run everything except for the calculation of the single ICAs. Try to leave out the problematic number in in the dimensionality range argument. For example (if you had 10 repetitions):
python '/Path_to_your_toolbox_installation/py/ic_corr /Path_to_folder_containing_the_samples_folder/ 10 2-10,12-20 1'
Let me know if that works.

Best,

Florian