span>
0079 [J2,Msp2,Nsp2]=mtfftpb(data2,tapers,nfft); <span class="comment">% fourier transform of point process data</span>
0080 zerosp(Nsp2==0)=1; <span class="comment">% set the trials where no spikes were found to have zerosp=1</span>
0081 J1=J1(findx,:,:);
0082 J2=J2(findx,:,:);
0083 S12=squeeze(mean(conj(J1).*J2,2));
0084 S1=squeeze(mean(conj(J1).*J1,2));
0085 S2=squeeze(mean(conj(J2).*J2,2));
0086 <span class="keyword">if</span> trialave; S12=squeeze(mean(S12,2)); S1=squeeze(mean(S1,2)); S2=squeeze(m