I'm not sure I understand, Gang. I have carried out the analysis with percentage change scores (a single score) and also with symptom severity scores (also a single figure) and this worked out. I also followed the example 1 in the 3dRegAna manual. Are the RT scores somehow different? I have multiple RT figures for each subject but the number I used in the script posted was an average of these scores.