Bayesian Quadrature Optimization for Probability Threshold Robustness Measure