我有一个二维矩阵,在第一列中包含实验条件的索引,在第二列中包含相应实验的索引,即[condition experiment] 。每行对应一个有趣的事件(一个实验可以产生一个或多个事件)。
计算条件和事件很容易。我想知道如何计算每种给定条件下有多少个独特的实验。
这是我现在使用
ACCUMARRAY的解决方案,但我认为应该有一个更简单或更优雅的解决方案:
idxList = [1 1;... %# There are two experiments for condition 1... 1 2;... 1 2;... 2 1;... %# ...and 1 experiment for condition 2. 2 1]; accumarray(idxList(:,1),idxList(:,2),[],@(x)length(unique(x))) ans = 2 1
回答:
以下是一些替代方案:
- 使用整个矩阵idxList为subs参数ACCUMARRAY (即指定用于累积行和列索引),然后求和横跨结果的行非零的数目:
experCounts = sum(accumarray(idxList,1) > 0,2);
- 首先在idxList上使用UNIQUE删除重复的行,从而简化了对ACCUMARRAY的调用:
idxList = unique(idxList,'rows'); experCounts = accumarray(idxList(:,1),1);
更多&回答...