Understanding 3dLME output for ICC

Hello!

I’ve just recently discovered what ICC is and 3dLME, so forgive my ignorance. I am trying to compute test-retest reliability for 3 runs of resting state. The file I input is a whole brain correlation map from a seed in the left inferior frontal gyrus. I am struggling to understand the output of the command. It seems that the only output is a 3D dataset and I’m not sure how to interpret the results in the GUI. Also, is there a way to add an option to quantify values in a text file or something? Thanks so much!! Here’s my code:

3dLME -prefix test_ICC
-jobs 12
-model “1”
-ranEff ‘Cond+Subj’
-ICC
-dataTable
Subj Cond InputFile
p1 EO1 p1_left_IFG_rest_run1_Tcorr+tlrc
p1 EC2 p1_left_IFG_rest_run2_Tcorr+tlrc
p1 EC3 p1_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p96 EO1 p96_left_IFG_rest_run1_Tcorr+tlrc
p96 EC2 p96_left_IFG_rest_run2_Tcorr+tlrc
p96 EC3 p96_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc

I assume that you have more than two subjects as input.

It seems that the only output is a 3D dataset and I’m not sure how to interpret the results in the GUI.

There should have more than one sub-brick in the output. The one with the label ‘Subj’ is the ICC you’re looking for.

Also, is there a way to add an option to quantify values in a text file or something?

This is a whole-brain voxel-wise ICC analysis. I’m not so sure what you’re asking. If you want to the ICC for a particular list of ROIs, you can extract them from the whole brain analysis results.

I assume that you have more than two subjects as input.

Yes, I just wanted to run a dummy before inputting the larger dataset.

There should have more than one sub-brick in the output. The one with the label ‘Subj’ is the ICC you’re looking for.

So the ICC map under ‘Subj’ sub-brick can be interpreted as what? What is the typical threshold for the map?

This is a whole-brain voxel-wise ICC analysis. I’m not so sure what you’re asking. If you want to the ICC for a particular list of ROIs, you can extract them from the whole brain analysis results.

how would you do that?

So the ICC map under ‘Subj’ sub-brick can be interpreted as what?

I’m not sure what you’re asking. That sub-brick is the voxel-wise ICC value.

What is the typical threshold for the map?

A loose rule of thumb has been adopted as the following: [0, 0.4), poor; [0.4, 0.6), fair; [0.6, 0.75), good; and [0.75, 1], strong. See more discussion here: https://onlinelibrary.wiley.com/doi/pdf/10.1002/hbm.23909

Thank you! I have a follow up question. Like I said earlier, I am trying to look at ICC across 3 runs of resting state fMRI data for 31 patients. When I run the 3dLME command, it gets killed early on. Is there a reason for this? Below the command and options I used, as well as the output:

3dLME -prefix test_31_patients_ICC_left_IFG
-jobs 120
-model “1”
-ranEff ‘Cond+Subj’
-ICC
-dataTable
Subj Cond InputFile
p3 EO1 p3_left_IFG_rest_run1_Tcorr+tlrc
p3 EC2 p3_left_IFG_rest_run2_Tcorr+tlrc
p3 EC3 p3_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p25 EO1 p25_left_IFG_rest_run1_Tcorr+tlrc
p25 EC2 p25_left_IFG_rest_run2_Tcorr+tlrc
p25 EC3 p25_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p32 EO1 p32_left_IFG_rest_run1_Tcorr+tlrc
p32 EC2 p32_left_IFG_rest_run2_Tcorr+tlrc
p32 EC3 p32_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p33 EO1 p33_left_IFG_rest_run1_Tcorr+tlrc
p33 EC2 p33_left_IFG_rest_run2_Tcorr+tlrc
p33 EC3 p33_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p43 EO1 p43_left_IFG_rest_run1_Tcorr+tlrc
p43 EC2 p43_left_IFG_rest_run2_Tcorr+tlrc
p43 EC3 p43_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p47 EO1 p47_left_IFG_rest_run1_Tcorr+tlrc
p47 EC2 p47_left_IFG_rest_run2_Tcorr+tlrc
p47 EC3 p47_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p48 EO1 p48_left_IFG_rest_run1_Tcorr+tlrc
p48 EC2 p48_left_IFG_rest_run2_Tcorr+tlrc
p48 EC3 p48_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p64 EO1 p64_left_IFG_rest_run1_Tcorr+tlrc
p64 EC2 p64_left_IFG_rest_run2_Tcorr+tlrc
p64 EC3 p64_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p65 EO1 p65_left_IFG_rest_run1_Tcorr+tlrc
p65 EC2 p65_left_IFG_rest_run2_Tcorr+tlrc
p65 EC3 p65_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p67 EO1 p67_left_IFG_rest_run1_Tcorr+tlrc
p67 EC2 p67_left_IFG_rest_run2_Tcorr+tlrc
p67 EC3 p67_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p71 EO1 p71_left_IFG_rest_run1_Tcorr+tlrc
p71 EC2 p71_left_IFG_rest_run2_Tcorr+tlrc
p71 EC3 p71_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p73 EO1 p73_left_IFG_rest_run1_Tcorr+tlrc
p73 EC2 p73_left_IFG_rest_run2_Tcorr+tlrc
p73 EC3 p73_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p76 EO1 p76_left_IFG_rest_run1_Tcorr+tlrc
p76 EC2 p76_left_IFG_rest_run2_Tcorr+tlrc
p76 EC3 p76_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p82 EO1 p82_left_IFG_rest_run1_Tcorr+tlrc
p82 EC2 p82_left_IFG_rest_run2_Tcorr+tlrc
p82 EC3 p82_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p83 EO1 p83_left_IFG_rest_run1_Tcorr+tlrc
p83 EC2 p83_left_IFG_rest_run2_Tcorr+tlrc
p83 EC3 p83_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p86 EO1 p86_left_IFG_rest_run1_Tcorr+tlrc
p86 EC2 p86_left_IFG_rest_run2_Tcorr+tlrc
p86 EC3 p86_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p92 EO1 p92_left_IFG_rest_run1_Tcorr+tlrc
p92 EC2 p92_left_IFG_rest_run2_Tcorr+tlrc
p92 EC3 p92_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p96 EO1 p96_left_IFG_rest_run1_Tcorr+tlrc
p96 EC2 p96_left_IFG_rest_run2_Tcorr+tlrc
p96 EC3 p96_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p5 EO1 p5_left_IFG_rest_run1_Tcorr+tlrc
p5 EC2 p5_left_IFG_rest_run2_Tcorr+tlrc
p5 EC3 p5_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p14 EO1 p14_left_IFG_rest_run1_Tcorr+tlrc
p14 EC2 p14_left_IFG_rest_run2_Tcorr+tlrc
p14 EC3 p14_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p17 EO1 p17_left_IFG_rest_run1_Tcorr+tlrc
p17 EC2 p17_left_IFG_rest_run2_Tcorr+tlrc
p17 EC3 p17_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p40 EO1 p40_left_IFG_rest_run1_Tcorr+tlrc
p40 EC2 p40_left_IFG_rest_run2_Tcorr+tlrc
p40 EC3 p40_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p78 EO1 p78_left_IFG_rest_run1_Tcorr+tlrc
p78 EC2 p78_left_IFG_rest_run2_Tcorr+tlrc
p78 EC3 p78_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p79 EO1 p79_left_IFG_rest_run1_Tcorr+tlrc
p79 EC2 p79_left_IFG_rest_run2_Tcorr+tlrc
p79 EC3 p79_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p80 EO1 p80_left_IFG_rest_run1_Tcorr+tlrc
p80 EC2 p80_left_IFG_rest_run2_Tcorr+tlrc
p80 EC3 p80_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p85 EO1 p85_left_IFG_rest_run1_Tcorr+tlrc
p85 EC2 p85_left_IFG_rest_run2_Tcorr+tlrc
p85 EC3 p85_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p97 EO1 p97_left_IFG_rest_run1_Tcorr+tlrc
p97 EC2 p97_left_IFG_rest_run2_Tcorr+tlrc
p97 EC3 p97_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p23 EO1 p23_left_IFG_rest_run1_Tcorr+tlrc
p23 EC2 p23_left_IFG_rest_run2_Tcorr+tlrc
p23 EC3 p23_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p29 EO1 p29_left_IFG_rest_run1_Tcorr+tlrc
p29 EC2 p29_left_IFG_rest_run2_Tcorr+tlrc
p29 EC3 p29_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p12 EO1 p12_left_IFG_rest_run1_Tcorr+tlrc
p12 EC2 p12_left_IFG_rest_run2_Tcorr+tlrc
p12 EC3 p12_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc
p22 EO1 p22_left_IFG_rest_run1_Tcorr+tlrc
p22 EC2 p22_left_IFG_rest_run2_Tcorr+tlrc
p22 EC3 p22_left_IFG_rest_run3_eyes_closed_Tcorr+tlrc

Output:

*+ WARNING: didn’t find any environment equations in ~/.afnirc
Loading required package: lme4
Loading required package: Matrix
Package lme4 loaded successfully!

++++++++++++++++++++++++++++++++++++++++++++++++++++
***** Summary information of data structure *****
31 subjects : p12 p14 p17 p22 p23 p25 p29 p3 p32 p33 p40 p43 p47 p48 p5 p64 p65 p67 p71 p73 p76 p78 p79 p80 p82 p83 p85 p86 p92 p96 p97
93 response values
3 levels for factor Cond : EC2 EC3 EO1
0 post hoc tests

Contingency tables of subject distributions among the categorical variables:

***** End of data structure information *****
++++++++++++++++++++++++++++++++++++++++++++++++++++

Reading input files now…

Reading input files: Done!

If the program hangs here for more than, for example, half an hour,
kill the process because the model specification or the GLT coding
is likely inappropriate.

[1] “Great, test run passed at voxel (32, 57, 48)!”
[1] “Start to compute 96 slices along Z axis. You can monitor the progress”
[1] “and estimate the total run time as shown below.”
[1] “09/18/18 15:39:36.897”
Loading required package: snow
Package snow loaded successfully!

Killed

Any thoughts?

Several possibilities might have killed the jobs. Just rerun the script.

thanks for all your help!