Dynamics of state price densities

State price densities (SPDs) are an important element in applied quantitative finance. In a Black–Scholes world they are lognormal distributions, but in practice volatility changes and the distribution deviates from log-normality. In order to study the degree of this deviation, we estimate SPDs using EUREX option data on the DAX index via a nonparametric estimator of the second derivative of the (European) call pricing function. The estimator is constrained so as to satisfy no-arbitrage constraints and corrects for the intraday covariance structure in option prices. In contrast to existing methods, we do not use any parametric or smoothness assumptions.     

Nonparametric state price density estimation using constrained least squares and the bootstrap

The economic theory of option pricing imposes constraints on the structure of call functions and state price densities. Except in a few polar cases, it does not prescribe functional forms. This paper proposes a nonparametric estimator of option pricing models which incorporates various restrictions (such as monotonicity and convexity) within a single least squares procedure. The bootstrap is used to produce confidence intervals for the call function and its first two derivatives and to calibrate a residual regression test of shape constraints. We apply the techniques to option pricing data on the DAX.     

Asymptotic normality of narrow-band least squares in the stationary fractional cointegration model and volatility forecasting

We consider semiparametric frequency domain analysis of cointegration between long memory processes, i.e. fractional cointegration, allowing derivation of useful long-run relations even among stationary processes. The approach is due to Robinson (1994b. Annals of Statistics 22, 515–539) and uses a degenerating part of the periodogram near the origin to form a narrow-band frequency domain least squares (FDLS) estimator of the cointegrating relation, which is consistent for arbitrary short-run dynamics. We derive the asymptotic distribution theory for the FDLS estimator of the cointegration vector in the stationary long memory case, thus complementing Robinson's consistency result. An application to the relation between the volatility realized in the stock market and the associated implicit volatility derived from option prices is offered.     

A two-stage procedure for partially identified models

This paper studies a two-stage procedure for estimating partially identified models, based on Chernozhukov, Hong, and Tamer’s (2007) theory of set estimation and inference. We consider the case where a sub-vector of parameters or their identified set can be estimated separately from the rest, possibly subject to a priori restrictions. Our procedure constructs the second-stage set estimator and confidence set by taking appropriate level sets of a criterion function, using a first-stage estimator to impose restrictions on the parameter of interest. We give conditions under which the two-stage set estimator is a set-valued random element that is measurable in an appropriate sense. We also establish the consistency of the two-stage set estimator.     

Testing and imposing Slutsky symmetry in nonparametric demand systems

Maximization of utility implies that consumer demand systems have a Slutsky matrix which is everywhere symmetric. However, previous non- and semi-parametric approaches to the estimation of consumer demand systems do not give estimators that are restricted to satisfy this condition, nor do they offer powerful tests of this restriction. We use nonparametric modeling to test and impose Slutsky symmetry in a system of expenditure share equations over prices and expenditure. In this context, Slutsky symmetry is a set of nonlinear cross-equation restrictions on levels and derivatives of consumer demand equations. The key insight is that due to the differing convergence rates of levels and derivatives and due to the fact that the symmetry restrictions are linear in derivatives, both the test and the symmetry restricted estimator behave asymptotically as if these restrictions were (locally) linear. We establish large and finite sample properties of our methods, and show that our test has advantages over the only other comparable test. All methods we propose are implemented with Canadian micro-data. We find that our nonparametric analysis yields statistically significantly and qualitatively different results from traditional parametric estimators and tests.     

Instrumental variable estimation of nonseparable models

There are many environments where knowledge of a structural relationship is required to answer questions of interest. Also, nonseparability of a structural disturbance is a key feature of many models. Here, we consider nonparametric identification and estimation of a model that is monotonic in a nonseparable scalar disturbance, which disturbance is independent of instruments. This model leads to conditional quantile restrictions. We give local identification conditions for the structural equations from those quantile restrictions. We find that a modified completeness condition is sufficient for local identification. We also consider estimation via a nonparametric minimum distance estimator. The estimator minimizes the sum of squares of predicted values from a nonparametric regression of the quantile residual on the instruments. We show consistency of this estimator.     

On implied volatility for options—Some reasons to smile and more to correct

We analyze the properties of the implied volatility, the commonly used volatility estimator by direct option price inversion. It is found that the implied volatility is subject to a systematic bias in the presence of pricing errors, which makes it inconsistent to the underlying volatility. We propose an estimator of the underlying volatility by first estimating nonparametrically the option price function, followed by inverting the nonparametrically estimated price. It is shown that the approach removes the adverse impacts of the pricing errors and produces a consistent volatility estimator for a wide range of option price models. We demonstrate the effectiveness of the proposed approach by numerical simulation and empirical analysis on S&P 500 option data.     

Local GMM estimation of time series models with conditional moment restrictions

This paper investigates statistical properties of the local generalized method of moments (LGMM) estimator for some time series models defined by conditional moment restrictions. First, we consider Markov processes with possible conditional heteroskedasticity of unknown forms and establish the consistency, asymptotic normality, and semi-parametric efficiency of the LGMM estimator. Second, we undertake a higher-order asymptotic expansion and demonstrate that the LGMM estimator possesses some appealing bias reduction properties for positively autocorrelated processes. Our analysis of the asymptotic expansion of the LGMM estimator reveals an interesting contrast with the OLS estimator that helps to shed light on the nature of the bias correction performed by the LGMM estimator. The practical importance of these findings is evaluated in terms of a bond and option pricing exercise based on a diffusion model for spot interest rate.     

Estimation of partial differential equations with applications in finance

Linear parabolic partial differential equations (PDE’s) and diffusion models are closely linked through the celebrated Feynman–Kac representation of solutions to PDE’s. In asset pricing theory, this leads to the representation of derivative prices as solutions to PDE’s. Very often implied derivative prices are calculated given preliminary estimates of the diffusion model for the underlying variable. We demonstrate that the implied derivative prices are consistent and derive their asymptotic distribution under general conditions. We apply this result to three leading cases of preliminary estimators: Nonparametric, semiparametric and fully parametric ones. In all three cases, the asymptotic distribution of the solution is derived. We demonstrate the use of these results in obtaining confidence bands and standard errors for implied prices of bonds, options and other derivatives. Our general results also are of interest for the estimation of diffusion models using either historical data of the underlying process or option prices; these issues are also discussed.     

Semi-nonparametric estimation of the call-option price surface under strike and time-to-expiry no-arbitrage constraints

We suggest a semi-nonparametric estimator for the call-option price surface. The estimator is a bivariate tensor-product B-spline. To enforce no-arbitrage constraints across strikes and expiry dates, we establish sufficient no-arbitrage conditions on the control net of the B-spline surface. The conditions are linear and therefore allow for an implementation of the estimator by means of standard quadratic programming techniques. The consistency of the estimator is proved. By means of simulations, we explore the statistical efficiency benefits that are associated with estimating option price surfaces and state-price densities under the full set of no-arbitrage constraints. We estimate a call-option price surface, families of first-order strike derivatives, and state-price densities for S&P 500 option data.     

Identification and information in monotone binary models

This paper considers binary response models where errors are uncorrelated with a set of instrumental variables and are independent of a continuous regressor v$v$, conditional on all other variables. It is shown that these exclusion restrictions are not sufficient for identification and that additional identifying assumptions are needed. Such an assumption, introduced by Lewbel [Semiparametric qualitative response model estimation with unknown heteroskedasticity or instrumental variables. Journal of Econometrics 97, 145–177], is that the support of the continuous regressor is large, but we show that it significantly restricts the class of binary phenomena which can be analysed. We propose an alternative additional assumption under which β$\beta$ remains just identified and the estimation unchanged. This alternative assumption does not impose specific restrictions on the data, which broadens the scope of the estimation method in empirical work. The semiparametric efficiency bound of the model is also established and an existing estimator is shown to achieve that bound. The efficient estimator uses a plug-in density estimate. It is shown that plugging in the true density rather than an estimate is inefficient. Extensions to ordered choice models are provided.     

Pre-averaging estimators of the ex-post covariance matrix in noisy diffusion models with non-synchronous data

We show how pre-averaging can be applied to the problem of measuring the ex-post covariance of financial asset returns under microstructure noise and non-synchronous trading. A pre-averaged realised covariance is proposed, and we present an asymptotic theory for this new estimator, which can be configured to possess an optimal convergence rate or to ensure positive semi-definite covariance matrix estimates. We also derive a noise-robust Hayashi–Yoshida estimator that can be implemented on the original data without prior alignment of prices. We uncover the finite sample properties of our estimators with simulations and illustrate their practical use on high-frequency equity data.     

Estimation of conditional moment restrictions without assuming parameter identifiability in the implied unconditional moments

A well-known difficulty in estimating conditional moment restrictions is that the parameters of interest need not be globally identified by the implied unconditional moments. In this paper, we propose an approach to constructing a continuum of unconditional moments that can ensure parameter identifiability. These unconditional moments depend on the “instruments” generated from a “generically comprehensively revealing” function, and they are further projected along the exponential Fourier series. The objective function is based on the resulting Fourier coefficients, from which an estimator can be easily computed. A novel feature of our method is that the full continuum of unconditional moments is incorporated into each Fourier coefficient. We show that, when the number of Fourier coefficients in the objective function grows at a proper rate, the proposed estimator is consistent and asymptotically normally distributed. An efficient estimator is also readily obtained via the conventional two-step GMM method. Our simulations confirm that the proposed estimator compares favorably with that of Domínguez and Lobato (2004, Econometrica) in terms of bias, standard error, and mean squared error.     

Sparse estimators and the oracle property,or the return of Hodges’ estimator

We point out some pitfalls related to the concept of an oracle property as used in Fan and Li [2001. Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Statistical Association 96, 1348–1360; 2002. Variable selection for Cox's proportional hazards model and frailty model. Annals of Statistics 30, 74–99; 2004. New estimation and model selection procedures for semiparametric modeling in longitudinal data analysis. Journal of the American Statistical Association 99, 710–723] which are reminiscent of the well-known pitfalls related to Hodges’ estimator. The oracle property is often a consequence of sparsity of an estimator. We show that any estimator satisfying a sparsity property has maximal risk that converges to the supremum of the loss function; in particular, the maximal risk diverges to infinity whenever the loss function is unbounded. For ease of presentation the result is set in the framework of a linear regression model, but generalizes far beyond that setting. In a Monte Carlo study we also assess the extent of the problem in finite samples for the smoothly clipped absolute deviation (SCAD) estimator introduced in Fan and Li [2001. Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Statistical Association 96, 1348–1360]. We find that this estimator can perform rather poorly in finite samples and that its worst-case performance relative to maximum likelihood deteriorates with increasing sample size when the estimator is tuned to sparsity.     

Smooth minimum distance estimation and testing with conditional estimating equations: Uniform in bandwidth theory

To study the influence of a bandwidth parameter in inference with conditional moments, we propose a new class of estimators and establish an asymptotic representation of our estimator as a process indexed by a bandwidth, which can vary within a wide range including bandwidths independent of the sample size. We study its behavior under misspecification. We also propose an efficient version of our estimator. We develop a procedure based on a distance metric statistic for testing restrictions on parameters as well as a bootstrap technique to account for the bandwidth’s influence. Our new methods are simple to implement, apply to non-smooth problems, and perform well in our simulations.     

A semiparametric model for binary response and continuous outcomes under index heteroscedasticity

This paper formulates a likelihood‐based estimator for a double‐index, semiparametric binary response equation. A novel feature of this estimator is that it is based on density estimation under local smoothing. While the proofs differ from those based on alternative density estimators, the finite sample performance of the estimator is significantly improved. As binary responses often appear as endogenous regressors in continuous outcome equations, we also develop an optimal instrumental variables estimator in this context. For this purpose, we specialize the double‐index model for binary response to one with heteroscedasticity that depends on an index different from that underlying the ‘mean response’. We show that such (multiplicative) heteroscedasticity, whose form is not parametrically specified, effectively induces exclusion restrictions on the outcomes equation. The estimator developed exploits such identifying information. We provide simulation evidence on the favorable performance of the estimators and illustrate their use through an empirical application on the determinants, and affect, of attendance at a government‐financed school. Copyright © 2009 John Wiley & Sons, Ltd.     

Hermite polynomial based expansion of European option prices

We seek a closed-form series approximation of European option prices under a variety of diffusion models. The proposed convergent series are derived using the Hermite polynomial approach. Departing from the usual option pricing routine in the literature, our model assumptions have no requirements for affine dynamics or explicit characteristic functions. Moreover, convergent expansions provide a distinct insight into how and on which order the model parameters affect option prices, in contrast with small-time asymptotic expansions in the literature. With closed-form expansions, we explicitly translate model features into option prices, such as mean-reverting drift and self-exciting or skewed jumps. Numerical examples illustrate the accuracy of this approach and its advantage over alternative expansion methods.     

Structural vector autoregressions with Markov switching: Combining conventional with statistical identification of shocks

In structural vector autoregressive (SVAR) analysis a Markov regime switching (MS) property can be exploited to identify shocks if the reduced form error covariance matrix varies across regimes. Unfortunately, these shocks may not have a meaningful structural economic interpretation. It is discussed how statistical and conventional identifying information can be combined. The discussion is based on a VAR model for the US containing oil prices, output, consumer prices and a short-term interest rate. The system has been used for studying the causes of the early millennium economic slowdown based on traditional identification with zero and long-run restrictions and using sign restrictions. We find that previously drawn conclusions are questionable in our framework.     

Do interest rate options contain information about excess returns?

There is strong empirical evidence that long-term interest rates contain a time-varying risk premium. Options may contain valuable information about this risk premium because their prices are sensitive to the underlying interest rates. We use the joint time series of swap rates and interest rate option prices to estimate dynamic term structure models. The risk premiums that we estimate using option prices are better able to predict excess returns for long-term swaps over short-term swaps. Moreover, in contrast to the previous literature, the most successful models for predicting excess returns have risk factors with stochastic volatility. We also show that the stochastic volatility models we estimate using option prices match the failure of the expectations hypothesis.     

Estimation and inference of dynamic structural factor models with over-identifying restrictions

This paper develops a new estimator for the impulse response functions in structural factor models with a fixed number of over-identifying restrictions. The proposed identification scheme nests the conventional just-identified recursive scheme as a special case. We establish the asymptotic distributions of the new estimator and develop test statistics for the over-identifying restrictions. Simulation results show that adding a few more over-identifying restrictions can lead to a substantial improvement in estimation accuracy for impulse response functions at both zero and nonzero horizons. We estimate the effects of a monetary policy shock using a U.S. data set. The results show that our over-identified scheme can help to detect incorrect specifications that lead to spurious impulse responses.