Hello AFNI Experts,

I have a mixed design of Group (hc, sci) x Run (run1, run2) for which I wanted to use 3dLMEr (in case I don't have a balanced design in the end for 3dANOVA).

Below is the code that I have tested:

derivatives=/mnt/LABS/tbi/implicitmem_SCI/derivatives
outdir=$derivatives/3dLME
indir=/mnt/LABS/tbi/implicitmem_SCI/derivatives/feat_lev1_hp100
3dLMEr \
  -prefix ${outdir}/hp100 \
  -jobs 2 \
  -mask $FSLDIR/data/standard/MNI152_T1_2mm_brain.nii.gz \
  -model 'run*group+(1|Subj)' \
  -dataTable \
      Subj  run   group InputFile \
      sub-IPLHC201      run1  hc    $indir/sub-IPLHC201_run-1.feat/stats/cope4.nii.gz \
      sub-IPLHC201      run2  hc    $indir/sub-IPLHC201_run-2.feat/stats/cope4.nii.gz \
      sub-IPLHC202      run1  hc    $indir/sub-IPLHC202_run-1.feat/stats/cope4.nii.gz \
      sub-IPLHC202      run2  hc    $indir/sub-IPLHC202_run-2.feat/stats/cope4.nii.gz \
      sub-IPLHC203      run1  hc    $indir/sub-IPLHC203_run-1.feat/stats/cope4.nii.gz \
      sub-IPLHC203      run2  hc    $indir/sub-IPLHC203_run-2.feat/stats/cope4.nii.gz \
      sub-IPLHC206      run1  hc    $indir/sub-IPLHC206_run-1.feat/stats/cope4.nii.gz \
      sub-IPLHC206      run2  hc    $indir/sub-IPLHC206_run-2.feat/stats/cope4.nii.gz \
      sub-IPLSCI102     run1  sci   $indir/sub-IPLSCI102_run-1.feat/stats/cope4.nii.gz \
      sub-IPLSCI102     run2  sci   $indir/sub-IPLSCI102_run-2.feat/stats/cope4.nii.gz \
      sub-IPLSCI103     run1  sci   $indir/sub-IPLSCI103_run-1.feat/stats/cope4.nii.gz \
      sub-IPLSCI103     run2  sci   $indir/sub-IPLSCI103_run-2.feat/stats/cope4.nii.gz \
      sub-IPLSCI104     run1  sci   $indir/sub-IPLSCI104_run-1.feat/stats/cope4.nii.gz \
      sub-IPLSCI104     run2  sci   $indir/sub-IPLSCI104_run-2.feat/stats/cope4.nii.gz \
      sub-IPLSCI105     run1  sci   $indir/sub-IPLSCI105_run-1.feat/stats/cope4.nii.gz \
      sub-IPLSCI105     run2  sci   $indir/sub-IPLSCI105_run-2.feat/stats/cope4.nii.gz \

All of the required R packages seemed to be able to load without issue, but afterward, I received this error: missing value where TRUE/FALSE needed

++++++++++++++++++++++++++++++++++++++++++++++++++++
***** Summary information of data structure *****
16 response values
2 levels for factor group : hc sci

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!


Range of input data: [-5181.417, 2601.186]

If the program hangs here for more than, for example, half an hour,
kill the process because the model specification or something else
is likely inappropriate.

Loading required package: lmerTest
Loading required package: lme4
Loading required package: Matrix

Attaching package: \u2018lmerTest\u2019

The following object is masked from \u2018package:lme4\u2019:

  lmer

The following object is masked from \u2018package:stats\u2019:

  step

Loading required package: phia
Loading required package: car
Loading required package: carData

Error in if (ii < dimx) ii <- ii + 1 else if (jj < dimy) { :
  missing value where TRUE/FALSE needed
Execution halted

I saw an user reported the same issue in this thread, but wasn't able to follow up because the topic has been closed: [afni.nimh.nih.gov]

I'm using afni 22.1.14 and R 4.0.4. The only R package that I wasn't able to install is 'brms' because there were issues installing igraph in R 4.0. Please let me know if you need any other information.

Thank you,
Tien