Honorable AFNI gurus:

I'm currently wanting to look at relations between a single non-neural variable (blood cytokine levels) and a voxel-wise estimate of neural inflammation but I need to first account for effects of both a voxel-wise noise covariate (VBM Jacobians) and a single motion estimate. Currently I'm doing this by using 3dttest++ in "just -setA" mode where the neural inflammation volumes are the DV of interest and the covariates option is used with a covariate file containing three columns (in addition to the subject identifier) for 1) the cytokine; 2) motion; and 3) the file name of the VBM Jacobians for each subject. Is this the proper way to proceed? I think so but I wanted to check.

Now, assuming my approach above is correct, let's say I found a significant relation between my two variables of interest (cytokines and neural inflammation, accounting for the noise covariate effects) in the resulting statistical map AND I wanted to show this relation in a scatterplot accompanying the statistical map. Could you suggest a path for doing this using the input and output data?

Big thanks!!

Paul