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 https://afni.nimh.nih.gov/sscc/gangc/BF.html. 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 ‘TimeAge’
-qVars ‘Age’
-num_glt 8
-gltLabel 1 t0_Age -gltCode 1 'Time : 1t0 Age : 1* Age’
-gltLabel 2 t1_Age -gltCode 2 ‘Time : 1t1 Age : 1 Age’
-gltLabel 3 t2_Age -gltCode 3 ‘Time : 1t2 Age : 1 Age’
-gltLabel 4 t3_Age -gltCode 4 ‘Time : 1t3 Age : 1 Age’
-gltLabel 5 t4_Age -gltCode 5 ‘Time : 1t4 Age : 1 Age’
-gltLabel 6 t5_Age -gltCode 6 ‘Time : 1t5 Age : 1 Age’
-gltLabel 7 t6_Age -gltCode 7 ‘Time : 1t6 Age : 1 Age’
-gltLabel 8 t7_Age -gltCode 8 ‘Time : 1t7 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 ‘TimeAge’
-num_glt 8
-gltLabel 1 t0_Age -gltCode 1 'Time : 1t0 Age : 1* Age’
-gltLabel 2 t1_Age -gltCode 2 ‘Time : 1t1 Age : 1 Age’
-gltLabel 3 t2_Age -gltCode 3 ‘Time : 1t2 Age : 1 Age’
-gltLabel 4 t3_Age -gltCode 4 ‘Time : 1t3 Age : 1 Age’
-gltLabel 5 t4_Age -gltCode 5 ‘Time : 1t4 Age : 1 Age’
-gltLabel 6 t5_Age -gltCode 6 ‘Time : 1t5 Age : 1 Age’
-gltLabel 7 t6_Age -gltCode 7 ‘Time : 1t6 Age : 1 Age’
-gltLabel 8 t7_Age -gltCode 8 ‘Time : 1t7 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