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.

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.

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.

In this article, we investigate the properties of heteroskedasticity and autocorrelation robust (HAR) test statistics in time series regression settings when observations are missing. We primarily focus on the nonrandom missing process case where we treat the missing locations to be fixed as T → ∞ by mapping the missing and observed cutoff dates into points on [0,1] based on the proportion of time periods in the sample that occur up to those cutoff dates. We consider two models, the amplitude modulated series (Parzen, 1963) regression model, which amounts to plugging in zeros for missing observations, and the equal space regression model, which simply ignores the missing observations. When the amplitude modulated series regression model is used, the fixed-b limits of the HAR test statistics depend on the locations of missing observations but are otherwise pivotal. When the equal space regression model is used, the fixed-b limits of the HAR test statistics have the standard fixed-b limits as in Kiefer and Vogelsang (2005). We discuss methods for obtaining fixed-b critical values with a focus on bootstrap methods and find the naive i.i.d. bootstrap with missing dates fixed to be an effective and practical way to obtain the fixed-b critical values.

Let x be a transformation of y, whose distribution is unknown. We derive an expansion formulating the expectations of x in terms of the expectations of y. Apart from the intrinsic interest in such a fundamental relation, our results can be applied to calculating E(x) by the low-order moments of a transformation which can be chosen to give a good approximation for E(x). To do so, we generalize the approach of bounding the terms in expansions of characteristic functions, and use our result to derive an explicit and accurate bound for the remainder when a finite number of terms is taken. We illustrate one of the implications of our method by providing accurate naive bootstrap confidence intervals for the mean of any fat-tailed distribution with an infinite variance, in which case currently available bootstrap methods are asymptotically invalid or unreliable in finite samples.

We show boundedness in probability uniformly in sample size of a general M-estimator for multiple linear regression in time series. The positive criterion function for the M-estimator is assumed lower semicontinuous and sufficiently large for large argument. Particular cases are the Huber-skip and quantile regression. Boundedness requires an assumption on the frequency of small regressors. We show that this is satisfied for a variety of deterministic and stochastic regressors, including stationary and random walks regressors. The results are obtained using a detailed analysis of the condition on the regressors combined with some recent martingale results.