> the resulting individual average beta values from 3dmaskave plotted with the behavioral scores don't match the polarity of the clusters from the regression.
In what sense they don't match? Could you show (1) the scatter plot of contrast vs behavioral score at one ROI, and (2) the regression coefficient from 3dttest++ at that ROI?
Gang