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