significant-spike-field-coherenceMBL2005

0001 [fname,pname]=uigetfile('*');
0002 eval(['load ' pname '\' fname]);
0003 params.tapers=[5 9]; params.pad=2; params.Fs=1000; params.fpass=[0 100]; params.trialave=1;params.err=[1 0.05];
0004 movingwin=[0.5 0.1]; % duration of moving window used to evaluate spectrograms
0005 win=[2 2]; % window around events
0006 fignum=[2 3 6 9 8 7 4 1]; % figure numbers for particular directions
0007 for targ=0:7;
0008     E=targoff(find(targets==targ));
0009     datasp=createdatamatpt(dsp(1),E,win);
0010     datalfp=createdatamatc(dlfp(:,1),E,Fs,win);
0011     [C,phi,S12,S1,S2,t,f,zerosp,confC,phierr]=cohgramcpt(datalfp,datasp,movingwin,params); % Note cohgramcpt does give you spectra but not the errors.
0012     subplot(3,3,fignum(targ+1));plotsig(C,confC,t,f); axis xy; colorbar;
0013 end;