<div dir="ltr">Dear FT users,<div><br></div><div>I am writing to request how pairwise phase consistency (PPC) is  implemented in Fieldtrip. </div><div><br></div><div>I trace the code and find that the following lines are relevant for PPC computation</div><div><br></div><div>input    = input./abs(input);  % normalize the crosspectrum</div><div>outsum   = nansum(input);      % compute the sum; this is 1 x size(2:end)</div><div> c        = (outsum.*conj(outsum) - n)./(n*(n-1)); % do the pairwise thing in a handy way</div><div><br></div><div>However, when I try to verify the results by performing the equations indicated by Vinck et. al. (2010), the outcomes are totally different from those performed using ft_connectivity_ppc. </div><div><br></div><div>My code is as follows and many thanks for any help.</div><div><br></div><div>%%%%%%%%%%%%%%%%%%%%%%%%%</div><div>nTrial = size(input,1);</div><div><br></div><div>ppc_all = [];</div><div>for j = 1:(nTrial-1)</div><div>    for k = (j+1):nTrial</div><div>        ppc_temp = input(j,1,:,:).*conj(input(k,1,:,:));</div><div>        ppc_all = [ppc_all;ppc_temp];</div><div>     end</div><div>end</div><div>            </div><div>PPC = squeeze(abs(nansum(ppc_all,1)*2/(nTrial*(nTrial-1))));</div><div>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div><br></div><div>Best,</div><div>Shen-Mou Hsu</div></div>