Hello,

I am attempting to run a between subjects 3dMEMA with the -unequal_variance option. Other 3dMEMA script run with this data give output when the -unequal_variance option is omitted. Does this option conflict with other options i've passed into the script?

I am receiving the following message:

[1] "Analysis finished: 06/13/12 10:17:25.132"
[1] "#++++++++++++++++++++++++++++++++++++++++++++"
Error in ifelse(DF >= 1, qt(pt(tval, DF), fullDF), 0) :
object 'i' not found
Calls: tConvert -> ifelse
Execution halted

My script is as follows:

cd ${rootDir}/
gp1='ASD'
gp2='TD'
list1='S01A S01A/3dREMLfit/S01A_glm_hrf_stats_SVC_REML+tlrc'[2]' S01A/3dREMLfit/S01A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S02A S02A/3dREMLfit/S02A_glm_hrf_stats_SVC_REML+tlrc'[2]' S02A/3dREMLfit/S02A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S03A S03A/3dREMLfit/S03A_glm_hrf_stats_SVC_REML+tlrc'[2]' S03A/3dREMLfit/S03A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S04A S04A/3dREMLfit/S04A_glm_hrf_stats_SVC_REML+tlrc'[2]' S04A/3dREMLfit/S04A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S06A S06A/3dREMLfit/S06A_glm_hrf_stats_SVC_REML+tlrc'[2]' S06A/3dREMLfit/S06A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S07A S07A/3dREMLfit/S07A_glm_hrf_stats_SVC_REML+tlrc'[2]' S07A/3dREMLfit/S07A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S08A S08A/3dREMLfit/S08A_glm_hrf_stats_SVC_REML+tlrc'[2]' S08A/3dREMLfit/S08A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S09A S09A/3dREMLfit/S09A_glm_hrf_stats_SVC_REML+tlrc'[2]' S09A/3dREMLfit/S09A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S10A S10A/3dREMLfit/S10A_glm_hrf_stats_SVC_REML+tlrc'[2]' S10A/3dREMLfit/S10A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S11A S11A/3dREMLfit/S11A_glm_hrf_stats_SVC_REML+tlrc'[2]' S11A/3dREMLfit/S11A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S12A S12A/3dREMLfit/S12A_glm_hrf_stats_SVC_REML+tlrc'[2]' S12A/3dREMLfit/S12A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S13A S13A/3dREMLfit/S13A_glm_hrf_stats_SVC_REML+tlrc'[2]' S13A/3dREMLfit/S13A_glm_hrf_stats_SVC_REML+tlrc'[3]'
S14A S14A/3dREMLfit/S14A_glm_hrf_stats_SVC_REML+tlrc'[2]' S14A/3dREMLfit/S14A_glm_hrf_stats_SVC_REML+tlrc'[3]''
list2='S01 S01/3dREMLfit/S01_glm_hrf_stats_SVC_REML+tlrc'[2]' S01/3dREMLfit/S01_glm_hrf_stats_SVC_REML+tlrc'[3]'
S02 S02/3dREMLfit/S02_glm_hrf_stats_SVC_REML+tlrc'[2]' S02/3dREMLfit/S02_glm_hrf_stats_SVC_REML+tlrc'[3]'
S03 S03/3dREMLfit/S03_glm_hrf_stats_SVC_REML+tlrc'[2]' S03/3dREMLfit/S03_glm_hrf_stats_SVC_REML+tlrc'[3]'
S04 S04/3dREMLfit/S04_glm_hrf_stats_SVC_REML+tlrc'[2]' S04/3dREMLfit/S04_glm_hrf_stats_SVC_REML+tlrc'[3]'
S06 S06/3dREMLfit/S06_glm_hrf_stats_SVC_REML+tlrc'[2]' S06/3dREMLfit/S06_glm_hrf_stats_SVC_REML+tlrc'[3]'
S07 S07/3dREMLfit/S07_glm_hrf_stats_SVC_REML+tlrc'[2]' S07/3dREMLfit/S07_glm_hrf_stats_SVC_REML+tlrc'[3]'
S08 S08/3dREMLfit/S08_glm_hrf_stats_SVC_REML+tlrc'[2]' S08/3dREMLfit/S08_glm_hrf_stats_SVC_REML+tlrc'[3]'
S09 S09/3dREMLfit/S09_glm_hrf_stats_SVC_REML+tlrc'[2]' S09/3dREMLfit/S09_glm_hrf_stats_SVC_REML+tlrc'[3]'
S10 S10/3dREMLfit/S10_glm_hrf_stats_SVC_REML+tlrc'[2]' S10/3dREMLfit/S10_glm_hrf_stats_SVC_REML+tlrc'[3]'
S11 S11/3dREMLfit/S11_glm_hrf_stats_SVC_REML+tlrc'[2]' S11/3dREMLfit/S11_glm_hrf_stats_SVC_REML+tlrc'[3]'
S12 S12/3dREMLfit/S12_glm_hrf_stats_SVC_REML+tlrc'[2]' S12/3dREMLfit/S12_glm_hrf_stats_SVC_REML+tlrc'[3]'
S13 S13/3dREMLfit/S13_glm_hrf_stats_SVC_REML+tlrc'[2]' S13/3dREMLfit/S13_glm_hrf_stats_SVC_REML+tlrc'[3]'
S14 S14/3dREMLfit/S14_glm_hrf_stats_SVC_REML+tlrc'[2]' S14/3dREMLfit/S14_glm_hrf_stats_SVC_REML+tlrc'[3]''

3dMEMA \
-missing_data 0 \
-residual_Z \
-unequal_variance \
-prefix ./3dMEMA/data/${gp1}v${gp2}_123 \
-jobs 8 \
-groups ${gp1} ${gp2} \
-set ${gp1} \
${list1} \
-set ${gp2} \
${list2}