ancova on jacobian determinant maps

AFNI version info (afni -ver): AFNI_25.1.15

Dear experts
I have a set of jacobian determinant maps extracted from T1-to-MNI warps, one map per subject. I would like to perform an ANCOVA to test the effect of a quantitative variable (a blood biomarker) on morphometry, while correcting for the effects of sex (categorical) and age (quantitative). I would like to obtain p values for the quantitative variable, and if possible for its interaction with sex.
I see that there are quite a lot of functions to perform GLM / ANCOVA, what would you suggest?
This is the command I'm using at the moment

3dMVM                                                                        \
    -prefix       HCTEMP_reg                                                 \
    -jobs         10                                                         \
    -bsVars       'Sex*Age+tau'                                              \
    -qVars        'Age,tau'                                                  \
    -qVarCenters  '65,0'                                                     \
    -mask         HCTEMP_template0_mask.nii.gz                               \
    -dataTable    @HCTEMP.csv

I would like to obtain p values for the quantitative variable, and if possible for its interaction with sex.

If you want to make inference about the interaction between sex and tau, try the following model:

-bsVars       'Sex*Age+Sex*tau'   \

Gang Chen