In the two-dimensional space consisting of the failure rate and space time yield, we first When δ = 1, the entire DoS is compared. Otherwise, the distance is computed only with respect to a low failure rate region. Figure shows that EPA achieves better accuracy than random sampling in low failure rate regions, while random sampling is better in higher failure rate regions. EPA is preferred for discovering best configurations of a chemical process, because samples at low failure regions are more relevant than the others. Note that as increasing the number of samples, i.e., the number of MCMC steps, Hellinger distances for all thresholds become smaller (see Figure ) in EPA. For E-factor, material cost, and yield, same analysis is also performed, and the results are shown in Figures S2, S3 To understand high quality samples whose failure rate is below 0.01, they are classified into four clusters by applying the weighted k-means method to the nine-dimensional outcome space. Here, five-types of failure rates for total, compounds 1, 2, 4, and 5 are included.