A Robust Test for Weak Instruments

This article deals with the estimation of continuous-time stochastic volatility models of option pricing. We argue that option prices are much more informative about the parameters than are asset prices. This is confirmed in a Monte Carlo experiment that compares two very simple strategies based on the different information sets. Both approaches are based on indirect inference and avoid any discretization bias by simulating the continuous-time model. We assume an Ornstein-Uhlenbeck process for the log of the volatility, a zero-volatility risk premium, and no leverage effect. We do not pursue asymptotic efficiency or specification issues; rather, we stick to a framework with no overidentifying restrictions and show that, given our option-pricing model, estimation based on option prices is much more precise in samples of typical size, without increasing the computational burden.     

Volatility and GMM — Monte Carlo studies and empirical estimations

In this article we examine small sample properties of a generalized method of moments (GMM) estimation using Monte Carlo simulations. We assume that the generated time series describe the stochastic variance rate of a stock index. we use a mean reverting square-root process to simulate the dynamics of this instantaneous variance rate. The time series obtained are used to estimate the parameters of the assumed variance rate process by applying GMM. Our results are described and compared to estimates from empirical data which consist of volatility as well as daily volume data of the German stock market. One of our main findings is that estimates of the mean reverting parameter that are not significantly different from zero do not necessarily imply a rejection of the hypothesis of a mean reverting behavior of the underlying stochastic process.     

Enhancing Estimation for Interest Rate Diffusion Models With Bond Prices

We consider improving estimating parameters of diffusion processes for interest rates by incorporating information in bond prices. This is designed to improve the estimation of the drift parameters, which are known to be subject to large estimation errors. It is shown that having the bond prices together with the short rates leads to more efficient estimation of all parameters for the interest rate models. It enhances the estimation efficiency of the maximum likelihood estimation based on the interest rate dynamics alone. The combined estimation based on the bond prices and the interest rate dynamics can also provide inference to the risk premium parameter. Simulation experiments were conducted to confirm the theoretical properties of the estimators concerned. We analyze the overnight Fed fund rates together with the U.S. Treasury bond prices. Supplementary materials for this article are available online.     

Value at risk estimation under stochastic volatility models using adaptive PMCMC methods

In this paper, we propose a value-at-risk (VaR) estimation technique based on a new stochastic volatility model with leverage effect, nonconstant conditional mean and jump. In order to estimate the model parameters and latent state variables, we integrate the particle filter and adaptive Markov Chain Monte Carlo (MCMC) algorithms to develop a novel adaptive particle MCMC (A-PMCMC) algorithm. Comprehensive simulation experiments based on three stock indices and two foreign exchange time series show effectiveness of the proposed A-PMCMC algorithm and the VaR estimation technique.     

Multivariate stochastic volatility with Bayesian dynamic linear models

This paper develops a Bayesian procedure for estimation and forecasting of the volatility of multivariate time series. The foundation of this work is the matrix-variate dynamic linear model, for the volatility of which we adopt a multiplicative stochastic evolution, using Wishart and singular multivariate beta distributions. A diagonal matrix of discount factors is employed in order to discount the variances element by element and therefore allowing a flexible and pragmatic variance modelling approach. Diagnostic tests and sequential model monitoring are discussed in some detail. The proposed estimation theory is applied to a four-dimensional time series, comprising spot prices of aluminium, copper, lead and zinc of the London metal exchange. The empirical findings suggest that the proposed Bayesian procedure can be effectively applied to financial data, overcoming many of the disadvantages of existing volatility models.     

Analysis of a jump-diffusion option pricing model with serially correlated jump sizes

This paper extends the classical jump-diffusion option pricing model to incorporate serially correlated jump sizes which have been documented in recent empirical studies. We model the series of jump sizes by an autoregressive process and provide an analysis on the underlying stock return process. Based on this analysis, the European option price and the hedging parameters under the extended model are derived analytically. Through numerical examples, we investigate how the autocorrelation of jump sizes influences stock returns, option prices and hedging parameters, and demonstrate its effects on hedging portfolios and implied volatility smiles. A calibration example based on real market data is provided to show the advantage of incorporating the autocorrelation of jump sizes.     

Parameter Estimation Robust to Low-Frequency Contamination

We provide methods to robustly estimate the parameters of stationary ergodic short-memory time series models in the potential presence of additive low-frequency contamination. The types of contamination covered include level shifts (changes in mean) and monotone or smooth time trends, both of which have been shown to bias parameter estimates toward regions of persistence in a variety of contexts. The estimators presented here minimize trimmed frequency domain quasi-maximum likelihood (FDQML) objective functions without requiring specification of the low-frequency contaminating component. When proper sample size-dependent trimmings are used, the FDQML estimators are consistent and asymptotically normal, asymptotically eliminating the presence of any spurious persistence. These asymptotic results also hold in the absence of additive low-frequency contamination, enabling the practitioner to robustly estimate model parameters without prior knowledge of whether contamination is present. Popular time series models that fit into the framework of this article include autoregressive moving average (ARMA), stochastic volatility, generalized autoregressive conditional heteroscedasticity (GARCH), and autoregressive conditional heteroscedasticity (ARCH) models. We explore the finite sample properties of the trimmed FDQML estimators of the parameters of some of these models, providing practical guidance on trimming choice. Empirical estimation results suggest that a large portion of the apparent persistence in certain volatility time series may indeed be spurious. Supplementary materials for this article are available online.     

Monte Carlo Maximum Likelihood Estimation for Generalized Long-Memory Time Series Models

An exact maximum likelihood method is developed for the estimation of parameters in a non-Gaussian nonlinear density function that depends on a latent Gaussian dynamic process with long-memory properties. Our method relies on the method of importance sampling and on a linear Gaussian approximating model from which the latent process can be simulated. Given the presence of a latent long-memory process, we require a modification of the importance sampling technique. In particular, the long-memory process needs to be approximated by a finite dynamic linear process. Two possible approximations are discussed and are compared with each other. We show that an autoregression obtained from minimizing mean squared prediction errors leads to an effective and feasible method. In our empirical study, we analyze ten daily log-return series from the S&P 500 stock index by univariate and multivariate long-memory stochastic volatility models. We compare the in-sample and out-of-sample performance of a number of models within the class of long-memory stochastic volatility models.     

Estimating the Parameters of Stochastic Volatility Models Using Option Price Data

This article describes a maximum likelihood method for estimating the parameters of the standard square-root stochastic volatility model and a variant of the model that includes jumps in equity prices. The model is fitted to data on the S&P 500 Index and the prices of vanilla options written on the index, for the period 1990 to 2011. The method is able to estimate both the parameters of the physical measure (associated with the index) and the parameters of the risk-neutral measure (associated with the options), including the volatility and jump risk premia. The estimation is implemented using a particle filter whose efficacy is demonstrated under simulation. The computational load of this estimation method, which previously has been prohibitive, is managed by the effective use of parallel computing using graphics processing units (GPUs). The empirical results indicate that the parameters of the models are reliably estimated and consistent with values reported in previous work. In particular, both the volatility risk premium and the jump risk premium are found to be significant.     

Multivariate temporal disaggregation with cross-sectional constraints

Multivariate temporal disaggregation deals with the historical reconstruction and nowcasting of economic variables subject to temporal and contemporaneous aggregation constraints. The problem involves a system of time series that are related not only by a dynamic model but also by accounting constraints. The paper introduces two fundamental (and realistic) models that implement the multivariate best linear unbiased estimation approach that has potential application to the temporal disaggregation of the national accounts series. The multivariate regression model with random walk disturbances is most suitable to deal with the chained linked volumes (as the nature of the national accounts time series suggests); however, in this case the accounting constraints are not binding and the discrepancy has to be modeled by either a trend-stationary or an integrated process. The tiny, compared with other driving disturbances, size of the discrepancy prevents maximum-likelihood estimation to be carried out, and the parameters have to be estimated separately. The multivariate disaggregation with integrated random walk disturbances is suitable for the national accounts aggregates expressed at current prices, in which case the accounting constraints are binding.     

Comparing volatility forecasting models during the global financial crisis

Volatility estimation in financial markets has always been a challenge especially in time of crisis. Once asset prices and investment decisions are highly sensitive to such variable, many different models have been proposed in literature. This article estimates the volatility from a new family of stochastic volatility models called non-Gaussian State Space Models, a subclass of state space models where it is possible to compute exact likelihood. Volatilities of important Asian and Oceanian stock market indexes have been estimated and compared to APARCH model estimates. Results showed that non-Gaussian State Space Models outperformed significantly in both in-sample and forecasting cases.     

Bayesian analysis of stochastic volatility models with flexible tails

An alternative distributional assumption is proposed for the stochastic volatility model. This results in extremely flexible tail behaviour of the sampling distribution for the observables, as well as in the availability of a simple Markov Chain Monte Carlo strategy for posterior analysis. By allowing the tail behaviour to be determined by a separate parameter, we reserve the parameters of the volatility process to dictate the degree of volatility clustering. Treatment of a mean function is formally integrated in the analysis.

Some empirical examples on both stock prices and exchange rates clearly indicate the presence of fat tails, in combination with high levels of volatility clustering. In addition, predictive distributions indicate a good fit with these typical financial data sets.     

Bayesian analysis of stochastic volatility models with flexible tails

An alternative distributional assumption is proposed for the stochastic volatility model. This results in extremely flexible tail behaviour of the sampling distribution for the observables, as well as in the availability of a simple Markov Chain Monte Carlo strategy for posterior analysis. By allowing the tail behaviour to be determined by a separate parameter, we reserve the parameters of the volatility process to dictate the degree of volatility clustering. Treatment of a mean function is formally integrated in the analysis.

Some empirical examples on both stock prices and exchange rates clearly indicate the presence of fat tails, in combination with high levels of volatility clustering. In addition, predictive distributions indicate a good fit with these typical financial data sets.     

Time-varying autoregressive conditional duration model

The main goal of this work is to generalize the autoregressive conditional duration (ACD) model applied to times between trades to the case of time-varying parameters. The use of wavelets allows that parameters vary through time and makes possible the modeling of non-stationary processes without preliminary data transformations. The time-varying ACD model estimation was done by maximum-likelihood with standard exponential distributed errors. The properties of the estimators were assessed via bootstrap. We present a simulation exercise for a non-stationary process and an empirical application to a real series, namely the TELEMAR stock. Diagnostic and goodness of fit analysis suggest that the time-varying ACD model simultaneously modeled the dependence between durations, intra-day seasonality and volatility.     

Extreme quantile estimation based on financial time series

Estimation of market risk is an important problem in finance. Two well-known risk measures, viz., value at risk and median shortfall, turn out to be extreme quantiles of the marginal distribution of asset return. Time series on asset returns are known to exhibit certain stylized facts, such as heavy tails, skewness, volatility clustering, etc. Therefore, estimation of extreme quantiles in the presence of such features in the data seems to be of natural interest. It is difficult to capture most of these stylized facts using one specific time series model. This motivates nonparametric and extreme value theory-based estimation of extreme quantiles that do not require exact specification of the asset return model. We review these quantile estimators and compare their known properties. Their finite sample performance are compared using Monte Carlo simulation. We propose a new estimator that exhibits encouraging finite sample performance while estimating extreme quantile in the right tail region.     

Bounds for the expected value of the stochastic Divisia's price index

In this article, we estimate bounds for the expected value of the stochastic Divisia's price index, that is, we assume that prices and quantities of the given commodities are stochastic processes with continuous time. We consider some special case of the stochastic model in which prices and quantities are described by the geometric Brownian motion. It is shown that the precision of this estimation depends rather on the volatility of prices than quantities volatilities.     

Robust online scale estimation in time series: A model-free approach

This paper presents variance extraction procedures for univariate time series. The volatility of a times series is monitored allowing for non-linearities, jumps and outliers in the level. The volatility is measured using the height of triangles formed by consecutive observations of the time series. This idea was proposed by Rousseeuw and Hubert [1996. Regression-free and robust estimation of scale for bivariate data. Comput. Statist. Data Anal. 21, 67–85] in the bivariate setting. This paper extends their procedure to apply for online scale estimation in time series analysis. The statistical properties of the new methods are derived and finite sample properties are given. A financial and a medical application illustrate the use of the procedures.     

基于拟似然方法的股票收益与波动率关系及其应用研究

股票市场中收益与波动率的关系研究在金融证券领域起着很重要的作用,而随机波动率模型能够很好地拟合这种关系。本文将拟似然方法和渐近拟似然方法运用在随机波动率模型的参数估计方面,渐近拟似然方法可以避免因为人为的结构错误指定而造成的偏差,比较稳健。本文采用拟似然和渐近拟似然方法对随机波动率模型的参数估计进行了模拟探索,并和两种已有估计方法进行了对比,结果表明拟似然和渐近拟似然方法在模型的参数估计方面有着很好的估计结果。实证研究中,选取2000-2015年标普500指数作为研究对象,结果显示所选数据具有金融时间序列的常见特征。本文为金融证券领域中股票收益与波动率关系及其应用研究提供了一定的启示。     

中国股票市场波动率的多重分形分析与实证

我国股票市场波动表现出随时间变化的动态特征。文章采用多重消除趋势波动分析法(MFDFA),对沪深股市四个主要指数的日波动率时间序列进行了分析。结果表明,沪深股市四个主要指数的日波动率时间序列均表现出多重分形特征,且上证指数和中证500指数日波动率序列相对于其他两个指数日波动率序列表现出更强的多重分形特征。各指数日波动率时间序列的多重分形特征均是自身的长程相关性和波动的厚尾分布共同作用的结果,且波动的厚尾分布对原始序列的多重分形特征的影响比长程相关性大。     

A recursive approach to parameter estimation in regression and time series models

In this paper we discuss the recursive (or on line) estimation in (i) regression and (ii) autoregressive integrated moving average (ARIMA) time series models. The adopted approach uses Kalman filtering techniques to calculate estimates recursively. This approach is used for the estimation of constant as well as time varying parameters. In the first section of the paper we consider the linear regression model. We discuss recursive estimation both for constant and time varying parameters. For constant parameters, Kalman filtering specializes to recursive least squares. In general, we allow the parameters to vary according to an autoregressive integrated moving average process and update the parameter estimates recursively. Since the stochastic model for the parameter changes will "be rarely known, simplifying assumptions have to be made. In particular we assume a random walk model for the time varying parameters and show how to determine whether the parameters are changing over time. This is illustrated with an example.