
A number of economic, financial, and climatic time series exhibit persistent cycles which are characterized by time-dependence patterns and peaks in the spectrum. In this article, we introduce a class of semiparametric cyclical–memory processes which enable the modeling of random cyclical patterns in stationary and non-stationary time series. We develop a theoretical background and asymptotic estimation theory for the frequency of a cycle represented by the location of a peak in the spectrum. The estimation procedure is easy to implement and allows for the construction of narrow confidence intervals around the location point. Monte Carlo simulations confirm the good finite sample performance of our estimator. We illustrate our method with three empirical applications. We uncover (quasi-)periodic cycles in macroeconomic series, both nominal and real (U.S. nominal GDP and real industrial production), and CO 2 $\mathrm {CO}_{2}$ upper C upper O 2 concentration levels.
We explore how information from covariates can be incorporated into the CUSUM-based real-time monitoring procedure for explosive asset price bubbles developed in Homm and Breitung (2012, Journal of Financial Econometrics 10, 198-231). Where dynamic covariates are present in the data generating process (DGP), the false positive rate (FPR) of the basic CUSUM procedure, which is based on the assumption that prices follow a univariate DGP, under the null of no explosivity will not, in general, be properly controlled, even asymptotically. In contrast, accounting for these relevant covariates in the construction of the CUSUM statistics leads to a procedure whose FPR can be controlled using the same asymptotic crossing function as employed by Homm and Breitung (2012). Doing so is also shown to have the potential to significantly increase the chance of detecting an emerging bubble episode in finite samples. We additionally allow for time-varying volatility in the innovations driving the model through the use of a kernel-based variance estimator.
Linear factor models dominate the field of empirical asset pricing and are largely considered a successful framework even though testing for model adequacy is uncommon in empirical research. This article develops a probabilistic argument under the assumption that return data are generated by a multivariate L & eacute;vy process and shows that linear factor models are asymptotically adequate as the return frequency declines. Rate-of-convergence results are provided assuming finite or infinite second moments and the moment conditions that determine the convergence rate are analyzed. Separately, we consider the combination of cross-sectional and temporal aggregation, which concerns portfolio applications, and we find that the model error of the linear framework can be smaller for portfolios than individual stocks. Although the intended application is finance, the results in the article are valid in other areas where data aggregation is natural.
We consider the problem of estimating and deriving confidence intervals for change points in linear models with heteroscedastic errors. A CUSUM process-based estimator is proposed, and we establish its asymptotic properties when the linear regression model exhibits change points in both the regression parameters and the distribution of the errors. This theory motivates the construction of confidence sets for multiple change points by refining preliminary change point estimators and approximating their distribution in a way that is robust to heteroscedasticity. Monte Carlo experiments indicate that the proposed confidence intervals achieve accurate empirical coverage for change-point locations under both homoscedastic and heteroscedastic error structures. In two data applications, we apply the proposed confidence intervals to examine changes in the flattening of the New Keynesian Phillips curve and in cryptocurrency risk factors.
We propose a high-dimensional extension of the heteroscedasticity test proposed in Newey and Powell (1987). Our test is based on expectile regression in the proportional asymptotic regime where n / p → δ ∈ ( 0 , 1 ] . The asymptotic analysis of the test statistic uses the approximate message passing algorithm, from which we obtain the limiting distribution of the test and establish its asymptotic power. The numerical performance of the test is validated through an extensive simulation study. As real-data applications, we present the analysis based on "international economic growth" data (Belloni et al., 2013), which is found to be homoscedastic, and "supermarket" data (Lan et al., 2016), which is found to be heteroscedastic.
Marco Lippi was born in Rome in 1943. An indefatigable and inspiring pedagogue, he has been teaching mathematics, economics, the history of economic thought, and econometrics to generations of students at the Universities of Perugia, Rome (La Sapienza, Tor Vergata, and LUISS), Modena, the Scuola Superiore Sant'Anna in Pisa, and the European Center for Advanced Research in Economics and Statistics (ECARES) in Brussels. As a fellow of the Einaudi Institute for Economics and Finance (EIEF), he still teaches, with the indomitable enthusiasm that has become legendary among his students and colleagues, Master and Ph.D. courses offered by this renowned Roman institution.
This article proposes sequential randomized tests to locate the presence of jumps on the paths of efficient asset prices in a continuous-time model. The randomized statistics are generated by artificially adding randomness to the robust approximations of the locally averaged returns of the efficient price. In the case of finite activity jumps, we derive the asymptotic distribution of the maximum of all the local statistics unaffected by jumps, which makes it feasible to control the limiting probability of the global type I error and demonstrate the power of the test. We also present the theoretical results to illustrate the behaviors of the test statistics in the presence of infinite activity jumps. Simulation studies indicate the favorable performance of the proposed test in finite samples, and we also apply the test to the stock price data of Apple and Microsoft.
This article studies uniform inference on a function $g(\cdot )$ and its functionals in a nonparametric panel data model with fixed effects. The nonparametric panel model relaxes restrictions on time-series behavior by allowing for arbitrary types of stationary or nonstationary dependence (e.g., stationary mixingale, mildly stationary, or local-to-unity process). After removing the fixed effects via transformations, a sieve estimator is proposed, accompanied by Yurinskii's coupling principle of Gaussian processes and uniform confidence bands (UCBs) that rely on the sieve score bootstrap method to test for linear functionals of $g(\cdot )$ . Under the asymptotic framework of an increasing cross-sectional dimension and either a fixed or diverging time dimension, we prove that the bootstrapping Kolmogorov-Smirnov (sup-type) test has asymptotic uniform size controls. This article shows that our uniform inference procedure can be extended to the two-way fixed-effects nonparametric panel model with stationary mixingale regressors. Extensive simulations confirm that our sieve estimators and their UCBs work well in finite samples. The present article further applies the above methods to empirical settings and finds some interesting results in nonlinear patterns of consumption concerning income shocks and asset holdings.
In high-dimensional (HD) sparse linear regression, parameter selection and estimation are addressed using a constraint $l_0$ on the direction of the parameter vector. We begin by establishing a general result that identifies this direction through the leading generalized eigenspace of specific measurable matrices. Using this result, we propose a novel approach to the selection of the best subsets by solving an empirical generalized eigenvalue problem to estimate the direction of the HD parameter. We then introduce a new estimator based on the RIFLE algorithm, providing a non-asymptotic bound for the estimation risk, minimax convergence, and a central limit theorem. Simulations demonstrate the superiority of our method over existing $l_0$ -constrained estimators.
Pretesting for exogeneity has become routine in many empirical applications involving instrumental variables (IVs) to decide whether the ordinary least squares or IV-based method is appropriate. Guggenberger (2010a, Econometric Theory, 26, 369-382) shows that the second-stage test - based on the outcome of a Durbin-Wu-Hausman-type pretest in the first stage - exhibits extreme size distortion, with asymptotic size equal to 1 when the standard critical values are used. In this paper, we first show that both conditional and unconditional on the data, standard wild bootstrap procedures are invalid for two-stage testing. Second, we propose an identification-robust two-stage test statistic that switches between OLS-based and weak-IV-robust statistics. Third, we develop a size-adjusted wild bootstrap approach for our two-stage test that integrates specific wild bootstrap critical values with an appropriate size-adjustment method. We establish uniform validity of this procedure under conditional heteroskedasticity or clustering in the sense that the resulting tests achieve correct asymptotic size, regardless of whether the identification is strong or weak. Our procedure is especially valuable for empirical researchers facing potential weak identification. In such settings, its power advantage is notable: whereas weak-IV-robust methods maintain correct size but often suffer from relatively low power, our approach achieves better performance.
This article proposes and studies two Huber-type estimation approaches, namely, the Huber instrumental variable (IV) estimation and the Huber generalized method of moments (GMM) estimation, for a spatial autoregressive model. We establish the consistency, asymptotic distributions, finite sample breakdown points, and influence functions of these estimators. Simulation studies show that compared to the corresponding traditional estimators (the two-stage least squares estimator, the best IV estimator, and the GMM estimator), our estimators are more robust when the unknown disturbances are long-tailed, and our estimators only lose a little efficiency when the disturbances are short-tailed. Moreover, the Huber GMM estimator also outperforms several robust estimators in the literature. Finally, we apply our estimation method to investigate the impact of the urban heat island effect on housing prices. A package is published on GitHub for practitioners to use in their empirical studies.
This article studies the principal component analysis (PCA) estimation of weak factor models with sparse loadings. We uncover an intrinsic near-sparsity preservation property for the PCA estimators of loadings, which comes from the approximately (block) upper triangular structure of the rotation matrix. It suggests an asymmetric relationship among factors: the sparsity of the rotated loadings for a stronger factor can be contaminated by the loadings from weaker ones, but the sparsity of the rotated loadings of a weaker factor is almost unaffected by the loadings of stronger ones. Then, we propose a simple alternative to the existing penalized approaches to sparsify the loading estimators by screening out the small PCA loading estimators directly, and construct consistent estimators for factor strengths. The proposed estimators perform well in finite samples, as shown by a set of Monte Carlo simulations.
This article presents novel methods and theories for estimation and inference about parameters in statistical models using machine learning for nuisance parameter estimation when data are dyadic. We propose a dyadic cross-fitting method to remove over-fitting biases under arbitrary dyadic dependence. Together with the use of Neyman orthogonal scores, this novel cross-fitting method enables root-n consistent estimation and inference robustly against dyadic dependence. We demonstrate its versatility by applying it to high-dimensional network formation models and reexamine the determinants of free trade agreements.
A first-order Gaussian autoregressive model is considered. The exact finite-sample joint density of the minimal sufficient statistic is derived, for any value of the autoregressive parameter. This allows us to derive explicitly the exact density of the autocorrelation coefficient and its Studentized t-ratio, whose densities were available only in the asymptotic case and not for all values of the parameter and the statistic. This article also demonstrates how to solve a general problem in statistical distribution theory (well beyond the specific case of autoregressive models), that of inverting confluent characteristic functions in multiple variables.
In a recent paper, Juodis and Reese (2022, Journal of Business & Economic Statistics, 40, 1191-1203) (JR) show that the application of the CD test proposed by Pesaran (2004, General diagnostic tests for cross-sectional dependence in panels, CWPE 0435, Cambridge) to residuals from panels with latent factors results in over-rejection. They propose a randomized test statistic to correct for over-rejection, and add a screening component to achieve power. This article considers the same problem but from a different perspective and shows that the standard CD test remains valid if the latent factors are weak. A bias-corrected version, CD $<^>{\ast}$ , is proposed which is shown to be asymptotically standard normal under the null of error cross-sectional independence which has power against network-type alternatives. This result is shown to hold for pure latent factor models as well as for panel regression models with latent factors. The case where the errors are serially correlated is also considered. Small sample properties of the CD $<^>{\ast}$ test are investigated by Monte Carlo experiments and are shown to have satisfactory small sample properties. In an empirical application, using the CD $<^>{\ast}$ test, it is shown that there remains spatial error dependence in a panel data model for real house price changes across 377 Metropolitan Statistical Areas in the United States, even after the effects of latent factors are filtered out.
This article introduces and analyzes a framework that accommodates general heterogeneity in regression modeling. It demonstrates that regression models with fixed or time-varying parameters can be estimated using the ordinary least squares (OLS) and time-varying OLS methods, respectively, across a broad class of regressors and noise processes not covered by existing theory. The proposed setting facilitates the development of asymptotic theory and the estimation of robust standard errors. The robust confidence interval estimators accommodate substantial heterogeneity in both regressors and noise. The resulting robust standard error estimates coincide with White's (1980, Econometrica 48, 817-838) heteroskedasticity-consistent estimator but are applicable to a broader range of conditions, including models with missing data. They are computationally simple and perform well in Monte Carlo simulations. Their robustness, generality, and ease of implementation make them highly suitable for empirical applications. Finally, the article provides a brief empirical illustration.
This paper develops distribution theory and bootstrap-based inference methods for a broad class of convex pairwise difference estimators. These estimators minimize a kernel-weighted convex-in-parameter function over observation pairs that are similar in terms of certain covariates, where the similarity is governed by a localization (bandwidth) parameter. While classical results establish asymptotic normality under restrictive bandwidth conditions, we show that valid Gaussian and bootstrap-based inference remains possible under substantially weaker assumptions. First, we extend the theory of small bandwidth asymptotics to convex pairwise estimation settings, deriving robust Gaussian approximations even when a smaller than standard bandwidth is used. Second, we employ a debiasing procedure based on generalized jackknifing to enable inference with larger bandwidths, while preserving convexity of the objective function. Third, we construct a novel bootstrap method that adjusts for bandwidth-induced variance distortions, yielding valid inference across a wide range of bandwidth choices. Our proposed inference method enjoys demonstrable more robustness, while retaining the practical appeal of convex pairwise difference estimators.
A general asymptotic theory is established for sample cross moments of nonstationary time series, allowing for long-range dependence and local unit roots. The theory provides a substantial extension of earlier results on nonparametric regression that include near-cointegrated nonparametric regression as well as spurious nonparametric regression. Many new models are covered by the limit theory, among which are functional coefficient regressions in which both regressors and the functional covariate are nonstationary. Simulations show finite sample performance matching well with the asymptotic theory and having broad relevance to applications, while revealing how dual nonstationarity in regressors and covariates raises sensitivity to bandwidth choice and the impact of dimensionality in nonparametric regression. An empirical example is provided involving climate data regression to assess Earth's climate sensitivity to CO $_2$ , where nonstationarity is a prominent feature of both the regressors and covariates in the model. To our knowledge, this application is the first nonparametric empirical analysis to assess potential nonlinear impacts of CO $_2$ on Earth's climate while allowing for nonstationarity in both the regressors and covariates.
Panel data often contain stayers (units with no within-variations) and slow movers (units with little within-variations). In the presence of many slow movers, conventional econometric methods can fail to work. We propose a novel method of robust inference for the average partial effects in correlated random coefficient models robustly across various distributions of within-variations, including the cases with many stayers and/or many slow movers in a unified manner. In addition to this robustness property, our proposed method entails smaller biases and hence improves accuracy in inference compared to existing alternatives. Simulation studies demonstrate our theoretical claims about these properties: the conventional 95% confidence interval covers the true parameter value with 37-93% frequencies, whereas our proposed one achieves 93-96% coverage frequencies.