Forecasting crude oil prices with DSGE models

In this study, we conducted an oil prices forecasting competition among a set of structural models, including vector autoregression and dynamic stochastic general equilibrium (DSGE) models. Our results highlight two principles. First, forecasts should exploit the fact that real oil prices are mean reverting over long horizons. Second, models should not replicate the high volatility of the oil prices observed in samples. By following these principles, we show that an oil sector DSGE model performs much better at real oil price forecasting than random walk or vector autoregression.     

Forecasting UK GDP growth and inflation under structural change. A comparison of models with time-varying parameters

Evidence from a large and growing body of empirical literature strongly suggests that there have been changes in the inflation and output dynamics in the United Kingdom. The majority of these papers base their results on a class of econometric models that allows for time-variation in the coefficients and volatilities of shocks. While these models have been used extensively for studying evolving dynamics and for structural analysis, there has been little evidence that they are useful for forecasting UK output growth and inflation. This paper attempts to fill this gap by comparing the performances of a wide range of time-varying parameter models in forecasting output growth and inflation. We find that allowing for time-varying parameters can lead to large and statistically significant gains in forecast accuracy.     

Bayesian prior elicitation in DSGE models: Macro- vs micropriors

Bayesian approaches to the estimation of DSGE models are becoming increasingly popular. Prior knowledge is normally formalized either directly on deep parameters' values (‘microprior’) or indirectly, on macroeconomic indicators, e.g. moments of observable variables (‘macroprior’). We introduce a non-parametric macroprior which is elicited from impulse response functions and assess its performance in shaping posterior estimates. We find that using a macroprior can lead to substantially different posterior estimates. We probe into the details of our result, showing that model misspecification is likely to be responsible of that. In addition, we assess to what extent the use of macropriors is impaired by the need of calibrating some hyperparameters.     

Non-linear DSGE models and the optimized central difference particle filter

We improve the accuracy and speed of particle filtering for non-linear DSGE models with potentially non-normal shocks. This is done by introducing a new proposal distribution which (i) incorporates information from new observables and (ii) has a small optimization step that minimizes the distance to the optimal proposal distribution. A particle filter with this proposal distribution is shown to deliver a high level of accuracy even with relatively few particles, and it is therefore much more efficient than the standard particle filter.     

Second-order approximation of dynamic models without the use of tensors

Several approaches to finding the second-order approximation to a dynamic model have been proposed recently. This paper differs from the existing literature in that it makes use of the Magnus and Neudecker (1999) definition of the Hessian matrix. The key result is a linear system of equations that characterizes the second-order coefficients. No use is made of multi-dimensional arrays or tensors, a practical implication of which is that it is much easier to transcribe the mathematical representation of the solution into usable computer code. Matlab code is available from http://paulklein.se/newsite/codes/codes.php; Fortran 90 code is available from http://alcor.concordia.ca/∼pgomme/.     

Bayesian inference for nonlinear structural time series models

We consider efficient methods for likelihood inference applied to structural models. In particular, we introduce a particle filter method which concentrates upon disturbances in the Markov state of the approximating solution to the structural model. A particular feature of such models is that the conditional distribution of interest for the disturbances is often multimodal. We provide a fast and effective method for approximating such distributions. We estimate a neoclassical growth model using this approach. An asset pricing model with persistent habits is also considered. The methodology we employ allows many fewer particles to be used than alternative procedures for a given precision.     

Quasi-maximum likelihood estimation and testing for nonlinear models with endogenous explanatory variables

I propose a quasi-maximum likelihood framework for estimating nonlinear models with continuous or discrete endogenous explanatory variables. Joint and two-step estimation procedures are considered. The joint procedure is a quasi-limited information maximum likelihood procedure, as one or both of the log likelihoods may be misspecified. The two-step control function approach is computationally simple and leads to straightforward tests of endogeneity. In the case of discrete endogenous explanatory variables, I argue that the control function approach can be applied with generalized residuals to obtain average partial effects. I show how the results apply to nonlinear models for fractional and nonnegative responses.     

Tailored randomized block MCMC methods with application to DSGE models

In this paper we develop new Markov chain Monte Carlo schemes for the estimation of Bayesian models. One key feature of our method, which we call the tailored randomized block Metropolis–Hastings (TaRB-MH) method, is the random clustering of the parameters at every iteration into an arbitrary number of blocks. Then each block is sequentially updated through an M–H step. Another feature is that the proposal density for each block is tailored to the location and curvature of the target density based on the output of simulated annealing, following  and  and Chib and Ergashev (in press). We also provide an extended version of our method for sampling multi-modal distributions in which at a pre-specified mode jumping iteration, a single-block proposal is generated from one of the modal regions using a mixture proposal density, and this proposal is then accepted according to an M–H probability of move. At the non-mode jumping iterations, the draws are obtained by applying the TaRB-MH algorithm. We also discuss how the approaches of Chib (1995) and Chib and Jeliazkov (2001) can be adapted to these sampling schemes for estimating the model marginal likelihood. The methods are illustrated in several problems. In the DSGE model of Smets and Wouters (2007), for example, which involves a 36-dimensional posterior distribution, we show that the autocorrelations of the sampled draws from the TaRB-MH algorithm decay to zero within 30–40 lags for most parameters. In contrast, the sampled draws from the random-walk M–H method, the algorithm that has been used to date in the context of DSGE models, exhibit significant autocorrelations even at lags 2500 and beyond. Additionally, the RW-MH does not explore the same high density regions of the posterior distribution as the TaRB-MH algorithm. Another example concerns the model of An and Schorfheide (2007) where the posterior distribution is multi-modal. While the RW-MH algorithm is unable to jump from the low modal region to the high modal region, and vice-versa, we show that the extended TaRB-MH method explores the posterior distribution globally in an efficient manner.     

Forecasting sales response for multiple time horizons and temporally aggregated data: A comparison of constant and stochastic coefficient models

Past research on time-varying sales-response models emphasized the application of different estimation techniques in examining variation in advertising effectiveness over time. This study focuses on comparing sales forecasts using constant and stochastic coefficients sales-response models. Selected constant and stochastic coefficient models are applied to six sets of bimonthly and one set of annual advertising and sales data to assess forecasting accuracy for time horizons of various lengths. Results show improved forecasting accuracy for a first-order autoregressive stochastic coefficient model, particularly in short-run forecasting applications.     

Kuhn-Tucker,likelihood ratio and Wald tests for nonlinear models with inequality constraints on the parameters

This paper considers the problem of testing statistical hypothesis in nonlinear regression models with inequality constraints on the parameters. First, the Kuhn-Tucker test procedure is defined. Next, it is shown that the distribution of the Kuhn-Tucker, the likelihood ratio and the Wald test statistics converges to the same mixture of chi-square distributions under the null hypothesis. To illustrate these results two examples are considered: (1) the problem of testing that individual effects are missing in an error component model, and (2) the problem of testing equilibrium for a model of markets in disequilibrium.     

Maximum likelihood estimation for dynamic factor models with missing data

This paper concerns estimating parameters in a high-dimensional dynamic factor model by the method of maximum likelihood. To accommodate missing data in the analysis, we propose a new model representation for the dynamic factor model. It allows the Kalman filter and related smoothing methods to evaluate the likelihood function and to produce optimal factor estimates in a computationally efficient way when missing data is present. The implementation details of our methods for signal extraction and maximum likelihood estimation are discussed. The computational gains of the new devices are presented based on simulated data sets with varying numbers of missing entries.     

Testing for unit roots in time series models with non-stationary volatility

Many of the key macro-economic and financial variables in developed economies are characterized by permanent volatility shifts. It is known that conventional unit root tests are potentially unreliable in the presence of such behaviour, depending on a particular function (the variance profile) of the underlying volatility process. Somewhat surprisingly then, very little work has been undertaken to develop unit root tests which are robust to the presence of permanent volatility shifts. In this paper we fill this gap in the literature by proposing tests which are valid in the presence of a quite general class of permanent variance changes which includes single and multiple (abrupt and smooth-transition) volatility change processes as special cases. Our solution uses numerical methods to simulate the asymptotic null distribution of the statistics based on a consistent estimate of the variance profile which we also develop. The practitioner is not required to specify a parametric model for volatility. An empirical illustration using producer price inflation series from the Stock–Watson database is reported.     

Density estimation for nonlinear parametric models with conditional heteroscedasticity

This article studies density and parameter estimation problems for nonlinear parametric models with conditional heteroscedasticity. We propose a simple density estimate that is particularly useful for studying the stationary density of nonlinear time series models. Under a general dependence structure, we establish the root n$n$ consistency of the proposed density estimate. For parameter estimation, a Bahadur type representation is obtained for the conditional maximum likelihood estimate. The parameter estimate is shown to be asymptotically efficient in the sense that its limiting variance attains the Cramér–Rao lower bound. The performance of our density estimate is studied by simulations.     

Exact and asymptotic tests for possibly non-regular hypotheses on stochastic volatility models

We study the problem of testing hypotheses on the parameters of one- and two-factor stochastic volatility models (SV), allowing for the possible presence of non-regularities such as singular moment conditions and unidentified parameters, which can lead to non-standard asymptotic distributions. We focus on the development of simulation-based exact procedures–whose level can be controlled in finite samples–as well as on large-sample procedures which remain valid under non-regular conditions. We consider Wald-type, score-type and likelihood-ratio-type tests based on a simple moment estimator, which can be easily simulated. We also propose a C(α)-type test which is very easy to implement and exhibits relatively good size and power properties. Besides usual linear restrictions on the SV model coefficients, the problems studied include testing homoskedasticity against a SV alternative (which involves singular moment conditions under the null hypothesis) and testing the null hypothesis of one factor driving the dynamics of the volatility process against two factors (which raises identification difficulties). Three ways of implementing the tests based on alternative statistics are compared: asymptotic critical values (when available), a local Monte Carlo (or parametric bootstrap) test procedure, and a maximized Monte Carlo (MMC) procedure. The size and power properties of the proposed tests are examined in a simulation experiment. The results indicate that the C(α)-based tests (built upon the simple moment estimator available in closed form) have good size and power properties for regular hypotheses, while Monte Carlo tests are much more reliable than those based on asymptotic critical values. Further, in cases where the parametric bootstrap appears to fail (for example, in the presence of identification problems), the MMC procedure easily controls the level of the tests. Moreover, MMC-based tests exhibit relatively good power performance despite the conservative feature of the procedure. Finally, we present an application to a time series of returns on the Standard and Poor’s Composite Price Index.     

Efficient estimation and computation for the generalised additive models with unknown link function

The generalised additive models (GAM) are widely used in data analysis. In the application of the GAM, the link function involved is usually assumed to be a commonly used one without justification. Motivated by a real data example with binary response where the commonly used link function does not work, we propose a generalised additive models with unknown link function (GAMUL) for various types of data, including binary, continuous and ordinal. The proposed estimators are proved to be consistent and asymptotically normal. Semiparametric efficiency of the estimators is demonstrated in terms of their linear functionals. In addition, an iterative algorithm, where all estimators can be expressed explicitly as a linear function of $Y$, is proposed to overcome the computational hurdle for the GAM type model. Extensive simulation studies conducted in this paper show the proposed estimation procedure works very well. The proposed GAMUL are finally used to analyze a real dataset about loan repayment in China, which leads to some interesting findings.     

Nonparametric model validations for hidden Markov models with applications in financial econometrics

We address the nonparametric model validation problem for hidden Markov models with partially observable variables and hidden states. We achieve this goal by constructing a nonparametric simultaneous confidence envelope for transition density function of the observable variables and checking whether the parametric density estimate is contained within such an envelope. Our specification test procedure is motivated by a functional connection between the transition density of the observable variables and the Markov transition kernel of the hidden states. Our approach is applicable for continuous-time diffusion models, stochastic volatility models, nonlinear time series models, and models with market microstructure noise.     

Efficiencies for optimum designs when transforming the response in nonlinear models with nonconstant variance

Often the responses from mechanistic models have to be transformed to achieve error distributions that are symmetric and have constant variance. Because of the nature of the relationship between the response and the mechanistic model, it is necessary to transform both sides of the model. Expressions are given for the parameter sensitivities in the transformed model and examples given of optimum designs for particular values of λ, together with the efficiency of these designs as λ varies. Approaches to finding designs robust to variations in λ are indicated and exemplified.     

Convergence rates of recombining trees for pricing options on stocks under GBM and regime-switching models with known cash dividends

In the literature there appear various kinds of binomial trees for pricing options on stocks under geometric Brownian motions (GBMs) with known cash dividends. The aim of this paper is to compare the performance of the existing binomial trees in aspect of the convergence rates, which are usually used to measure precisely how fast the approximate values converge to the exact one, and to give a theoretical proof of the convergence rates for the interpolation binomial trees which are based on a model that excludes the arbitrage possibilities. Also the paper extends the studies to the regime-switching models with known cash dividend payment.     

Experiments with mixtures: Optimal allocations for becker’s models

Extending Scheffé’s simplex-centroid design for experiments with mixtures, we introduce aweighted simplex-centroid design for a class of mixture models. Becker’s homogeneous functions of degree one belong to this class. By applying optimal design theory, we obtainA-, D- andI-optimal allocations of observations for Becker’s models.