ET_PDToolkit/PDToolkit/@PDSession/mean.m

40 lines
1.1 KiB
Matlab

function[mean] = mean(session, varargin)
%% options : '''Trials', [list of trials]}, ...
% '''Type'', {raw|interpolated|filtered}'
% '''BaselineCorrected''0/1'
% '''Marker'', marker id'
% '''Settings''' , settings object
params = struct('Trials', [], ...
'Type', 'none', ...
'BaselineCorrected', 1, ...
'Marker', 'stimulus onset',...
'Settings', []);
params = parseparams(varargin);
if isempty(params.Trials)
params.Trials = 1:length(sessions.trials);
end
if isstruct(params.Settings)
threshold = params.Settings.QualityThreshold;
quality_values = [ session.trials(Trials).quality ];
trials_to_kill = find(quality_values < threshold);
%exclude trials that fail the quality measure
params.Trials = setdiff(params.Trials, trials_to_kill);
end
for t = 1:length(params.Trials)
trial_id = params.Trials(t);
trial = session.trials(trial_id);
[times{t} traces{t}] = trial.getTrialDataShiftedForMarker(marker_id);
end
if nargin == 1