Hi everybody,

I have still a problem with the ICA by using MELODIC program on data preprocessed with AFNI software.

The program stopped after the step "STATS".

Here is the error message: in the same way, if I run the program with command-lines, it fails

Main script

/bin/cp /tmp/feat_gI6Lcd.fsf design.fsf

mkdir .files;cp /usr/local/fsl/doc/fsl.css .files;ln -s /usr/local/fsl/doc/images .files/images

/usr/local/fsl/bin/fsl_sub -T 10 -l logs -N feat1a_init   /usr/local/fsl/bin/feat /Users/CaroleG/Desktop/ICA/CELIA/+.ica/design.fsf -D /Users/CaroleG/Desktop/ICA/CELIA/+.ica -I 1 -init
29121

/usr/local/fsl/bin/fsl_sub -T 197 -l logs -N feat2_pre -j 29121  /usr/local/fsl/bin/feat /Users/CaroleG/Desktop/ICA/CELIA/+.ica/design.fsf -D /Users/CaroleG/Desktop/ICA/CELIA/+.ica -I 1 -prestats
divide by zero
    while executing
"if { [ catch {



for { set argindex 1 } { $argindex < $argc } { incr argindex 1 } {
    switch -- [ lindex $argv $argindex ] {

	-I {
	    incr argin..."
    (file "/usr/local/fsl/bin/feat" line 119)

---

Initialisation

/usr/local/fsl/bin/fslmaths /Users/CaroleG/Desktop/ICA/CELIA/RSN.celia.AllRuns prefiltered_func_data -odt float
Total original volumes = 1975

/usr/local/fsl/bin/fslroi prefiltered_func_data example_func 987 1

---

Prestats

/usr/local/fsl/bin/fslstats prefiltered_func_data -p 2 -p 98
0.000000 3.401314 

/usr/local/fsl/bin/fslmaths prefiltered_func_data -thr 0.3401314 -Tmin -bin mask -odt char

/usr/local/fsl/bin/fslstats prefiltered_func_data -k mask -p 50
0.000000 

/usr/local/fsl/bin/fslmaths mask -dilF mask

/usr/local/fsl/bin/fslmaths prefiltered_func_data -mas mask prefiltered_func_data_thresh

/usr/local/fsl/bin/fslmaths prefiltered_func_data_thresh -Tmean mean_func

/usr/local/fsl/bin/susan prefiltered_func_data_thresh 0.0 2.12314225053 3 1 1 mean_func 0.0 prefiltered_func_data_smooth

/usr/local/fsl/bin/fslmaths prefiltered_func_data_smooth -mas mask prefiltered_func_data_smooth
Cannot open volume prefiltered_func_data_smooth for reading!
divide by zero
    while executing
"if { [ catch {



for { set argindex 1 } { $argindex < $argc } { incr argindex 1 } {
    switch -- [ lindex $argv $argindex ] {

	-I {
	    incr argin..."
    (file "/usr/local/fsl/bin/feat" line 119)

---

Stats

/usr/local/fsl/bin/fsl_sub -T 180 -l logs -N feat4_ica   /usr/local/fsl/bin/melodic -i filtered_func_data -o filtered_func_data.ica -v --nobet --bgthreshold=3 --tr=2.0000000000 --report --guireport=../../report.html -d 0 --mmthresh=0.5 --Ostats

Tank you in advance for your help!

Le 23 déc. 2012 à 09:22, Carole Guedj a écrit :

Hi Stephen and Matthew,

Thank you very much for your help. I will try that!!

Carole