Hello,
I'm seeking a straightforward method for generating peristimulus plots from a first level design run in FEAT. That is, I want to see the percent signal change within a particular ROI, time-locked to an event onset, with a separate line for each condition. I understand that featquery does this to some extent (using tsplot), but doesn't give you the ability to plot a separate line for each condition, and only generates static plots. This seems like a somewhat common thing people would do, but I've had a hard time finding a solution for FSL.
I've managed to use the extracted timecourses after running featquery (mean_mask_ts.txt) in combination with the ev*.txt files under custom_timing_files in order to do simple event-related averaging of the data. The problem is, we have a relatively quick event related design, and I don't think simple averaging is what we want. Everything I've read suggests estimating a FIR model in order to generate the peristimulus plots (e.g., Serences 2004, Neuroimage).
It looks like Marsbar can do this fairly easily, IF you've processed everything in SPM beforehand. I've tried my best to hack something together to import the FSL design into matlab and make it work with marsbar, but it's hard to determine the minimum amount of information Marsbar would need in order for this to work (since it only asks for an SPM.mat file).
It also looks like AFNI's 3dDeconvolve should be able to do this too, but I'm not familiar with it.
Rather than trying to reinvent the wheel, I thought I'd reach out to the FSL community and see if someone has come up with a more standard solution. Ideally we'd need something that could be implemented in Matlab or Python so we can have some flexibility with the data.
Thank you for any help you can offer!
Jason
|