SwitchingFCalongAtask/getERP_SUBJwindow.m

function [mWin,sdWin]=getERP_SUBJwindow(data,iXstg,win, onsets, x, y, subj)
nroi    = size(data,2);
nTr     = length(x);
nStg    = length(iXstg);
nSj     = max(subj);
N       = abs(win(1)) + win(2) + 1;

dataWin = zeros(N, nroi,nStg, nTr, 'single');
for tr = 1:nTr,
    dataTr = data(x(tr):y(tr),:);
    for iSt = 1:nStg        
        idXpos = onsets(tr,iXstg(iSt))-win(1) : onsets(tr,iXstg(iSt))+win(2);
        idXpos(idXpos<1) = 1;

        data_ = dataTr(idXpos, :) ;       
        dataWin(:,:,iSt,tr) = data_;
    end
end
dataSj_ = zeros(N, nroi, nStg, nSj, 'single');
for sj = 1:nSj,    
    dataSj_(:,:,:,sj) = squeeze( mean(dataWin(:,:,:, subj==sj),4));
end
mWin  = squeeze(mean(dataSj_,4));
sdWin = squeeze(std(abs(dataSj_),[],4)) ./ sqrt(nSj);