百分位数聚合编辑

一种 多值 指标聚合,用于计算从聚合文档中提取的数值的一个或多个百分位数。这些值可以从文档中的特定数值或 直方图字段 中提取。

百分位数显示观察值中一定百分比出现的点。例如,第 95 个百分位数是大于 95% 观察值的值。

百分位数通常用于查找异常值。在正态分布中,第 0.13 个和第 99.87 个百分位数表示与平均值相差三个标准差。任何超出三个标准差的数据通常被认为是异常值。

检索百分位数范围时,可以使用它们来估计数据分布并确定数据是否偏斜、双峰等。

假设您的数据由网站加载时间组成。平均加载时间和中位加载时间对管理员来说并不是很有用。最大值可能很有趣,但它很容易被单个缓慢的响应所扭曲。

让我们看一下代表加载时间的一系列百分位数

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    aggregations: {
      load_time_outlier: {
        percentiles: {
          field: 'load_time'
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "aggs": {
    "load_time_outlier": {
      "percentiles": {
        "field": "load_time" 
      }
    }
  }
}

字段 load_time 必须是数字字段

默认情况下,percentile 指标将生成一系列百分位数:[ 1, 5, 25, 50, 75, 95, 99 ]。响应将如下所示

{
  ...

 "aggregations": {
    "load_time_outlier": {
      "values": {
        "1.0": 10.0,
        "5.0": 30.0,
        "25.0": 170.0,
        "50.0": 445.0,
        "75.0": 720.0,
        "95.0": 940.0,
        "99.0": 980.0
      }
    }
  }
}

如您所见,聚合将返回默认范围内每个百分位数的计算值。如果我们假设响应时间以毫秒为单位,则可以立即看出网页通常在 10-725 毫秒内加载,但偶尔会飙升至 945-985 毫秒。

通常,管理员只对异常值(极端百分位数)感兴趣。我们可以仅指定我们感兴趣的百分比(请求的百分比必须是 0-100 之间的值,包括 0 和 100)

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    aggregations: {
      load_time_outlier: {
        percentiles: {
          field: 'load_time',
          percents: [
            95,
            99,
            99.9
          ]
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "aggs": {
    "load_time_outlier": {
      "percentiles": {
        "field": "load_time",
        "percents": [ 95, 99, 99.9 ] 
      }
    }
  }
}

使用 percents 参数指定要计算的特定百分位数

键控响应编辑

默认情况下,keyed 标志设置为 true,它将唯一的字符串键与每个存储桶相关联,并将范围作为哈希而不是数组返回。将 keyed 标志设置为 false 将禁用此行为

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    aggregations: {
      load_time_outlier: {
        percentiles: {
          field: 'load_time',
          keyed: false
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "aggs": {
    "load_time_outlier": {
      "percentiles": {
        "field": "load_time",
        "keyed": false
      }
    }
  }
}

响应

{
  ...

  "aggregations": {
    "load_time_outlier": {
      "values": [
        {
          "key": 1.0,
          "value": 10.0
        },
        {
          "key": 5.0,
          "value": 30.0
        },
        {
          "key": 25.0,
          "value": 170.0
        },
        {
          "key": 50.0,
          "value": 445.0
        },
        {
          "key": 75.0,
          "value": 720.0
        },
        {
          "key": 95.0,
          "value": 940.0
        },
        {
          "key": 99.0,
          "value": 980.0
        }
      ]
    }
  }
}

脚本编辑

如果需要对未编入索引的值运行聚合,请使用 运行时字段。例如,如果我们的加载时间以毫秒为单位,但您希望以秒为单位计算百分位数

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    runtime_mappings: {
      'load_time.seconds' => {
        type: 'long',
        script: {
          source: "emit(doc['load_time'].value / params.timeUnit)",
          params: {
            "timeUnit": 1000
          }
        }
      }
    },
    aggregations: {
      load_time_outlier: {
        percentiles: {
          field: 'load_time.seconds'
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "runtime_mappings": {
    "load_time.seconds": {
      "type": "long",
      "script": {
        "source": "emit(doc['load_time'].value / params.timeUnit)",
        "params": {
          "timeUnit": 1000
        }
      }
    }
  },
  "aggs": {
    "load_time_outlier": {
      "percentiles": {
        "field": "load_time.seconds"
      }
    }
  }
}

百分位数(通常)是近似值编辑

有许多不同的算法可以计算百分位数。朴素的实现只是将所有值存储在一个排序数组中。要找到第 50 个百分位数,您只需找到位于 my_array[count(my_array) * 0.5] 处的值。

显然,朴素的实现无法扩展——排序数组的大小与数据集中值的數量成线性增长。为了计算 Elasticsearch 集群中可能数十亿个值的百分位数,将计算*近似*百分位数。

percentile 指标使用的算法称为 TDigest(由 Ted Dunning 在 使用 T-Digest 计算精确分位数 中引入)。

使用此指标时,请记住以下一些准则

  • 准确度与 q(1-q) 成正比。这意味着极端百分位数(例如 99%)比不太极端的百分位数(例如中位数)更准确
  • 对于少量值,百分位数非常准确(如果数据足够小,则可能 100% 准确)。
  • 随着存储桶中值的數量增加,算法开始近似百分位数。它有效地以牺牲准确性来节省内存。确切的不准确程度很难概括,因为它取决于您的数据分布和正在聚合的数据量

下图显示了均匀分布上的相对误差,具体取决于收集的值的数量和请求的百分位数

percentiles error

它显示了极端百分位数的精度如何更高。误差随着值的数量增加而减少的原因是大数定律使值的分布越来越均匀,而 t-digest 树可以更好地对其进行汇总。在更偏斜的分布上,情况并非如此。

百分位数聚合也是 非确定性的。这意味着您可以使用相同的数据获得略有不同的结果。

压缩编辑

近似算法必须在内存使用量和估计准确性之间取得平衡。可以使用 compression 参数控制此平衡

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    aggregations: {
      load_time_outlier: {
        percentiles: {
          field: 'load_time',
          tdigest: {
            compression: 200
          }
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "aggs": {
    "load_time_outlier": {
      "percentiles": {
        "field": "load_time",
        "tdigest": {
          "compression": 200    
        }
      }
    }
  }
}

压缩控制内存使用量和近似误差

TDigest 算法使用多个“节点”来近似百分位数——可用的节点越多,准确性越高(以及更大的内存占用),与数据量成正比。compression 参数将最大节点数限制为 20 * compression

因此,通过增加压缩值,您可以提高百分位数的准确性,但会占用更多内存。较大的压缩值也会降低算法速度,因为底层树数据结构的大小会增加,从而导致操作成本更高。默认压缩值为 100

一个“节点”大约使用 32 字节的内存,因此在最坏的情况下(大量数据按顺序到达并按顺序排列),默认设置将产生大约 64KB 大小的 TDigest。在实践中,数据往往更加随机,TDigest 将使用更少的内存。

执行提示编辑

TDigest 的默认实现针对性能进行了优化,可以扩展到数百万甚至数十亿个样本值,同时保持可接受的准确度水平(在某些情况下,数百万个样本的相对误差接近 1%)。可以通过将参数 execution_hint 设置为值 high_accuracy 来选择使用针对准确性优化的实现

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    aggregations: {
      load_time_outlier: {
        percentiles: {
          field: 'load_time',
          tdigest: {
            execution_hint: 'high_accuracy'
          }
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "aggs": {
    "load_time_outlier": {
      "percentiles": {
        "field": "load_time",
        "tdigest": {
          "execution_hint": "high_accuracy"    
        }
      }
    }
  }
}

以牺牲性能为代价优化 TDigest 的准确性

此选项可以提高准确性(在某些情况下,数百万个样本的相对误差接近 0.01%),但百分位数查询的完成时间会延长 2 到 10 倍。

HDR 直方图编辑

HDR 直方图(高动态范围直方图)是一种替代实现,在计算延迟测量的百分位数时非常有用,因为它比 t-digest 实现更快,但代价是内存占用更大。此实现保持固定的最坏情况百分比误差(指定为有效数字位数)。这意味着,如果在设置为 3 位有效数字的直方图中记录的数据值从 1 微秒到 1 小时(3,600,000,000 微秒)不等,则对于高达 1 毫秒的值,它将保持 1 微秒的值分辨率,对于最大跟踪值(1 小时),它将保持 3.6 秒(或更好)的值分辨率。

可以通过在请求中指定 hdr 参数来使用 HDR 直方图

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    aggregations: {
      load_time_outlier: {
        percentiles: {
          field: 'load_time',
          percents: [
            95,
            99,
            99.9
          ],
          hdr: {
            number_of_significant_value_digits: 3
          }
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "aggs": {
    "load_time_outlier": {
      "percentiles": {
        "field": "load_time",
        "percents": [ 95, 99, 99.9 ],
        "hdr": {                                  
          "number_of_significant_value_digits": 3 
        }
      }
    }
  }
}

hdr 对象指示应使用 HDR 直方图来计算百分位数,并且可以在对象内部指定此算法的特定设置

number_of_significant_value_digits 指定直方图的值分辨率(以有效数字位数表示)

HDRHistogram 仅支持正值,如果传递负值,则会出错。如果值的范围未知,则也不建议使用 HDRHistogram,因为这可能会导致内存使用量过高。

缺失值编辑

missing 参数定义了如何处理缺少值的文档。默认情况下,它们将被忽略,但也可以将它们视为具有值。

response = client.search(
  index: 'latency',
  body: {
    size: 0,
    aggregations: {
      grade_percentiles: {
        percentiles: {
          field: 'grade',
          missing: 10
        }
      }
    }
  }
)
puts response
GET latency/_search
{
  "size": 0,
  "aggs": {
    "grade_percentiles": {
      "percentiles": {
        "field": "grade",
        "missing": 10       
      }
    }
  }
}

grade 字段中没有值的文档将与值为 10 的文档属于同一个存储桶。