Robust cointegration testing in the presence of weak trends,with an application to the human origin of global warming

Standard tests for the rank of cointegration of a vector autoregressive process present distributions that are affected by the presence of deterministic trends. We consider the recent approach of Demetrescu et al. (2009 Demetrescu, M., Lütkepohl, H., Saikkonen, P. (2009). Testing for the cointegrating rank of a vector autoregressive process with an uncertain deterministic trend term. Econometrics Journal 12:414–35.[Crossref], [Web of Science ®] , [Google Scholar]) who recommend testing a composite null. We assess this methodology in the presence of trends (linear or broken) whose magnitude is small enough not to be always detectable at conventional significance levels. We model them using local asymptotics and derive the properties of the test statistics. We show that whether the trend is orthogonal to the cointegrating vector has a major impact on the distributions but that the test combination approach remains valid. We apply of the methodology to the study of cointegration properties between global temperatures and the radiative forcing of human gas emissions. We find new evidence of Granger Causality.     

On the Properties of the Likelihood Function of Spanos' Conditional t Heteroskedastic Model

This article investigates the properties of the likelihood function of Spanos’ conditional t heteroskedastic model (Spanos, 1994 Spanos , A. ( 1994 ). On modeling heteroskedasticity: The student's t and elliptical linear regression models . Econometric Theor. 10 : 286 – 315 .[Crossref], [Web of Science ®] , [Google Scholar]) On modeling heteroskedasticity: the student's t and elliptical linear regression models. It is shown that estimability of the degrees of freedom of t distribution and the block-diagonality of the information matrix of the joint likelihood function with respect to conditional mean parameters and remaining parameters hold for the model. The joint maximum likelihood estimator and its inference based on the t-statistic and χ²-statistic are examined in finite samples by simulation when the degrees of freedom is known and unknown.     

Estimation of Long Memory in Integrated Variance

A stylized fact is that realized variance has long memory. We show that, when the instantaneous volatility is a long memory process of order d, the integrated variance is characterized by the same long-range dependence. We prove that the spectral density of realized variance is given by the sum of the spectral density of the integrated variance plus that of a measurement error, due to the sparse sampling and market microstructure noise. Hence, the realized volatility has the same degree of long memory as the integrated variance. The additional term in the spectral density induces a finite-sample bias in the semiparametric estimates of the long memory. A Monte Carlo simulation provides evidence that the corrected local Whittle estimator of Hurvich et al. (2005 Hurvich , C. M. , Moulines , E. , Soulier , P. ( 2005 ). Estimating long memory in volatility . Econometrica 73 ( 4 ): 1283 – 1328 .[Crossref], [Web of Science ®] , [Google Scholar]) is much less biased than the standard local Whittle estimator and the empirical application shows that it is robust to the choice of the sampling frequency used to compute the realized variance. Finally, the empirical results suggest that the volatility series are more likely to be generated by a nonstationary fractional process.     

On the Rounded Integer-Valued Autoregressive Process

In recent years, there has been a growing interest in modelling integred-valued time series. In this article, we propose a modified and generalized version of the first order rounded integer-valued autoregressive RINAR(1) model, originally introduced by Kachour and Yao (2009 Kachour , M. , Yao , J. F. ( 2009 ). First-order rounded integer-valued autoregressive (RINAR(1)) process . Journal of Time Series Analysis 30 ( 4 ): 417 – 448 .[Crossref], [Web of Science ®] , [Google Scholar]). Indeed, this class can be considered as an alternative of classical models based on the thinning operators. Using a Markov chain method, conditions for stationarity and the existence of moments are investigated. Least squares estimator of the model parameters is considered and its consistence is established. Finally, we describe the price change data using a model of the new class.     

Testing Parameter Constancy in Unit Root Autoregressive Models Against Multiple Continuous Structural Changes

This article considers tests for logistic smooth transition autoregressive (LSTAR) models accommodating multiple time dependent transitions between regimes when the data generating process is a random walk. The asymptotic null distributions of the tests, in contrast to the standard results in Lin and Teräsvirta (1994 Lin , C. F. J. , Teräsvirta , T. ( 1994 ). Testing the constancy of regression parameters against continuous structural change . Journal of Econometrics 62 : 211 – 228 .[Crossref], [Web of Science ®] , [Google Scholar]), are nonstandard. Monte Carlo experiments reveal that the tests have modest size distortions and satisfactory power against LSTAR models with multiple smooth breaks. The tests are applied to Swedish unemployment rates and the hysteresis hypothesis is over-turned in favour of an LSTAR model with two transitions between extreme regimes.     

First difference transformation in panel VAR models: Robustness,estimation, and inference

This article considers estimation of Panel Vector Autoregressive Models of order 1 (PVAR(1)) with focus on fixed T consistent estimation methods in First Differences (FD) with additional strictly exogenous regressors. Additional results for the Panel FD ordinary least squares (OLS) estimator and the FDLS type estimator of Han and Phillips (2010 Han, C., Phillips, P. C. B. (2010). Gmm estimation for dynamic panels with fixed effects and strong instruments at unity. Econometric Theory 26:119–151.[Crossref], [Web of Science ®] , [Google Scholar]) are provided. Furthermore, we simplify the analysis of Binder et al. (2005 Binder, M., Hsiao, C., Pesaran, M. H. (2005). Estimation and inference in short panel vector autoregressions with unit root and cointegration. Econometric Theory 21:795–837.[Crossref], [Web of Science ®] , [Google Scholar]) by providing additional analytical results and extend the original model by taking into account possible cross-sectional heteroscedasticity and presence of strictly exogenous regressors. We show that in the three wave panel the log-likelihood function of the unrestricted Transformed Maximum Likelihood (TML) estimator might violate the global identification assumption. The finite-sample performance of the analyzed methods is investigated in a Monte Carlo study.     

Using Implied Probabilities to Improve the Estimation of Unconditional Moment Restrictions for Weakly Dependent Data

In this article, we investigate the use of implied probabilities (Back and Brown, 1993) to improve estimation in unconditional moment conditions models. Using the seminal contributions of Bonnal and Renault (2001 Bonnal, H., Renault, E. (2001). Minimal Chi-Square Estimation with Conditional Moment Restrictions, Document de Travail, CESG, September 2001. [Google Scholar]) and Antoine et al. (2007 Antoine, B., Bonnal, H., Renault, E. (2007). On the efficient use of the informational content of estimating equations: Implied probabilities and euclidean empirical likelihood. Journal of Econometrics 138(2):461–487.[Crossref], [Web of Science ®] , [Google Scholar]), we propose two three-step Euclidian empirical likelihood (3S-EEL) estimators for weakly dependent data. Both estimators make use of a control variates principle that can be interpreted in terms of implied probabilities in order to achieve higher-order improvements relative to the traditional two-step GMM estimator. A Monte Carlo study reveals that the finite and large sample properties of the three-step estimators compare favorably to the existing approaches: the two-step GMM and the continuous updating estimator.     

Likelihood-based panel cointegration test in the presence of a linear time trend and cross-sectional dependence

This article proposes a new likelihood-based panel cointegration rank test which extends the test of Örsal and Droge (2014 Örsal, D. D. K., Droge, B. (2014). Panel cointegration testing in the presence of a time trend. Computational Statistics and Data Analysis 76:377–390.[Crossref], [Web of Science ®] , [Google Scholar]) (henceforth panel SL test) to dependent panels. The dependence is modelled by unobserved common factors which affect the variables in each cross-section through heterogeneous loadings. The data are defactored following the panel analysis of nonstationarity in idiosyncratic and common components (PANIC) approach of Bai and Ng (2004 Bai, J., Ng, S. (2004). A PANIC attack on unit roots and cointegration. Econometrica 72(4):1127–1177.[Crossref], [Web of Science ®] , [Google Scholar]) and the cointegrating rank of the defactored data is then tested by the panel SL test. A Monte Carlo study demonstrates that the proposed testing procedure has reasonable size and power properties in finite samples.     

Inference on locally ordered breaks in multiple regressions

We consider issues related to inference about locally ordered breaks in a system of equations, as originally proposed by Qu and Perron (2007 Qu, Z., Perron, P. (2007). Estimating and testing structural changes in multivariate regressions. Econometrica 75:459–502.[Crossref], [Web of Science ®] , [Google Scholar]). These apply when break dates in different equations within the system are not separated by a positive fraction of the sample size. This allows constructing joint confidence intervals of all such locally ordered break dates. We extend the results of Qu and Perron (2007 Qu, Z., Perron, P. (2007). Estimating and testing structural changes in multivariate regressions. Econometrica 75:459–502.[Crossref], [Web of Science ®] , [Google Scholar]) in several directions. First, we allow the covariates to be any mix of trends and stationary or integrated regressors. Second, we allow for breaks in the variance-covariance matrix of the errors. Third, we allow for multiple locally ordered breaks, each occurring in a different equation within a subset of equations in the system. Via some simulation experiments, we show first that the limit distributions derived provide good approximations to the finite sample distributions. Second, we show that forming confidence intervals in such a joint fashion allows more precision (tighter intervals) compared to the standard approach of forming confidence intervals using the method of Bai and Perron (1998 Bai, J., Perron, P. (1998). Estimating and testing linear models with multiple structural changes. Econometrica 66:47–78.[Crossref], [Web of Science ®] , [Google Scholar]) applied to a single equation. Simulations also indicate that using the locally ordered break confidence intervals yields better coverage rates than using the framework for globally distinct breaks when the break dates are separated by roughly 10% of the total sample size.     

Empirical likelihood ratio-based goodness-of-fit test for the logistic distribution

The logistic distribution has been used to model growth curves in survival analysis and biological studies. In this article, we propose a goodness-of-fit test for the logistic distribution based on the empirical likelihood ratio. The test is constructed based on the methodology introduced by Vexler and Gurevich [17 A. Vexler and G. Gurevich, Empirical likelihood ratios applied to goodness-of-fit tests based on sample entropy, Comput. Stat. Data Anal. 54 (2010), pp. 531–545. doi: 10.1016/j.csda.2009.09.025[Crossref], [Web of Science ®] , [Google Scholar]]. In order to compute the test statistic, parameters of the distribution are estimated by the method of maximum likelihood. Power comparisons of the proposed test with some known competing tests are carried out via simulations. Finally, an illustrative example is presented and analyzed.     

A Prediction-Oriented Criterion for Choosing the Biasing Parameter in Liu Estimation

Distribution-free tests have been proposed in the literature for comparing the hazard rates of two probability distributions when the available samples are complete. In this article, we generalize the test of Kochar (1981 Kochar , S. C. ( 1981 ). A new distribution-free test for the equality of two failure rates . Biometrika 68 : 423 – 426 .[Crossref], [Web of Science ®] , [Google Scholar]) to the case when the available sample is Type-II censored, and then examine its power properties.     

Longitudinal profiles of bounded outcome scores as predictors for disease activity in rheumatoid arthritis patients: a joint modeling approach

A variety of statistical approaches have been suggested in the literature for the analysis of bounded outcome scores (BOS). In this paper, we suggest a statistical approach when BOSs are repeatedly measured over time and used as predictors in a regression model. Instead of directly using the BOS as a predictor, we propose to extend the approaches suggested in [16 E. Lesaffre, D. Rizopoulos, and R. Tsonaka, The logistics-transform for bounded outcome scores, Biostatistics 8 (2007), pp. 72–85. doi: 10.1093/biostatistics/kxj034[Crossref], [PubMed], [Web of Science ®] , [Google Scholar],21 M. Molas and E. Lesaffre, A comparison of the three random effects approaches to analyse repeated bounded outcome scores with an application in a stroke revalidation study, Stat. Med. 27 (2008), pp. 6612–6633. doi: 10.1002/sim.3432[Crossref], [PubMed], [Web of Science ®] , [Google Scholar],28 R. Tsonaka, D. Rizopoulos, and E. Lesaffre, Power and sample size calculations for discrete bounded outcome scores, Stat. Med. 25 (2006), pp. 4241–4252. doi: 10.1002/sim.2679[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]] to a joint modeling setting. Our approach is illustrated on longitudinal profiles of multiple patients’ reported outcomes to predict the current clinical status of rheumatoid arthritis patients by a disease activities score of 28 joints (DAS28). Both a maximum likelihood as well as a Bayesian approach is developed.     

Testing for sphericity in a two-way error components panel data model

This article is concerned with sphericity test for the two-way error components panel data model. It is found that the John statistic and the bias-corrected LM statistic recently developed by Baltagi et al. (2011 Baltagi, B. H., Feng, Q., Kao, C. (2011). Testing for sphericity in a fixed effects panel data model. Econometrics Journal 14:25–47.[Crossref], [Web of Science ®] , [Google Scholar])Baltagi et al. (2012 Baltagi, B. H., Feng, Q., Kao, C. (2012). A Lagrange multiplier test for cross-sectional dependence in a fixed effects panel data model. Journal of Econometrics 170:164–177.[Crossref], [Web of Science ®] , [Google Scholar], which are based on the within residuals, are not helpful under the present circumstances even though they are in the one-way fixed effects model. However, we prove that when the within residuals are properly transformed, the resulting residuals can serve to construct useful statistics that are similar to those of Baltagi et al. (2011 Baltagi, B. H., Feng, Q., Kao, C. (2011). Testing for sphericity in a fixed effects panel data model. Econometrics Journal 14:25–47.[Crossref], [Web of Science ®] , [Google Scholar])Baltagi et al. (2012 Baltagi, B. H., Feng, Q., Kao, C. (2012). A Lagrange multiplier test for cross-sectional dependence in a fixed effects panel data model. Journal of Econometrics 170:164–177.[Crossref], [Web of Science ®] , [Google Scholar]). Simulation results show that the newly proposed statistics perform well under the null hypothesis and several typical alternatives.     

Bootstrapping unit root tests with covariates

We consider the bootstrap method for the covariates augmented Dickey–Fuller (CADF) unit root test suggested in Hansen (1995 Hansen, B. E. (1995). Rethinking the univariate approach to unit root testing: Using covariates to increase power. \ Econometric Theory 11:1148–1171.[Crossref], [Web of Science ®] , [Google Scholar]) which uses related variables to improve the power of univariate unit root tests. It is shown that there are substantial power gains from including correlated covariates. The limit distribution of the CADF test, however, depends on the nuisance parameter that represents the correlation between the equation error and the covariates. Hence, inference based directly on the CADF test is not possible. To provide a valid inferential basis for the CADF test, we propose to use the parametric bootstrap procedure to obtain critical values, and establish the asymptotic validity of the bootstrap CADF test. Simulations show that the bootstrap CADF test significantly improves the asymptotic and the finite sample size performances of the CADF test, especially when the covariates are highly correlated with the error. Indeed, the bootstrap CADF test offers drastic power gains over the conventional unit root tests. Our testing procedures are applied to the extended Nelson and Plosser data set.     

Weighted simulated integrated conditional moment tests for parametric conditional distributions of stationary time series processes

In this article, we propose a weighted simulated integrated conditional moment (WSICM) test of the validity of parametric specifications of conditional distribution models for stationary time series data, by combining the weighted integrated conditional moment (ICM) test of Bierens (1984 Bierens, H. J. (1984). Model specification testing of time series regressions. Journal of Econometrics 26:323–353.[Crossref], [Web of Science ®] , [Google Scholar]) for time series regression models with the simulated ICM test of Bierens and Wang (2012 Bierens, H. J., Wang, L. (2012). Integrated conditional moment tests for parametric conditional distributions. Econometric Theory 28:328–362.[Crossref], [Web of Science ®] , [Google Scholar]) of conditional distribution models for cross-section data. To the best of our knowledge, no other consistent test for parametric conditional time series distributions has been proposed yet in the literature, despite consistency claims made by some authors.     

A comparison of alternative methods to construct confidence intervals for the estimate of a break date in linear regression models

This article considers constructing confidence intervals for the date of a structural break in linear regression models. Using extensive simulations, we compare the performance of various procedures in terms of exact coverage rates and lengths of the confidence intervals. These include the procedures of Bai (1997 Bai, J. (1997). Estimation of a change point in multiple regressions. Review of Economics and Statistics 79:551–563.[Crossref], [Web of Science ®] , [Google Scholar]) based on the asymptotic distribution under a shrinking shift framework, Elliott and Müller (2007 Elliott, G., Müller, U. (2007). Confidence sets for the date of a single break in linear time series regressions. Journal of Econometrics 141:1196–1218.[Crossref], [Web of Science ®] , [Google Scholar]) based on inverting a test locally invariant to the magnitude of break, Eo and Morley (2015 Eo, Y., Morley, J. (2015). Likelihood-ratio-based confidence sets for the timing of structural breaks. Quantitative Economics 6:463–497.[Crossref], [Web of Science ®] , [Google Scholar]) based on inverting a likelihood ratio test, and various bootstrap procedures. On the basis of achieving an exact coverage rate that is closest to the nominal level, Elliott and Müller's (2007 Elliott, G., Müller, U. (2007). Confidence sets for the date of a single break in linear time series regressions. Journal of Econometrics 141:1196–1218.[Crossref], [Web of Science ®] , [Google Scholar]) approach is by far the best one. However, this comes with a very high cost in terms of the length of the confidence intervals. When the errors are serially correlated and dealing with a change in intercept or a change in the coefficient of a stationary regressor with a high signal-to-noise ratio, the length of the confidence interval increases and approaches the whole sample as the magnitude of the change increases. The same problem occurs in models with a lagged dependent variable, a common case in practice. This drawback is not present for the other methods, which have similar properties. Theoretical results are provided to explain the drawbacks of Elliott and Müller's (2007 Elliott, G., Müller, U. (2007). Confidence sets for the date of a single break in linear time series regressions. Journal of Econometrics 141:1196–1218.[Crossref], [Web of Science ®] , [Google Scholar]) method.     

Negative variance components for non-negative hierarchical data with correlation,over-, and/or underdispersion

The concept of negative variance components in linear mixed-effects models, while confusing at first sight, has received considerable attention in the literature, for well over half a century, following the early work of Chernoff [7 H. Chernoff, On the distribution of the likelihood ratio, Ann. Math. Statist. 25 (1954), pp. 573–578.[Crossref] , [Google Scholar]] and Nelder [21 J.A. Nelder, The interpretation of negative components of variance, Biometrika 41 (1954), pp. 544–548.[Crossref], [Web of Science ®] , [Google Scholar]]. Broadly, negative variance components in linear mixed models are allowable if inferences are restricted to the implied marginal model. When a hierarchical view-point is adopted, in the sense that outcomes are specified conditionally upon random effects, the variance–covariance matrix of the random effects must be positive-definite (positive-semi-definite is also possible, but raises issues of degenerate distributions). Many contemporary software packages allow for this distinction. Less work has been done for generalized linear mixed models. Here, we study such models, with extension to allow for overdispersion, for non-negative outcomes (counts). Using a study of trichomes counts on tomato plants, it is illustrated how such negative variance components play a natural role in modeling both the correlation between repeated measures on the same experimental unit and over- or underdispersion.     

Signal extraction goodness-of-fit diagnostic tests under model parameter uncertainty: Formulations and empirical evaluation

We present a time-domain goodness-of-fit (gof) diagnostic test that is based on signal-extraction variances for nonstationary time series. This diagnostic test extends the time-domain gof statistic of Maravall (2003 Maravall, A. (2003). A class of diagnostics in the ARIMA-model-based decomposition of a time series. Memorandum, Bank of Spain. Available at http://www.bde.es/servicio/software/tramo/diagnosticsamb.pdf [Google Scholar]) by taking into account the effects of model parameter uncertainty, utilizing theoretical results of McElroy and Holan (2009 McElroy, T., Holan, S. (2009). A local spectral approach for assessing time series model misspeci?cation. Journal of Multivariate Analysis 100:604–621.[Crossref], [Web of Science ®] , [Google Scholar]). We demonstrate that omitting this correction results in a severely undersized statistic. Adequate size and power are obtained in Monte Carlo studies for fairly short time series (10 to 15 years of monthly data). Our Monte Carlo studies of finite sample size and power consider different combinations of both signal and noise components using seasonal, trend, and irregular component models obtained via canonical decomposition. Details of the implementation appropriate for SARIMA models are given. We apply the gof diagnostic test statistics to several U.S. Census Bureau time series. The results generally corroborate the output of the automatic model selection procedure of the X-12-ARIMA software, which in contrast to our diagnostic test statistic does not involve hypothesis testing. We conclude that these diagnostic test statistics are a useful supplementary model-checking tool for practitioners engaged in the task of model-based seasonal adjustment.     

A modified confidence set for the structural break date in linear regression models

Elliott and Müller (EM) (2007 Elliott, G., Müller, U. K. (2007). Confidence sets for the date of a single break in linear time series regressions. Journal of Econometrics 141:1196–1218.[Crossref], [Web of Science ®] , [Google Scholar]) provide a method for constructing a confidence set for the structural break date by inverting a variant of the locally best test statistic. Previous studies have shown that the EM method produces a set with an accurate coverage ratio even for a small break; however, the set is often overly lengthy. This study proposes a simple modification to rehabilitate their method through the long-run variance estimation. Following the literature, we provide an asymptotic justification for the improvement of the modified method over the original method under a nonlocal asymptotic framework. A Monte Carlo simulation shows that the modified method achieves a shorter confidence set than the EM method, especially when the break is large or the HAC correction is conducted. The modified method may exhibit minor errors in the coverage rate when the break is small; however, the coverage is more stable than alternative methods when the break is large. We apply our method to a level shift in post-1980s Japanese inflation data.     

Permutation Matrix with Application to Optimal Invariant Quadratic Prediction in Finite Populations

In this article, we primarily aim to apply the permutation matrix techniques to the problem of the optimal invariant quadratic prediction in a finite population. An alternative to the work of Liu and Rong (2007 Liu , X. , Rong , J. ( 2007 ). Quadratic prediction problems in finite populations . Statist. Probab. Lett. 77 : 483 – 489 .[Crossref], [Web of Science ®] , [Google Scholar]) is offered. In addition, we derive the OIQP for the population variance and show that it has less PMSE than the ordinary optimal unbiased predictor.