Hi Chrisine,
You could use the same method that we apply in
AFNI_data6/ht03/@stim_analyze_modern.
Note that the script repeatedly calls make_random_timing.py to generate the timing files, then runs cmd.3dd.$iter.txt (which was created by the prior call to MRT.py) to run "3dDeconvolve -nodata" on them. You could just swap these commands out with 3dDeconvolve -nodata run on your own sets of timing files to generate a similar ranking. The script then tells you how to rank them.
How does that seem?
- rick