0

0

如何用 Python 迭代生成单变量敏感性分析配置(P10/P90)

碧海醫心

碧海醫心

发布时间:2026-01-08 19:54:08

|

230人浏览过

|

来源于php中文网

原创

如何用 Python 迭代生成单变量敏感性分析配置(P10/P90)

本文介绍如何利用 `itertools.product` 和动态字段替换,高效生成单变量扰动的敏感性分析配置列表,支持任意数量输入参数,避免硬编码,适用于模型输入的 p10/p90 单因素敏感性分析。

在进行模型敏感性分析时,一个常见且严谨的做法是:每次仅改变一个输入变量(如设为 P10 或 P90),其余变量保持基准值(P50)不变。这种“单变量扰动”策略能清晰归因输出变化,但手动构造配置极易出错,尤其当输入维度升高(如 10+ 个参数)时。

Python 提供了优雅的解决方案——结合 itertools.product 与命名元组(namedtuple)的字段反射能力,可全自动、可扩展地生成全部所需配置。

首先,定义你的敏感性用例结构(如题中所示):

from collections import namedtuple
from itertools import product

Sensitivity_Case = namedtuple('Sensitivity_Case', ['a', 'b', 'c'])

p50 = Sensitivity_Case(5, 50, 'medium')
p10 = Sensitivity_Case(1, 5, 'low')
p90 = Sensitivity_Case(10, 500, 'high')

关键在于:我们不预先写死 'a.p10' 这类字符串,而是将每个字段的三种取值(P10/P50/P90)组织为独立列表,并明确指定“仅允许一个字段偏离 P50”。

立即学习Python免费学习笔记(深入)”;

OneAI
OneAI

将生成式AI技术打包为API,整合到企业产品和服务中

下载

以下函数 generate_sensitivity_configs 实现该逻辑:

def generate_sensitivity_configs(base_case, cases: list):
    """
    生成单变量敏感性配置列表。

    Args:
        base_case: 基准用例(如 p50),用于提供字段名和默认值
        cases: 包含多个 Sensitivity_Case 实例的列表,顺序应为 [p10, p50, p90]

    Returns:
        list of tuples: 每个 tuple 是一个完整配置(对应模型一次输入)
    """
    field_names = base_case._fields
    n_fields = len(field_names)

    # 对每个字段,提取其在所有 cases 中的取值(如 a: [p10.a, p50.a, p90.a])
    field_values = [
        [case._asdict()[field] for case in cases]
        for field in field_names
    ]

    configs = []
    # 遍历每个字段索引(即哪个字段被扰动)
    for i in range(n_fields):
        # 当前扰动字段:取 P10 和 P90(跳过 P50)
        for perturb_val in [field_values[i][0], field_values[i][2]]:  # p10, p90
            # 其余字段固定为 P50(即 cases[1])
            config_vals = []
            for j, field in enumerate(field_names):
                if j == i:
                    config_vals.append(perturb_val)
                else:
                    config_vals.append(cases[1]._asdict()[field])
            configs.append(tuple(config_vals))

    return configs

# 使用示例
configs = generate_sensitivity_configs(p50, [p10, p50, p90])
for i, cfg in enumerate(configs, 1):
    print(f"config{i} = {cfg}")

输出结果完全匹配需求:

config1 = (1, 50, 'medium')     # a.p10, b.p50, c.p50
config2 = (10, 50, 'medium')    # a.p90, b.p50, c.p50
config3 = (5, 5, 'medium')      # a.p50, b.p10, c.p50
config4 = (5, 500, 'medium')    # a.p50, b.p90, c.p50
config5 = (5, 50, 'low')        # a.p50, b.p50, c.p10
config6 = (5, 50, 'high')       # a.p50, b.p50, c.p90

优势说明

  • 自动适配字段数:无论 Sensitivity_Case 含 3 个还是 100 个字段,逻辑不变;
  • 清晰分离关注点:cases 列表定义取值集(P10/P50/P90),函数负责组合逻辑;
  • 类型安全 & 可读性强:输出为 tuple,可直接解包传入模型函数(如 model(*cfg));
  • 无冗余配置:不生成全 P50 基准配置(若需,可单独添加 configs.append(tuple(p50)))。

⚠️ 注意事项

  • 确保 cases 参数中 p50 严格位于索引 1(即 [p10, p50, p90]),否则 cases[1] 将取错基准值;
  • 若未来需支持更多分位点(如 P5/P25/P75/P95),只需扩展 cases 列表并调整内层循环逻辑;
  • 如需返回 namedtuple 而非 tuple,可在最后用 Sensitivity_Case(*cfg) 封装。

通过此方法,你彻底告别手动枚举,以不到 20 行核心代码,构建出健壮、可维护、可扩展的敏感性配置生成器——这才是 Python “优雅解决实际问题”的典型范式。

热门AI工具

更多
DeepSeek
DeepSeek

幻方量化公司旗下的开源大模型平台

豆包大模型
豆包大模型

字节跳动自主研发的一系列大型语言模型

通义千问
通义千问

阿里巴巴推出的全能AI助手

腾讯元宝
腾讯元宝

腾讯混元平台推出的AI助手

文心一言
文心一言

文心一言是百度开发的AI聊天机器人,通过对话可以生成各种形式的内容。

讯飞写作
讯飞写作

基于讯飞星火大模型的AI写作工具,可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

即梦AI
即梦AI

一站式AI创作平台,免费AI图片和视频生成。

ChatGPT
ChatGPT

最最强大的AI聊天机器人程序,ChatGPT不单是聊天机器人,还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

相关专题

更多
js 字符串转数组
js 字符串转数组

js字符串转数组的方法:1、使用“split()”方法;2、使用“Array.from()”方法;3、使用for循环遍历;4、使用“Array.split()”方法。本专题为大家提供js字符串转数组的相关的文章、下载、课程内容,供大家免费下载体验。

760

2023.08.03

js截取字符串的方法
js截取字符串的方法

js截取字符串的方法有substring()方法、substr()方法、slice()方法、split()方法和slice()方法。本专题为大家提供字符串相关的文章、下载、课程内容,供大家免费下载体验。

220

2023.09.04

java基础知识汇总
java基础知识汇总

java基础知识有Java的历史和特点、Java的开发环境、Java的基本数据类型、变量和常量、运算符和表达式、控制语句、数组和字符串等等知识点。想要知道更多关于java基础知识的朋友,请阅读本专题下面的的有关文章,欢迎大家来php中文网学习。

1564

2023.10.24

字符串介绍
字符串介绍

字符串是一种数据类型,它可以是任何文本,包括字母、数字、符号等。字符串可以由不同的字符组成,例如空格、标点符号、数字等。在编程中,字符串通常用引号括起来,如单引号、双引号或反引号。想了解更多字符串的相关内容,可以阅读本专题下面的文章。

649

2023.11.24

java读取文件转成字符串的方法
java读取文件转成字符串的方法

Java8引入了新的文件I/O API,使用java.nio.file.Files类读取文件内容更加方便。对于较旧版本的Java,可以使用java.io.FileReader和java.io.BufferedReader来读取文件。在这些方法中,你需要将文件路径替换为你的实际文件路径,并且可能需要处理可能的IOException异常。想了解更多java的相关内容,可以阅读本专题下面的文章。

1208

2024.03.22

php中定义字符串的方式
php中定义字符串的方式

php中定义字符串的方式:单引号;双引号;heredoc语法等等。想了解更多字符串的相关内容,可以阅读本专题下面的文章。

1184

2024.04.29

go语言字符串相关教程
go语言字符串相关教程

本专题整合了go语言字符串相关教程,阅读专题下面的文章了解更多详细内容。

192

2025.07.29

c++字符串相关教程
c++字符串相关教程

本专题整合了c++字符串相关教程,阅读专题下面的文章了解更多详细内容。

111

2025.08.07

Go高并发任务调度与Goroutine池化实践
Go高并发任务调度与Goroutine池化实践

本专题围绕 Go 语言在高并发任务处理场景中的实践展开,系统讲解 Goroutine 调度模型、Channel 通信机制以及并发控制策略。内容包括任务队列设计、Goroutine 池化管理、资源限制控制以及并发任务的性能优化方法。通过实际案例演示,帮助开发者构建稳定高效的 Go 并发任务处理系统,提高系统在高负载环境下的处理能力与稳定性。

4

2026.03.10

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
最新Python教程 从入门到精通
最新Python教程 从入门到精通

共4课时 | 22.5万人学习

Django 教程
Django 教程

共28课时 | 4.9万人学习

SciPy 教程
SciPy 教程

共10课时 | 1.9万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号