bayesopt.py 文件源码

def _poi(x, gp, y_max, xi):
        """
        Probability of Improvement, (Kushner, 1964)
        """
        mean, var = gp.predict(x, eval_MSE=True)
        var = np.maximum(var, 1e-9 + 0 * var)   # ????????
        z = (mean - y_max - xi) / np.sqrt(var)
        return norm.cdf(z)