Hi all,
I need to do linear mixed-effects (LME) modeling analysis of two within-subject factors in AFNI.
Factor1 has two levels: level1OfFactor1, level2OfFactor1.
And Factor2 has three levels : level1OfFactor2,level2OfFactor2, level3OfFactor2.
There are no demos of two within-subject factors in the help files.So I wrote my analysis script by rewriting the demo of one within-subject factor and one between-subject factor in the help file of 3dLME.Maybe there is some difference between
'one between-subject &one within-subject factors' and
'two within-subject factors' in the analysis code that I didn't know .An error came out when my script processed the 28th slice.
The feedback of error is:
#########################################################################
.
.
.
Z slice 25 done: 03/15/16 00:05:33.344
Z slice 26 done: 03/15/16 00:06:01.827
Z slice 27 done: 03/15/16 00:06:27.911
Error in checkForRemoteErrors(val) :
one node produced an error: 更换参数长度为零
Calls: aperm ... clusterApply -> staticClusterApply -> checkForRemoteErrors
停止执行
#########################################################################
My script is:
#########################################################################
3dLME -prefix $output/results.LME \
-jobs 10 \
-model "Factor1*Factor2" \
-ranEff '~1' \
-SS_type 3 \
-num_glt 4 \
-gltLabel 1 'level1OfFactor1-level2OfFactor1' -gltCode 1 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1' \
-gltLabel 2 'level1OfFactor2_Difference' -gltCode 2 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level1OfFactor2' \
-gltLabel 3 'level2OfFactor2_Difference' -gltCode 3 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level2OfFactor2' \
-gltLabel 4 'level3OfFactor2_Difference' -gltCode 4 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level3OfFactor2' \
-num_glf 6 \
-glfLabel 1 'Flevel1OfFactor1-level2OfFactor1' -glfCode 1 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level1OfFactor2 & 1*level2OfFactor2 & 1*level3OfFactor2' \
-glfLabel 2 'Flevel1OfFactor2Difference' -glfCode 2 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level1OfFactor2' \
-glfLabel 3 'Flevel2OfFactor2Difference' -glfCode 3 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level2OfFactor2' \
-glfLabel 4 'Flevel3OfFactor2Difference' -glfCode 4 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level3OfFactor2' \
##Interaction Effect
-glfLabel 5 'Interaction_Factor1onFactor2' -glfCode 5 'Factor1 : 1*level1OfFactor1 -1*level2OfFactor1 Factor2 : 1*level1OfFactor2-1*level2OfFactor2 & 1*level1OfFactor2-1*level3OfFactor2'\
-glfLabel 6 'Interaction_Factor2onFactor1' -glfCode 6 'Factor2 : 1*level1OfFactor2-1*level2OfFactor2 & 1*level1OfFactor2-1*level3OfFactor2 Factor1 : 1*level1OfFactor1 -1*level2OfFactor1'\
-dataTable @$output/datatable
#########################################################################
And the datatable is:
#########################################################################
Subj Factor1 Factor2 InputFile \
subject1 level1OfFactor1 level1OfFactor2 sub001.results/stats.sub001+tlrc[7] \
subject1 level1OfFactor1 level2OfFactor2 sub001.results/stats.sub001+tlrc[13] \
subject1 level1OfFactor1 level3OfFactor2 sub001.results/stats.sub001+tlrc[1] \
subject1 level2OfFactor1 level1OfFactor2 sub001.results/stats.sub001+tlrc[10] \
subject1 level2OfFactor1 level2OfFactor2 sub001.results/stats.sub001+tlrc[16] \
subject1 level2OfFactor1 level3OfFactor2 sub001.results/stats.sub001+tlrc[4] \
.
.
.
subject6 level1OfFactor1 level1OfFactor2 sub006.results/stats.sub006+tlrc[7] \
subject6 level1OfFactor1 level2OfFactor2 sub006.results/stats.sub006+tlrc[13] \
subject6 level1OfFactor1 level3OfFactor2 sub006.results/stats.sub006+tlrc[1] \
subject6 level2OfFactor1 level1OfFactor2 sub006.results/stats.sub006+tlrc[10] \
subject6 level2OfFactor1 level2OfFactor2 sub006.results/stats.sub006+tlrc[16] \
subject6 level2OfFactor1 level3OfFactor2 sub006.results/stats.sub006+tlrc[4]
#########################################################################
Can you tell me what's wrong with my script and how to improve it?
Thank you!