Accumarray 如何设置“有趣”以仅使用每个 bin 的最后一个 obs

Question

有没有办法让accumarray 放弃所有观察，但每组的最后一个？

我想到的是类似的东西：

lastobs=accumarray(bin,x,[],@(x){pick the observation with the max index in each group};  

举个例子，假设我有以下内容：

bin=[1 2 3  3 3 4 4];  %#The bin where the observations should be put
x= [21 3 12 5 6 8 31]; %#The vector of observations
%#The output I would like is as follow  
lastobs=[21 3 6 31];  

我实际上只是在考虑accumarray，因为我只是用它来计算每个 bin 的观察值的平均值。所以每一个能成功的函数对我来说都很好。

Answer 1

当然，您可以使用accumarray 执行此操作。 x(end) 是数组中的最后一个观察值。请注意，bin 需要对其进行排序才能使其正常工作，因此如果不是，请运行 [bin,sortIdx]=sort(bin);x = x(sortIdx); 首先。

lastobs = accumarray(bin(:),x(:),[],@(x)x(end)); %# bin, x, should be n-by-1

Answer 2

您已经得到了accumarray 的答案，但是由于您正在寻找可以完成这项工作的任何解决方案，请考虑以下unique 的应用。

使用unique 和'legacy' 选项可以根据需要给出每个值最后次出现的索引：

>> [~,ia] = unique(bin,'legacy')
ia =
     1     2     5     7
>> lastobs = x(ia)
lastobs =
    21     3     6    31

现在，我喜欢 accumarray，这里的很多人都知道，但我实际上更喜欢这个解决方案。