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.     

Bootstrapping realized multivariate volatility measures

We propose a bootstrap method for statistics that are a function of multivariate high frequency returns such as realized regression, covariance and correlation coefficients. We show that the finite sample performance of the bootstrap is superior to the existing first-order asymptotic theory. Nevertheless, and contrary to the existing results in the bootstrap literature for regression models subject to error heteroskedasticity, the Edgeworth expansion for the pairs bootstrap that we develop here shows that this method is not second-order accurate. We argue that this is due to the fact that the conditional mean parameters of realized regression models are heterogeneous under stochastic volatility.     

Testing multivariate distributions in GARCH models

In this paper, we consider testing distributional assumptions in multivariate GARCH models based on empirical processes. Using the fact that joint distribution carries the same amount of information as the marginal together with conditional distributions, we first transform the multivariate data into univariate independent data based on the marginal and conditional cumulative distribution functions. We then apply the Khmaladze's martingale transformation (K-transformation) to the empirical process in the presence of estimated parameters. The K-transformation eliminates the effect of parameter estimation, allowing a distribution-free test statistic to be constructed. We show that the K-transformation takes a very simple form for testing multivariate normal and multivariate t-distributions. The procedure is applied to a multivariate financial time series data set.     

Efficient estimation of a multivariate multiplicative volatility model

We propose a multivariate generalization of the multiplicative volatility model of Engle and Rangel (2008), which has a nonparametric long run component and a unit multivariate GARCH short run dynamic component. We suggest various kernel-based estimation procedures for the parametric and nonparametric components, and derive the asymptotic properties thereof. For the parametric part of the model, we obtain the semiparametric efficiency bound. Our method is applied to a bivariate stock index series. We find that the univariate model of Engle and Rangel (2008) appears to be violated in the data whereas our multivariate model is more consistent with the data.     

Forecasting multivariate realized stock market volatility

We present a new matrix-logarithm model of the realized covariance matrix of stock returns. The model uses latent factors which are functions of lagged volatility, lagged returns and other forecasting variables. The model has several advantages: it is parsimonious; it does not require imposing parameter restrictions; and, it results in a positive-definite estimated covariance matrix. We apply the model to the covariance matrix of size-sorted stock returns and find that two factors are sufficient to capture most of the dynamics.     

Local rank tests in a multivariate nonparametric relationship

Consider a multivariate nonparametric model where the unknown vector of functions depends on two sets of explanatory variables. For a fixed level of one set of explanatory variables, we provide consistent statistical tests, called local rank tests, to determine whether the multivariate relationship can be explained by a smaller number of functions. We also provide estimators for the smallest number of functions, called local rank, explaining the relationship. The local rank tests and the estimators of local rank are defined in terms of the eigenvalues of a kernel-based estimator of some matrix. The asymptotics of the eigenvalues is established by using the so-called Fujikoshi expansion along with some techniques of the theory of U-statistics. We present a simulation study which examines the small sample properties of local rank tests. We also apply the local rank tests and the local rank estimators to a demand system given by a newly constructed data set. This work can be viewed as a “local” extension of the tests for a number of factors in a nonparametric relationship introduced by Stephen Donald.     

Estimation and model selection of semiparametric multivariate survival functions under general censorship

We study estimation and model selection of semiparametric models of multivariate survival functions for censored data, which are characterized by possibly misspecified parametric copulas and nonparametric marginal survivals. We obtain the consistency and root-n$n$ asymptotic normality of a two-step copula estimator to the pseudo-true copula parameter value according to KLIC, and provide a simple consistent estimator of its asymptotic variance, allowing for a first-step nonparametric estimation of the marginal survivals. We establish the asymptotic distribution of the penalized pseudo-likelihood ratio statistic for comparing multiple semiparametric multivariate survival functions subject to copula misspecification and general censorship. An empirical application is provided.     

Testing for Seasonal Unit Roots in Monthly Panels of Time Series*

We consider the problem of testing for seasonal unit roots in monthly panel data. To this aim, we generalize the quarterly cross‐sectionally augmented Hylleberg–Engle–Granger–Yoo (CHEGY) test to the monthly case. This parametric test is contrasted with a new non‐parametric test, which is the panel counterpart to the univariate record unit–root seasonal (RURS) test that relies on counting extrema in time series. All methods are applied to an empirical data set on tourism in Austrian provinces. The power properties of the tests are evaluated in simulation experiments that are tuned to the tourism data.     

Measuring Conditional Persistence in Nonlinear Time Series*

The persistence properties of economic time series have been a primary object of investigation in a variety of guises since the early days of econometrics. Recently, work on nonlinear modelling for time series has introduced the idea that persistence of a shock at a point in time may vary depending on the state of the process at that point in time. This article suggests investigating the persistence of processes conditioning on their history as a tool that may aid parametric nonlinear modelling. In particular, we suggest that examining the nonparametrically estimated derivatives of the conditional expectation of a variable with respect to its lag(s) may be a useful indicator of the variation in persistence with respect to its past history. We discuss in detail the implementation of the measure and present a Monte Carlo investigation. We further apply the persistence analysis to real exchange rates.     

A Bayesian approach to bandwidth selection for multivariate kernel regression with an application to state-price density estimation

This paper presents a Bayesian approach to bandwidth selection for multivariate kernel regression. A Monte Carlo study shows that under the average squared error criterion, the Bayesian bandwidth selector is comparable to the cross-validation method and clearly outperforms the bootstrapping and rule-of-thumb bandwidth selectors. The Bayesian bandwidth selector is applied to a multivariate kernel regression model that is often used to estimate the state-price density of Arrow–Debreu securities with the S&P 500 index options data and the DAX index options data. The proposed Bayesian bandwidth selector represents a data-driven solution to the problem of choosing bandwidths for the multivariate kernel regression involved in the nonparametric estimation of the state-price density pioneered by Aït-Sahalia and Lo [Aït-Sahalia, Y., Lo, A.W., 1998. Nonparametric estimation of state-price densities implicit in financial asset prices. The Journal of Finance, 53, 499, 547.]     

Practical estimation of multivariate densities using wavelet methods

This paper describes a practical method for estimating multivariate densities using wavelets. As in kernel methods, wavelet methods depend on two types of parameters. On the one hand we have a functional parameter: the wavelet Ø (comparable to the kernel K ) and on the other hand we have a smoothing parameter: the resolution index (comparable to the bandwidth h ). Classically, we determine the resolution index with a cross-validation method. The advantage of wavelet methods compared to kernel methods is that we have a technique for choosing the wavelet Ø among a fixed family. Moreover, the wavelets method simplifies significantly both the theoretical and the practical computations.     

Nonlinear causality testing with stepwise multivariate filtering: Evidence from stock and currency markets

We examine the spillovers of the US subprime crisis to Asian and European economies and in particular to what extent currency and stock markets have been affected by the crisis. Linear and nonlinear dependencies are detected after pairwise and system-wise causality analysis. A new stepwise multivariate filtering approach is implemented after controlling for conditional heteroskedasticity in the raw data and in VAR/VECM residuals using multivariate GARCH models. Significant nonlinear causal linkages persisted even after the application of GARCH-BEKK, CCC-GARCH and DCC-GARCH modelling. This indicates that volatility effects might partly induce nonlinear causality. Perhaps new short-term asset-pricing models could be developed to explain this stylized fact. These results might also have important implications for hedging, trading strategies and financial market regulation.     

On the Geometry of the Instrumental Variable Estimator*

I derive the exact distribution of the exact determined instrumental variable estimator using a geometric approach. The approach provides a decomposition of the exact estimator. The results show that by geometric reasoning one may efficiently derive the distribution of the estimation error. The often striking non‐normal shape of the instrumental variable estimator, in the case of weak instruments and small samples, follows intuitively by the geometry of the problem. The method allows for intuitive interpretations of how the shape of the distribution is determined by instrument quality and endogeneity. The approach can also be used when deriving the exact distribution of any ratio of stochastic variables.     

Studying co-movements in large multivariate data prior to multivariate modelling

For non-stationary vector autoregressive models (VAR hereafter, or VAR with moving average, VARMA hereafter), we show that the presence of common cyclical features or cointegration leads to a reduction of the order of the implied univariate autoregressive-integrated-moving average (ARIMA hereafter) models. This finding can explain why we identify parsimonious univariate ARIMA models in applied research although VAR models of typical order and dimension used in macroeconometrics imply non-parsimonious univariate ARIMA representations.     

Time-varying combinations of predictive densities using nonlinear filtering

We propose a Bayesian combination approach for multivariate predictive densities which relies upon a distributional state space representation of the combination weights. Several specifications of multivariate time-varying weights are introduced with a particular focus on weight dynamics driven by the past performance of the predictive densities and the use of learning mechanisms. In the proposed approach the model set can be incomplete, meaning that all models can be individually misspecified. A Sequential Monte Carlo method is proposed to approximate the filtering and predictive densities. The combination approach is assessed using statistical and utility-based performance measures for evaluating density forecasts of simulated data, US macroeconomic time series and surveys of stock market prices. Simulation results indicate that, for a set of linear autoregressive models, the combination strategy is successful in selecting, with probability close to one, the true model when the model set is complete and it is able to detect parameter instability when the model set includes the true model that has generated subsamples of data. Also, substantial uncertainty appears in the weights when predictors are similar; residual uncertainty reduces when the model set is complete; and learning reduces this uncertainty. For the macro series we find that incompleteness of the models is relatively large in the 1970’s, the beginning of the 1980’s and during the recent financial crisis, and lower during the Great Moderation; the predicted probabilities of recession accurately compare with the NBER business cycle dating; model weights have substantial uncertainty attached. With respect to returns of the S&P 500 series, we find that an investment strategy using a combination of predictions from professional forecasters and from a white noise model puts more weight on the white noise model in the beginning of the 1990’s and switches to giving more weight to the professional forecasts over time. Information on the complete predictive distribution and not just on some moments turns out to be very important, above all during turbulent times such as the recent financial crisis. More generally, the proposed distributional state space representation offers great flexibility in combining densities.     

Combining VAR and DSGE forecast densities

A popular macroeconomic forecasting strategy utilizes many models to hedge against instabilities of unknown timing; see (among others) Stock and Watson (2004), Clark and McCracken (2010), and Jore et al. (2010). Existing studies of this forecasting strategy exclude dynamic stochastic general equilibrium (DSGE) models, despite the widespread use of these models by monetary policymakers. In this paper, we use the linear opinion pool to combine inflation forecast densities from many vector autoregressions (VARs) and a policymaking DSGE model. The DSGE receives a substantial weight in the pool (at short horizons) provided the VAR components exclude structural breaks. In this case, the inflation forecast densities exhibit calibration failure. Allowing for structural breaks in the VARs reduces the weight on the DSGE considerably, but produces well-calibrated forecast densities for inflation.     

Nonparametric/semiparametric estimation and testing of econometric models with data dependent smoothing parameters

We consider nonparametric/semiparametric estimation and testing of econometric models with data dependent smoothing parameters. Most of the existing works on asymptotic distributions of a nonparametric/semiparametric estimator or a test statistic are based on some deterministic smoothing parameters, while in practice it is important to use data-driven methods to select the smoothing parameters. In this paper we give a simple sufficient condition that can be used to establish the first order asymptotic equivalence of a nonparametric estimator or a test statistic with stochastic smoothing parameters to those using deterministic smoothing parameters. We also allow for general weakly dependent data.     

Manipulation of the running variable in the regression discontinuity design: A density test

Standard sufficient conditions for identification in the regression discontinuity design are continuity of the conditional expectation of counterfactual outcomes in the running variable. These continuity assumptions may not be plausible if agents are able to manipulate the running variable. This paper develops a test of manipulation related to continuity of the running variable density function. The methodology is applied to popular elections to the House of Representatives, where sorting is neither expected nor found, and to roll call voting in the House, where sorting is both expected and found.     

Temporal aggregation of multivariate GARCH processes

This paper derives results for the temporal aggregation of multivariate GARCH(1,1) processes in the general vector specification. It is shown that the class of weak multivariate GARCH(1,1) processes is closed under temporal aggregation. Fourth moment characteristics turn out to be crucial for the low frequency dynamics for both stock and flow variables. In some aspects, the aggregation characteristics of multivariate GARCH processes are shown to be different from those of vector autoregressive moving average processes. A numerical example illustrates some of the results.     

Jumps in equilibrium prices and market microstructure noise

Asset prices observed in financial markets combine equilibrium prices and market microstructure noise. In this paper, we study how to tell apart large shifts in equilibrium prices from noise using high frequency data. We propose a new nonparametric test which allows us to asymptotically remove the noise from observable price data and to discover jumps in fundamental asset values. We provide its asymptotic distribution to decide when such jumps occur. In finite samples, our test offers reasonable power for distinguishing between noise and jumps. Empirical evidence indicates that it is necessary to incorporate the presence of jumps in equilibrium prices.