timecourse extraction for single trials
Joshua Buckholtz <joshua.buckholtz <at> vanderbilt.edu>
2009-07-14 07:28:49 GMT
Dear Matthew and others,
We have a slow event-related task with 60 trials, divided into 3 conditions. We’ve had luck using Marsbar to extract, for each subject, timecourses for each condition. However, we’re also interested in extracting, for each subject, timecourses for each of the 60 single trials. In thinking about how to do this, we’ve had the idea of creating a design matrix with 60 predictors (i.e. each trial onset has its own vector), and using the following code to extract FIRs for each event.
% Get definitions of all events in model
[e_specs, e_names] = event_specs(E);
n_events = size(e_specs, 2);
% Bin size in seconds for FIR
bin_size = tr(E);
% Length of FIR in seconds
fir_length = 24;
% Number of FIR time bins to cover length of FIR
bin_no = fir_length / bin_size;
% Options - here 'single' FIR model, return estimated
opts = struct('single', 1, 'percent', 1);
% Return time courses for all events in fir_tc matrix
for e_s = 1:n_events
fir_tc(:, e_s) = event_fitted_fir(E, e_specs(:,e_s), bin_size, ...
bin_no, opts);
end
In your opinion(s), is this a valid means of obtaining our desired result?
Many thanks,
Josh
------------------------------------------------------------------------------
Enter the BlackBerry Developer Challenge
This is your chance to win up to $100,000 in prizes! For a limited time,
vendors submitting new applications to BlackBerry App World(TM) will have
the opportunity to enter the BlackBerry Developer Challenge. See full prize
details at: http://p.sf.net/sfu/Challenge
_______________________________________________
Marsbar-users mailing list
Marsbar-users <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/marsbar-users