Difficulties with 3dICC

AFNI Message Board
1

Hi folks,

I am attempting to use the new 3dICC program but running into the following error:


Error in if (!is.na(lop$dataStr$tStat)) inDataV <- array(apply(inDataV,  : 
  argument is of length zero
Execution halted

It apears to start off correctly, here is the full output - though is it missing the contingency tables?


Read 163 items

++++++++++++++++++++++++++++++++++++++++++++++++++++
***** Summary information of data structure *****
20 subjects :  GM1001 GM1002 GM1003 GM1004 GM1005 GM1006 GM1007 GM1009 GM1010 GM1011 GM1012 GM1013 GM1014 GM1015 GM1016 GM1017 GM1018 GM1019 GM1020 GM1021 
40 response values
2 levels for factor Session : V1 V2 

Contingency tables of subject distributions among the categorical variables:

***** End of data structure information *****
++++++++++++++++++++++++++++++++++++++++++++++++++++

Reading input files now...

Reading input files for effect estimates: Done!

Error in if (!is.na(lop$dataStr$tStat)) inDataV <- array(apply(inDataV,  : 
  argument is of length zero
Execution halted

3dICC reads the files and reports all other values correctly. This run is using a data table that worked with the ICC implementation in 3dLME previously. Here is the bulk of the code


3dICC -prefix RightVRest_ICC_test \
-jobs 4        \
-model  "1+(1|Session)+(1|Subj)" \
-mask /Volumes/alexandria/hypnos/morpheus/gait_mate/models/afni/group_mask.nii.gz \
-dataTable \
Subj Session InputFile \
GM1001 V1 /Volumes/alexandria/hypnos/morpheus/gait_mate/models/afni/sub-GM1001_model-icorLVe_ses-one.results/stats.GM1001_ses-one_REML_nods+tlrc.HEAD'[25]' \
GM1001 V2 /Volumes/alexandria/hypnos/morpheus/gait_mate/models/afni/sub-GM1001_model-icorLVe_ses-two.results/stats.GM1001_ses-two_REML_nods+tlrc.HEAD'[25]' \
[...]
GM1021 V2 /Volumes/alexandria/hypnos/morpheus/gait_mate/models/afni/sub-GM1021_model-icorLVe_ses-two.results/stats.GM1021_ses-two_REML_nods+tlrc.HEAD'[25]'

AFNI version details, just updated this am.


afni --version
Precompiled binary macos_10.12_local: Jul 26 2019 (Version AFNI_19.2.06 'Claudius')
2

Thanks for reporting the problem!

There is indeed a bug in the program. Open the file 3dICC.R under your AFNI binary directory (e.g., under ~/abin/), and replace the following string

lop$dataStr$tStat

with

lop$tStat

Save it, and run your script again. I’ll fix it in the next AFNI release.

3

Works perfectly now, thanks for the quick fix!