Dear experts,
I'm trying to set up a model to examine whether Age might affect the BOLD signal (8 tent functions) over time and space (i.e., when and where is the BOLD signal associated with Age). This is within one group (N=30) currently and kind of following Gang's Example 3 at [
afni.nimh.nih.gov]. However, I can't get it to work.
Here are two of my attempts:
3dMVM -prefix $output_dir/mvm_time_age \
-mask $output_dir/some_mask+tlrc \
-bsVars 1 \
-wsVars 'Time*Age' \
-qVars 'Age' \
-num_glt 8 \
-gltLabel 1 t0_Age -gltCode 1 'Time : 1*t0 Age : 1* Age' \
-gltLabel 2 t1_Age -gltCode 2 'Time : 1*t1 Age : 1* Age' \
-gltLabel 3 t2_Age -gltCode 3 'Time : 1*t2 Age : 1* Age' \
-gltLabel 4 t3_Age -gltCode 4 'Time : 1*t3 Age : 1* Age' \
-gltLabel 5 t4_Age -gltCode 5 'Time : 1*t4 Age : 1* Age' \
-gltLabel 6 t5_Age -gltCode 6 'Time : 1*t5 Age : 1* Age' \
-gltLabel 7 t6_Age -gltCode 7 'Time : 1*t6 Age : 1* Age' \
-gltLabel 8 t7_Age -gltCode 8 'Time : 1*t7 Age : 1* Age' \
-dataTable \
Subj Time Group Age Gender InputFile \
01 t0 control 67 M $input_dir/01.results/stats.01_REML+tlrc'[bh#0_Coef]' \
01 t1 control 67 M $input_dir/01.results/stats.01_REML+tlrc'[bh#1_Coef]' \
01 t2 control 67 M $input_dir/01.results/stats.01_REML+tlrc'[bh#2_Coef]' \
...
02 t0 control 70 M $input_dir/02.results/stats.02_REML+tlrc'[bh#0_Coef]' \
02 t1 control 70 M $input_dir/02.results/stats.02_REML+tlrc'[bh#1_Coef]' \
02 t2 control 70 M $input_dir/02.results/stats.02_REML+tlrc'[bh#2_Coef]' \
...
### Error message:
Error:
Quitting due to model test failure...
There were 50 or more warnings (use warnings() to see the first 50)
3dMVM -prefix $output_dir/mvm_time_age \
-mask $output_dir/some_mask+tlrc \
-bsVars 1 \
-wsVars 'Time*Age' \
-num_glt 8 \
-gltLabel 1 t0_Age -gltCode 1 'Time : 1*t0 Age : 1* Age' \
-gltLabel 2 t1_Age -gltCode 2 'Time : 1*t1 Age : 1* Age' \
-gltLabel 3 t2_Age -gltCode 3 'Time : 1*t2 Age : 1* Age' \
-gltLabel 4 t3_Age -gltCode 4 'Time : 1*t3 Age : 1* Age' \
-gltLabel 5 t4_Age -gltCode 5 'Time : 1*t4 Age : 1* Age' \
-gltLabel 6 t5_Age -gltCode 6 'Time : 1*t5 Age : 1* Age' \
-gltLabel 7 t6_Age -gltCode 7 'Time : 1*t6 Age : 1* Age' \
-gltLabel 8 t7_Age -gltCode 8 'Time : 1*t7 Age : 1* Age' \
-dataTable \
Subj Time Group Age Gender InputFile \
01 t0 control 67 M $input_dir/01.results/stats.01_REML+tlrc'[bh#0_Coef]' \
01 t1 control 67 M $input_dir/01.results/stats.01_REML+tlrc'[bh#1_Coef]' \
01 t2 control 67 M $input_dir/01.results/stats.01_REML+tlrc'[bh#2_Coef]' \
...
02 t0 control 70 M $input_dir/02.results/stats.02_REML+tlrc'[bh#0_Coef]' \
02 t1 control 70 M $input_dir/02.results/stats.02_REML+tlrc'[bh#1_Coef]' \
02 t2 control 70 M $input_dir/02.results/stats.02_REML+tlrc'[bh#2_Coef]' \
...
### Error message:
Inappropriate coding in test No.1!
I'd be very grateful about pointers where I may have gone wrong.
Thanks a lot in advance!
Nic