Hi Genevieve,

This can be achieved with a for-loop.

Here is how to do it. Note that I have not tested this but it should work:

global nbs; N=nbs.STATS.N;

ind_upper=find(triu(ones(N,N),1));
cross_ref=zeros(N,N);
cross_ref(ind_upper)=1:length(ind_upper);
cross_ref=cross_ref+cross_ref';

[i,j]=find(nbs.NBS.con_mat{1});
for n=1:length(i)
     ind=cross_ref(i(n),j(n));
     fprintf('%0.3f\n',nbs.GLM.y(:,ind));
end


Originally posted by gmcphilemy: