In this paper we derive the optimal linear shrinkage estimator for the large-dimensional mean vector using random matrix theory. The results are obtained under the assumption that both the dimension p and the sample size n tend to infinity such that p/n \to c\in(0, + \infty). Under weak conditions imposed on the underlying data generating process , we find the asymptotic equivalents to the optimal shrinkage intensities and estimate them consistently. The obtained non-parametric estimator for the high-dimensional mean vector has a simple structure and is proven to minimize asymptotically  with probability 1  the quadratic loss in the case of c\in(0,1). For c\in(1, + \infty) we modify the suggested estimator by using a feasible estimator for the precision covariance matrix. To this end, an exhaustive simulation study and an application to real data are provided where the proposed estimator is compared with known benchmarks from the literature. It turns out that the existent estimators of the mean vector including the suggested one converge to the sample mean vector when the true mean vector possesses unbounded Euclidean norm.
In this paper we derive the optimal linear shrinkage estimator for the high-dimensional mean vector using random matrix theory. The results are obtained under the assumption that both the dimension p and the sample size n tend to infinity in such a way that p/n \to c\in(0, \infty). Under weak conditions imposed on the underlying data generating mechanism , we find the asymptotic equivalents to the optimal shrinkage intensities and estimate them consistently. The proposed nonparametric estimator for the high-dimensional mean vector has a simple structure and is proven to minimize asymptotically , with probability 1 , the quadratic loss when c\in(0,1). When c\in(1, \infty) we modify the estimator by using a feasible estimator for the precision covariance matrix. To this end, an exhaustive simulation study and an application to real data are provided where the proposed estimator is compared with known benchmarks from the literature. It turns out that the existing estimators of the mean vector , including the new proposal, converge to the sample mean vector when the true mean vector has an unbounded Euclidean norm.
R
large-dimensional
high-dimensional
R
such
in such a way
D
+
R
process
mechanism
R
obtained non-parametric
proposed nonparametric
A
,
A
,
R
in the case of
when
R
For
When
D
+
D
suggested
R
existent
existing
R
including the suggested one
, including the new proposal,
R
possesses
has an