[讨论] 分组求平均问题

楼主: jenny33377 (澄)   2018-12-07 23:14:55
我的问题是我有一组data如下
Record_ProEst_P_and_ASR = [0.2 0.5
0.6 0.6
0.5 0.8
0.9 0.5]
想利用data第一行资料分组
例如分组区间为interval = [0, 0.01, 0.21, 0.41, 0.61, 0.81, 1]
之后判断在data(1,:)落于
0~0.01、0.01~0.21、0.21~0.41、0.41~0.61、0.61~0.81内的data第二行资料平均
也就是会得到[0, 0.5, 0, 0.7, 0, 0,5]
我目前的写法如下
average_Record_ASR = zeros(length(interval)-1,1);
sum_Record_ASR = [];
for i = 1:length(interval)-1
for j = 1:length(Record_ProEst_P_and_ASR)
if Record_ProEst_P_and_ASR(j,1) >= interval(i) & Record_ProEst_P_and_ASR(j,1) < interval(i+1);
sum_Record_ASR = [sum_Record_ASR, Record_ProEst_P_and_ASR(j,2)];
average_Record_ASR(i,1) = sum(sum_Record_ASR)/length(sum_Record_ASR);
end
end
sum_Record_ASR = [];
end
但这样写法等于是每个区间都要判断一次数列,
当我的data有10^8笔的时候,
就跑超久的
不知道有没有更好的办法
有想说先将data进行排序,然后判断过的不要在判断
但实在想不出来QQ
作者: sunev (Veritas)   2018-12-08 04:12:00
data=Record_ProEst_P_and_ASR;[~,~,subs]=histcounts(data(:,1),interval);accumarray(subs,data(:,2),[length(interval)-1 1],@mean)

Links booklink

Contact Us: admin [ a t ] ucptt.com