Combining the MGHyp distribution with nonlinear shrinkage in modeling financial asset returns

The present paper combines nonlinear shrinkage with the multivariate generalized hyperbolic (MGHyp) distribution, thereby extending a flexible parametric model to 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.