Testing whether the underlying continuous-time process follows a diffusion: An infinitesimal operator-based approach

We develop a nonparametric test to check whether a process can be represented by a stochastic differential equation driven only by a Brownian motion. Our testing procedure utilizes the infinitesimal operator-based martingale characterization combined with a generalized spectral approach. Such a testing procedure is feasible and convenient because the infinitesimal operator of the diffusion process has a closed-form expression. The proposed test is applicable to both univariate and multivariate processes and has an N(0,1)$N(0,1)$ limit distribution under the diffusion hypothesis. Simulation and empirical studies show that the proposed test has reasonable performance in small samples.     

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.     

A quasi-maximum likelihood method for estimating the parameters of multivariate diffusions

A quasi-maximum likelihood procedure for estimating the parameters of multi-dimensional diffusions is developed in which the transitional density is a multivariate Gaussian density with first and second moments approximating the true moments of the unknown density. For affine drift and diffusion functions, the moments are exactly those of the true transitional density and for nonlinear drift and diffusion functions the approximation is extremely good and is as effective as alternative methods based on likelihood approximations. The estimation procedure generalises to models with latent factors. A conditioning procedure is developed that allows parameter estimation in the absence of proxies.     

Parameter estimation and bias correction for diffusion processes

This paper considers parameter estimation for continuous-time diffusion processes which are commonly used to model dynamics of financial securities including interest rates. To understand why the drift parameters are more difficult to estimate than the diffusion parameter, as observed in previous studies, we first develop expansions for the bias and variance of parameter estimators for two of the most employed interest rate processes, Vasicek and CIR processes. Then, we study the first order approximate maximum likelihood estimator for linear drift processes. A parametric bootstrap procedure is proposed to correct bias for general diffusion processes with a theoretical justification. Simulation studies confirm the theoretical findings and show that the bootstrap proposal can effectively reduce both the bias and the mean square error of parameter estimates, for both univariate and multivariate processes. The advantages of using more accurate parameter estimators when calculating various option prices in finance are demonstrated by an empirical study.     

Empirical likelihood-based inference for nonparametric recurrent diffusions

This paper provides a new approach to constructing confidence intervals for nonparametric drift and diffusion functions in the continuous-time diffusion model via empirical likelihood (EL). The log EL ratios are constructed through the estimating equations satisfied by the local linear estimators. Limit theories are developed by means of increasing time span and shrinking observational intervals. The results apply to both stationary and nonstationary recurrent diffusion processes. Simulations show that for both drift and diffusion functions, the new procedure performs remarkably well in finite samples and clearly dominates the conventional method in constructing confidence intervals based on asymptotic normality. An empirical example is provided to illustrate the usefulness of the proposed method.     

Asymptotically distribution-free tests for the volatility function of a diffusion

This paper develops two tests for parametric volatility function of a diffusion model based on Khmaladze (1981)’s martingale transformation. The tests impose no restrictions on the functional form of the drift function and are shown to be asymptotically distribution-free. The tests are consistent against a large class of fixed alternatives and have nontrivial power against a class of root-n$n$ local alternatives. The paper also extends the tests of volatility to testing for joint specifications of drift and volatility. Monte Carlo simulations show that the tests perform well in finite samples. The proposed tests are then applied to testing models of short-term interest, using data of Treasury bill rate and Eurodollar deposit rate. The empirical results show that the commonly used CKLS volatility function of Chan et al. (1992) fits volatility function poorly and none of the parametric interest rate models considered in the paper fit data well.     

Extended Neyman smooth goodness-of-fit tests,applied to competing heavy-tailed distributions

A simplified version of the Neyman (1937) “Smooth” goodness-of-fit test is extended to account for the presence of estimated model parameters, thereby removing overfitting bias. Using a Lagrange Multiplier approach rather than the Likelihood Ratio statistic proposed by Neyman greatly simplifies the calculations. Polynomials, splines, and the step function of Pearson’s test are compared as alternative perturbations to the theoretical uniform distribution. The extended tests have negligible size distortion and more power than standard tests. The tests are applied to competing symmetric leptokurtic distributions with US stock return data. These are generally rejected, primarily because of the presence of skewness.     

Evaluation of dynamic stochastic general equilibrium models based on distributional comparison of simulated and historical data

We take as a starting point the existence of a joint distribution implied by different dynamic stochastic general equilibrium (DSGE) models, all of which are potentially misspecified. Our objective is to compare “true” joint distributions with ones generated by given DSGEs. This is accomplished via comparison of the empirical joint distributions (or confidence intervals) of historical and simulated time series. The tool draws on recent advances in the theory of the bootstrap, Kolmogorov type testing, and other work on the evaluation of DSGEs, aimed at comparing the second order properties of historical and simulated time series. We begin by fixing a given model as the “benchmark” model, against which all “alternative” models are to be compared. We then test whether at least one of the alternative models provides a more “accurate” approximation to the true cumulative distribution than does the benchmark model, where accuracy is measured in terms of distributional square error. Bootstrap critical values are discussed, and an illustrative example is given, in which it is shown that alternative versions of a standard DSGE model in which calibrated parameters are allowed to vary slightly perform equally well. On the other hand, there are stark differences between models when the shocks driving the models are assigned non-plausible variances and/or distributional assumptions.     

Estimation of market power in the presence of firm level inefficiencies

“The quiet life hypothesis” (QLH) by Hicks (1935) argues that, due to management’s subjective cost of reaching optimal profits, firms use their market power to allow inefficient allocation of resources. Increasing competitive pressure is therefore likely to force management to work harder to reach optimal profits. Another hypothesis, which also relates market power to efficiency is “the efficient structure hypothesis” (ESH) by Demsetz (1973). ESH argues that firms with superior efficiencies or technologies have lower costs and therefore higher profits. These firms are assumed to gain larger market shares which lead to higher concentration. Ignoring the efficiency levels of the firms in a market power model might cause both estimation and interpretation problems. Unfortunately, the literature on market power measurement largely ignores this relationship. In the context of a dynamic setting, we estimate the market power of US airlines in two city-pairs by both allowing inefficiencies of the firms and not allowing inefficiencies of the firms. Using industry level cost data, we estimate the cost function parameters and time-varying efficiencies. An instrumental variables version of the square root Kalman filter is used to estimate time-varying conduct parameters.     

A Wald test for the cointegration rank in nonstationary fractional systems

This paper develops new methods for determining the cointegration rank in a nonstationary fractionally integrated system, extending univariate optimal methods for testing the degree of integration. We propose a simple Wald test based on the singular value decomposition of the unrestricted estimate of the long run multiplier matrix. When the “strength” of the cointegrating relationship is less than 1/2, the test statistic has a standard asymptotic distribution, like Lagrange Multiplier tests exploiting local properties. We consider the behavior of our test under estimation of short run parameters and local alternatives. We compare our procedure with other cointegration tests based on different principles and find that the new method has better properties in a range of situations by using information on the alternative obtained through a preliminary estimate of the cointegration strength.     

Trend/cycle decomposition of regime-switching processes

We present a new approach to trend/cycle decomposition of time series that follow regime-switching processes. The proposed approach, which we label the “regime-dependent steady-state” (RDSS) decomposition, is motivated as the appropriate generalization of the Beveridge and Nelson decomposition [Beveridge, S., Nelson, C.R., 1981. A new approach to decomposition of economic time series into permanent and transitory components with particular attention to measurement of the business cycle. Journal of Monetary Economics 7, 151–174] to the setting where the reduced-form dynamics of a given series can be captured by a regime-switching forecasting model. For processes in which the underlying trend component follows a random walk with possibly regime-switching drift, the RDSS decomposition is optimal in a minimum mean-squared-error sense and is more broadly applicable than directly employing an Unobserved Components model.     

Semi-nonparametric estimation and misspecification testing of diffusion models

Novel transition-based misspecification tests of semiparametric and fully parametric univariate diffusion models based on the estimators developed in [Kristensen, D., 2010. Pseudo-maximum likelihood estimation in two classes of semiparametric diffusion models. Journal of Econometrics 156, 239-259] are proposed. It is demonstrated that transition-based tests in general lack power in detecting certain departures from the null since they integrate out local features of the drift and volatility. As a solution to this, tests that directly compare drift and volatility estimators under the relevant null and alternative are also developed which exhibit better power against local alternatives.     

Testing for structural stability of factor augmented forecasting models

Mild factor loading instability, particularly if sufficiently independent across the different constituent variables, does not affect the estimation of the number of factors, nor subsequent estimation of the factors themselves (see e.g.  Stock and Watson (2009)). This result does not hold in the presence of large common breaks in the factor loadings, however. In this case, information criteria overestimate the number of breaks. Additionally, estimated factors are no longer consistent estimators of “true” factors. Hence, various recent research papers in the diffusion index literature focus on testing the constancy of factor loadings. However, forecast failure of factor augmented models can be due to either factor loading instability, regression coefficient instability, or both. To address this issue, we develop a test for the joint hypothesis of structural stability of both factor loadings and factor augmented forecasting model regression coefficients. Our proposed test statistic has a chi-squared limiting distribution, and we are able to establish the first order validity of (block) bootstrap critical values. Empirical evidence is also presented for 11 US macroeconomic indicators.     

A Bayesian mixed logit–probit model for multinomial choice

In this paper, we introduce a new flexible mixed model for multinomial discrete choice where the key individual- and alternative-specific parameters of interest are allowed to follow an assumption-free nonparametric density specification, while other alternative-specific coefficients are assumed to be drawn from a multivariate Normal distribution, which eliminates the independence of irrelevant alternatives assumption at the individual level. A hierarchical specification of our model allows us to break down a complex data structure into a set of submodels with the desired features that are naturally assembled in the original system. We estimate the model, using a Bayesian Markov Chain Monte Carlo technique with a multivariate Dirichlet Process (DP) prior on the coefficients with nonparametrically estimated density. We employ a “latent class” sampling algorithm, which is applicable to a general class of models, including non-conjugate DP base priors. The model is applied to supermarket choices of a panel of Houston households whose shopping behavior was observed over a 24-month period in years 2004–2005. We estimate the nonparametric density of two key variables of interest: the price of a basket of goods based on scanner data, and driving distance to the supermarket based on their respective locations. Our semi-parametric approach allows us to identify a complex multi-modal preference distribution, which distinguishes between inframarginal consumers and consumers who strongly value either lower prices or shopping convenience.     

Statistical estimation of multivariate Ornstein–Uhlenbeck processes and applications to co-integration

Ornstein–Uhlenbeck models are continuous-time processes which have broad applications in finance as, e.g., volatility processes in stochastic volatility models or spread models in spread options and pairs trading. The paper presents a least squares estimator for the model parameter in a multivariate Ornstein–Uhlenbeck model driven by a multivariate regularly varying Lévy process with infinite variance. We show that the estimator is consistent. Moreover, we derive its asymptotic behavior and test statistics. The results are compared to the finite variance case. For the proof we require some new results on multivariate regular variation of products of random vectors and central limit theorems. Furthermore, we embed this model in the setup of a co-integrated model in continuous time.     

Bilateral matching with Latin squares

We develop a general procedure to construct pairwise meeting processes characterized by two features. First, in each period the process maximizes the number of matches in the population. Second, over time agents meet everybody else exactly once. We call this type of meetings “absolute strangers.” Our methodological contribution to economics is to offer a simple procedure to construct a type of decentralized trading environments usually employed in both theoretical and experimental economics. In particular, we demonstrate how to make use of the mathematics of Latin squares to enrich the modeling of matching economies.     

Regression discontinuity inference with specification error

A regression discontinuity (RD) research design is appropriate for program evaluation problems in which treatment status (or the probability of treatment) depends on whether an observed covariate exceeds a fixed threshold. In many applications the treatment-determining covariate is discrete. This makes it impossible to compare outcomes for observations “just above” and “just below” the treatment threshold, and requires the researcher to choose a functional form for the relationship between the treatment variable and the outcomes of interest. We propose a simple econometric procedure to account for uncertainty in the choice of functional form for RD designs with discrete support. In particular, we model deviations of the true regression function from a given approximating function—the specification errors—as random. Conventional standard errors ignore the group structure induced by specification errors and tend to overstate the precision of the estimated program impacts. The proposed inference procedure that allows for specification error also has a natural interpretation within a Bayesian framework.     

A robust version of the KPSS test based on indicators

This paper proposes a test of the null hypothesis of stationarity that is robust to the presence of fat-tailed errors. The test statistic is a modified version of the so-called KPSS statistic. The modified statistic uses the “sign” of the data minus the sample median, whereas KPSS used deviations from means. This “indicator” KPSS statistic has the same limit distribution as the standard KPSS statistic under the null, without relying on assumptions about moments, but a different limit distribution under unit root alternatives. The indicator test has lower power than standard KPSS when tails are thin, but higher power when tails are fat.     

Quadratic hedging schemes for non-Gaussian GARCH models

We propose different schemes for option hedging when asset returns are modeled using a general class of GARCH models. More specifically, we implement local risk minimization and a minimum variance hedge approximation based on an extended Girsanov principle that generalizes Duan׳s (1995) delta hedge. Since the minimal martingale measure fails to produce a probability measure in this setting, we construct local risk minimization hedging strategies with respect to a pricing kernel. These approaches are investigated in the context of non-Gaussian driven models. Furthermore, we analyze these methods for non-Gaussian GARCH diffusion limit processes and link them to the corresponding discrete time counterparts. A detailed numerical analysis based on S&P 500 European call options is provided to assess the empirical performance of the proposed schemes. We also test the sensitivity of the hedging strategies with respect to the risk neutral measure used by recomputing some of our results with an exponential affine pricing kernel.     

Processes with volatility-induced stationarity: an application for interest rates

In this paper we propose a refinement of the existing definition of volatility-induced stationarity that allows us to distinguish between processes with drift and diffusion induced stationarity and processes with pure volatility-induced stationarity. We also propose a classification of stationary processes with volatility-induced stationarity according to the volatility that is needed to inject stationarity. Processes with volatility-induced stationarity are potentially applicable to interest rate time-series since, as has been acknowledged, mean-reversion effects occur mainly in periods of high volatility. As such, we provide evidence that the logarithm of the Fed funds rate can be modelled as a local martingale with volatility-induced stationarity.