inconsistent 3dMEMA errors

Hi Gang - we have been getting some weird and inconsistent errors with 3dMEMA recently. These errors happen for new analyses and old scripts that previously worked (and nothing was changed). Sometimes the program seems to work fine, however, many more times we are halted by random errors. I have tried to find a pattern and cannot. I can run the same code on the same data multiple times in a row and sometimes it will give a random error (possibly different errors even though the same code was run) and sometimes it will work.

Here are some examples of the errors we have received:


Error in gzfile(file, "rb") : invalid 'compress' argument
Execution halted


Error in dim(rs2$brk) <- ddb : dims do not match the length of object
Calls: process.MEMA.opts -> lapply -> FUN -> read.c.AFNI
Execution halted


Error in routines[[103]] : subscript out of bounds
Execution halted


Warning in read.MEMA.opts.batch(args, verb = 0) :
  You must have either one or two groups
Error: Error parsing input
Execution halted


During startup - Warning message:
package "methods" in options("defaultPackages") was not found

** Error:
   Expecting 101 parameters for option "-prefix".
   Have 1 parameter(s) in string "VSppi.agree.peer-comp" instead.
NULL options

** Error:
   Error parsing arguments. See 3dMEMA -help for details.


During startup - Warning messages:
1: package "methods" in options("defaultPackages") was not found
2: package 'datasets' in options("defaultPackages") was not found
3: package 'utils' in options("defaultPackages") was not found
4: package 'grDevices' in options("defaultPackages") was not found
5: package 'graphics' in options("defaultPackages") was not found
6: package 'stats' in options("defaultPackages") was not found
7: package 'methods' in options("defaultPackages") was not found
Error in strsplit(Sys.getenv("PATH"), ":")[[101]] :
  subscript out of bounds
Calls: source -> first.in.path -> paste
Execution halted


Error in lop$subjLab[[ii]] : subscript out of bounds
Calls: process.MEMA.opts
Execution halted

Here are some other errors we have received from another bit of code that used to work:


Writing xMat to copyASDvTD_11_Age+meanFD_P-C_G.RxMat
Error in array(apply(comArr[, , kk, ], 3, function(x) x * (abs(lop$maskData[,  :
  negative length vectors are not allowed
In addition: Warning message:
In array(apply(comArr[, , kk, ], 3, function(x) x * (abs(lop$maskData[,  :
  NAs introduced by coercion to integer range
Execution halted


Error in dim(rs2$brk) <- ddb : the dims contain missing values
Calls: process.MEMA.opts -> lapply -> FUN -> read.c.AFNI
In addition: Warning messages:
1: In get.c.AFNI.attribute(hatr, "DATASET_DIMENSIONS") :
  NAs introduced by coercion to integer range
2: In get.c.AFNI.attribute(hatr, "DATASET_RANK") :
  NAs introduced by coercion to integer range
Execution halted

We are using Precompiled binary linux_xorg7_64: Feb 21 2018 (Version AFNI_18.0.18) and R 3.2.2.

Running afni_system_check is all fine except that I don’t have the R library ‘psych’ installed. However, I haven’t needed it for MEMA or LME previously.

We are also running a lot of analysis on our server: multiple 3dDeconvolve and 3dMEMA commands, each of which with multiple -jobs. Do you think this could get in the way?

I never saw such drunken behavior of 3dMEMA before! And I cannot think of any reason for that either… If you run out of options, upload the data and script, I’ll see if I have a better luck.

Hello,

I’m am seeing some of the same errors and was wondering if anything else been figured out with regard to this?

Errors -
** Error:
Expecting more than 4 parameters for option “-set”.
Have 1 parameter(s) in string “pos” instead.
**Error:
Error parsing arguments. See 3dMEMA -help for details
021: Command not found.
-model_outliers: Command not found.

I’m attempting to run the following script:

#run from group_analysis_scripts

3dMEMA -prefix …/…/group_analysis/craving_between_groups
-groups neutral positive
-set neut
001 …/…/preproc_data/001_glm/stats.001+tlrc’[43]’ …/…/preproc_data/001_glm/stats.001+tlrc’[44]’
002 …/…/preproc_data/002_glm/stats.002+tlrc’[43]’ …/…/preproc_data/002_glm/stats.002+tlrc’[44]’
003 …/…/preproc_data/003_glm/stats.003+tlrc’[43]’ …/…/preproc_data/003_glm/stats.003+tlrc’[44]’
004 …/…/preproc_data/004_glm/stats.004+tlrc’[43]’ …/…/preproc_data/004_glm/stats.004+tlrc’[44]’
005 …/…/preproc_data/005_glm/stats.005+tlrc’[43]’ …/…/preproc_data/005_glm/stats.005+tlrc’[44]’
006 …/…/preproc_data/006_glm/stats.006+tlrc’[43]’ …/…/preproc_data/006_glm/stats.006+tlrc’[44]’
007 …/…/preproc_data/007_glm/stats.007+tlrc’[43]’ …/…/preproc_data/007_glm/stats.007+tlrc’[44]’
008 …/…/preproc_data/008_glm/stats.008+tlrc’[43]’ …/…/preproc_data/008_glm/stats.008+tlrc’[44]’
009 …/…/preproc_data/009_glm/stats.009+tlrc’[43]’ …/…/preproc_data/009_glm/stats.009+tlrc’[44]’
010 …/…/preproc_data/010_glm/stats.010+tlrc’[43]’ …/…/preproc_data/010_glm/stats.010+tlrc’[44]’
101 …/…/preproc_data/101_glm/stats.101+tlrc’[43]’ …/…/preproc_data/101_glm/stats.101+tlrc’[44]’
102 …/…/preproc_data/102_glm/stats.102+tlrc’[43]’ …/…/preproc_data/102_glm/stats.102+tlrc’[44]’
011 …/…/preproc_data/011_glm/stats.011+tlrc’[43]’ …/…/preproc_data/011_glm/stats.011+tlrc’[44]’
012 …/…/preproc_data/012_glm/stats.012+tlrc’[43]’ …/…/preproc_data/012_glm/stats.012+tlrc’[44]’
013 …/…/preproc_data/013_glm/stats.013+tlrc’[43]’ …/…/preproc_data/013_glm/stats.013+tlrc’[44]’
014 …/…/preproc_data/014_glm/stats.014+tlrc’[43]’ …/…/preproc_data/014_glm/stats.014+tlrc’[44]’
015 …/…/preproc_data/015_glm/stats.015+tlrc’[43]’ …/…/preproc_data/015_glm/stats.015+tlrc’[44]’
016 …/…/preproc_data/016_glm/stats.016+tlrc’[43]’ …/…/preproc_data/016_glm/stats.016+tlrc’[44]’
017 …/…/preproc_data/017_glm/stats.017+tlrc’[43]’ …/…/preproc_data/017_glm/stats.017+tlrc’[44]’
018 …/…/preproc_data/018_glm/stats.018+tlrc’[43]’ …/…/preproc_data/018_glm/stats.018+tlrc’[44]’
-set pos \
021 …/…/preproc_data/021_glm/stats.021+tlrc’[43]’ …/…/preproc_data/021_glm/stats.021+tlrc’[44]’
022 …/…/preproc_data/022_glm/stats.022+tlrc’[43]’ …/…/preproc_data/022_glm/stats.022+tlrc’[44]’
023 …/…/preproc_data/023_glm/stats.023+tlrc’[43]’ …/…/preproc_data/023_glm/stats.023+tlrc’[44]’
024 …/…/preproc_data/024_glm/stats.024+tlrc’[43]’ …/…/preproc_data/024_glm/stats.024+tlrc’[44]’
025 …/…/preproc_data/025_glm/stats.025+tlrc’[43]’ …/…/preproc_data/025_glm/stats.025+tlrc’[44]’
026 …/…/preproc_data/026_glm/stats.026+tlrc’[43]’ …/…/preproc_data/026_glm/stats.026+tlrc’[44]’
027 …/…/preproc_data/027_glm/stats.027+tlrc’[43]’ …/…/preproc_data/027_glm/stats.027+tlrc’[44]’
028 …/…/preproc_data/028_glm/stats.028+tlrc’[43]’ …/…/preproc_data/028_glm/stats.028+tlrc’[44]’
029 …/…/preproc_data/029_glm/stats.029+tlrc’[43]’ …/…/preproc_data/029_glm/stats.029+tlrc’[44]’
030 …/…/preproc_data/030_glm/stats.030+tlrc’[43]’ …/…/preproc_data/030_glm/stats.030+tlrc’[44]’
121 …/…/preproc_data/121_glm/stats.121+tlrc’[43]’ …/…/preproc_data/121_glm/stats.121+tlrc’[44]’
031 …/…/preproc_data/031_glm/stats.031+tlrc’[43]’ …/…/preproc_data/031_glm/stats.031+tlrc’[44]’
032 …/…/preproc_data/032_glm/stats.032+tlrc’[43]’ …/…/preproc_data/032_glm/stats.032+tlrc’[44]’
033 …/…/preproc_data/033_glm/stats.033+tlrc’[43]’ …/…/preproc_data/033_glm/stats.033+tlrc’[44]’
034 …/…/preproc_data/034_glm/stats.034+tlrc’[43]’ …/…/preproc_data/034_glm/stats.034+tlrc’[44]’
035 …/…/preproc_data/035_glm/stats.035+tlrc’[43]’ …/…/preproc_data/035_glm/stats.035+tlrc’[44]’
036 …/…/preproc_data/036_glm/stats.036+tlrc’[43]’ …/…/preproc_data/036_glm/stats.036+tlrc’[44]’
037 …/…/preproc_data/037_glm/stats.037+tlrc’[43]’ …/…/preproc_data/037_glm/stats.037+tlrc’[44]’
038 …/…/preproc_data/038_glm/stats.038+tlrc’[43]’ …/…/preproc_data/038_glm/stats.038+tlrc’[44]’
039 …/…/preproc_data/039_glm/stats.039+tlrc’[43]’ …/…/preproc_data/039_glm/stats.039+tlrc’[44]’
-mask …/…/group_analysis/group_mask_frac0.8+tlrc
-model_outliers
-unequal_variance
-residual_Z

If I add a return between my two sets, I get the following error -
Error in lop$subjLab[[ii]] : subscript out of bounds
Calls: read.MEMA.opts.batch
Execution halter
-set: Command not found.
021: command not found.
-model_outliers: Command not found.

Thanks,
Lea

Lea,

There are a couple of extra spaces after those backslashes. Use the following command on your 3dMEMA script to figure out where those troubling spaces are:

file_tool -test -infile …

Hi Gang,

You are completely right, thank you so much.

Lea