Lieberman and Phillips (2017; LP) introduced a multivariate stochastic unit root (STUR) model, which allows for random, time varying local departures from a unit root (UR) model, where nonlinear least squares (NLLS) may be used for estimation and inference on the STUR coefficient. In a structural version of this model where the driver variables of the STUR coefficient are endogenous, the NLLS estimate of the STUR parameter is inconsistent, as are the corresponding estimates of the associated covariance parameters. This paper develops a nonlinear instrumental variable (NLIV) as well as GMM estimators of the STUR parameter which conveniently addresses endogeneity. We derive the asymptotic distributions of the NLIV and GMM estimators and establish consistency under similar orthogonality and relevance conditions to those used in the linear model. An overidentification test and its asymptotic distribution are also developed. The results enable inference about structural STUR models and a mechanism for testing the local STUR model against a simple UR null, which complements usual UR tests. Simulations reveal that the asymptotic distributions of the NLIV and GMM estimators of the STUR parameter as well as the test for overidentifying restrictions perform well in small samples and that the distribution of the NLIV estimator is heavily leptokurtic with a limit theory which has Cauchy-like tails. Comparisons of STUR coefficient and standard UR coefficient tests show that the one-sided UR test performs poorly against the one-sided STUR coefficient test both as the sample size and departures from the null rise. The results are applied to study the relationships between stock returns and bond spread changes.

In this paper we study doubly robust estimators of various average and quantile treatment effects under unconfoundedness; we also consider an application to a setting with an instrumental variable. We unify and extend much of the recent literature by providing a very general identification result which covers binary and multi-valued treatments; unnormalized and normalized weighting; and both inverse-probability weighted (IPW) and doubly robust estimators. We also allow for subpopulation-specific average treatment effects where subpopulations can be based on covariate values in an arbitrary way. Similar to Wooldridge (2007), we then discuss estimation of the conditional mean using quasi-log likelihoods (QLL) from the linear exponential family.

The contribution of this paper is two-fold. First, we derive the asymptotic null distribution of the familiar augmented Dickey-Fuller [ADF] statistics in the case where the shocks follow a linear process driven by infinite variance innovations. We show that these distributions are free of serial correlation nuisance parameters but depend on the tail index of the infinite variance process. These distributions are shown to coincide with the corresponding results for the case where the shocks follow a finite autoregression, provided the lag length in the ADF regression satisfies the same o(T1/3) rate condition as is required in the finite variance case. In addition, we establish the rates of consistency and (where they exist) the asymptotic distributions of the ordinary least squares sieve estimates from the ADF regression. Given the dependence of their null distributions on the unknown tail index, our second contribution is to explore sieve wild bootstrap implementations of the ADF tests. Under the assumption of symmetry, we demonstrate the asymptotic validity (bootstrap consistency) of the wild bootstrap ADF tests. This is done by establishing that (conditional on the data) the wild bootstrap ADF statistics attain the same limiting distribution as that of the original ADF statistics taken conditional on the magnitude of the innovations.

Estimation of the intensity of a point process is considered within a nonparametric framework. The intensity measure is unknown and depends on covariates, possibly many more than the observed number of jumps. Only a single trajectory of the counting process is observed. Interest lies in estimating the intensity conditional on the covariates. The impact of the covariates is modelled by an additive model where each component can be written as a linear combination of possibly unknown functions. The focus is on prediction as opposed to variable screening. Conditions are imposed on the coefficients of this linear combination in order to control the estimation error. The rates of convergence are optimal when the number of active covariates is large. As an application, the intensity of the buy and sell trades of the New Zealand Dollar futures is estimated and a test for forecast evaluation is presented. A simulation is included to provide some finite sample intuition on the model and asymptotic properties.

Statistical models of unobserved heterogeneity are typically formalized as mixtures of simple parametric models and interest naturally focuses on testing for homogeneity versus general mixture alternatives. Many tests of this type can be interpreted as C(α) tests, as in Neyman (1959), and shown to be locally asymptotically optimal. These C(α) tests will be contrasted with a new approach to likelihood ratio testing for general mixture models. The latter tests are based on estimation of general nonparametric mixing distribution with the Kiefer and Wolfowitz (1956) maximum likelihood estimator. Recent developments in convex optimization have dramatically improved upon earlier EM methods for computation of these estimators, and recent results on the large sample behavior of likelihood ratios involving such estimators yield a tractable form of asymptotic inference. Improvement in computation efficiency also facilitates the use of a bootstrap method to determine critical values that are shown to work better than the asymptotic critical values in finite samples. Consistency of the bootstrap procedure is also formally established. We compare performance of the two approaches identifying circumstances in which each is preferred.

This article studies two-step sieve M estimation of general semi/nonparametric models, where the second step involves sieve estimation of unknown functions that may use the nonparametric estimates from the first step as inputs, and the parameters of interest are functionals of unknown functions estimated in both steps. We establish the asymptotic normality of the plug-in two-step sieve M estimate of a functional that could be root-n estimable. The asymptotic variance may not have a closed form expression, but can be approximated by a sieve variance that characterizes the effect of the first-step estimation on the second-step estimates. We provide a simple consistent estimate of the sieve variance, thereby facilitating Wald type inferences based on the Gaussian approximation. The finite sample performance of the two-step estimator and the proposed inference procedure are investigated in a simulation study.

This paper considers linear rational expectations models from the linear systems point of view. Using a generalization of the Wiener-Hopf factorization, the linear systems approach is able to furnish very simple conditions for existence and uniqueness of both particular and generic linear rational expectations models. To illustrate the applicability of this approach, the paper characterizes the structure of stationary and cointegrated solutions, including a generalization of Granger’s representation theorem.

We investigate the asymptotic and finite sample properties of a number of methods for estimating the cointegration rank in integrated vector autoregressive systems of unknown autoregressive order driven by heteroskedastic shocks. We allow for both conditional and unconditional heteroskedasticity of a very general form. We establish the conditions required on the penalty functions such that standard information criterion-based methods, such as the Bayesian information criterion [BIC], when employed either sequentially or jointly, can be used to consistently estimate both the cointegration rank and the autoregressive lag order. In doing so we also correct errors which appear in the proofs provided for the consistency of information-based estimators in the homoskedastic case by Aznar and Salvador (2002, Econometric Theory 18, 926–947). We also extend the corpus of available large sample theory for the conventional sequential approach of Johansen (1995, Likelihood-Based Inference in Cointegrated Vector Autoregressive Models. Oxford University Press) and the associated wild bootstrap implementation thereof of Cavaliere, Rahbek, and Taylor (2014, Econometric Reviews 33, 606–650) to the case where the lag order is unknown. In particular, we show that these methods remain valid under heteroskedasticity and an unknown lag length provided the lag length is first chosen by a consistent method, again such as the BIC. The relative finite sample properties of the different methods discussed are investigated in a Monte Carlo simulation study. The two best performing methods in this study are a wild bootstrap implementation of the Johansen (1995, Likelihood-Based Inference in Cointegrated Vector Autoregressive Models. Oxford University Press) procedure implemented with BIC selection of the lag length and joint IC approach (cf. Phillips, 1996, Econometrica 64, 763–812) which uses the BIC to jointly select the lag order and the cointegration rank.

The current article examines the limit distribution of the quasi-maximum likelihood estimator obtained from a directionally differentiable quasi-likelihood function and represents its limit distribution as a functional of a Gaussian stochastic process indexed by direction. In this way, the standard analysis that assumes a differentiable quasi-likelihood function is treated as a special case of our analysis. We also examine and redefine the standard quasi-likelihood ratio, Wald, and Lagrange multiplier test statistics so that their null limit behaviors are regular under our model framework.

This paper studies linear factor models that have arbitrarily dependent factors. Assuming that the coefficients are known and that their matrix representation satisfies rank conditions, we identify the nonparametric joint distribution of the unobserved factors using first and then second-order partial derivatives of the log characteristic function of the observed variables. In conjunction with these identification strategies the mean and variance of the vector of factors are identified. The main result provides necessary and sufficient conditions for identification of the joint distribution of the factors. In an illustrative example, we show identification of an earnings dynamics model with a subset of arbitrarily dependent income shocks. Closed-form formulas lead to estimators that converge uniformly and despite being based on inverse Fourier transforms have tight confidence bands around their theoretical counterparts in Monte Carlo simulations.

We propose a nonparametric estimation theory for the occupation density, the drift vector, and the diffusion matrix of multivariate diffusion processes. The estimators are sample analogues to infinitesimal conditional expectations constructed as Nadaraya-Watson kernel averages. Mild assumptions are imposed on the statistical properties of the multivariate system to obtain limiting results. Harris recurrence is all that we require to show consistency and asymptotic (mixed) normality of the proposed functional estimators. The identification method and asymptotic theory apply to both stationary and nonstationary multivariate diffusion processes of the recurrent type.

Limit theory with stochastic integrals plays a major role in time series econometrics. In earlier contributions on weak convergence to stochastic integrals, the literature commonly uses martingale and semi-martingale structures. Liang, Phillips, Wang, and Wang (2016) (see also Wang (2015), Chap. 4.5) currently extended weak convergence to stochastic integrals by allowing for a linear process or a α-mixing sequence in innovations. While these martingale, linear process and α-mixing structures have wide relevance, they are not sufficiently general to cover many econometric applications that have endogeneity and nonlinearity. This paper provides new conditions for weak convergence to stochastic integrals. Our frameworks allow for long memory processes, causal processes, and near-epoch dependence in innovations, which have applications in a wide range of econometric areas such as TAR, bilinear, and other nonlinear models.

This paper provides sufficient conditions for the nonparametric identification of the regression function $m\left( \cdot \right)$ in a regression model with an endogenous regressor x and an instrumental variable z. It has been shown that the identification of the regression function from the conditional expectation of the dependent variable on the instrument relies on the completeness of the distribution of the endogenous regressor conditional on the instrument, i.e., $f\left( {x|z} \right)$. We show that (1) if the deviation of the conditional density $f\left( {x|{z_k}} \right)$ from a known complete sequence of functions is less than a sequence of values determined by the complete sequence in some distinct sequence $\left\{ {{z_k}:k = 1,2,3, \ldots } \right\}$ converging to ${z_0}$, then $f\left( {x|z} \right)$ itself is complete, and (2) if the conditional density $f\left( {x|z} \right)$ can form a linearly independent sequence $\{ f( \cdot |{z_k}):k = 1,2, \ldots \}$ in some distinct sequence $\left\{ {{z_k}:k = 1,2,3, \ldots } \right\}$ converging to ${z_0}$ and its relative deviation from a known complete sequence of functions under some norm is finite then $f\left( {x|z} \right)$ itself is complete. We use these general results to provide specific sufficient conditions for completeness in three different specifications of the relationship between the endogenous regressor x and the instrumental variable $z.$

The paper studies spatial autoregressive models with group interaction structure, focussing on estimation and inference for the spatial autoregressive parameter λ. The quasi-maximum likelihood estimator for λ usually cannot be written in closed form, but using an exact result obtained earlier by the authors for its distribution function, we are able to provide a complete analysis of the properties of the estimator, and exact inference that can be based on it, in models that are balanced. This is presented first for the so-called pure model, with no regression component, but is also extended to some special cases of the more general model. We then study the much more difficult case of unbalanced models, giving analogues of some, but by no means all, of the results obtained for the balanced case earlier. In both balanced and unbalanced models, results obtained for the pure model generalize immediately to the model with group-specific regression components.

Proxies for unobserved skills and technologies are increasingly available in empirical data. For dynamic discrete choice models of forward-looking agents where a continuous state variable is unobserved but its proxy is available, we derive closed-form identification of the structure by explicitly solving integral equations. In the first step, we derive closed-form identification of Markov components, including the conditional choice probabilities and the law of state transition. In the second step, we plug in these first-step identifying formulas to obtain primitive structural parameters of dynamically optimizing agents.

We consider the method of moments estimation of a structural equation in a panel dynamic simultaneous equations model under different sample size combinations of cross-sectional dimension, N, and time series dimension, T. Two types of linear transformation to remove the individual-specific effects from the model, first difference and forward orthogonal demeaning, are considered. We show that the Alvarez and Arellano (2003) type GMM estimator under both transformations is consistent only if ${T \over N} \to 0$ as $\left( {N,T} \right) \to \infty $. However, it is asymptotically biased if ${{{T^3}} \over N} \to \kappa \ne 0 < \infty$ as $\left( {N,T} \right) \to \infty $. Since the validity of statistical inference depends critically on whether an estimator is asymptotically unbiased, we suggest a jackknife bias reduction method and derive its limiting distribution. Monte Carlo studies are conducted to demonstrate the importance of using an asymptotically unbiased estimator to obtain valid statistical inference.

Many economic models yield conditional moment inequalities that can be used for inference on parameters of these models. In this paper, I construct new tests of parameter hypotheses in conditional moment inequality models based on studentized kernel estimates of moment functions. The tests automatically adapt to the unknown smoothness of the moment functions, have uniformly correct asymptotic size, and are rate-optimal against certain classes of alternatives. Some existing tests have nontrivial power against n−1/2-local alternatives of a certain type whereas my methods only allow for nontrivial testing against (n / log n)−1/2-local alternatives of this type. There exist, however, large classes of sequences of well-behaved alternatives against which the tests developed in this paper are consistent and those tests are not.

Barndorff-Nielsen’s celebrated p*-formula and variations thereof have amongst their various attractions the ability to approximate bimodal distributions. In this paper we show that in general this requires a crucial adjustment to the basic formula. The adjustment is based on a simple idea and straightforward to implement, yet delivers important improvements. It is based on recognizing that certain outcomes are theoretically impossible and the density of the MLE should then equal zero, rather than the positive density that a straight application of p* would suggest. This has implications for inference and we show how to use the new p**-formula to construct improved confidence regions. These can be disjoint as a consequence of the bimodality. The degree of bimodality depends heavily on the value of an approximate ancillary statistic and conditioning on the observed value of this statistic is therefore desirable. The p**-formula naturally delivers the relevant conditional distribution. We illustrate these results in small and large samples using a simple nonlinear regression model and errors in variables model where the measurement errors in dependent and explanatory variables are correlated and allow for weak proxies.

This paper studies the weak convergence of renorming volatilities in a family of GARCH(1,1) models from a functional point of view. After suitable renormalization, it is shown that the limiting distribution is a geometric Brownian motion when the associated top Lyapunov exponent γ > 0 and is an exponential functional of the maximum process of a Brownian motion when γ = 0. This indicates that the volatility of the GARCH(1,1)-type model has a completely different random structure according to the sign of γ. The obtained results further strengthen our understanding of volatilities in GARCH-type models. Simulation studies are carried out to assess our findings.