123456789101112131415161718192021222324 |
- 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);
|