We consider inference about coefficients on a small number of variables of interest in a linear panel data model with additive unobserved individual and time specific effects and a large number of additional time-varying confounding variables. We suppose that, in addition to unrestricted time and individual specific effects, these confounding variables are generated by a small number of common factors and high-dimensional weakly dependent disturbances. We allow that both the factors and the disturbances are related to the outcome variable and other variables of interest. To make informative inference feasible, we impose that the contribution of the part of the confounding variables not captured by time specific effects, individual specific effects, or the common factors can be captured by a relatively small number of terms whose identities are unknown. Within this framework, we provide a convenient inferential procedure based on factor extraction followed by lasso regression and show that the procedure has good asymptotic properties. We also provide a simple k-step bootstrap procedure that may be used to construct inferential statements about the low-dimensional parameters of interest and prove its asymptotic validity. We provide simulation evidence about the performance of our procedure and illustrate its use in an empirical application.

Estimating a treatment’s effect on an outcome conditional on covariates is a primary goal of many empirical investigations. Accurate estimation of the treatment effect given covariates can enable the optimal treatment to be applied to each unit or guide the deployment of limited treatment resources for maximum program benefit. Applications of conditional treatment effect estimation are found in direct marketing, economic policy, and personalized medicine. When estimating conditional treatment effects, the typical practice is to select a statistical model or procedure based on sample data. However, combining estimates from the candidate procedures often provides a more accurate estimate than the selection of a single procedure. This article proposes a method of model combination that targets accurate estimation of the treatment effect conditional on covariates. We provide a risk bound for the resulting estimator under squared error loss and illustrate the method using data from a labor skills training program.

We develop parametric inference procedures for large panels of noisy option data in a setting, where the underlying process is of pure-jump type, i.e., evolves only through a sequence of jumps. The panel consists of options written on the underlying asset with a (different) set of strikes and maturities available across the observation times. We consider an asymptotic setting in which the cross-sectional dimension of the panel increases to infinity, while the time span remains fixed. The information set is augmented with high-frequency data on the underlying asset. Given a parametric specification for the risk-neutral asset return dynamics, the option prices are nonlinear functions of a time-invariant parameter vector and a time-varying latent state vector (or factors). Furthermore, no-arbitrage restrictions impose a direct link between some of the quantities that may be identified from the return and option data. These include the so-called jump activity index as well as the time-varying jump intensity. We propose penalized least squares estimation in which we minimize the L2 distance between observed and model-implied options. In addition, we penalize for the deviation of the model-implied quantities from their model-free counterparts, obtained from the high-frequency returns. We derive the joint asymptotic distribution of the parameters, factor realizations and high-frequency measures, which is mixed Gaussian. The different components of the parameter and state vector exhibit different rates of convergence, depending on the relative (asymptotic) informativeness of the high-frequency return data and the option panel.

This article investigates the statistical inference problem of whether a measurement equation is self-consistent in the logarithmic realized GARCH model (log-RealGARCH). First, we provide the sufficient and necessary conditions for the strict stationarity of both the log-RealGARCH model and the log-GARCH-X model. Under these conditions, strong consistency and asymptotic normality of the quasi-maximum likelihood estimators of these two models are obtained. Then, based on the asymptotic results, we propose a Hausman-type self-consistency test for diagnosing the suitability of the measurement equation in the log-RealGARCH model. Finally, the results of simulations and an empirical study are found to accord with the theoretical results.

The asymptotic distribution of the Gaussian quasi-maximum likelihood estimator (QMLE) is obtained for a wide class of asymmetric GARCH models with exogenous covariates. The true value of the parameter is not restricted to belong to the interior of the parameter space, which allows us to derive tests for the significance of the parameters. In particular, the relevance of the exogenous variables can be assessed. The results are obtained without assuming that the innovations are independent, which allows conditioning on different information sets. Monte Carlo experiments and applications to financial series illustrate the asymptotic results. In particular, an empirical study demonstrates that the realized volatility can be a helpful covariate for predicting squared returns.

Monotonicity is a key qualitative prediction of a wide array of economic models derived via robust comparative statics. It is therefore important to design effective and practical econometric methods for testing this prediction in empirical analysis. This article develops a general nonparametric framework for testing monotonicity of a regression function. Using this framework, a broad class of new tests is introduced, which gives an empirical researcher a lot of flexibility to incorporate ex ante information she might have. The article also develops new methods for simulating critical values, which are based on the combination of a bootstrap procedure and new selection algorithms. These methods yield tests that have correct asymptotic size and are asymptotically nonconservative. It is also shown how to obtain an adaptive and rate optimal test that has the best attainable rate of uniform consistency against models whose regression function has Lipschitz-continuous first-order derivatives and that automatically adapts to the unknown smoothness of the regression function. Simulations show that the power of the new tests in many cases significantly exceeds that of some prior tests, e.g., that of Ghosal, Sen, and Van der Vaart (2000).

We provide novel characterizations of multivariate normality that incorporate both the characteristic function and the moment generating function, and we employ these results to construct a class of affine invariant, consistent and easy-to-use goodness-of-fit tests for normality. The test statistics are suitably weighted L2-statistics, and we provide their asymptotic behavior both for i.i.d. observations as well as in the context of testing that the innovation distribution of a multivariate GARCH model is Gaussian. We also study the finite-sample behavior of the new tests and compare the new criteria with alternative existing tests.

This article considers testing for independence in a time series of small counts within an Integer Autoregressive (INAR) model, taking a semiparametric approach that avoids any distributional assumption on the arrivals process of the model. The nature of the testing problem is shown to differ depending on whether or not the support of the arrivals distribution is the full set of natural numbers (as would be the case for Poisson or Negative Binomial distributions for example) or some strict subset of the natural numbers (such as for a Binomial or Uniform distribution). The theory for these two cases is studied separately.

For the case where the arrivals have support on the natural numbers, a new asymptotically efficient semiparametric test, the effective score (Neyman-Rao) test, is derived. The semiparametric Likelihood-Ratio, Wald and score tests are shown to be asymptotically equivalent to the effective score test, and hence also asymptotically efficient. Asymptotic relative efficiency calculations demonstrate that the semiparametric effective score test can provide substantial power advantages over the first order autocorrelation coefficient, which is most commonly applied in practice.

For the case where the arrivals have support that is a strict subset of the natural numbers, the theory is considerably altered because the support of the observations becomes different under the null and alternative hypotheses. The semiparametric Likelihood-Ratio, Wald and score tests become asymptotically degenerate in this case, while the effective score test remains valid. Remarkably, in this case the effective score test is also found to have power against local alternatives that shrink to the null at the rate T−1. In rare cases where the arrival support is partly or totally known, additional tests exploiting this information are considered.

Finite sample properties of the tests in these various cases demonstrate the semiparametric effective score test can provide substantial power advantages over the first order autocorrelation test implied by a parametric Poisson specification. The simulations also reveal situations in which the first order autocorrelation is preferable in finite samples, so a hybrid of the effective score and autocorrelation tests is proposed to capture most of the benefits of each test.

We establish oracle inequalities for a version of the Lasso in high-dimensional fixed effects dynamic panel data models. The inequalities are valid for the coefficients of the dynamic and exogenous regressors. Separate oracle inequalities are derived for the fixed effects. Next, we show how one can conduct uniformly valid inference on the parameters of the model and construct a uniformly valid estimator of the asymptotic covariance matrix which is robust to conditional heteroskedasticity in the error terms. Allowing for conditional heteroskedasticity is important in dynamic models as the conditional error variance may be nonconstant over time and depend on the covariates. Furthermore, our procedure allows for inference on high-dimensional subsets of the parameter vector of an increasing cardinality. We show that the confidence bands resulting from our procedure are asymptotically honest and contract at the optimal rate. This rate is different for the fixed effects than for the remaining parts of the parameter vector.

The present article discusses panel unit root tests for two classes of models. One is the moving average (MA) model and the other is the error components model. We conduct score type tests for these models, allowing for various types of regressors, and examine the cross-sectional effect explicitly by presenting an efficient way of computing limiting local powers. It is found that the existence of common regressors does not affect the asymptotic behavior of tests, although heterogeneous regressors do. We also derive the limiting power envelopes for some simple panel models, which shows that the score type panel tests are asymptotically efficient, unlike the time series case.

We consider a semiparametric transformation model, in which the regression function has an additive nonparametric structure and the transformation of the response is assumed to belong to some parametric family. We suppose that endogeneity is present in the explanatory variables. Using a control function approach, we show that the proposed model is identified under suitable assumptions, and propose a profile estimation method for the transformation. The proposed estimator is shown to be asymptotically normal under certain regularity conditions. A simulation study shows that the estimator behaves well in practice. Finally, we give an empirical example using the U.K. Family Expenditure Survey.

This article introduces a local Gaussian bootstrap method useful for the estimation of the asymptotic distribution of high-frequency data-based statistics such as functions of realized multivariate volatility measures as well as their asymptotic variances. The new approach consists of dividing the original data into nonoverlapping blocks of M consecutive returns sampled at frequency h (where h−1 denotes the sample size) and then generating the bootstrap observations at each frequency within a block by drawing them randomly from a mean zero Gaussian distribution with a variance given by the realized variance computed over the corresponding block.

Our main contributions are as follows. First, we show that the local Gaussian bootstrap is first-order consistent when used to estimate the distributions of realized volatility and realized betas under assumptions on the log-price process which follows a continuous Brownian semimartingale process. Second, we show that the local Gaussian bootstrap matches accurately the first four cumulants of realized volatility up to o(h), implying that this method provides third-order refinements. This is in contrast with the wild bootstrap of Gonçalves and Meddahi (2009, Econometrica 77(1), 283–306), which is only second-order correct. Third, we show that the local Gaussian bootstrap is able to provide second-order refinements for the realized beta, which is also an improvement of the existing bootstrap results in Dovonon, Gonçalves, and Meddahi (2013, Journal of Econometrics 172, 49–65) (where the pairs bootstrap was shown not to be second-order correct under general stochastic volatility). In addition, we highlight the connection between the local Gaussian bootstrap and the local Gaussianity approximation of continuous semimartingales established by Mykland and Zhang (2009, Econometrica 77, 1403–1455) and show the suitability of this bootstrap method to deal with the new class of estimators introduced in that article. Lastly, we provide Monte Carlo simulations and use empirical data to compare the finite sample accuracy of our new bootstrap confidence intervals for integrated volatility with the existing results.

We examine a test for weak stationarity against alternatives that covers both local-stationarity and break point models. A key feature of the test is that its asymptotic distribution is a functional of the standard Brownian bridge sheet in [0,1]2, so that it does not depend on any unknown quantity. The test has nontrivial power against local alternatives converging to the null hypothesis at a T−1/2 rate, where T is the sample size. We also examine an easy-to-implement bootstrap analogue and present the finite sample performance in a Monte Carlo experiment. Finally, we implement the methodology to assess the stability of inflation dynamics in the United States and on a set of neuroscience tremor data.

We propose a test for a generalized regression monotonicity (GRM) hypothesis. The GRM hypothesis is the sharp testable implication of the monotonicity of certain latent structures, as we show in this article. Examples include the monotonicity of the conditional mean function when only interval data are available for the dependent variable and the monotone instrumental variable assumption of Manski and Pepper (2000). These instances of latent monotonicity can be tested using our test. Moreover, the GRM hypothesis includes regression monotonicity and stochastic monotonicity as special cases. Thus, our test also serves as an alternative to existing tests for those hypotheses. We show that our test controls the size uniformly over a broad set of data generating processes asymptotically, is consistent against fixed alternatives, and has nontrivial power against some ${n^{ - 1/2}}$ local alternatives.

A mixture copula is a linear combination of several individual copulas that can be used to generate dependence structures not belonging to existing copula families. Because different pairs of markets may exhibit quite different dependence structures in empirical studies, mixture copulas are useful in modeling the dependence in financial data. Therefore, rather than selecting a single copula based on certain criteria, we propose using a model averaging approach to estimate financial data dependence structures in a mixture copula framework. We select weights (for averaging) by a J-fold Cross-Validation procedure. We prove that the model averaging estimator is asymptotically optimal in the sense that it minimizes the squared estimation loss. Our simulation results show that the model averaging approach outperforms some competing methods when the working mixture model is misspecified. Using 12 years of data on daily returns from four developed economies’ stock indexes, we show that the model averaging approach more accurately estimates their dependence structures than some competing methods.

We examine a higher-order spatial autoregressive model with stochastic, but exogenous, spatial weight matrices. Allowing a general spatial linear process form for the disturbances that permits many common types of error specifications as well as potential ‘long memory’, we provide sufficient conditions for consistency and asymptotic normality of instrumental variables, ordinary least squares, and pseudo maximum likelihood estimates. The implications of popular weight matrix normalizations and structures for our theoretical conditions are discussed. A set of Monte Carlo simulations examines the behaviour of the estimates in a variety of situations. Our results are especially pertinent in situations where spatial weights are functions of stochastic economic variables, and this type of setting is also studied in our simulations.

We present novel theory for testing for reduction of GARCH-X type models with an exogenous (X) covariate to standard GARCH type models. To deal with the problems of potential nuisance parameters on the boundary of the parameter space as well as lack of identification under the null, we exploit a noticeable property of specific zero-entries in the inverse information of the GARCH-X type models. Specifically, we consider sequential testing based on two likelihood ratio tests and as demonstrated the structure of the inverse information implies that the proposed test neither depends on whether the nuisance parameters lie on the boundary of the parameter space, nor on lack of identification. Asymptotic theory is derived essentially under stationarity and ergodicity, coupled with a regularity assumption on the exogenous covariate X. Our general results on GARCH-X type models are applied to Gaussian based GARCH-X models, GARCH-X models with Student’s t-distributed innovations as well as integer-valued GARCH-X (PAR-X) models.