Nonparametric estimation of a maximum of quantiles

A simulation model of a complex system is considered, for which the outcome is described by m(p, X), where p is a parameter of the system, X is a random input of the system and m is a real-valued function. The maximum (with respect to p) of the quantiles of m(p, X) is estimated. The quantiles of m(p, X) of a given level are estimated for various values of p from an order statistic of values m(pi, Xi) where X, X1, X2, . . . are independent and identically distributed and where pi is close to p, and the maximal quantile is estimated by the maximum of these quantile estimates. Under assumptions on the smoothness of the function describing the dependency of the values of the quantiles on the parameter p the rate of convergence of this estimate is analyzed. The finite sample size behavior of the estimate is illustrated by simulated data and by applying it in a simulation model of a real mechanical system.