在矩阵中查找唯一对

Question

我们假设我有以下矩阵：

A = [1 1 2 1;1 2 2 1;2 1 3 0;2 2 2 0;3 1 2 1]

第一列是索引而下两个是交互，最后一个是逻辑，表示是或否。 所以我知道我想根据相互作用生成以下热图。 “X”轴表示相互作用，“Y”轴表示指数。

   1-2  1-3  2-2
1   1   NaN   1
2  NaN   0    0
3   1   NaN  NaN

我目前的做法：

B = sortrows(A,[2,3]);

然后我分别为每行和每列应用find。

是否有一个类似于unique的函数可以同时检查两列？

Answer 1

这是一种使用unique(...,'rows') ：

A = [1 1 2 1; 1 2 2 1; 2 1 3 0; 2 2 2 0; 3 1 2 1]; % data
[~, ~, jj] = unique(A(:,[2 3]),'rows'); % get interaction identifiers
B = accumarray([A(:,1) jj], A(:,4), [], @sum, NaN); % build result, with NaN as fill value

这给了

B =
     1   NaN     1
   NaN     0     0
     1   NaN   NaN

Answer 2

>> A

A =

     1     1     2     1
     1     2     2     1
     2     1     3     0
     2     2     2     0
     3     1     2     1

>> [C, IA, IC] = unique(A(:, [2, 3]), 'rows')

C =

     1     2
     1     3
     2     2

IA =

     1
     3
     2

IC =

     1
     3
     2
     3
     1

C是一组独特的对。 IA是C的对应索引（即， C == A(IA, [2, 3]) ）。 IC是每行的相应索引（即A(:, [2, 3]) == C(IC, :) ）。

Answer 3

借助@Jeon的答案（ 更新 ），这是一个可能的解决方案：

A = [1 1 2 1;1 2 2 1;2 1 3 0;2 2 2 0;3 1 2 1]
[~,IA,idx] = unique(A(:, [2, 3]), 'rows');
r = max(A(:,1));
c = numel(IA);
out= NaN(r,c );
out(sub2ind([r ,c], A(:,1),idx)) = A(:,4)