Dear all,

I would like to use @stim_analyze or @stim_analyze_modern to iteratively create & evaluate (via 3dDeconvolve -nodata) random event-related timing parameters in order to maximize efficiency & power of our design. I have 2 questions -- 1) what are the differences between the two scripts? Am I correct to think the only differences are in the make_random_timing.py arguments, with the _modern version taking advantage of the timing classes functionality? and 2) in our experiments we use a MION contrast agent that significantly alters the HRF by lengthening time to return to baseline (recovery is on the order of 30-40, rather than 10-20 seconds). How can I incorporate this aspect of the response function into the pipeline?

thanks much.
best,
Geena