ECM Algorithm for Estimating Vector ARMA Model with Variance Gamma Distribution and Possible Unbounded Density

Abstract

The simultaneous analysis of several financial time series is important in portfolio setting and risk management. This paper proposes a novel alternating Expectation conditional Maximisation (AECM) algorithm to estimate the vector autoregressive moving average (VARMA) model with variance gamma (VG) error distribution in the multivariate skewed setting. We explain why VARMA-VG model is suitable for high frequency returns (HFRs) because VG distribution provides thick tails to capture the high kurtosis in the data and the unbounded central density further captures the majority of near zero HFRs. The distribution can also be expressed in normal mean-variance mixtures which facilitate model implementation using Bayesian and expectation maximisation (EM) approaches. We adopt the EM approach which avoids the time consuming Markov chain Monto Carlo sampling and solve the unbounded density problem in the classical maximum likelihood estimation. We discuss some properties of VARMA model, conduct extensive simulation studies to evaluate the accuracy of the proposed AECM estimator and apply the models to analyse the dependency between several HFR series.