Moving average stochastic volatility models with application to inflation forecast

We introduce a new class of models that has both stochastic volatility and moving average errors, where the conditional mean has a state space representation. Having a moving average component, however, means that the errors in the measurement equation are no longer serially independent, and estimation becomes more difficult. We develop a posterior simulator that builds upon recent advances in precision-based algorithms for estimating these new models. In an empirical application involving US inflation we find that these moving average stochastic volatility models provide better in-sample fitness and out-of-sample forecast performance than the standard variants with only stochastic volatility.     

Econometric estimation in long-range dependent volatility models: Theory and practice

It is commonly accepted that some financial data may exhibit long-range dependence, while other financial data exhibit intermediate-range dependence or short-range dependence. These behaviours may be fitted to a continuous-time fractional stochastic model. The estimation procedure proposed in this paper is based on a continuous-time version of the Gauss–Whittle objective function to find the parameter estimates that minimize the discrepancy between the spectral density and the data periodogram. As a special case, the proposed estimation procedure is applied to a class of fractional stochastic volatility models to estimate the drift, standard deviation and memory parameters of the volatility process under consideration. As an application, the volatility of the Dow Jones, S&P 500, CAC 40, DAX 30, FTSE 100 and NIKKEI 225 is estimated.     

Estimation of affine term structure models with spanned or unspanned stochastic volatility

We develop new procedures for maximum likelihood estimation of affine term structure models with spanned or unspanned stochastic volatility. Our approach uses linear regression to reduce the dimension of the numerical optimization problem yet it produces the same estimator as maximizing the likelihood. It improves the numerical behavior of estimation by eliminating parameters from the objective function that cause problems for conventional methods. We find that spanned models capture the cross-section of yields well but not volatility while unspanned models fit volatility at the expense of fitting the cross-section.     

Option Pricing under the Double Exponential Jump-Diffusion Model with Stochastic Volatility and Interest Rate

This paper proposes an efficient option pricing model that incorporates stochastic interest rate (SIR), stochastic volatility (SV), and double exponential jump into the jump-diffusion settings. The model comprehensively considers the leptokurtosis and heteroscedasticity of the underlying asset’s returns, rare events, and an SIR. Using the model, we deduce the pricing characteristic function and pricing formula of a European option. Then, we develop the Markov chain Monte Carlo method with latent variable to solve the problem of parameter estimation under the double exponential jump-diffusion model with SIR and SV. For verification purposes, we conduct time efficiency analysis, goodness of fit analysis, and jump/drift term analysis of the proposed model. In addition, we compare the pricing accuracy of the proposed model with those of the Black–Scholes and the Kou (2002) models. The empirical results show that the proposed option pricing model has high time efficiency, and the goodness of fit and pricing accuracy are significantly higher than those of the other two models.     

Estimation of objective and risk-neutral distributions based on moments of integrated volatility

In this paper, we present an estimation procedure which uses both option prices and high-frequency spot price feeds to estimate jointly the objective and risk-neutral parameters of stochastic volatility models. The procedure is based on a method of moments that uses analytical expressions for the moments of the integrated volatility and series expansions of option prices and implied volatilities. This results in an easily implementable and rapid estimation technique. An extensive Monte Carlo study compares various procedures and shows the efficiency of our approach. Empirical applications to the Deutsche mark–US dollar exchange rate futures and the S&P 500 index provide evidence that the method delivers results that are in line with the ones obtained in previous studies where much more involved estimation procedures were used.     

Econometric analysis of jump-driven stochastic volatility models

This paper introduces and studies the econometric properties of a general new class of models, which I refer to as jump-driven stochastic volatility models, in which the volatility is a moving average of past jumps. I focus attention on two particular semiparametric classes of jump-driven stochastic volatility models. In the first, the price has a continuous component with time-varying volatility and time-homogeneous jumps. The second jump-driven stochastic volatility model analyzed here has only jumps in the price, which have time-varying size. In the empirical application I model the memory of the stochastic variance with a CARMA(2,1) kernel and set the jumps in the variance to be proportional to the squared price jumps. The estimation, which is based on matching moments of certain realized power variation statistics calculated from high-frequency foreign exchange data, shows that the jump-driven stochastic volatility model containing continuous component in the price performs best. It outperforms a standard two-factor affine jump–diffusion model, but also the pure-jump jump-driven stochastic volatility model for the particular jump specification.     

Analysis of high dimensional multivariate stochastic volatility models

This paper is concerned with the Bayesian estimation and comparison of flexible, high dimensional multivariate time series models with time varying correlations. The model proposed and considered here combines features of the classical factor model with that of the heavy tailed univariate stochastic volatility model. A unified analysis of the model, and its special cases, is developed that encompasses estimation, filtering and model choice. The centerpieces of the estimation algorithm (which relies on MCMC methods) are: (1) a reduced blocking scheme for sampling the free elements of the loading matrix and the factors and (2) a special method for sampling the parameters of the univariate SV process. The resulting algorithm is scalable in terms of series and factors and simulation-efficient. Methods for estimating the log-likelihood function and the filtered values of the time-varying volatilities and correlations are also provided. The performance and effectiveness of the inferential methods are extensively tested using simulated data where models up to 50 dimensions and 688 parameters are fit and studied. The performance of our model, in relation to various multivariate GARCH models, is also evaluated using a real data set of weekly returns on a set of 10 international stock indices. We consider the performance along two dimensions: the ability to correctly estimate the conditional covariance matrix of future returns and the unconditional and conditional coverage of the 5% and 1% value-at-risk (VaR) measures of four pre-defined portfolios.     

Dynamic correlation multivariate stochastic volatility with latent factors

Modeling the correlation structure of returns is essential in many financial applications. Considerable evidence from empirical studies has shown that the correlation among asset returns is not stable over time. A recent development in the multivariate stochastic volatility literature is the application of inverse Wishart processes to characterize the evolution of return correlation matrices. Within the inverse Wishart multivariate stochastic volatility framework, we propose a flexible correlated latent factor model to achieve dimension reduction and capture the stylized fact of ‘correlation breakdown’ simultaneously. The parameter estimation is based on existing Markov chain Monte Carlo methods. We illustrate the proposed model with several empirical studies. In particular, we use high‐dimensional stock return data to compare our model with competing models based on multiple performance metrics and tests. The results show that the proposed model not only describes historic stylized facts reasonably but also provides the best overall performance.     

A discrete-time model for daily S & P500 returns and realized variations: Jumps and leverage effects

We develop an empirically highly accurate discrete-time daily stochastic volatility model that explicitly distinguishes between the jump and continuous-time components of price movements using nonparametric realized variation and Bipower variation measures constructed from high-frequency intraday data. The model setup allows us to directly assess the structural inter-dependencies among the shocks to returns and the two different volatility components. The model estimates suggest that the leverage effect, or asymmetry between returns and volatility, works primarily through the continuous volatility component. The excellent fit of the model makes it an ideal candidate for an easy-to-implement auxiliary model in the context of indirect estimation of empirically more realistic continuous-time jump diffusion and Lévy-driven stochastic volatility models, effectively incorporating the interdaily dependencies inherent in the high-frequency intraday data.     

Efficient estimation of Bayesian VARMAs with time‐varying coefficients

Empirical work in macroeconometrics has been mostly restricted to using vector autoregressions (VARs), even though there are strong theoretical reasons to consider general vector autoregressive moving averages (VARMAs). A number of articles in the last two decades have conjectured that this is because estimation of VARMAs is perceived to be challenging and proposed various ways to simplify it. Nevertheless, VARMAs continue to be largely dominated by VARs, particularly in terms of developing useful extensions. We address these computational challenges with a Bayesian approach. Specifically, we develop a Gibbs sampler for the basic VARMA, and demonstrate how it can be extended to models with time‐varying vector moving average (VMA) coefficients and stochastic volatility. We illustrate the methodology through a macroeconomic forecasting exercise. We show that in a class of models with stochastic volatility, VARMAs produce better density forecasts than VARs, particularly for short forecast horizons.     

Bayesian semiparametric stochastic volatility modeling

This paper extends the existing fully parametric Bayesian literature on stochastic volatility to allow for more general return distributions. Instead of specifying a particular distribution for the return innovation, nonparametric Bayesian methods are used to flexibly model the skewness and kurtosis of the distribution while the dynamics of volatility continue to be modeled with a parametric structure. Our semiparametric Bayesian approach provides a full characterization of parametric and distributional uncertainty. A Markov chain Monte Carlo sampling approach to estimation is presented with theoretical and computational issues for simulation from the posterior predictive distributions. An empirical example compares the new model to standard parametric stochastic volatility models.     

Modeling spot rate using a realized stochastic volatility model with level effect and dynamic drift☆

This paper proposes a class of realized stochastic volatility model based on both various realized volatility measures and spot rate. It applies the realized stochastic volatility model (Takahashi, Omori, & Watanabe, 2009, and Koopman & Scharth, 2013) to the spot rate model with dynamic drift and level effect setups (RSVL). A jointly approximated maximum likelihood procedure is used to estimate this model. The simulation results show that the RSVL model can be consistently estimated and noise-and-jump-robust realized volatility measures improve the accuracy of the estimation. This study empirically investigates the Chinese interbank repo market with RSVL model, which manifested the advantage of taking the level effect and nonlinear drift into consideration. The noise-and-jump-robust realized volatility measures (e.g. subsample realized volatility and threshold pre-average realized volatility) decrease the volatility fitting error. The nonparametric testing suggests that the RSVL model with noise-and-jump-robust realized volatility measures has more power on forecasting excess kurtosis and fat tails and predicting dynamics of higher order autocorrelations.     

Indirect inference with time series observed with error

We propose the indirect inference estimator as a consistent method to estimate the parameters of a structural model when the observed series are contaminated by measurement error by considering the noise as a structural feature. We show that the indirect inference estimates are asymptotically biased if the error is neglected. When the condition for identification is satisfied, the structural and measurement error parameters can be consistently estimated. The issues of identification and misspecification of measurement error are discussed in detail. We illustrate the reliability of this procedure in the estimation of stochastic volatility models based on realized volatility measures contaminated by microstructure noise.     

Kalman filtering estimation of unobserved rational expectations with an application to the German hyperinflation

The assumption that rational expectations always lie on a convergent path is subject to an empirical test using the German hyperinflation data. The estimation technique employs a Kalman filtering algorithm. After presenting a brief background for the convergent expectations problem and a derivation of the various model specifications, a generalized expectations model and its attendant Kalman filtering estimation technique are discussed. Additional estimation details and empirical results are then presented. Based on an assumption of normally distributed errors, the null hypothesis of convergent paths is rejected in all situations involving a deterministic specification of the evolution of the unobserved parameter which characterizes the convergent path. The same null hypothesis is rejected in four of the six cases corresponding to a stochastic specification of the evolution of the unobserved parameter which characterizes the convergent path. A discussion of these findings, their economic significance, and suggestions for further research concludes the paper.     

Approximate analytic solution for Asian options with stochastic volatility

The valuation of Asian options is complicated because the arithmetic average of lognormal random variables is no longer lognormal. Furthermore, the stochastic volatility inherent in financial asset prices is easily observed. However, few academic studies consider the pricing and hedging of Asian options with stochastic volatility, despite the popularity of such options. This study extends the work of Hull and White (1987) and integrates the Taylor series expansion technique to derive an approximate analytic solution for Asian options with stochastic volatility. Numerical experiments show that the proposed approximate analytic solution performs favorably and is computationally efficient compared with large-sample simulations. The approximate analytic solution provides a practical approach for pricing and hedging Asian options with stochastic volatility and is both easy to implement and desirable in terms of computing speed.     

Dynamic stochastic copula models: estimation,inference and applications

We propose a new dynamic copula model in which the parameter characterizing dependence follows an autoregressive process. As this model class includes the Gaussian copula with stochastic correlation process, it can be viewed as a generalization of multivariate stochastic volatility models. Despite the complexity of the model, the decoupling of marginals and dependence parameters facilitates estimation. We propose estimation in two steps, where first the parameters of the marginal distributions are estimated, and then those of the copula. Parameters of the latent processes (volatilities and dependence) are estimated using efficient importance sampling. We discuss goodness‐of‐fit tests and ways to forecast the dependence parameter. For two bivariate stock index series, we show that the proposed model outperforms standard competing models. Copyright © 2010 John Wiley & Sons, Ltd.     

The estimation of a state space model by estimating functions with an application

Conventionally the parameters of a linear state space model are estimated by maximizing a Gaussian likelihood function, even when the input errors are not Gaussian. In this paper we propose estimation by estimating functions fulfilling Godambe's optimality criterion. We discuss the issue of an unknown starting state vector, and we also develop recursive relations for the third- and fourth-order moments of the state predictors required for the calculations. We conclude with a simulation study demonstrating the proposed procedure on the estimation of the stochastic volatility model. The results suggest that the new estimators outperform the Gaussian likelihood.     

Forecasting volatility: A reality check based on option pricing, utility function, value-at-risk, and predictive likelihood

We analyze the predictive performance of various volatility models for stock returns. To compare their performance, we choose loss functions for which volatility estimation is of paramount importance. We deal with two economic loss functions (an option pricing function and an utility function) and two statistical loss functions (a goodness-of-fit measure for a value-at-risk (VaR) calculation and a predictive likelihood function). We implement the tests for superior predictive ability of White [Econometrica 68 (5) (2000) 1097] and Hansen [Hansen, P. R. (2001). An unbiased and powerful test for superior predictive ability. Brown University]. We find that, for option pricing, simple models like the Riskmetrics exponentially weighted moving average (EWMA) or a simple moving average, which do not require estimation, perform as well as other more sophisticated specifications. For a utility-based loss function, an asymmetric quadratic GARCH seems to dominate, and this result is robust to different degrees of risk aversion. For a VaR-based loss function, a stochastic volatility model is preferred. Interestingly, the Riskmetrics EWMA model, proposed to calculate VaR, seems to be the worst performer. For the predictive likelihood-based loss function, modeling the conditional standard deviation instead of the variance seems to be a dominant modeling strategy.     

Nowcasting GDP with a pool of factor models and a fast estimation algorithm

We propose a novel mixed-frequency dynamic factor model with time-varying parameters and stochastic volatility for macroeconomic nowcasting and develop a fast estimation algorithm. This enables us to generate forecast densities based on a large space of factor models. We apply our framework to nowcast US GDP growth in real time. Our results reveal that stochastic volatility seems to improve the accuracy of point forecasts the most, compared to the constant-parameter factor model. These gains are most prominent during unstable periods such as the Covid-19 pandemic. Finally, we highlight indicators driving the US GDP growth forecasts and associated downside risks in real time.     

An empirical application of stochastic volatility models

This paper studies the empirical performance of stochastic volatility models for twenty years of weekly exchange rate data for four major currencies. We concentrate on the effects of the distribution of the exchange rate innovations for both parameter estimates and for estimates of the latent volatility series. The density of the log of squared exchange rate innovations is modelled as a flexible mixture of normals. We use three different estimation techniques: quasi-maximum likelihood, simulated EM, and a Bayesian procedure. The estimated models are applied for pricing currency options. The major findings of the paper are that: (1) explicitly incorporating fat-tailed innovations increases the estimates of the persistence of volatility dynamics; (2) the estimation error of the volatility time series is very large; (3) this in turn causes standard errors on calculated option prices to be so large that these prices are rarely significantly different from a model with constant volatility. © 1998 John Wiley & Sons, Ltd.