MORRIS 敏感性分析

听风阁主

已于 2022-08-04 22:04:56 修改

阅读量3.4k

点赞数 1

分类专栏：编码日常数据分析_py 文章标签：大数据

于 2022-07-09 23:50:22 首次发布

本文链接：https://blog.csdn.net/weixin_43388193/article/details/125700734

版权

数据分析_py 同时被 2 个专栏收录

8 篇文章 0 订阅

订阅专栏

编码日常

1 篇文章 1 订阅

订阅专栏

from SALib.analyze import morris
from SALib.sample.morris import sample
from SALib.test_functions import Sobol_G
from SALib.util import read_param_file
from SALib.plotting.morris import horizontal_bar_plot, covariance_plot, \
    sample_histograms
import matplotlib.pyplot as plt


problem = {
 'num_vars': 8,
 'names': ['x1', 'x2', 'x3', 'x4', 'x5', 'x6', 'x7', 'x8'],
 'groups': None,
 'bounds': [[0.0, 1.0],
            [0.0, 1.0],
            [0.0, 1.0],
            [0.0, 1.0],
            [0.0, 1.0],
            [0.0, 1.0],
            [0.0, 1.0],
            [0.0, 1.0]]
}

# 生成样本
param_values = sample(problem, N=1000, num_levels=4,
                      optimal_trajectories=None)

# 生成结果，可以自己套一个模型
Y = Sobol_G.evaluate(param_values)

# 计算morris
Si = morris.analyze(problem, param_values, Y, conf_level=0.95,
                    print_to_console=True,
                    num_levels=4, num_resamples=100)

# 画图呗
fig, (ax1, ax2) = plt.subplots(1, 2)
horizontal_bar_plot(ax1, Si, {}, sortby='mu_star', unit=r"tCO$_2$/year")
covariance_plot(ax2, Si, {}, unit=r"tCO$_2$/year")

fig2 = plt.figure()
sample_histograms(fig2, param_values, problem, {'color': 'y'})
plt.show()