Comrades--
We'd like to do some multi-voxel pattern analysis using data from a simple block-design study with baseline epochs long enough for the HRF to settle. One fork in the decision tree for implementing this approach is whether to use for training and testing a series of betas from individual blocks of each condition (using, e.g., 3dLSS) or to use the preprocessed BOLD time-series. We're opting for the latter since we have lots of TRs and not-so-many blocks per condition. With that decided, we need for training and testing to wed the actual task timing to the BOLD data which means deconvolving the HRF out of neural activity. To do this, would you recommend using 3dTfitter?
All best,
Paul