Simulation-based estimation of quantiles

This paper discusses the implementation of a sequential quantile-estimation algorithm for highly correlated steady-state simulation output. The primary focus is on issues related to computational and storage requirements of order statistics. The algorithm can compute exact sample quantiles and process sample sizes up to several billion without storing and sorting the whole sequence. The algorithm dynamically increases the sample size so that the quantile estimated satisfies a pre-specified precision requirement.