scipy.stats.qmc.

update_discrepancy#

scipy.stats.qmc.update_discrepancy(x_new, sample, initial_disc)[源代码]#

用新的样本更新中心偏差。

参数:

x_newarray_like (1, d): 要添加到 sample 中的新样本。
samplearray_like (n, d): 初始样本。
initial_discfloat: sample 的中心偏差。

返回:

discrepancyfloat: 由 x_new 和 sample 组成的样本的中心偏差。

示例

我们也可以通过使用 iterative=True 迭代地计算偏差。

>>> import numpy as np
>>> from scipy.stats import qmc
>>> space = np.array([[1, 3], [2, 6], [3, 2], [4, 5], [5, 1], [6, 4]])
>>> l_bounds = [0.5, 0.5]
>>> u_bounds = [6.5, 6.5]
>>> space = qmc.scale(space, l_bounds, u_bounds, reverse=True)
>>> disc_init = qmc.discrepancy(space[:-1], iterative=True)
>>> disc_init
0.04769081147119336
>>> qmc.update_discrepancy(space[-1], space[:-1], disc_init)
0.008142039609053513