Numerical implementation of the QuEST function

Certain estimation problems involving the covariance matrix in large dimensions are considered. Due to the breakdown of finite-dimensional asymptotic theory when the dimension is not negligible with respect to the sample size, it is necessary to resort to an alternative framework known as large-dimensional asymptotics. Recently, an estimator of the eigenvalues of the population covariance matrix has been proposed that is consistent according to a mean-squared criterion under large-dimensional asymptotics. It requires numerical inversion of a multivariate nonrandom function called the QuEST function. The numerical implementation of this QuEST function in practice is explained through a series of six successive steps. An algorithm is provided in order to compute the Jacobian of the QuEST function analytically, which is necessary for numerical inversion via a nonlinear optimizer. Monte Carlo simulations document the effectiveness of the code.