Abstract
The present paper combines nonlinear shrinkage with the Multivariate Generalized Hyperbolic (MGHyp) distribution to account for heavy tails in estimating the first and second moments in high dimensions. An Expectation-Maximization (EM) algorithm is developed that is fast, stable, and applicable in high dimensions. Theoretical arguments for the monotonicity of the proposed algorithm are provided and it is shown in simulations that it is able to accurately retrieve parameter estimates. Finally, in an extensive Markowitz portfolio optimization analysis, the approach is compared to state-of-the-art benchmark models. The proposed model excels with a strong out-of-sample portfolio performance combined with a comparably low turnover.