afni.nimh.nih.gov

You can use 3dRegAna or 3dttest++.

If you use 3dRegAna, you can run a simple regression with one predictor, then take the square root of the R^2 term and take the sign of the b-value coefficient.

With 3dttest++, you can use the covariate file.