A quantile-based method to efficiently determine model response at boundary probabilities