Hello,

I have 2 groups, kids and adults in my surface-based analysis stream. For each group I have 6 stats files per subject. These stats files correspond to 6 seed-to-whole-brain functional connectivity maps. I would like to compare whole brain connectivity of each seed regions across groups. Thus, I would like the results of this step to be 6 stats files corresponding to kid-adult for each seed region. I have set up my 3dMEMA command as follows:

#stats sub brik 2 = coef and 3 = t stat
set i = 1
while ($i <= 6)
cd $dir/group/mema/all/sub${i}
foreach h (lh rh)
3dMEMA \
-prefix KidminusAdult_roi${i}.$h.niml.dset \
-cio -max_zeros 2 -jobs 6 -unequal_variance -verb 2 \
-groups adult kid \
-set adult \
sub002 sub002.roi${i}.$h.stats.niml.dset'[2]' sub002.roi${i}.$h.stats.niml.dset'[3]' \
sub004 sub004.roi${i}.$h.stats.niml.dset'[2]' sub004.roi${i}.$h.stats.niml.dset'[3]' \
sub005 sub005.roi${i}.$h.stats.niml.dset'[2]' sub005.roi${i}.$h.stats.niml.dset'[3]' \
sub006 sub006.roi${i}.$h.stats.niml.dset'[2]' sub006.roi${i}.$h.stats.niml.dset'[3]' \
sub007 sub007.roi${i}.$h.stats.niml.dset'[2]' sub007.roi${i}.$h.stats.niml.dset'[3]' \
sub008 sub008.roi${i}.$h.stats.niml.dset'[2]' sub008.roi${i}.$h.stats.niml.dset'[3]' \
sub009 sub009.roi${i}.$h.stats.niml.dset'[2]' sub009.roi${i}.$h.stats.niml.dset'[3]' \
sub010 sub010.roi${i}.$h.stats.niml.dset'[2]' sub010.roi${i}.$h.stats.niml.dset'[3]' \
sub011 sub011.roi${i}.$h.stats.niml.dset'[2]' sub011.roi${i}.$h.stats.niml.dset'[3]' \
sub012 sub012.roi${i}.$h.stats.niml.dset'[2]' sub012.roi${i}.$h.stats.niml.dset'[3]' \
sub013 sub013.roi${i}.$h.stats.niml.dset'[2]' sub013.roi${i}.$h.stats.niml.dset'[3]' \
sub014 sub014.roi${i}.$h.stats.niml.dset'[2]' sub014.roi${i}.$h.stats.niml.dset'[3]' \
sub015 sub015.roi${i}.$h.stats.niml.dset'[2]' sub015.roi${i}.$h.stats.niml.dset'[3]' \
sub016 sub016.roi${i}.$h.stats.niml.dset'[2]' sub016.roi${i}.$h.stats.niml.dset'[3]' \
sub017 sub017.roi${i}.$h.stats.niml.dset'[2]' sub017.roi${i}.$h.stats.niml.dset'[3]' \
sub018 sub018.roi${i}.$h.stats.niml.dset'[2]' sub018.roi${i}.$h.stats.niml.dset'[3]' \
sub019 sub019.roi${i}.$h.stats.niml.dset'[2]' sub019.roi${i}.$h.stats.niml.dset'[3]' \
sub020 sub020.roi${i}.$h.stats.niml.dset'[2]' sub020.roi${i}.$h.stats.niml.dset'[3]' \
sub021 sub021.roi${i}.$h.stats.niml.dset'[2]' sub021.roi${i}.$h.stats.niml.dset'[3]' \
sub022 sub022.roi${i}.$h.stats.niml.dset'[2]' sub022.roi${i}.$h.stats.niml.dset'[3]' \
sub023 sub023.roi${i}.$h.stats.niml.dset'[2]' sub023.roi${i}.$h.stats.niml.dset'[3]' \
-set kid \
sub104 sub104.roi${i}.$h.stats.niml.dset'[2]' sub104.roi${i}.$h.stats.niml.dset'[3]' \
sub105 sub105.roi${i}.$h.stats.niml.dset'[2]' sub105.roi${i}.$h.stats.niml.dset'[3]' \
sub109 sub109.roi${i}.$h.stats.niml.dset'[2]' sub109.roi${i}.$h.stats.niml.dset'[3]' \
sub110 sub110.roi${i}.$h.stats.niml.dset'[2]' sub110.roi${i}.$h.stats.niml.dset'[3]' \
sub115 sub115.roi${i}.$h.stats.niml.dset'[2]' sub115.roi${i}.$h.stats.niml.dset'[3]' \
sub116 sub116.roi${i}.$h.stats.niml.dset'[2]' sub116.roi${i}.$h.stats.niml.dset'[3]' \
sub120 sub120.roi${i}.$h.stats.niml.dset'[2]' sub120.roi${i}.$h.stats.niml.dset'[3]' \
sub121 sub121.roi${i}.$h.stats.niml.dset'[2]' sub121.roi${i}.$h.stats.niml.dset'[3]' \
sub125 sub125.roi${i}.$h.stats.niml.dset'[2]' sub125.roi${i}.$h.stats.niml.dset'[3]' \
sub126 sub126.roi${i}.$h.stats.niml.dset'[2]' sub126.roi${i}.$h.stats.niml.dset'[3]' \
sub130 sub130.roi${i}.$h.stats.niml.dset'[2]' sub130.roi${i}.$h.stats.niml.dset'[3]' \
sub133 sub133.roi${i}.$h.stats.niml.dset'[2]' sub133.roi${i}.$h.stats.niml.dset'[3]' \
sub134 sub134.roi${i}.$h.stats.niml.dset'[2]' sub134.roi${i}.$h.stats.niml.dset'[3]' \
sub140 sub140.roi${i}.$h.stats.niml.dset'[2]' sub140.roi${i}.$h.stats.niml.dset'[3]' \
sub146 sub146.roi${i}.$h.stats.niml.dset'[2]' sub146.roi${i}.$h.stats.niml.dset'[3]' \
sub154 sub154.roi${i}.$h.stats.niml.dset'[2]' sub154.roi${i}.$h.stats.niml.dset'[3]' \
sub157 sub157.roi${i}.$h.stats.niml.dset'[2]' sub157.roi${i}.$h.stats.niml.dset'[3]' \
sub164 sub164.roi${i}.$h.stats.niml.dset'[2]' sub164.roi${i}.$h.stats.niml.dset'[3]' \
sub170 sub170.roi${i}.$h.stats.niml.dset'[2]' sub170.roi${i}.$h.stats.niml.dset'[3]' \
sub173 sub173.roi${i}.$h.stats.niml.dset'[2]' sub173.roi${i}.$h.stats.niml.dset'[3]' \
sub174 sub174.roi${i}.$h.stats.niml.dset'[2]' sub174.roi${i}.$h.stats.niml.dset'[3]' \
sub177 sub177.roi${i}.$h.stats.niml.dset'[2]' sub177.roi${i}.$h.stats.niml.dset'[3]' \
sub178 sub178.roi${i}.$h.stats.niml.dset'[2]' sub178.roi${i}.$h.stats.niml.dset'[3]' \
sub181 sub181.roi${i}.$h.stats.niml.dset'[2]' sub181.roi${i}.$h.stats.niml.dset'[3]' \
sub184 sub184.roi${i}.$h.stats.niml.dset'[2]' sub184.roi${i}.$h.stats.niml.dset'[3]' \
sub108 sub108.roi${i}.$h.stats.niml.dset'[2]' sub108.roi${i}.$h.stats.niml.dset'[3]' \
sub111 sub111.roi${i}.$h.stats.niml.dset'[2]' sub111.roi${i}.$h.stats.niml.dset'[3]' \
sub123 sub123.roi${i}.$h.stats.niml.dset'[2]' sub123.roi${i}.$h.stats.niml.dset'[3]' \
sub129 sub129.roi${i}.$h.stats.niml.dset'[2]' sub129.roi${i}.$h.stats.niml.dset'[3]' \
sub131 sub131.roi${i}.$h.stats.niml.dset'[2]' sub131.roi${i}.$h.stats.niml.dset'[3]' \
sub132 sub132.roi${i}.$h.stats.niml.dset'[2]' sub132.roi${i}.$h.stats.niml.dset'[3]' \
sub141 sub141.roi${i}.$h.stats.niml.dset'[2]' sub141.roi${i}.$h.stats.niml.dset'[3]' \
sub143 sub143.roi${i}.$h.stats.niml.dset'[2]' sub143.roi${i}.$h.stats.niml.dset'[3]' \
sub151 sub151.roi${i}.$h.stats.niml.dset'[2]' sub151.roi${i}.$h.stats.niml.dset'[3]' \
sub155 sub155.roi${i}.$h.stats.niml.dset'[2]' sub155.roi${i}.$h.stats.niml.dset'[3]' \
sub156 sub156.roi${i}.$h.stats.niml.dset'[2]' sub156.roi${i}.$h.stats.niml.dset'[3]' \
sub158 sub158.roi${i}.$h.stats.niml.dset'[2]' sub158.roi${i}.$h.stats.niml.dset'[3]' \
sub159 sub159.roi${i}.$h.stats.niml.dset'[2]' sub159.roi${i}.$h.stats.niml.dset'[3]' \
sub160 sub160.roi${i}.$h.stats.niml.dset'[2]' sub160.roi${i}.$h.stats.niml.dset'[3]' \
sub165 sub165.roi${i}.$h.stats.niml.dset'[2]' sub165.roi${i}.$h.stats.niml.dset'[3]' \
sub167 sub167.roi${i}.$h.stats.niml.dset'[2]' sub167.roi${i}.$h.stats.niml.dset'[3]' \
sub169 sub169.roi${i}.$h.stats.niml.dset'[2]' sub169.roi${i}.$h.stats.niml.dset'[3]' \
sub172 sub172.roi${i}.$h.stats.niml.dset'[2]' sub172.roi${i}.$h.stats.niml.dset'[3]' \
sub183 sub183.roi${i}.$h.stats.niml.dset'[2]' sub183.roi${i}.$h.stats.niml.dset'[3]' \
sub185 sub185.roi${i}.$h.stats.niml.dset'[2]' sub185.roi${i}.$h.stats.niml.dset'[3]' \
sub187 sub187.roi${i}.$h.stats.niml.dset'[2]' sub187.roi${i}.$h.stats.niml.dset'[3]'
end
@ i++
end


For some reason only the .RxMat is written, no niml.dset stats file and I see this error:

Error in statsym[[1]]$sb : $ operator is invalid for atomic vectors
Calls: write.AFNI ... set.c.AFNI.attribute -> r.NI_new_element -> statsym.list2code
Execution halted

Please let me know if I set something up wrong. Thanks!



Edited 1 time(s). Last edit at 08/13/2015 03:38PM by dmoracze.