怎么在Python中实现带权重的无偏随机采样结构

落枫大大_2183

落枫大大_2183

2026-10-19

220人浏览

原创

random.choices()是python标准库中用于带权重、可放回、无偏采样的最直接方法,支持非归一化非负权重,k参数控制采样数量,权重为0表示永不被选中。

怎么在python中实现带权重的无偏随机采样结构

用 random.choices() 做带权重的无偏采样最直接

Python 标准库的 random.choices() 就是为此设计的:它支持重复采样(可放回),按权重比例决定每个元素被选中的概率,且每次采样相互独立——这正是“无偏”的核心:每个样本的生成不依赖历史结果,权重只影响单次选择的概率分布。

注意它和 random.sample() 的根本区别:sample() 是等概率、不可重复的;而 choices() 是加权、可重复的,且不改变原序列结构。

  • 权重可以是任意非负数,不要求归一化(内部自动处理)
  • 若权重全为 0,会抛出 ValueError: sequence too short
  • 采样数量由 k 参数控制,省略时默认返回单个元素(非列表)
import random
data = ['a', 'b', 'c']
weights = [1, 2, 10]  # 'c' 被选中概率约 10/13 ≈ 77%
sample = random.choices(data, weights=weights, k=5)
# 示例输出:['c', 'c', 'b', 'c', 'a']

权重为浮点数或零值时要注意精度与边界

浮点权重本身没问题,但极小值(如 1e-15)可能因浮点累积误差导致归一化后总和偏离 1.0,不过 random.choices() 内部用的是前缀和 + 二分查找,对精度鲁棒性较强。真正要警惕的是含零权重的情况:

  • 权重为 0 表示该元素**永不被选中**,但不会报错(只要至少有一个正权重)
  • 如果所有权重都是 0 或负数,会触发 ValueError: weights sum to zero
  • 空序列或空权重列表会抛出 IndexError: Cannot choose from empty sequence

实际使用中,建议显式过滤掉权重 ≤ 0 的项,避免隐式跳过引发逻辑偏差:

valid_pairs = [(x, w) for x, w in zip(data, weights) if w > 0]
if not valid_pairs:
    raise ValueError("No positive weights available")
filtered_data, filtered_weights = zip(*valid_pairs)
result = random.choices(filtered_data, weights=filtered_weights, k=10)

需要不放回的加权采样?标准库不支持,得自己实现

random.choices() 天然是可放回的;若你明确要求“不放回”+“按权重比例初始概率”,那就不是标准意义上的“无偏随机采样”,而是带约束的加权抽样(例如类似轮盘赌但逐次移除)。Python 没有内置函数,常见做法是用别名法(Alias Method)或重采样+拒绝法。

361SCAN security check your skill python scripts
361SCAN security check your skill python scripts

扫描并分析已安装技能,用于查看指定技能目录的名称、描述和详情,或列出所有已安装的技能。

下载

简单场景下可用“权重动态衰减”模拟,但严格来说已非无偏——因为后续选择依赖前面的结果。如果你真需要数学上严格的不放回加权抽样(即第 i 次选中 x 的概率 = wₓ / Σ(剩余权重)),推荐用 numpy.random.Generator.choice() 并设 replace=False:

import numpy as np
rng = np.random.default_rng()
data = np.array(['a', 'b', 'c'])
weights = np.array([1, 2, 10])
sample = rng.choice(data, size=2, replace=False, p=weights/weights.sum())

注意:这要求 weights 必须是正数且和为 1(或让 NumPy 归一化),且 size 不能超过 len(data)。

性能敏感时别在循环里反复调用 random.choices()

如果要在高频循环中做大量小规模采样(比如每轮选 1 个,跑百万次),random.choices() 每次都会重新计算权重前缀和——虽然单次很快,但累积开销可观。此时应预计算前缀和,再手动二分查找:

  • 用 itertools.accumulate(weights) 一次性算好前缀和
  • 用 bisect.bisect_left(prefix_sum, r * total) 查找索引,其中 r = random.random()
  • 这样把 O(n) 权重处理降为 O(1) 初始化 + O(log n) 每次查询

不过绝大多数业务场景完全没必要这么干——先用 choices(),profile 真发现它是瓶颈再说。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

python

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1751

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

4544

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1729

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

26797

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3247

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3267

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1183

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

596

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2523

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Uvicorn官方安装与入门手册
Uvicorn官方安装与入门手册

共0课时 | 0人学习

Python uv官方项目管理指南
Python uv官方项目管理指南

共0课时 | 0人学习

Python uv官方安装与入门文档
Python uv官方安装与入门文档

共0课时 | 0人学习