A note on the Estimation of a Gamma-Variance Process: Learning from a Failure

This paper confirms that, as originally reported in Seneta (Journal of Applied Probability 41:177–187, 2004, p. 183), it is impossible to replicate Madan et al. (European Finance Review 2:135–156, 1998) results using log daily returns on S&P 500 Index from January 1992 to September 1994. This failure leads to a close investigation of the computational problems associated with finding maximum likelihood estimates of the parameters of the popular VG model. Both standard econometric software, such as R, low level programming languages, such as Matlab®, and non-standard optimization software, such as Ezgrad described in Tucci (Journal of Economic Dynamics and Control 26:1739–1764, 2002), are used. The complexity of the log-likelihood function is studied. It is shown that it looks very complicated, with many local optima, and may be incredibly sensitive to very small changes in the sample used. Adding or removing a single observation may cause huge changes both in the maximum of the log-likelihood function and in the estimated parameter values. An intuitive procedure which works nicely both when implemented in R and in Matlab® is presented