Plotting trial-related blood flow (FSL)

You need to extract the peristimulus timing data. This is stored in the ‘tsplot’ folder output into your .feat directory. A great overview of the output from FEAT is stored here.

Specifically, this folder contains:

Reference: What the tsplots are:

tsplot/tsplot_zstat1.gif the full model vs data plot for the maximum Z statistic voxel from contrast 1.

tsplot/tsplot_zstat1p.gif the plot of reduced data vs cope partial model fit – i.e. data-full_model+partial_model vs partial_model.

tsplot/tsplot_zstat1.txt text file of values used for the above plots, for the maximum Z statistic voxel from contrast 1. The first column is the data, the second column is the partial model fit for contrast 1, the third column is the full model fit and the fourth column is the reduced data for contrast 1.

tsplot/tsplotc_zstat1.* plots and text file as above, but instead of using the peak Z stat voxel, here the data is averaged over all voxels in all significant clusters for contrast 1.

ps_tsplot_zstat1_ev1.gif (etc); peristimulus text files and plots, showing model fits averaged over multiple stimulus repeats and “data scatter” over these repeats. The format for the text files is the same as above, except for the insertion of an extra first column which encodes the peristimulus timing (i.e., time relative to the start time within the peristimulus window). The _ev1 part of the filename means that this file/plot relates to EV1; separate plots and text files are generated for each EV for every contrast.
You can extract and reconstruct the peristimulus data in excel. For description, see Benjamin & Gaab, 2012.

I’m sure there is an easier way to extract this data using awk recursively, even if it is just to pull all the relevant tsplot_zstat columns into a single file, but I have not yet done this.

Leave a comment