AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

|
nic
December 19, 2016 03:55PM
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
Subject Author Posted

3dMVM one group two within subj variables

nic December 19, 2016 03:55PM

Re: 3dMVM one group two within subj variables

gang December 20, 2016 08:15AM

Re: 3dMVM one group two within subj variables

nic December 20, 2016 11:20AM

Re: 3dMVM one group two within subj variables

gang December 20, 2016 12:30PM