Hi Gang thanks for answering me, I feel very confused with the design of my 3dMVM script, would you mind to help me a bit?
I want to see the relationship between the CBF_maps (voxel wise, .nii file per subject) with a quantitative clinical measure (loneliness, which is one single number per subject). I need to control for age, gender, depression scores and gray matter concentration (the gray matter concentration is a voxel wise covariate).

Dependent variable: voxel wise CBF
Independent variable: Loneliness (clinical measure)
Covariates: age, gender, depression and voxel wise gray matter concentration.

I have done this part of the script but I feel lost with the other part

3dMVM -prefix 3dMVM_test -jobs 4 \
-bsVars 'loneliness' \
-mask GM2mm_WFU_Mask.nii.gz \
-qVars "age,gender,hamd" \
-vVars 'gmc' \


I also create the txt file of the data and it looks something like this:

Subj age gender hamd hama gmc loneliness group InputFile
172_120916 68 1 21 9 c1w172_120916.nii 37 1 ASL_CBF_2MNI_6mm_172_120916.nii
222_103015 77 -1 19 13 c1w222_103015.nii 36 1 ASL_CBF_2MNI_6mm_222_103015.nii
225_010416 65 1 2 1 c1w225_010416.nii 23 1 ASL_CBF_2MNI_6mm_225_010416nii

Would you mind to help me with some guidance to build the script?

Thank you