24

这里的这个问题似乎没有帮助：Calculating Percentiles (Ruby)

我想从一组数字中计算第 95 个百分位（或者，实际上，任何其他所需的百分位）。最终，这将在 Rails 中应用，以针对大量记录计算分布。

但是，如果我能确定如何从一组数字中准确地确定给定的百分位数，我就可以从那里得到它。

坦率地说，我很惊讶我还没有找到某种具有这种功能的宝石——我还没有找到。

非常感谢您的帮助。

4

4 回答 4

40

如果要复制 Excel 的 PERCENTILE 函数，请尝试以下操作：

def percentile(values, percentile)
    values_sorted = values.sort
    k = (percentile*(values_sorted.length-1)+1).floor - 1
    f = (percentile*(values_sorted.length-1)+1).modulo(1)
    
    return values_sorted[k] + (f * (values_sorted[k+1] - values_sorted[k]))
end

values = [1, 2, 3, 4]
p = 0.95
puts percentile(values, p)
#=> 3.85

该公式基于 QUARTILE 方法，它实际上只是一个特定的百分位数 - https://support.microsoft.com/en-us/office/quartile-inc-function-1bbacc80-5075-42f1-aed6-47d735c4819d。

于 2012-08-02T20:41:05.827 回答
18

如果您对现有的 gem 感兴趣，descriptive_statistics那么到目前为止我发现 gem 是最好的percentile功能。

内部审查委员会会议

> require 'descriptive_statistics'
=> true
irb(main):009:0> data = [1, 2, 3, 4]
=> [1, 2, 3, 4]
irb(main):010:0> data.percentile(95)
=> 3.8499999999999996
irb(main):011:0> data.percentile(95).round(2)
=> 3.85

gem 的优点在于它描述“我想要 95% 的数据”的优雅方式。

于 2014-04-04T10:05:11.813 回答
7

基于项目数的百分位数

a = [1,2,3,4,5,6,10,11,12,13,14,15,20,30,40,50,60,61,91,99,120]

def percentile_by_count(array,percentile)
  count = (array.length * (1.0-percentile)).floor
  array.sort[-count..-1]
end

# 80th percentile (21 items*80% == 16.8 items are below; pick the top 4)
p percentile_by_count(a,0.8) #=> [61, 91, 99, 120]

基于值范围的百分位数

def percentile_by_value(array,percentile)
  min, max = array.minmax
  range = max - min
  min_value = (max-min)*percentile + min
  array.select{ |v| v >= min_value }
end

# 80th percentile (119 * 80% = 95.2; pick values above this)
p percentile_by_value(a,0.8) #=> [99, 120]

有趣的是，Excel 的PERCENTILE函数返回60第 80 个百分位数的第一个值。如果你想要这个结果——如果你想要包含一个落在限制尖端的项目——那么将.floor上面的内容更改为.ceil.

于 2012-08-02T20:18:16.923 回答
0

这是我在自己的统计库中开发的方法：

def quantiles(data, probs=[0.25, 0.50, 0.75])
  values = data.sort

  probs.map do |prob|
    h = 1 + (values.count - 1) * prob
    mod = h % 1
    (1 - mod) * values[h.floor - 1] + (mod) * values[h.ceil - 1]
  end
end

如果您只想要一个分位数，请执行quantiles(data, [0.95]).

于 2020-10-05T16:04:58.273 回答