> I expect these behavioral scores to differ a little between the three
You may consider centering the covariate per each difficulty so that it would not ruin the interpretation of cross-difficulty effects.
> Is 3dMVM the best way to do this?
Use 3dLME, and check out example 2 in the help:
3dLME -help | less
Gang