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免费学习笔记(深入)”;

以下函数 generate_sensitivity_configs 实现该逻辑:

会译·对照式翻译
会译·对照式翻译

会译是一款AI智能翻译浏览器插件,支持多语种对照式翻译

下载
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

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

豆包大模型
豆包大模型

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

WorkBuddy
WorkBuddy

腾讯云推出的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字符串转数组的相关的文章、下载、课程内容,供大家免费下载体验。

761

2023.08.03

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

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

221

2023.09.04

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

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

1570

2023.10.24

字符串介绍
字符串介绍

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

651

2023.11.24

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

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

1269

2024.03.22

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

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

1206

2024.04.29

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

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

194

2025.07.29

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

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

131

2025.08.07

c++ 字符处理
c++ 字符处理

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

0

2026.03.17

热门下载

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

精品课程

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

共4课时 | 22.5万人学习

Django 教程
Django 教程

共28课时 | 5.1万人学习

SciPy 教程
SciPy 教程

共10课时 | 2万人学习

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

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