在 MATLAB 中查找每个单元格的阈值以上元素的频率

Question

我有一个 4-D 矩阵。维度是经度、纬度、天、年，如 [17,14,122,16]。 我必须找出每个单元格的值高于 98% 的频率，以便最终输出以 17x14 数组的形式出现，其中包含高于 98% 阈值的值的出现次数。

我做了一些事情，它为我提供了一个 17x14 的矩阵，该矩阵包含与每个单元格的 98 个百分位数相关的值，但我无法确定出现的频率。

k=0;
p=cell(1,238);
r=cell(1,238);

for i=1:17
   for j=1:14
      n=m(i,j,[1:122],[1:16]);
      n=squeeze(n);
      k=k+1;
      q=prctile(n(:),98);
      r{k}=nansum(nansum(n>=q));
      p{k}=q;
   end
end

此代码给出矩阵p 很好，但矩阵r 包含所有单元格的相同值。这怎么可能？我做错了什么？请帮忙。

Answer 1

根据定义，第 98 个百分位以上的值的频率为 2%。

我猜你得到的 r 的值是 39； 122x16 矩阵前 2% 的元素数量（即 1952 个元素）。

r = 0.02*1952; 

r = 
      39.040 

您的代码正在验证理论值。也许您正在考虑一个不同的问题？

这是一个模拟示例，使用从 0 到 100 的随机生成（均匀分布）数据 (n)。

p=cell(1,238);
r=cell(1,238);
for i=1:17
    for j=1:14
        %         n=m(i,j,[1:122],[1:16]);
        %         n=squeeze(n);

        % After you do n=squeeze(n), it gives 2-D matrix of 122x16
        % dimensions.
        n = rand(122,16)*100;  % simulation for your 2-D matrix
        k=k+1;
        q=prctile(n(:),98);
        r{k}=nansum(nansum(n>=q));
        p{k}=q;
    end
end