HI Gang,
I was using 3 subjects in each group simply to test the Script. I have 3 groups all with more than 20 subjects. I have acquired a single resting scan for each subject and computed Z-score coverted pearson correlation maps from a single ROI for each subject (input files). The goal is to compare all 3 groups using all combinations (g1 vs g2, g2 vs g3, g1 vs g3) while covarying for age and bmi across subjects.
Here is an example of the script I created.
3dMVM -prefix MVMtest1 -jobs 1 \
-mask /Users/jessedesimone/Desktop/MVM_analysis/Templates/network_masks_thr005/DAN_12mm_ntwk+tlrc \
-bsVars "group*age+bmi" \
-qVars age,bmi \
-num_glt 3 \
-gltLabel 1 g1_vs_g2 -gltCode 1 'group : 1*g1 -1*g2' \
-gltLabel 2 g1_vs_g3 -gltCode 2 'group : 1*g1 -1*g3' \
-gltLabel 3 g2_vs_g3 -gltCode 3 'group : 1*g2 -1*g3' \
-dataTable @/Users/jessedesimone/Desktop/MVM_analysis/flists/datatab.txt
Here is the input for datatab.txt. These are manufactured values for age and bmi - simply for testing purposes before input of the actual values.
Subj group age bmi InputFile \
subj1 g1 11 20 subj1.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj2 g1 12 21 subj2.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj3 g1 13 22 subj3.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj4 g2 11 20 subj4.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj5 g2 12 21 subj5.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj6 g2 13 22 subj6.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj7 g3 11 20 subj7.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj8 g3 12 21 subj8.DAN_12mm_RS+tlrc.HEAD[PearCorr#0] \
subj9 g3 13 22 subj9.DAN_12mm_RS+tlrc.HEAD[PearCorr#0]
