wolffd@0: function r = corrcoef(x,y)
wolffd@0: 
wolffd@0: dx = get(x,'Data');
wolffd@0: dy = get(y,'Data');
wolffd@0: fx = get(x,'FramePos');
wolffd@0: fy = get(y,'FramePos');
wolffd@0: dr = cell(1,length(fx));
wolffd@0: for h = 1:length(fx)
wolffd@0:     dr{h} = cell(1,length(fx{h}));
wolffd@0:     for i = 1:length(fx{h})
wolffd@0:         [f,ix,iy] = intersect(fx{h}{i}',fy{h}{i}','rows');
wolffd@0:         dxi = dx{h}{i}(:,ix,:,:);
wolffd@0:         dyi = dy{h}{i}(:,iy,:,:);
wolffd@0:         dr{h}{i} = corrcoef(dxi,dyi);
wolffd@0:     end
wolffd@0: end
wolffd@0: r = mirmatrix(x,'Data',dr,'Title','Correlation coefficients','Unit','');