x'*x;</span>
0061 <span class="comment">%</span>
0062 <span class="comment">% y=zeros(size(x));</span>
0063 <span class="comment">% for k=1:n</span>
0064 <span class="comment">%     idx=[1:k-1,k+1:n];</span>
0065 <span class="comment">%     [topcs,eigenvalues]=nt_pcarot(cc(idx,idx)); % PCA to orthogonalize the other channels</span>
0066 <span class="comment">%     topcs=topcs(:,1:60);</span>
0067 <span class="comment">%     %topcs=topcs(:,find(eigenvalues/max(eigenvalues)&gt;thresh));</span>
0068 <span clas