Nonparametric estimation of a quantile qm(X),αof a random variable m(X) is considered, where m : ℝd→ ℝ is a function, which is costly to compute and X is an ℝd-valued random variable with known distribution. Monte Carlo surrogate quantile estimates are considered, where in a first step, the function m is estimated by some estimate (surrogate) mn, and then, the quantile qm(X),αis estimated by a Monte Carlo estimate of the quantile qmn(X),α. A general error bound on the error of this quantile estimate is derived, which depends on the local error of the function estimate mn, and the rates of convergence of the corresponding Monte Carlo surrogate quantile estimates are analyzed for two different function estimates. The finite sample size behavior of the estimates is investigated in simulations.