Heteroskedasticity Robust Panel Unit Root Testing Under Variance Breaks in Pooled Regressions

Noting that many economic variables display occasional shifts in their second order moments, we investigate the performance of homogenous panel unit root tests in the presence of permanent volatility shifts. It is shown that in this case the test statistic proposed by Herwartz and Siedenburg (2008 Herwartz, H., Siedenburg, F. (2008). Homogenous panel unit root tests under cross-sectional dependence: Finite sample modifications and the wild bootstrap. Computational Statistics and Data Analysis 53(1):137–150.[Crossref], [Web of Science ®] , [Google Scholar]) is asymptotically standard Gaussian. By means of a simulation study we illustrate the performance of first and second generation panel unit root tests and undertake a more detailed comparison of the test in Herwartz and Siedenburg (2008 Herwartz, H., Siedenburg, F. (2008). Homogenous panel unit root tests under cross-sectional dependence: Finite sample modifications and the wild bootstrap. Computational Statistics and Data Analysis 53(1):137–150.[Crossref], [Web of Science ®] , [Google Scholar]) and its heteroskedasticity consistent Cauchy counterpart introduced in Demetrescu and Hanck (2012a Demetrescu, M., Hanck, C. (2012a). A simple nonstationary-volatility robust panel unit root test. Economics Letters 117(2):10–13.[Crossref], [Web of Science ®] , [Google Scholar]). As an empirical illustration, we reassess evidence on the Fisher hypothesis with data from nine countries over the period 1961Q2–2011Q2. Empirical evidence supports panel stationarity of the real interest rate for the entire subperiod. With regard to the most recent two decades, the test results cast doubts on market integration, since the real interest rate is diagnosed nonstationary.     

Random Effects,Fixed Effects and Hausman's Test for the Generalized Mixed Regressive Spatial Autoregressive Panel Data Model

This article suggests random and fixed effects spatial two-stage least squares estimators for the generalized mixed regressive spatial autoregressive panel data model. This extends the generalized spatial panel model of Baltagi et al. (2013 Baltagi, B. H., Egger, P., Pfaffermayr, M. (2013). A generalized spatial panel data model with random effects. Econometric Reviews 32:650–685.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) by the inclusion of a spatial lag term. The estimation method utilizes the Generalized Moments method suggested by Kapoor et al. (2007 Kapoor, M., Kelejian, H. H., Prucha, I. R. (2007). Panel data models with spatially correlated error components. Journal of Econometrics 127(1):97–130.[Crossref], [Web of Science ®] , [Google Scholar]) for a spatial autoregressive panel data model. We derive the asymptotic distributions of these estimators and suggest a Hausman test a la Mutl and Pfaffermayr (2011 Mutl, J., Pfaffermayr, M. (2011). The Hausman test in a Cliff and Ord panel model. Econometrics Journal 14:48–76.[Crossref], [Web of Science ®] , [Google Scholar]) based on the difference between these estimators. Monte Carlo experiments are performed to investigate the performance of these estimators as well as the corresponding Hausman test.     

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.     

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.     

Implementing Box–Cox Quantile Regression

The Box–Cox quantile regression model introduced by Powell (1991 Powell , J. ( 1991 ). Estimation of monotonic regression models under quantile restrictions . In: Barnett , W. , Powell , J. , Tauchen , G. , eds. Nonparametric and Semiparametric Methods in Econometrics . New York , NY : Cambridge University Press , pp. 357 – 384 . [Google Scholar]) is a flexible and numerically attractive extension of linear quantile regression techniques. Chamberlain (1994 Chamberlain , G. ( 1994 ). Quantile regression, censoring, and the structure of wages . In: Sims , C. , ed. Advances in Econometrics: Sixth World Congress . Vol. 1 . Econometric Society Monograph . Cambridge : Cambridge University Press . [Google Scholar]) and Buchinsky (1995 Buchinsky , M. ( 1995 ). Quantile regression, Box–Cox transformation model, and the U.S. wage structure, 1963–1987 . Journal of Econometrics 65 : 109 – 154 .[Crossref], [Web of Science ®] , [Google Scholar]) suggest a two stage estimator for this model but the objective function in stage two of their method may not be defined in an application. We suggest a modification of the estimator which is easy to implement. A simulation study demonstrates that the modified estimator works well in situations, where the original estimator is not well defined.     

Sparse group variable selection based on quantile hierarchical Lasso

The group Lasso is a penalized regression method, used in regression problems where the covariates are partitioned into groups to promote sparsity at the group level [27 M. Yuan and Y. Lin, Model selection and estimation in regression with grouped variables, J. R. Stat. Soc. Ser. B 68 (2006), pp. 49–67. doi: 10.1111/j.1467-9868.2005.00532.x[Crossref] , [Google Scholar]]. Quantile group Lasso, a natural extension of quantile Lasso [25 Y. Wu and Y. Liu, Variable selection in quantile regression, Statist. Sinica 19 (2009), pp. 801–817.[Web of Science ®] , [Google Scholar]], is a good alternative when the data has group information and has many outliers and/or heavy tails. How to discover important features that are correlated with interest of outcomes and immune to outliers has been paid much attention. In many applications, however, we may also want to keep the flexibility of selecting variables within a group. In this paper, we develop a sparse group variable selection based on quantile methods which select important covariates at both the group level and within the group level, which penalizes the empirical check loss function by the sum of square root group-wise L₁-norm penalties. The oracle properties are established where the number of parameters diverges. We also apply our new method to varying coefficient model with categorial effect modifiers. Simulations and real data example show that the newly proposed method has robust and superior performance.     

Symmetric quantiles and their applications

To develop estimators with stronger efficiencies than the trimmed means which use the empirical quantile, Kim (1992) Kim, S. J. 1992. The metrically trimmed means as a robust estimator of location. Annals of Statistics, 20: 1534–1547. [Crossref], [Web of Science ®] [Google Scholar] and Chen & Chiang (1996) Chen, L. A. and Chiang, Y. C. 1996. Symmetric type quantile and trimmed means for location and linear regression model. Journal of Nonparametric Statistics, 7: 171–185. [Taylor & Francis Online] [Google Scholar], implicitly or explicitly used the symmetric quantile, and thus introduced new trimmed means for location and linear regression models, respectively. This study further investigates the properties of the symmetric quantile and extends its application in several aspects. (a) The symmetric quantile is more efficient than the empirical quantiles in asymptotic variances when quantile percentage α is either small or large. This reveals that for any proposal involving the α th quantile of small or large α s, the symmetric quantile is the right choice; (b) a trimmed mean based on it has asymptotic variance achieving a Cramer-Rao lower bound in one heavy tail distribution; (c) an improvement of the quantiles-based control chart by Grimshaw & Alt (1997) Grimshaw, S. D. and Alt, F. B. 1997. Control charts for quantile function values. Journal of Quality Technology, 29: 1–7. [Taylor & Francis Online] [Google Scholar] is discussed; (d) Monte Carlo simulations of two new scale estimators based on symmetric quantiles also support this new quantile.     

MSE Performance of a Heterogeneous Pre-Test Ridge Regression Estimator

This article describes how diagnostic procedures were derived for symmetrical nonlinear regression models, continuing the work carried out by Cysneiros and Vanegas (2008 Cysneiros , F. J. A. , Vanegas , L. H. ( 2008 ). Residuals and their statistical properties in symmetrical nonlinear models . Statist. Probab. Lett. 78 : 3269 – 3273 .[Crossref], [Web of Science ®] , [Google Scholar]) and Vanegas and Cysneiros (2010 Vanegas , L. H. , Cysneiros , F. J. A. ( 2010 ). Assesment of diagnostic procedures in symmetrical nonlinear regression models . Computat. Statist. Data Anal. 54 : 1002 – 1016 .[Crossref], [Web of Science ®] , [Google Scholar]), who showed that the parameters estimates in nonlinear models are more robust with heavy-tailed than with normal errors. In this article, we focus on assessing if the robustness of this kind of models is also observed in the inference process (i.e., partial F-test). Symmetrical nonlinear regression models includes all symmetric continuous distributions for errors covering both light- and heavy-tailed distributions such as Student-t, logistic-I and -II, power exponential, generalized Student-t, generalized logistic, and contaminated normal. Firstly, a statistical test is shown to evaluating the assumption that the error terms all have equal variance. The results of simulation studies which describe the behavior of the test for heteroscedasticity proposed in the presence of outliers are then given. To assess the robustness of inference process, we present the results of a simulation study which described the behavior of partial F-test in the presence of outliers. Also, some diagnostic procedures are derived to identify influential observations on the partial F-test. As ilustration, a dataset described in Venables and Ripley (2002 Venables , W. N. , Ripley , B. D. ( 2002 ). Modern Applied with S. , 4th ed. New York : Springer .[Crossref] , [Google Scholar]), is also analyzed.     

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.     

Some Theoretical and Simulation Results on the Frequency Domain Causality Test

Breitung and Candelon (2006 Breitung , J. , Candelon , B. ( 2006 ). Testing for short- and long-run causality: A frequency-domain approach . Journal of Econometrics 132 : 363 – 378 .[Crossref], [Web of Science ®] , [Google Scholar]) in Journal of Econometrics proposed a simple statistical testing procedure for the noncausality hypothesis at a given frequency. In their paper, however, they reported some theoretical results indicating that their test severely suffers from quite low power when the noncausality hypothesis is tested at a frequency close to 0 or pi. This paper examines whether or not these results indicate their procedure is useless at such frequencies.     

Testing for Unit Roots and the Impact of Quadratic Trends,with an Application to Relative Primary Commodity Prices

In practice a degree of uncertainty will always exist concerning what specification to adopt for the deterministic trend function when running unit root tests. While most macroeconomic time series appear to display an underlying trend, it is often far from clear whether this component is best modeled as a simple linear trend (so that long-run growth rates are constant) or by a more complicated nonlinear trend function which may, for instance, allow the deterministic trend component to evolve gradually over time. In this article, we consider the effects on unit root testing of allowing for a local quadratic trend, a simple yet very flexible example of the latter. Where a local quadratic trend is present but not modeled, we show that the quasi-differenced detrended Dickey–Fuller-type test of Elliott et al. (1996 Elliott , G. , Rothenberg , T. J. , Stock , J. H. ( 1996 ). Efficient tests for an autoregressive unit root . Econometrica 64 : 813 – 836 .[Crossref], [Web of Science ®] , [Google Scholar]) has both size and power which tend to zero asymptotically. An extension of the Elliott et al. (1996 Elliott , G. , Rothenberg , T. J. , Stock , J. H. ( 1996 ). Efficient tests for an autoregressive unit root . Econometrica 64 : 813 – 836 .[Crossref], [Web of Science ®] , [Google Scholar]) approach to allow for a quadratic trend resolves this problem but is shown to result in large power losses relative to the standard detrended test when no quadratic trend is present. We consequently propose a simple and practical approach to dealing with this form of uncertainty based on a union of rejections-based decision rule whereby the unit root is rejected whenever either of the detrended or quadratic detrended unit root tests rejects. A modification of this basic strategy is also suggested which further improves on the properties of the procedure. An application to relative primary commodity price data highlights the empirical relevance of the methods outlined in this article. A by-product of our analysis is the development of a test for the presence of a quadratic trend which is robust to whether the data admit a unit root.     

Weighting variables in Kohonen competitive learning algorithms

This paper presents a new variable weight method, called the singular value decomposition (SVD) approach, for Kohonen competitive learning (KCL) algorithms based on the concept of Varshavsky et al. [18 R. Varshavsky, A. Gottlieb, M. Linial, and D. Horn, Novel unsupervised feature filtering of bilogical data, Bioinformatics 22 (2006), pp. 507–513.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]]. Integrating the weighted fuzzy c-means (FCM) algorithm with KCL, in this paper, we propose a weighted fuzzy KCL (WFKCL) algorithm. The goal of the proposed WFKCL algorithm is to reduce the clustering error rate when data contain some noise variables. Compared with the k-means, FCM and KCL with existing variable-weight methods, the proposed WFKCL algorithm with the proposed SVD's weight method provides a better clustering performance based on the error rate criterion. Furthermore, the complexity of the proposed SVD's approach is less than Pal et al. [17 S.K. Pal, R.K. De, and J. Basak, Unsupervised feature evaluation: a neuro-fuzzy approach, IEEE. Trans. Neural Netw. 11 (2000), pp. 366–376.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]], Wang et al. [19 X.Z. Wang, Y.D. Wang, and L.J. Wang, Improving fuzzy c-means clustering based on feature-weight learning, Pattern Recognit. Lett. 25 (2004), pp. 1123–1132.[Crossref], [Web of Science ®] , [Google Scholar]] and Hung et al. [9 W. -L. Hung, M. -S. Yang, and D. -H. Chen, Bootstrapping approach to feature-weight selection in fuzzy c-means algorithms with an application in color image segmentation, Pattern Recognit. Lett. 29 (2008), pp. 1317–1325.[Crossref], [Web of Science ®] , [Google Scholar]].     

Spurious Rejections with Endogenous Break Unit Root Tests in the Presence of Outliers and Breaks

Unit root tests with structural break developed by Zivot and Andrews (1992 Zivot , E. , Andrews , D. W. K. ( 1992 ). Further evidence on the great crash, the oil price shock and the unit root hypothesis . Journal of Business and Economic Statistics 10 : 251 – 270 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Perron and Rodriguez (2003 Perron , P. , Rodriguez , G. ( 2003 ). GLS detrending, efficient unit root tests and structural change . Journal of Econometrics 115 : 1 – 27 .[Crossref], [Web of Science ®] , [Google Scholar]) in the presence of additive outliers and breaks are studied by simulation experiments. The results show that the Zivot–Andrews test appears to have size distortions due to the additive outliers whereas the Perron–Rodriguez test exhibits good properties of size and power. However, the two tests are biased when a second break is present but not taken into account. Furthermore, these endogenous break unit root tests tend to determine the break point incorrectly at one period behind the true break point, leading to spurious rejections of the unit root null hypothesis.     

Asymptotic Properties of Conditional Quantile Estimator Under Left-Truncated and α-Mixing Conditions

In this article, we establish strong uniform convergence and asymptotic normality of estimators of conditional quantile and conditional distribution function for a left truncated model when the data exhibit some kind of dependence. It is assumed that the observations form a stationary α-mixing sequence. The results of Lemdani et al. (2009 Lemdani , M. , Ould-Saïd , E. , Poulin , N. ( 2009 ). Asymptotic properties of a conditional quantile estimator with randomly truncated data . J. Multivariate Anal. 100 : 546 – 559 .[Crossref], [Web of Science ®] , [Google Scholar]) are relaxed from the i.i.d. assumption to α-mixing setting. Finite sample behavior of the estimators is investigated via simulations as well.     

Unified multivariate survival model with a surviving fraction: an application to a Brazilian customer churn data

In this paper we propose a new lifetime model for multivariate survival data in presence of surviving fractions and examine some of its properties. Its genesis is based on situations in which there are m types of unobservable competing causes, where each cause is related to a time of occurrence of an event of interest. Our model is a multivariate extension of the univariate survival cure rate model proposed by Rodrigues et al. [37 J. Rodrigues, V.G. Cancho, M. de Castro, and F. Louzada-Neto, On the unification of long-term survival models, Statist. Probab. Lett. 79 (2009), pp. 753–759. doi: 10.1016/j.spl.2008.10.029[Crossref], [Web of Science ®] , [Google Scholar]]. The inferential approach exploits the maximum likelihood tools. We perform a simulation study in order to verify the asymptotic properties of the maximum likelihood estimators. The simulation study also focus on size and power of the likelihood ratio test. The methodology is illustrated on a real data set on customer churn data.     

Evaluating the relative merits of competing models based on empirical likelihood ratio test

Competing models arise naturally in many research fields, such as survival analysis and economics, when the same phenomenon of interest is explained by different researcher using different theories or according to different experiences. The model selection problem is therefore remarkably important because of its great importance to the subsequent inference; Inference under a misspecified or inappropriate model will be risky. Existing model selection tests such as Vuong's tests [26 Q.H. Vuong, Likelihood ratio test for model selection and non-nested hypothesis, Econometrica 57 (1989), pp. 307–333. doi: 10.2307/1912557[Crossref], [Web of Science ®] , [Google Scholar]] and Shi's non-degenerate tests [21 X. Shi, A non-degenerate Vuong test, Quant. Econ. 6 (2015), pp. 85–121. doi: 10.3982/QE382[Crossref], [Web of Science ®] , [Google Scholar]] suffer from the variance estimation and the departure of the normality of the likelihood ratios. To circumvent these dilemmas, we propose in this paper an empirical likelihood ratio (ELR) tests for model selection. Following Shi [21 X. Shi, A non-degenerate Vuong test, Quant. Econ. 6 (2015), pp. 85–121. doi: 10.3982/QE382[Crossref], [Web of Science ®] , [Google Scholar]], a bias correction method is proposed for the ELR tests to enhance its performance. A simulation study and a real-data analysis are provided to illustrate the performance of the proposed ELR tests.     

Hip fracture prediction from a new classification algorithm based on recursive partitioning methods

Classification and regression tree has been useful in medical research to construct algorithms for disease diagnosis or prognostic prediction. Jin et al. 7 Jin, H., Lu, Y., Harris, R. T., Black, D., Stone, K., Hochberg, M. and Genant, H. 2004. Classification algorithms for hip fracture prediction base on recursive partitioning methods. Med. Decis. Mak., 24: 386–398. (doi:10.1177/0272989X04267009)[Crossref], [PubMed], [Web of Science ®] , [Google Scholar] developed a robust and cost-saving tree (RACT) algorithm with application in classification of hip fracture risk after 5-year follow-up based on the data from the Study of Osteoporotic Fractures (SOF). Although conventional recursive partitioning algorithms have been well developed, they still have some limitations. Binary splits may generate a big tree with many layers, but trinary splits may produce too many nodes. In this paper, we propose a classification approach combining trinary splits and binary splits to generate a trinary–binary tree. A new non-inferiority test of entropy is used to select the binary or trinary splits. We apply the modified method in SOF to construct a trinary–binary classification rule for predicting risk of osteoporotic hip fracture. Our new classification tree has good statistical utility: it is statistically non-inferior to the optimum binary tree and the RACT based on the testing sample and is also cost-saving. It may be useful in clinical applications: femoral neck bone mineral density, age, height loss and weight gain since age 25 can identify subjects with elevated 5-year hip fracture risk without loss of statistical efficiency.     

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.     

A modified Hosmer–Lemeshow test for large data sets

The Hosmer–Lemeshow test is a widely used method for evaluating the goodness of fit of logistic regression models. But its power is much influenced by the sample size, like other chi-square tests. Paul, Pennell, and Lemeshow (2013 Paul, P., M. L. Pennell, and S. Lemeshow. 2013. Standardizing the power of the Hosmer–Lemeshow goodness of fit test in large data sets. Statistics in Medicine 32:67–80.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) considered using a large number of groups for large data sets to standardize the power. But simulations show that their method performs poorly for some models. In addition, it does not work when the sample size is larger than 25,000. In the present paper, we propose a modified Hosmer–Lemeshow test that is based on estimation and standardization of the distribution parameter of the Hosmer–Lemeshow statistic. We provide a mathematical derivation for obtaining the critical value and power of our test. Through simulations, we can see that our method satisfactorily standardizes the power of the Hosmer–Lemeshow test. It is especially recommendable for enough large data sets, as the power is rather stable. A bank marketing data set is also analyzed for comparison with existing methods.     

Inference for impulse response coefficients from multivariate fractionally integrated processes

This article considers a multivariate system of fractionally integrated time series and investigates the most appropriate way for estimating Impulse Response (IR) coefficients and their associated confidence intervals. The article extends the univariate analysis recently provided by Baillie and Kapetanios (2013 Baillie, R. T., Kapetanios, G. (2013). Estimation and inference for impulse response functions form univariate strongly persistent processes. Econometrics Journal 16:373–399.[Crossref], [Web of Science ®] , [Google Scholar]), and uses a semiparametric, time domain estimator, based on a vector autoregression (VAR) approximation. Results are also derived for the orthogonalized estimated IRs which are generally more practically relevant. Simulation evidence strongly indicates the desirability of applying the Kilian small sample bias correction, which is found to improve the coverage accuracy of confidence intervals for IRs. The most appropriate order of the VAR turns out to be relevant for the lag length of the IR being estimated.