Error in if (len%%wd != 0) errex.AFNI(paste("The content under -dataTable is not rectangular !", :

defenderfer · August 31, 2018, 3:01pm

Thanks for any help. I am using a mac and tcsh shell.

My MVM file is titled MVM1.txt file using the command: tcsh -x MVM1.txt &. I receive this message:

3dMVM -prefix Accuracy -jobs 4
Error in if (len%%wd != 0) errex.AFNI(paste(“The content under -dataTable is not rectangular !”, :
argument is of length zero
Calls: process.MVM.opts
Execution halted

Everything else is followed by “Command not found.” Ive checked my syntax, but I don’t know what needs to be fixed. Here is my script; Ive only included 3 subjects for the example.

#!/bin/tcsh
3dMVM -prefix Accuracy -jobs 4
-model 1
-wsVars “oxyaccnoiseType”
-num_glt 3
-gltLabel 1 oxyEff -gltCode 1 ‘oxy : 1HbO -1HbR’
-gltLabel 2 accEff -gltCode 2 ‘acc : 1right -1wrong’
-gltLabel 3 noiseTypeEff -gltCode 3 ‘noiseType : 1vocod -1noise’
-dataTable
Subj oxy acc noiseType InputFile
S1001 HbO wrong vocod /p2process/S1001/viewer/Subject/1001_cond7_unmaskedOxy.nii
S1001 HbO wrong noise /p2process/S1001/viewer/Subject/1001_cond6_unmaskedOxy.nii
S1001 HbO right vocod /p2process/S1001/viewer/Subject/1001_cond5_unmaskedOxy.nii
S1001 HbO right noise /p2process/S1001/viewer/Subject/1001_cond4_unmaskedOxy.nii
S1001 HbR wrong vocod /p2process/S1001/viewer/Subject/1001_cond7_unmaskedDeoxy.nii
S1001 HbR wrong noise /p2process/S1001/viewer/Subject/1001_cond6_unmaskedDeoxy.nii
S1001 HbR right vocod /p2process/S1001/viewer/Subject/1001_cond5_unmaskedDeoxy.nii
S1001 HbR right noise /p2process/S1001/viewer/Subject/1001_cond4_unmaskedDeoxy.nii
S1003 HbO wrong vocod /p2process/S1003/viewer/Subject/1003_cond7_unmaskedOxy.nii
S1003 HbO wrong noise /p2process/S1003/viewer/Subject/1003_cond6_unmaskedOxy.nii
S1003 HbO right vocod /p2process/S1003/viewer/Subject/1003_cond5_unmaskedOxy.nii
S1003 HbO right noise /p2process/S1003/viewer/Subject/1003_cond4_unmaskedOxy.nii
S1003 HbR wrong vocod /p2process/S1003/viewer/Subject/1003_cond7_unmaskedDeoxy.nii
S1003 HbR wrong noise /p2process/S1003/viewer/Subject/1003_cond6_unmaskedDeoxy.nii
S1003 HbR right vocod /p2process/S1003/viewer/Subject/1003_cond5_unmaskedDeoxy.nii
S1003 HbR right noise /p2process/S1003/viewer/Subject/1003_cond4_unmaskedDeoxy.nii
S1004 HbO wrong vocod /p2process/S1004/viewer/Subject/1004_cond7_unmaskedOxy.nii
S1004 HbO wrong noise /p2process/S1004/viewer/Subject/1004_cond6_unmaskedOxy.nii
S1004 HbO right vocod /p2process/S1004/viewer/Subject/1004_cond5_unmaskedOxy.nii
S1004 HbO right noise /p2process/S1004/viewer/Subject/1004_cond4_unmaskedOxy.nii
S1004 HbR wrong vocod /p2process/S1004/viewer/Subject/1004_cond7_unmaskedDeoxy.nii
S1004 HbR wrong noise /p2process/S1004/viewer/Subject/1004_cond6_unmaskedDeoxy.nii
S1004 HbR right vocod p2process/S1004/viewer/Subject/1004_cond5_unmaskedDeoxy.nii
S1004 HbR right noise /p2process/S1004/viewer/Subject/1004_cond4_unmaskedDeoxy.nii \

Gang · August 31, 2018, 7:00pm

I suspect that there is some extra space after one of the backslashes that is causing the grief. Use the following command

file_tool -test -infile …

with your 3dMVM script and see if you can find the problem.

defenderfer · September 3, 2018, 12:40pm

Thank you for your response. I used the command to fix the script. Now, it will run until it says “Reading input files: Done!” but quickly followed by this error message:

** Error:
Failed at GLT No. 1! Make sure that the model or GLT specification syntax is correct.
Warning message:
In summary.Anova.mlm(object$Anova, multivariate = FALSE) :
Singular error SSP matrix:
non-sphericity test and corrections not available

Gang · September 4, 2018, 2:30pm

The problem is that there are only 3 subjects with the current script, which would not be enough to solve the system. You need at least 10 subjects in this case.

defenderfer · September 4, 2018, 5:24pm

Thanks again for replying. I’m sorry for the confusion- I only listed three because I have 38 subjects and didn’t think it was necessary for the post. I’ve copied and pasted below the output with the final error message. I’m still missing something, and I cant get past that first GLT 1 error.

***** Summary information of data structure *****
38 subjects : S1001 S1003 S1004 S1005 S1006 S1007 S1008 S1009 S1010 S1011 S1012 S1013 S1014 S1016 S1017 S1018 S1019 S1020 S1021 S1022 S1023 S1024 S1025 S1026 S1027 S1028 S1029 S1030 S1031 S1032 S1033 S1034 S1035 S1036 S1037 S1038 S1039 S1040
304 response values
2 levels for factor oxy : HbO HbR
2 levels for factor acc : right wrong
2 levels for factor noiseType : noise vocod
3 post hoc tests

Contingency tables of subject distributions among the categorical variables:

, , noiseType = noise

acc

oxy right wrong
HbO 38 38
HbR 38 38

, , noiseType = vocod

acc

oxy right wrong
HbO 38 38
HbR 38 38

Tabulation of subjects against each of the categorical variables:

lop$nSubj vs oxy:
       
        HbO HbR
  S1001   4   4
  S1003   4   4
  S1004   4   4
  S1005   4   4
  S1006   4   4
  S1007   4   4
  S1008   4   4
  S1009   4   4
  S1010   4   4
  S1011   4   4
  S1012   4   4
  S1013   4   4
  S1014   4   4
  S1016   4   4
  S1017   4   4
  S1018   4   4
  S1019   4   4
  S1020   4   4
  S1021   4   4
  S1022   4   4
  S1023   4   4
  S1024   4   4
  S1025   4   4
  S1026   4   4
  S1027   4   4
  S1028   4   4
  S1029   4   4
  S1030   4   4
  S1031   4   4
  S1032   4   4
  S1033   4   4
  S1034   4   4
  S1035   4   4
  S1036   4   4
  S1037   4   4
  S1038   4   4
  S1039   4   4
  S1040   4   4

lop$nSubj vs acc:

    right wrong

S1001 4 4
S1003 4 4
S1004 4 4
S1005 4 4
S1006 4 4
S1007 4 4
S1008 4 4
S1009 4 4
S1010 4 4
S1011 4 4
S1012 4 4
S1013 4 4
S1014 4 4
S1016 4 4
S1017 4 4
S1018 4 4
S1019 4 4
S1020 4 4
S1021 4 4
S1022 4 4
S1023 4 4
S1024 4 4
S1025 4 4
S1026 4 4
S1027 4 4
S1028 4 4
S1029 4 4
S1030 4 4
S1031 4 4
S1032 4 4
S1033 4 4
S1034 4 4
S1035 4 4
S1036 4 4
S1037 4 4
S1038 4 4
S1039 4 4
S1040 4 4

lop$nSubj vs noiseType:
       
        noise vocod
  S1001     4     4
  S1003     4     4
  S1004     4     4
  S1005     4     4
  S1006     4     4
  S1007     4     4
  S1008     4     4
  S1009     4     4
  S1010     4     4
  S1011     4     4
  S1012     4     4
  S1013     4     4
  S1014     4     4
  S1016     4     4
  S1017     4     4
  S1018     4     4
  S1019     4     4
  S1020     4     4
  S1021     4     4
  S1022     4     4
  S1023     4     4
  S1024     4     4
  S1025     4     4
  S1026     4     4
  S1027     4     4
  S1028     4     4
  S1029     4     4
  S1030     4     4
  S1031     4     4
  S1032     4     4
  S1033     4     4
  S1034     4     4
  S1035     4     4
  S1036     4     4
  S1037     4     4
  S1038     4     4
  S1039     4     4
  S1040     4     4

***** 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.


** Error: 
   Failed at GLT No. 1! Make sure that the model or GLT specification syntax is correct.
Warning message:
In summary.Anova.mlm(object$Anova, multivariate = FALSE) :
  Singular error SSP matrix:
non-sphericity test and corrections not available