Next Article in Journal
A Systems Perspective in Examining Industry Clusters: Case Studies of Clusters in Russia and India
Next Article in Special Issue
Contrasting Cryptocurrencies with Other Assets: Full Distributions and the COVID Impact
Previous Article in Journal
What Information in Financial Statements Could Be Used to Predict the Risk of Equity Investment?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spurious Relationships for Nearly Non-Stationary Series

1
School of Mathematics and Statistics, Xi’an Jiaotong University, Xi’an 710049, China
2
Department of Finance, College of Management, Asia University, Taichung City 41354, Taiwan
3
Department of Bioinformatics and Medical Engineering, College of Information and Electrical Engineering, Asia University, Taichung City 41354, Taiwan
4
Discipline of Business Analytics, University of Sydney Business School, Darlington, NSW 2006, Australia
5
Econometric Institute, Erasmus School of Economics, Erasmus University Rotterdam, 3062 PA Rotterdam, The Netherlands
6
Department of Economic Analysis and ICAE, Complutense University of Madrid, 28040 Madrid, Spain
7
Department of Finance, Fintech and Blockchain Research Center, and Big Data Research Center, Asia University, 500, Lioufeng Road, Wufeng, Taichung City 41354, Taiwan
8
Department of Medical Research, China Medical University Hospital, Taichung City 41354, Taiwan
9
Department of Economics and Finance, The Hang Seng University of Hong Kong, Siu Lek Yuen, Hong Kong, China
*
Author to whom correspondence should be addressed.
Shortly after finishing the first draft of this paper, Michael McAleer passed away. Deceased 8 July 2021.
J. Risk Financial Manag. 2021, 14(8), 366; https://doi.org/10.3390/jrfm14080366
Submission received: 27 May 2021 / Revised: 18 July 2021 / Accepted: 28 July 2021 / Published: 9 August 2021

Abstract

:
Literature shows that the regression of independent and (nearly) nonstationary time series could result in spurious outcomes. In this paper, we conjecture that under some situations, the regression of two independent and nearly non-stationary series does not have any spurious problem at all. To check whether our conjecture holds, we set up several situations and conduct simulations to justify our conjecture. Our simulations show that under some situations, the chance that the regressions being spurious is very high for all the cases simulated in our paper. Nonetheless, under some other situations, our simulation shows that the rejection rates are much smaller than the 5% level of significance for all the cases simulated in our paper, implying that our conjecture could hold under some situations that regression of two independent and nearly non-stationary series does not have any spurious problem at all.
JEL Classification:
C01; C15; C22; C58; C60

1. Introduction

Granger and Newbold (1974) and others show that regression of independent (nearly) nonstationary time series could result in spurious outcomes, Pesaran et al. (1999) and others find that a mixed integration of orders; that is, I(0) or I(1), could be cointegrated and the residual is stationary, and Westerlund (2008) documents that many studies commit a Type 1 error by failing to reject the no-cointegration hypothesis. On the other hand, Engle and Granger (1987) establish the relationship between cointegration and error correction models that first suggested in Granger (1981) and develop estimation procedures and tests for the cointegration model. In addition, Phillips (1986) develops an asymptotic theory for regressions of integrated random processes, including the spurious regressions discovered by Granger and Newbold (1974) and the cointegrating regressions developed by Engle and Granger (1987). Entorf (1997) analyses the regression of two independent random walks with drifts and shows that the convergence to pseudo true values applies to the estimation of spurious fixed-effects models. Readers may refer to Ventosa-Santaulária (2009) for an overview of spurious regression.
Is it possible that the regression of two independent and nearly non-stationary series does not have any spurious problem? In this paper, we explore the issue. To explore the problem, we first conjecture that under some situations, regression of two independent and nearly non-stationary series does not have any spurious problem at all. To check whether the conjecture we set holds, we first generate two independent and nearly nonstationary AR(1) processes, X t = α 1 X t 1 + ε t and Y t = α 2 Y t 1 + e t with 0.9 < | α 1 | , | α 2 | < 1 . We then regress Y t on the independent X t to get Y t = α + β X t + u t and check the proportion of rejecting the null hypothesis that the beta ( β ) is zero. We first find that under some situations, consistent with the literature, regressing two independent and (nearly) nonstationary time series could be spurious. Nonetheless, we also find that under some other situations, different from the literature, our results show that the rejection rates are much smaller than the 5% level of significance for all the cases simulated in our paper, implying that under some other situations, regressing nearly nonstationary Y t on independent and nearly nonstationary X t will not get any spurious problem at all as shown in all the cases being simulated in our paper.
The rest of the paper is organized as follows. In Section 2, we state the basic models for the regression and the regression with a spurious problem. In Section 3, we state our model setup and construct the algorithm for the simulation. In Section 4, we discuss our findings from our simulation and the last section concludes.

2. The Model

In this section, we state the basic models for a simple regression and the regression with a spurious problem. We first state the basic simple regression model.

2.1. Linear Regression Model

In this paper, we consider the following simple regression model:
Y t = α + β X t + u t ,
where u t is a random component denoted as error term assumed to be independent and identically distributed (iid) with mean 0 and variance σ u 2 , t = 1 , , T in which T is the sample size, α is the intercept parameter, and β is the slope parameter.
The most important hypothesis for the simple linear regression stated in (1) is
H 0 : β = 0 versus H 1 : β 0 ,
which is used to detect the linear relationship between two variables Y and X. If the null hypothesis H 0 in (2) is true, then for the model in (1), the population mean of Y t is always equal to α for any value of X t , concluding that Y t does not depend on the value of X t and there is no linear relationship between X t and Y t . On the other hand, if the alternative hypothesis H 1 is true, then one could conclude that a change in X t is associated with a change in Y t linearly.
To test whether the null hypothesis H 0 in (2) is true, one could use the following T test:
T = β ^ β S E ( β ^ ) = β ^ S E ( β ^ ) ,
where
β ^ = t = 1 N ( X t X ¯ ) ( Y t Y ¯ ) t = 1 N ( X t X ¯ ) 2 ,
is the estimate of β in which X ¯ = t = 1 T X t / T , Y ¯ = t = 1 T Y t / T , and
S E ( β ^ ) = t = 1 T ( u ^ t 2 ) / ( T 2 ) t = 1 T ( X t X ¯ ) 2
is the standard error of the estimate measuring the accuracy of prediction with u ^ t = ( Y t α ^ β ^ X t ) 2 . It is well known that the test statistic T follows a t-distribution with T 2 degrees of freedom if the null hypothesis H 0 is true. We note that Goodness of fit can refer to specific individual factors (or variables) in the regression equation, whereas R-squared refers to a specific numerical value of the entire regression model.

2.2. Spurious Regression

A spurious relationship is a relationship that does not make sense. In this situation, two or more independent variables could appear to be correlated with the effect of an unseen factor (“confounding factor” or “lurking variable”). That is, goodness-of-fit indicators from regression such as R 2 are likely to be large, implying a valid fit even if the underlying variables are not truly related.
A variable, say, Y t , is said to be integrated of order d, and denoted as Y t = I ( d ) , if it has a stationary, invertible, and stochastic ARMA representation after differencing d times, it is stationary if d < 1 , is non-stationary if d 1 , and is nearly non-stationary if d < 1 but close to 1. Granger and Newbold (1974) and others find that if
Y t I ( 1 ) and X t I ( 1 ) ,
then regression in (1) could be spurious. It means that even when Y t and X t are known to be independent, applying the test statistic T in (3) to an ordinary least square regression could misleadingly indicate a good fit when non-stationary time series data are involved. Granger and Newbold (1974) exhibit an example in which an equation has R 2 = 0.997 and the value of the Durbin–Watson (DW) statistic is 0.093 , implying that when the residuals are strongly autocorrelated in time series regression, the interpretability of the coefficients could be questionable.
Since Granger and Newbold (1974) and others show that regression of independent and nonstationary time series could result in spurious outcomes, it is common to believe that regression of independent and nearly nonstationary time series could also get spurious outcomes. In this paper, we believe that there are some cases in which regression of independent and nearly nonstationary time series may not be spurious. Thus, we set up the following conjecture:
Conjecture 1.
Under some situations, the regression of two independent and nearly non-stationary series does not have any spurious problem at all.
To examine whether the above conjecture could hold in some situations, we will discuss it in the next section.

3. Model Setup and Algorithm

In this section, we first state the model setting of generating two purely independent and nearly nonstationary time series, regressing one of them onto the other, and examining whether the corresponding regression is spurious. We then construct the algorithm for the simulation and discuss our simulation result in the next section.

3.1. Model Setup

We consider the simple linear regression in (1) between two unrelated nearly nonstationary AR(1) series X t and Y t such that
X t = α 1 X t 1 + ε t and Y t = α 2 Y t 1 + e t , with ε t i i d ( 0 , σ ε 2 ) and e t i i d ( 0 , σ e 2 )
in which 0.9 < | α i | < 1 ( i = 1 , 2 ) . For simplicity, we assume that both ε t and e t follow:
f ( a ; p ) 1 σ 1 + a 2 k σ 2 p ( < a < ) ,
where k = 2 p 1 and p 2 . We note that E ( a ) = 0 , V ( a ) = σ 2 , and t = ( k / ν ) ( a / σ ) follows a Student’s t distribution with ν = 2 p 1 degrees of freedom (df). For 1 p < 2 , k is equated to 1 and in this case, σ in (8) is simply a scale parameter. When ν = 1 , it becomes a Cauchy distribution and when ν = , it becomes a normal distribution. Readers may refer to Pötzelberger (1990), Tiku and Wong (1998), Tiku et al. (1999, 2000), Wong and Bian (2005), Fu and Fu (2015), and others to know more properties of AR(1) series.
To simulate X t and Y t properly, without loss of generality, we will consider different factors that could affect the behavior of the time series. First, we consider the distribution of the error terms. We choose a time series that follows the following four different iid error distributions in our study:
Situation 1.
We assume that the distribution of the error terms ε t and e t defined in (7) follow the following situations:
1. 
a standard normal distribution: that is, both ε t and e t ∼ N(0,1);
2. 
a t-distribution with df = 5: that is, both ε t and e t ∼ t(5);
3. 
a t-distribution with df = 2: that is, both ε t and e t ∼ t(2); and
4. 
a t-distribution with df = 1: that is, both ε t and e t follow the standard Cauchy distribution.
Second, we vary the lengths of the times series and simulate a time series with the following four different lengths in our study as stated in the following situations:
Situation 2.
We consider that the lengths of the times series Y t and X t defined in (7) to be: (i) T = 100; (ii) T = 200; (iii) T = 400; and (iv) T = 800.
After deciding the error distribution and the lengths of the AR(1) processes, we now consider the different values of α 1 and α 2 . In our model, since both X t and Y t are nearly nonstationary, we choose 0.9 < | α i | < 1 ( i = 1 , 2 ) and, in particular, we define A = { 0.99 , 0.97 , 0.95 , 0.92 , 0.9 } and A + = { 0.99 , 0.97 , 0.95 , 0.92 , 0.9 } 1 and consider the following values for both α 1 and α 2 as stated in Situations 3 and 4:
Situation 3.
We consider that the values of both α 1 and α 2 such that α 1 A + and α 2 A .
Situation 4.
We consider that the values of both α 1 and α 2 such that α 1 A and α 2 A + .
We note that in this paper, we consider Situations 3 and 4 because when two autoregressive processes in which one is associated to the zero frequency; that is, the AR(1) with a positive coefficient in our paper, and the other is associated to the Nyquist frequency ( π ); that is, the AR(1) with a negative coefficient in our paper that has power at frequency π and completes a cycle every 2 observations, are independent or even asymptotically orthogonal. Readers may refer to Johansen and Schaumburg (1999), Ghysels and Osborn (2001), and del Barrio Castro et al. (2018, 2019) for more information. Readers may also refer to seasonal unit root tests, see, for example, del Barrio Castro et al. (2012) and Smith et al. (2009), and cointegration for processes integrated at different frequencies, see, for example, del Barrio Castro et al. (2020) with properties that are related to the series we are using in our paper.2
With four different error distributions, four different time series lengths, and the above 50 combinations of α 1 and α 2 values as stated in Assumptions 1, 2, 3, and 4, there are in total 800 cases of simulation in our study for the cases when autoregressive coefficients α 1 and α 2 have different signs.
Nevertheless, in this paper, we also study the cases when both autoregressive coefficients α 1 and α 2 are of the same signs, either positive or negative. Thus, we include the following situations in our study:
Situation 5.
We consider that the values of both α 1 and α 2 such that α 1 A + and α 2 A + .
Situation 6.
We consider that the values of both α 1 and α 2 such that α 1 A and α 2 A .

3.2. Algorithm

The two series X t and Y t are generated from independent error terms, and thus, they are expected not to be related. However, Granger and Newbold (1974) and others have shown that regression of independent nonstationary time series could result in spurious outcomes. In this paper, we believe that it is possible that when regressinng independent and nearly nonstationary Y t and X t as shown in Equation (1) may not be spurious under some situations as we stated in Conjecture 1. To check whether Conjecture 1 could hold under some situations, we set the following algorithm for each situation (different error distributions, different time series lengths, different combinations of α 1 and α 2 ) as described in Section 3.1:
Algorithm 1: For each situation (different error distributions, different time series lengths, different combinations of α 1 and α 2 ) as described in Section 3.1, we will conduct the following steps in our simulation:
  • Simulate 10,000 pairs of X t and Y t defined in (7) with coefficients described in Section 3.1.
  • For each pair of simulated X t and Y t , fit model in (1). Thus, in each subcase, we will obtain 10,000 β ^ ’s and 10,000 corresponding p-values.
  • Plot the distribution of β ^ ’s and record the standard errors and t-statistics of β ^ .
  • Use the T test defined in Equation (3) to test whether the null hypothesis H 0 in (2) hold.  H 0 is rejected if p-value for the T test is less than 0.05. Calculate the proportion of significant β ^ ’s or proportion of p-values that are less than 0.05 among the 10,000 fitted linear regression models in each subcase. This proportion is denoted as the rejection rate in this paper.
For each situation (different error distributions, different time series lengths, different combinations of α 1 and α 2 ) as described in Section 3.1, we will conduct simulation as described in Algorithm 1 and discuss the results in the next section.

4. Simulation

We follow Algorithm 1 to conduct simulation for each situation (different error distributions, different time series lengths, different combinations of α 1 and α 2 ) as described in Section 3.1. The simulation helps us to examine whether the T statistic as shown in Equation (3) for the model as shown in Equation (1) follow a Student t-distribution. If X t and Y t are unrelated, the true null hypothesis that all β coefficients are zero should be rejected around 5% of the time at the significance level of 5%. If the T test is good, that is, β ^ ’s follow student t-distribution, the rejection rate should be close to 5%. If the rejection rate is significantly greater than 5%, then we conclude that there exists the spurious problem. In addition, we believe that it is possible that when regressing independent and nearly nonstationary Y t and X t as shown in (1) may not be spurious under some situations as we hypothesized in Conjecture 1. To check whether Conjecture 1 could hold under some situations, we discuss it in this section. We first discuss the results of the simulation for the cases when α 1 and α 2 are of different signs in the next subsection.

4.1. Simulation for the Cases When α 1 and α 2 Are of Different Signs

We first analyze cases as stated in Situation 3 and exhibit the results in Table A1, Table A2, Table A3 and Table A4 displaying in Appendix A that report the rejecting frequency of the T test when ( α 1 , α 2 ) ( A + , A ) . From Table A1, Table A2, Table A3 and Table A4, one can observe that when choosing the values of both α 1 and α 2 as stated in Situation 3 are from | 0.9 | to | 0.99 | , the rejection rate is about 0.0000 for any n and for any error distribution studied in our paper, except the situation when the error term follows a t ( 1 ) in which the rejection rates are close to 0.0004.
We then analyze the cases as stated in Situation 4 and show the results in Table A5, Table A6, Table A7 and Table A8 displaying in Appendix B that report the rejecting frequency of the T test when ( α 1 , α 2 ) ( A , A + ) . Similarly, from Table A5, Table A6, Table A7 and Table A8, one can observe that when choosing values of both α 1 and α 2 as stated in Situation 4 are between | 0.9 | and | 0.99 | , the rejection rate is zero or close to zero for any n and any error distribution studied in our paper.
Our analysis shows that for all the cases when choosing values for ( α 1 , α 2 ) as stated in Situations 3 and 4 and when choosing values of both | α 1 | and | α 2 | are between 0.9 and 0.99 , the rejection rates are much smaller than the 5% level of significance, implying that when ( α 1 , α 2 ) follow Situations 3 and 4 and when both | α 1 | and | α 2 | are between 0.9 and 0.99 , all the corresponding regressions do not encounter any spurious problem for all the cases simulated in our paper, confirming that Conjecture 1 holds. In other words, our analysis shows that when independent Y t and X t follow nearly nonstationary AR(1) model and the autoregressive coefficients α 1 and α 2 have opposite signs, there is no spurious problem in the regression stated in Equation (1) and Conjecture 1 holds.

4.2. Simulation for the Cases When α 1 and α 2 Are of the Same Sign

We turn to examine whether the regression shown in Equation (1) is spurious for the cases when both α 1 and α 2 are of the same signs; that is, both α 1 and α 2 are positive or both are negative. To do so, we follow Algorithm 1 to conduct simulations for the cases when both α 1 and α 2 are positive and both are negative as displayed in Situations 5 and 6 and exhibit the results in Table A9, Table A10, Table A11, Table A12, Table A13, Table A14, Table A15 and Table A16 displaying in Appendices Appendix C and Appendix D, respectively.
We first discuss the cases when both α 1 and α 2 are positive as stated in Situation 5. Compared with the results in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7 and Table A8, all of the rejection rates in Table A9, Table A10, Table A11 and Table A12 are significantly higher than 5% and the rejecting frequency of the T test is higher than 49% for any n and any error distribution studied in our paper, except the situation when the error term follows t ( 1 ) in which the rejection rates is higher than 32%. In addition, as n increases, or either α 1 or α 2 increases, or as the error distributions are further away from normal distribution, the rejecting rate increases even further.
We turn to discuss the cases when both α 1 and α 2 are negative as stated in Situation 6. Similar to the cases when both α 1 and α 2 are positive, when both α 1 and α 2 are negative, the rejecting frequency of the T test is higher than 50% for any n and for any error distribution studied in our paper, except the situation when the error term follows t ( 1 ) in which the rejection rates is higher than 31%. In addition, Similar to the cases when both α 1 and α 2 are positive, as n increases, or either α 1 or α 2 increases, or as the error distributions are further away from normal distribution, the rejecting rate increases even further.
Our analysis shows that, different from all the cases when for α 1 and α 2 are of different signs, for all the cases when α 1 and α 2 are of the same signs, either positive or negative, as stated in Situations 5 and 6, respectively, and when both | α 1 | and | α 2 | are between 0.9 and 0.99 , the rejection rate is much higher than the 5% level of significance for all the cases studied in our paper and it could be higher than 49%, implying that when ( α 1 , α 2 ) follow Situations 5 and 6 and when both | α 1 | and | α 2 | are between 0.9 and 0.99 , the chance that the regressions being spurious is very high for all the cases simulated in our paper, which, in turn, rejects Conjecture 1 for all the cases in Situations 5 and 6.

5. Concluding Remarks

In this paper, we conjecture that under some situations, the regression of two independent and nearly non-stationary series does not have any spurious problem at all. To check whether our conjecture holds, we first generate two independent and nearly nonstationary AR(1) processes, X t = α 1 X t 1 + ε t and Y t = α 2 Y t 1 + e t in which 0.9 < | α 1 | , | α 2 | < 1 . We then regress Y t on independent X t to get Y t = α + β X t + u t and check whether the proportion of rejecting the null hypothesis of the beta ( β ) to be zero. We first find that consistent with the literature that supports the hypothesis of regressing two independent and (nearly) nonstationary time series could be spurious, when both α 1 and α 2 are of the same signs, either positive or negative, and when the values of both | α 1 | and | α 2 | are between 0.9 and 0.99 , the rejection rate is much bigger than the 5% level of significance in all the cases examined in our simulation and it could be higher than 49% in many cases, implying that the chance that the regressions being spurious is very high for all the cases when both α 1 and α 2 are of the same signs.
Nonetheless, for all the cases when for α 1 and α 2 are of different signs, then different from the literature, our results show that when both | α 1 | and | α 2 | are between 0.9 and 0.99 , the rejection rates are much smaller than the 5% level of significance for all the cases studied in our paper, implying that when α 1 and α 2 are of different signs, regressing nearly nonstationary Y t on independent and nearly nonstationary X t will not get any spurious problem at all for all the cases being simulated in our paper.
We note that the literature shows that the regression of independent and (nearly) nonstationary time series could result in spurious outcomes. In this paper, we conjecture that under some situations, regression of two independent and nearly non-stationary series does not have any spurious problem at all, and in this paper, we aim to find some situations that our conjecture could hold. In this paper, we find that when ( α 1 , α 2 ) ( A + , A ) or ( A , A + ) , then our conjecture holds. We note that when ( α 1 , α 2 ) ( A + , A ) or ( A , A + ) , our conjecture holds which does not imply that these are only situations that our conjecture holds. There could have other situations that our conjecture could hold. We leave it to future studies to find other situations that our conjecture could hold. The purpose of our paper is to tell readers that when one finds regression of any two or more time series that do not have any spurious problem, this does not necessarily imply that the series are not independent. Thus, academics and practitioners should conduct some proper tests to show whether the series are independent.
Some academics may wonder whether there are some financial or economic time series that exhibit extreme negative autocorrelations. We believe there could have some financial or economic time series exhibit positive autocorrelations and some exhibit negative autocorrelations. We note that the time period used in our paper may not be daily or monthly, it should be set to fit the nature of the time series. It is well-known that stock returns could be overreacted or underreacted, this means that it could be positively auto-correlated or negatively auto-correlated and the true unobserved stock returns are positively auto-correlated or negatively auto-correlated. Whether they are extreme positively auto-correlated or negatively auto-correlated will depend on particular stocks. In addition, as we have mentioned before, when ( α 1 , α 2 ) ( A + , A ) or ( A , A + ) , our conjecture holds which does not imply that these are the only situations that our conjecture could hold. There may have other situations that our conjecture could hold. Some financial or economic time series could follow other situations that yet to be discovered, and thus, the conjecture could be important not only for statistics, but also for economics and finance. We also note that in our paper, we only consider ( A + , A ) to cover nearly non-stationary series but do not cover the situations A + = 1 and A = 1 . We do not cover the situation A + = 1 because this has been well-studied in the literature. On the other hand, we do not cover the situation A = 1 because this situation, we believe, is of no practice relevance.
We note that as far as we know, this paper is the first paper to discover that under some situations, the regression of two independent and nearly non-stationary series does not have any spurious problem at all. We follow Granger and Newbold (1974) and others to provide simulation results to show our discovery. Academics could follow Phillips (1986), Johansen and Schaumburg (1999), and others to provide formal proof of the finding in our paper to replace Brownian motions by using the OU processes with e x p ( c / T ) to approximate ( 1 + c / T ) .3 We will leave it to further research to develop the theoretical results to explain the phenomena discovered in this paper. We also note that in this paper, we get very good results by using 0.9 < | α | < 1 . One may get good results by using the near-integrated approach. We will leave this to future studies.4 Another problem in our study is that there is a serious problem with under-rejection. Further study could expose this problem and correct the test properly.

Author Contributions

Conceptualization and methodology, Y.C., Y.H. and W.-K.W.; formal analysis, Y.C., Y.H. and W.-K.W.; writing—original draft preparation, Y.C., Y.H., M.M., and W.-K.W.; writing—review and editing, Y.C., Y.H., M.M. and W.-K.W. project administration and funding, Y.C., Y.H., M.M. and W.-K.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been supported by the Asia University, China Medical University Hospital, The Hang Seng University of Hong Kong, Research Grants Council (RGC) of Hong Kong (project number 12500915), and Ministry of Science and Technology (MOST, Project Numbers 106-2410-H-468-002 and 107-2410-H-468-002-MY3), Taiwan.

Acknowledgments

The authors are grateful to the anonymous referees for substantive comments that have significantly improved this manuscript. For financial support, the third author wishes to thank the Australian Research Council and the National Science Council, Ministry of Science and Technology (MOST), Taiwan. The fourth author acknowledges the Research Grants Council of Hong Kong (Project Numbers 12500915), Ministry of Science and Technology, Taiwan (MOST, Project Numbers 106-2410-H-468-002 and 107-2410-H-468-002-MY3), Asia University, China Medical University Hospital, and The Hang Seng University of Hong Kong. The fourth author would also like to thank Robert B. Miller and Howard E. Thompson for their continuous guidance and encouragement.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. (α1, α2) ∈ (A+, A)

Table A1. Rejection rate, error ∼ N(0,1).
Table A1. Rejection rate, error ∼ N(0,1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.90.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.920.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.950.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.970.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.990.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A2. Rejection rate, error ∼ t(5).
Table A2. Rejection rate, error ∼ t(5).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.90.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.920.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.950.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
-0.970.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.990.950.00000.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A3. Rejection rate, error ∼ t(2).
Table A3. Rejection rate, error ∼ t(2).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.00000.00010.00000.00010.0001
0.920.00010.00000.00000.00000.0001
−0.90.950.00000.00010.00000.00000.00000.0000
0.970.00010.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.920.950.00010.00000.00000.00000.00000.0000
0.970.00010.00000.00000.00000.0000
0.990.00010.00000.00000.00000.0000
Caverage0.00010.00000.00000.00000.00000.0000
0.90.00020.00000.00000.00000.0001
0.920.00000.00000.00000.00000.0000
−0.950.950.00000.00000.00000.00000.00000.0000
0.970.00010.00010.00000.00000.0001
0.990.00000.00010.00000.00000.0000
Caverage0.00010.00000.00000.00000.00000.0000
0.90.00010.00000.00000.00000.0000
0.920.00010.00010.00000.00000.0001
−0.970.950.00010.00000.00000.00000.00000.0000
0.970.00030.00010.00010.00000.0001
0.990.00000.00000.00000.00000.0000
Caverage0.00010.00000.00000.00000.00000.0000
0.90.00000.00000.00000.00000.0000
0.920.00000.00000.00000.00000.0000
−0.990.950.00010.00000.00000.00000.00000.0000
0.970.00000.00000.00000.00000.0000
0.990.00000.00000.00000.00000.0000
Caverage0.00010.00000.00000.00000.00000.0000
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A4. Rejection Rate, Error ∼ t(1).
Table A4. Rejection Rate, Error ∼ t(1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.00200.00050.00070.00270.0015
0.920.00100.00060.00060.00130.0009
−0.90.950.00090.00060.00020.00070.00060.0008
0.970.00090.00070.00020.00010.0005
0.990.00090.00010.00000.00030.0003
Caverage0.00110.00050.00030.00100.00080.0008
0.90.00080.00050.00020.00060.0005
0.920.00120.00100.00020.00060.0008
−0.920.950.00040.00040.00030.00030.00040.0004
0.970.00050.00040.00000.00010.0003
0.990.00050.00020.00010.00000.0002
Caverage0.00070.00050.00020.00030.00040.0004
0.90.00080.00070.00020.00050.0006
0.920.00100.00020.00010.00030.0004
−0.950.950.00110.00030.00020.00010.00040.0004
0.970.00070.00060.00000.00010.0004
0.990.00070.00010.00000.00010.0002
Caverage0.00090.00040.00010.00020.00040.0004
0.90.00080.00060.00030.00030.0005
0.920.00090.00070.00010.00010.0005
−0.970.950.00140.00010.00010.00020.00050.0004
0.970.00020.00020.00010.00020.0002
0.990.00100.00020.00020.00000.0004
Caverage0.00090.00040.00020.00020.00040.0004
0.90.00110.00030.00040.00000.0005
0.920.00050.00020.00010.00000.0002
−0.990.950.00030.00030.00000.00010.00020.0003
0.970.00090.00030.00010.00000.0003
0.990.00050.00000.00000.00000.0001
Caverage0.00070.00020.00010.00000.00030.0003
Oaverage0.00080.00040.00020.00030.00040.0004
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.

Appendix B. (α1, α2) ∈ (A, A+)

Table A5. Rejection rate, error ∼ N(0,1).
Table A5. Rejection rate, error ∼ N(0,1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.9−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.92−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.95−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.97−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.99−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
Oaverage0.00000.00000.00000.00000.00000.0000
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A6. Rejection rate, error ∼ t(5).
Table A6. Rejection rate, error ∼ t(5).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.9−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00010.00000.00000.00000.0000
0.92−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.95−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.97−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00010.00000.00000.00000.0000
−0.920.00000.00000.00000.00000.0000
0.99−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
Oaverage0.00000.00000.00000.00000.00000.0000
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A7. Rejection rate, error ∼ t(2).
Table A7. Rejection rate, error ∼ t(2).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.00000.00020.00000.00010.0001
−0.920.00000.00000.00000.00000.0000
0.9−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00010.00000.00000.0000
Caverage0.00000.00010.00000.00000.00000.0000
−0.90.00010.00000.00000.00010.0001
−0.920.00000.00010.00000.00000.0000
0.92−0.950.00010.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00000.00000.00000.0000
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00000.00000.00000.00000.0000
−0.920.00000.00010.00000.00000.0000
0.95−0.950.00020.00000.00000.00000.00010.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00000.00010.00000.00010.0001
Caverage0.00000.00000.00000.00000.00000.0000
−0.90.00010.00000.00000.00000.0000
−0.920.00020.00010.00000.00000.0001
0.97−0.950.00000.00000.00000.00000.00000.0000
−0.970.00010.00010.00000.00000.0001
−0.990.00010.00000.00000.00000.0000
Caverage0.00010.00000.00000.00000.00000.0000
−0.90.00020.00010.00000.00000.0001
−0.920.00010.00000.00000.00000.0000
0.99−0.950.00000.00000.00000.00000.00000.0000
−0.970.00000.00000.00000.00000.0000
−0.990.00010.00000.00000.00000.0000
Caverage0.00010.00000.00000.00000.00000.0000
Oaverage0.00010.00000.00000.00000.00000.0000
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A8. Rejection rate, error ∼ t(1).
Table A8. Rejection rate, error ∼ t(1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.00140.00060.00090.00190.0012
−0.920.00100.00050.00040.00160.0009
0.9−0.950.00110.00030.00040.00080.00070.0007
−0.970.00080.00060.00000.00010.0004
−0.990.00080.00000.00010.00040.0003
Caverage0.00100.00040.00040.00100.00070.0007
−0.90.00090.00080.00020.00070.0007
−0.920.00120.00100.00020.00070.0008
0.92−0.950.00040.00030.00030.00010.00030.0004
−0.970.00040.00050.00000.00020.0003
−0.990.00050.00030.00010.00000.0002
Caverage0.00070.00060.00020.00030.00040.0004
−0.90.00090.00060.00020.00050.0006
−0.920.00170.00040.00020.00010.0006
0.95−0.950.00090.00030.00030.00010.00040.0004
−0.970.00080.00040.00010.00010.0004
−0.990.00040.00040.00000.00010.0002
Caverage0.00090.00040.00020.00020.00040.0004
−0.90.00080.00070.00050.00020.0006
−0.920.00090.00060.00000.00010.0004
0.97−0.950.00140.00040.00020.00010.00050.0004
−0.970.00040.00020.00000.00020.0002
−0.990.00120.00040.00010.00000.0004
Caverage0.00090.00050.00020.00010.00040.0004
−0.90.00100.00070.00020.00010.0005
−0.920.00080.00060.00020.00000.0004
0.99−0.950.00030.00020.00030.00010.00020.0003
−0.970.00070.00040.00020.00000.0003
−0.990.00030.00010.00020.00010.0002
Caverage0.00060.00040.00020.00010.00030.0003
Oaverage0.00080.00050.00020.00030.00050.0005
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.

Appendix C. (α1, α2) ∈ (A+, A+)

Table A9. Rejection rate, error ∼ N(0,1).
Table A9. Rejection rate, error ∼ N(0,1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.51050.51730.51910.52020.5168
0.920.51280.53790.54460.55520.5376
0.90.950.55280.57040.58000.57640.56990.5669
0.970.56400.58870.60540.60580.5910
0.990.58470.62030.63580.63700.6195
Caverage0.54500.56690.57700.57890.56690.5669
0.90.51930.53910.54220.54450.5363
0.920.53920.55960.57560.57440.5622
0.920.950.56330.59950.61200.61020.59630.5924
0.970.58800.62390.63680.64580.6236
0.990.60280.64680.65840.66600.6435
Caverage0.56250.59380.60500.60820.59240.5924
0.90.54510.58050.57950.58770.5732
0.920.56600.59700.61310.61640.5981
0.950.950.61200.63910.65540.65040.63920.6354
0.970.62520.67190.68530.68890.6678
0.990.64860.69810.72210.72480.6984
Caverage0.59940.63730.65110.65360.63540.6354
0.90.55910.59650.61130.60060.5919
0.920.58330.61070.63670.63500.6164
0.970.950.62820.66250.68970.68690.66680.6623
0.970.64330.69840.72860.73000.7001
0.990.67670.72880.76140.77850.7364
Caverage0.61810.65940.68550.68620.66230.6623
0.90.57750.60510.62740.63850.6121
0.920.61060.64620.66250.66770.6468
0.990.950.64540.69570.72250.72400.69690.6937
0.970.66690.72510.76320.76870.7310
0.990.71030.77300.80710.83700.7819
Caverage0.64210.68900.71650.72720.69370.6937
Oaverage0.59340.62930.64700.65080.63010.6301
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A10. Rejection rate, error ∼ t(5).
Table A10. Rejection rate, error ∼ t(5).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.51040.51520.52330.51510.5160
0.920.52280.54200.54470.54420.5384
0.90.950.55140.55880.57520.58770.56830.5667
0.970.56080.59790.60920.60650.5936
0.990.58050.61340.63220.64350.6174
Caverage0.54520.56550.57690.57940.56670.5667
0.90.52070.54260.54900.54900.5403
0.920.53840.55980.56970.56900.5592
0.920.950.57570.59390.60420.60330.59430.5922
0.970.58340.63360.63490.63190.6210
0.990.61020.64540.66070.66900.6463
Caverage0.56570.59510.60370.60440.59220.5922
0.90.55390.57240.59250.59020.5773
0.920.57290.59840.61140.61200.5987
0.950.950.60880.64130.64890.64840.63690.6438
0.970.61980.66650.68270.69000.6648
0.990.64270.70390.71590.72350.6965
Caverage0.59960.63650.65030.65280.63480.6438
0.90.57230.59620.60930.60420.5955
0.920.59500.63040.63600.63920.6252
0.970.950.62910.67530.69310.68590.67090.6656
0.970.64480.70610.72730.72950.7019
0.990.66690.73280.76390.77390.7344
Caverage0.62160.66820.68590.68650.66560.6656
0.90.57520.61140.62520.63050.6106
0.920.60560.63940.66010.66950.6437
0.990.950.64780.70000.70930.72030.69440.6941
0.970.68290.73360.76570.77840.7402
0.990.70920.77470.80990.83400.7820
Caverage0.64410.69180.71400.72650.69410.6941
Oaverage0.59520.63140.64620.64990.63070.6307
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A11. Rejection rate, error ∼ t(2).
Table A11. Rejection rate, error ∼ t(2).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.50750.51340.50770.49650.5063
0.920.51980.52900.53100.53370.5284
0.90.950.55020.58100.57910.57980.57250.5631
0.970.57050.59410.60000.59990.5911
0.990.58090.62370.63120.63380.6174
Caverage0.54580.56820.56980.56870.56310.5631
0.90.51600.53210.54580.53140.5313
0.920.54140.55570.56040.55910.5542
0.920.950.57300.59410.60580.60050.59340.5894
0.970.59760.61950.62630.63810.6204
0.990.62020.64600.66180.66270.6477
Caverage0.56960.58950.60000.59840.58940.5894
0.90.55930.57970.56960.57360.5706
0.920.58310.59270.59590.59820.5925
0.950.950.60640.64300.64750.64490.63550.6341
0.970.63540.67090.68580.68820.6701
0.990.65120.70450.72050.73210.7021
Caverage0.60710.63820.64390.64740.63410.6341
0.90.57360.59610.59930.60700.5940
0.920.59080.62600.63750.63080.6213
0.970.950.64070.66500.68680.69080.67080.6660
0.970.66130.70480.72840.73400.7071
0.990.67610.73920.75880.77350.7369
Caverage0.62850.66620.68220.68720.66600.6660
0.90.58660.62900.63150.63860.6214
0.920.61750.65080.66650.66760.6506
0.990.950.64940.69860.72360.73240.70100.6997
0.970.68910.73370.76660.77980.7423
0.990.71130.76870.81460.83800.7832
Caverage0.65080.69620.72060.73130.69970.6997
Oaverage0.60040.63170.64330.64660.63050.6305
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A12. Rejection rate, error ∼ t(1).
Table A12. Rejection rate, error ∼ t(1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
0.90.48040.44000.38740.32510.4082
0.920.50240.47030.41260.35880.4360
0.90.950.55450.52870.48060.42590.49740.5000
0.970.58830.57640.54890.48020.5485
0.990.60140.62500.62800.58580.6101
Caverage0.54540.52810.49150.43520.50000.5000
0.90.51480.48430.41740.35110.4419
0.920.54030.51290.45790.38360.4737
0.920.950.57030.55920.52330.45880.52790.5344
0.970.60720.62530.58280.52540.5852
0.990.62560.65970.66020.62810.6434
Caverage0.57160.56830.52830.46940.53440.5344
0.90.54660.52940.47870.41280.4919
0.920.57420.57620.51680.45730.5311
0.950.950.61910.62130.58840.53110.59000.5906
0.970.64060.66380.65190.59840.6387
0.990.65720.70630.72910.71280.7014
Caverage0.60750.61940.59300.54250.59060.5906
0.90.58020.57920.54150.48670.5469
0.920.60420.62310.58530.52400.5842
0.970.950.64350.67010.65160.60390.64230.6401
0.970.65240.71120.70700.67610.6867
0.990.69070.74150.76810.76250.7407
Caverage0.63420.66500.65070.61060.64010.6401
0.90.60170.63460.62520.60240.6160
0.920.62820.65110.66430.64140.6463
0.990.950.65620.69830.71580.69900.69230.6951
0.970.68320.74250.76510.76070.7379
0.990.71500.76730.82130.82950.7833
Caverage0.65690.69880.71830.70660.69510.6951
Oaverage0.60310.61590.59640.55290.59210.5921
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A+ = {0.99, 0.97, 0.95, 0.92, 0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.

Appendix D. (α1, α2) ∈ (A, A)

Table A13. Rejection rate, error ∼ N(0,1).
Table A13. Rejection rate, error ∼ N(0,1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.51730.51960.52280.52190.5204
−0.920.54130.54690.54440.54190.5436
−0.9−0.950.57590.58520.58600.57830.58140.5761
−0.970.60100.60060.59980.60840.6025
−0.990.62340.63600.63610.63480.6326
Caverage0.57180.57770.57780.57710.57610.5761
−0.90.54070.54050.54350.54340.5420
−0.920.56370.56910.56870.57300.5686
−0.92−0.950.59450.60850.61010.60990.60580.6034
−0.970.62550.63450.63300.63880.6330
−0.990.66000.66350.67460.67240.6676
Caverage0.59690.60320.60600.60750.60340.6034
−0.90.57010.57550.58110.57550.5756
−0.920.59480.61220.60340.60770.6045
−0.95−0.950.64450.65820.65780.65680.65430.6478
−0.970.67300.68680.69090.68670.6844
−0.990.70850.72170.72390.72670.7202
Caverage0.63820.65090.65140.65070.64780.6478
−0.90.60380.60320.61490.61070.6082
−0.920.62840.63940.64030.65240.6401
−0.97−0.950.68500.68620.68960.69390.68870.6868
−0.970.71720.73050.72680.72740.7255
−0.990.75460.76620.77810.78660.7714
Caverage0.67780.68510.68990.69420.68680.6868
−0.90.62540.62290.64090.64040.6324
−0.920.65880.66070.65950.67410.6633
−0.99−0.950.70590.72160.73280.72970.72250.7237
−0.970.75210.76570.77900.78670.7709
−0.990.80100.83280.83650.84800.8296
Caverage0.70860.72070.72970.73580.72370.7237
Oaverage0.63870.64750.65100.65300.64760.6476
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A14. Rejection rate, error ∼ t(5).
Table A14. Rejection rate, error ∼ t(5).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.52250.52480.52060.52010.5220
−0.920.54220.54760.54070.55360.5460
−0.9−0.950.57780.58810.57660.57530.57950.5783
−0.970.60670.60460.61020.60620.6069
−0.990.62920.63590.63620.64730.6372
Caverage0.57570.58020.57690.58050.57830.5783
−0.90.54040.54350.54340.54100.5421
−0.920.56740.57090.57170.57930.5723
−0.92−0.950.59390.59930.59880.61460.60170.6045
−0.970.63540.64460.63520.64760.6407
−0.990.65500.66700.67040.66950.6655
Caverage0.59840.60510.60390.61040.60450.6045
−0.90.57220.57340.58520.58190.5782
−0.920.60480.60940.60970.61230.6091
−0.95−0.950.64340.64990.65980.64390.64930.6501
−0.970.68410.68880.69330.68980.6890
−0.990.71030.72630.73650.72650.7249
Caverage0.64300.64960.65690.65090.65010.6501
−0.90.60330.60250.60840.60870.6057
−0.920.63630.62470.63390.63960.6336
−0.97−0.950.68010.69140.69170.69160.68870.6845
−0.970.71610.72730.72940.72950.7256
−0.990.75020.76910.77510.78110.7689
Caverage0.67720.68300.68770.69010.68450.6845
−0.90.63340.63240.64170.63490.6356
−0.920.66240.66380.66880.67200.6668
−0.99−0.950.70570.72180.72960.72990.72180.7248
−0.970.75790.76880.77980.77760.7710
−0.990.79670.83200.84020.84570.8287
Caverage0.71120.72380.73200.73200.72480.7248
Oaverage0.64110.64830.65150.65280.64840.6484
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A15. Rejection rate, error ∼ t(2).
Table A15. Rejection rate, error ∼ t(2).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.50840.52310.51000.50330.5112
−0.920.54750.53190.53500.53160.5365
−0.9−0.950.57800.57370.57930.56910.57500.5721
−0.970.60370.61190.59920.59840.6033
−0.990.61820.64340.64060.63540.6344
Caverage0.57120.57680.57280.56760.57210.5721
−0.90.53710.53200.53410.52350.5317
−0.920.57050.56490.56680.54580.5620
−0.92−0.950.61310.59700.60310.59870.60300.6009
−0.970.62990.64360.63540.62680.6339
−0.990.67330.66920.67110.68140.6738
Caverage0.60480.60130.60210.59520.60090.6009
−0.90.57740.57810.57040.57290.5747
−0.920.59830.60620.59380.60240.6002
−0.95−0.950.65200.65510.65880.65210.65450.6495
−0.970.68400.68880.69170.69020.6887
−0.990.72190.73210.73770.72700.7297
Caverage0.64670.65210.65050.64890.64950.6495
−0.90.59380.60010.60340.59630.5984
−0.920.63370.63960.63280.62120.6318
−0.97−0.950.69300.68710.68990.68730.68930.6846
−0.970.72650.73460.73490.72370.7299
−0.990.76350.77390.78080.77570.7735
Caverage0.68210.68710.68840.68080.68460.6846
−0.90.63590.63910.63610.63340.6361
−0.920.65810.67560.66750.67290.6685
−0.99−0.950.71820.73090.72990.71980.72470.7268
−0.970.76060.77100.78220.77040.7711
−0.990.81470.83370.84240.84280.8334
Caverage0.71750.73010.73160.72790.72680.7268
Oaverage0.64450.64950.64910.64410.64680.6468
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.
Table A16. Rejection rate, error ∼ t(1).
Table A16. Rejection rate, error ∼ t(1).
α 2 α 1 n = 100n = 200n = 400n = 800RaverageSaverage
−0.90.46950.42710.36740.31620.3951
−0.920.49810.44940.39630.34220.4215
−0.9−0.950.54770.51580.45380.40700.48110.4882
−0.970.60060.56260.51760.45920.5350
−0.990.64760.62770.59660.56130.6083
Caverage0.55270.51650.46630.41720.48820.4882
−0.90.49820.45730.40000.33200.4219
−0.920.53490.48500.43000.36410.4535
−0.92−0.950.58570.53570.49220.43300.51170.5174
−0.970.62680.59630.54190.48700.5630
−0.990.66420.66780.62970.58570.6369
Caverage0.58200.54840.49880.44040.51740.5174
−0.90.55490.50450.45310.39360.4765
−0.920.58940.54740.49040.43840.5164
−0.95−0.950.63710.60960.54870.50070.57400.5802
−0.970.68330.65880.62240.54700.6279
−0.990.73140.72800.70550.66050.7064
Caverage0.63920.60970.56400.50800.58020.5802
−0.90.59820.56160.50630.45460.5302
−0.920.62620.59160.54570.49650.5650
−0.97−0.950.68670.65790.61990.55630.63020.6343
−0.970.73890.71060.66810.62320.6852
−0.990.78510.77940.75570.72350.7609
Caverage0.68700.66020.61910.57080.63430.6343
−0.90.64370.63860.59230.55970.6086
−0.920.67210.66430.63930.60280.6446
−0.99−0.950.72950.73230.70470.66440.70770.7121
−0.970.78470.77970.75630.73010.7627
−0.990.83330.85110.83840.82460.8369
Caverage0.73270.73320.70620.67630.71210.7121
Oaverage0.63870.61360.57090.52250.58640.5864
Note: ‘Raverage’ (stands for row average) is the average rejection rate for different values of n for the same case; ‘Caverage’ (stands for column average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} for the same n for the same situation; ‘Saverage’ (stands for situation average) is the average rejection rate for different cases in which α1A = {−0.99, −0.97, −0.95, −0.92, −0.9} and different n for the same situation and ‘overall average’ in the 3 to 6 columns is the overall average for each n for all situations and all cases while ‘overall average’ in the last two column is the overall average for all the cases in the entire table.

Notes

1
We would like to thank the anonymous referee for introducing the notations of A + and A .
2
We thank the anonymous referee for providing us this information.
3
We would like to thank the anonymous referee for giving us this information.
4
We would like to thank the anonymous referee for giving us this information.

References

  1. del Barrio Castro, Tomás, Gianluca Cubadda, and Denise R. Osborn. 2020. On cointegration for processes integrated at different frequencies. In CEIS Working Paper No. 502. Glasgow: CEIS. [Google Scholar]
  2. del Barrio Castro, Tomás, Denise R. Osborn, and A. M. Robert Taylor. 2012. On Augmented Hegy Tests For Seasonal Unit Roots. Econometric Theory 28: 1121–43. [Google Scholar] [CrossRef] [Green Version]
  3. del Barrio Castro, Tomás, Paulo M. M. Rodrigues, and A. M. Robert Taylor. 2018. Semi-Parametric Seasonal Unit Root Tests. Econometric Theory 34: 447–76. [Google Scholar] [CrossRef] [Green Version]
  4. del Barrio Castro, Tomás, Paulo M. M. Rodrigues, and A. M. Robert Taylor. 2019. Temporal Aggregation of Seasonally Near-Integrated Processes. Journal of Time Series Analysis 40: 872–86. [Google Scholar] [CrossRef] [Green Version]
  5. Engle, Robert F., and Clive W. J. Granger. 1987. Co-integration and error correction: Representation, estimation and testing. Econometrica 55: 251–76. [Google Scholar] [CrossRef]
  6. Entorf, Horst. 1997. Random walks with drifts: Nonsense regression and spurious fixed-effect estimation. Journal of Econometrics 80: 287–96. [Google Scholar] [CrossRef]
  7. Fu, Ke-Ang, and Xiaoyong Fu. 2015. Asymptotics for the random coefficient first-order autoregressive model with possibly heavy-tailed innovations. Journal of Computational and Applied Mathematics 285: 116–24. [Google Scholar] [CrossRef]
  8. Ghysels, Eric, and Denise R. Osborn. 2001. The Econometric Analysis of Seasonal Time Series. Cambridge: Cambridge University Press. [Google Scholar]
  9. Granger, Clive. 1981. Some properties of time series data and their use in econometric model specification. Journal of Econometrics 16: 121–30. [Google Scholar] [CrossRef]
  10. Granger, Clive W. J., and Paul Newbold. 1974. Spurious Regressions in Econometrics. Journal of Econometrics 2: 111–20. [Google Scholar] [CrossRef] [Green Version]
  11. Johansen, Søren, and Ernst Schaumburg. 1999. Likelihood analysis of seasonal cointegration. Journal of Econometrics 88: 301–39. [Google Scholar] [CrossRef] [Green Version]
  12. Pesaran, M. Hashem, Yongcheol Shin, and Ron J. Smith. 1999. Pooled mean group estimation of dynamic heterogeneous panels. Journal of the American Statistical Association 94: 621–34. [Google Scholar] [CrossRef]
  13. Phillips, Peter C. B. 1986. Understanding spurious regressions in econometrics. Journal of Econometrics 33: 311–40. [Google Scholar] [CrossRef] [Green Version]
  14. Westerlund, Joakim. 2008. Panel cointegration tests of the Fisher effect. Journal of Applied Econometrics 23: 193–233. [Google Scholar] [CrossRef]
  15. Pötzelberger, Klaus. 1990. A characterization of random-coefficient AR(1) models. Stochastic Processes and Their Applications 34: 171–80. [Google Scholar] [CrossRef] [Green Version]
  16. Smith, Richard J., A. M. Robert Taylor, and Tomas del Barrio Castro. 2009. Regression-Based Seasonal Unit Root Tests. Econometric Theory 25: 527–60. [Google Scholar] [CrossRef] [Green Version]
  17. Tiku, Moti L., and Wing-Keung Wong. 1998. Testing for unit root in AR(1) model using three and four moment approximations. Communications in Statistics: Simulation and Computation 27: 185–98. [Google Scholar] [CrossRef] [Green Version]
  18. Tiku, Moti L., Wing Keung Wong, and Guorui Bian. 1999. Time series models with asymmetric innovations. Communications in Statistics: Theory and Methods 28: 1331–60. [Google Scholar] [CrossRef]
  19. Tiku, Moti L., Wing-Keung Wong, David C. Vaughan, and Guorui Bian. 2000. Time series models in non-normal situations: Symmetric innovations. Journal of Time Series Analysis 21: 571–96. [Google Scholar] [CrossRef]
  20. Ventosa-Santaulária, Daniel. 2009. Spurious Regression. Journal of Probability and Statistics 1: 1–27. [Google Scholar] [CrossRef]
  21. Wong, Wing-Keung, and Guorui Bian. 2005. Estimating Parameters in Autoregressive Models with asymmetric innovations. Statistics and Probability Letters 71: 61–70. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cheng, Y.; Hui, Y.; McAleer, M.; Wong, W.-K. Spurious Relationships for Nearly Non-Stationary Series. J. Risk Financial Manag. 2021, 14, 366. https://doi.org/10.3390/jrfm14080366

AMA Style

Cheng Y, Hui Y, McAleer M, Wong W-K. Spurious Relationships for Nearly Non-Stationary Series. Journal of Risk and Financial Management. 2021; 14(8):366. https://doi.org/10.3390/jrfm14080366

Chicago/Turabian Style

Cheng, Yushan, Yongchang Hui, Michael McAleer, and Wing-Keung Wong. 2021. "Spurious Relationships for Nearly Non-Stationary Series" Journal of Risk and Financial Management 14, no. 8: 366. https://doi.org/10.3390/jrfm14080366

APA Style

Cheng, Y., Hui, Y., McAleer, M., & Wong, W. -K. (2021). Spurious Relationships for Nearly Non-Stationary Series. Journal of Risk and Financial Management, 14(8), 366. https://doi.org/10.3390/jrfm14080366

Article Metrics

Back to TopTop