Hello!

I'm having trouble with 3dLMEr. I am working with a dataset in which each participant watched a movie, and for each movie I have a regressor representing the extracted "sentiment." I have .nii.gz outputs from 3dDeconvolve for each subject, and now I'd like to look at the effect of that sentiment arc across all movies where I specify subjects and movies as random factors.

Here is the script I'm working with so far with a subset of participants...

#!/bin/bash

3dLMEr -prefix LME \
-model '(1|Subj)+(1|movie)' \
-dataTable \
Subj movie InputFile \
sub02 500daysofsummer stats_500daysofsummer.sub-02.nii.gz'[audio_arc#0_Coef]' \
sub04 500daysofsummer stats_500daysofsummer.sub-04.nii.gz'[audio_arc#0_Coef]' \
sub39 usualsuspects stats_usualsuspects.sub-39.nii.gz'[audio_arc#0_Coef]' \
sub41 usualsuspects stats_usualsuspects.sub-41.nii.gz'[audio_arc#0_Coef]' \
sub48 pulpfiction stats_pulpfiction.sub-48.nii.gz'[audio_arc#0_Coef]' \
sub50 pulpfiction stats_pulpfiction.sub-50.nii.gz'[audio_arc#0_Coef]' \
sub51 shawshankredemption stats_shawshankredemption.sub-51.nii.gz'[audio_arc#0_Coef]' \
sub52 shawshankredemption stats_shawshankredemption.sub-52.nii.gz'[audio_arc#0_Coef]' \
sub60 theprestige stats_theprestige.sub-60.nii.gz'[audio_arc#0_Coef]' \
sub62 theprestige stats_theprestige.sub-62.nii.gz'[audio_arc#0_Coef]' \
sub63 backtothefuture stats_backtothefuture.sub-63.nii.gz'[audio_arc#0_Coef]' \
sub68 backtothefuture stats_backtothefuture.sub-68.nii.gz'[audio_arc#0_Coef]' \
sub81 12yearsaslave stats_12yearsaslave.sub-81.nii.gz'[audio_arc#0_Coef]' \
sub86 12yearsaslave stats_12yearsaslave.sub-86.nii.gz'[audio_arc#0_Coef]'

When running this script, I get the following output that shows a few noticeable things... There is nothing listed on "contingency tables...among the categorical variables" and the range of input data shows -0.0, 0.001. It then gets stuck trying to load R packages.

++++++++++++++++++++++++++++++++++++++++++++++++++++
***** Summary information of data structure *****
38 response values

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: [-0.000, 0.001]

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: ‘lmerTest’

The following object is masked from ‘package:lme4’:

lmer

The following object is masked from ‘package:stats’:

step

Loading required package: phia
Loading required package: car
Loading required package: carData
Registered S3 methods overwritten by 'car':
method from
influence.merMod lme4
cooks.distance.influence.merMod lme4
dfbeta.influence.merMod lme4
dfbetas.influence.merMod lme4


I've looked at my individual subject images, and everything looks normal. I've looked into my R installation and everything seems to be installed appropriately (R version 3.6.3 as recommended on the afni installation page). I've tested out a few different versions of quotes around the InputFiles to no avail. Ultimately, I'm not sure if I'm specifying something wrong in the script, there is an R problem, or it's something else entirely! (I'm new to afni). Can you help?

Thank you!
Clare