[R] What is the CADF test criterion="BIC" report?

p99323005 at ntu.edu.tw p99323005 at ntu.edu.tw
Mon Nov 14 04:34:55 CET 2011


Hello:
   I am a rookie in using R. When I used the unit root test in  
"CADFtest", I got the different t-test statistics between using  
criterion="BIC" and no using criterion. But when I checked the result  
with eviews, I find out that no using criterion is correct. Why after  
using criterion="BIC", I got the different result?


Paul


> data(Canada)

> ADFt <- CADFtest(Canada[,1], max.lag.y = 14, criterion="BIC")

> summary(ADFt)
Augmented DF test
                                             ADF test
t-test statistic:                          -1.389086
p-value:                                    0.855681
Max lag of the diff. dependent variable:    1.000000

Call:
dynlm(formula = formula(model), start = obs.1, end = obs.T)

Residuals:
      Min       1Q   Median       3Q      Max
-0.79726 -0.20587 -0.03332  0.23840  0.70460

Coefficients:
              Estimate Std. Error t value Pr(>|t|)
(Intercept) 24.342321  17.435476   1.396    0.167
trnd         0.009959   0.006941   1.435    0.156
L(y, 1)     -0.026068   0.018767  -1.389    0.856
L(d(y), 1)   0.615983   0.092632   6.650 7.18e-09 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.3533 on 65 degrees of freedom
Multiple R-squared: 0.413,	Adjusted R-squared: 0.3859
F-statistic:    NA on NA and NA DF,  p-value: NA

> ADFt1 <- CADFtest(Canada[,1], max.lag.y =1)

> summary(ADFt1)
Augmented DF test
                                              ADF test
t-test statistic:                          -2.7285715
p-value:                                    0.2282588
Max lag of the diff. dependent variable:    1.0000000

Call:
dynlm(formula = formula(model), start = obs.1, end = obs.T)

Residuals:
      Min       1Q   Median       3Q      Max
-0.84769 -0.24745 -0.02081  0.24187  0.82344

Coefficients:
              Estimate Std. Error t value Pr(>|t|)
(Intercept) 47.661910  17.439021   2.733  0.00776 **
trnd         0.019217   0.007005   2.743  0.00754 **
L(y, 1)     -0.051256   0.018785  -2.729  0.22826
L(d(y), 1)   0.753011   0.075724   9.944 1.61e-15 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.3937 on 78 degrees of freedom
Multiple R-squared: 0.5674,	Adjusted R-squared: 0.5508
F-statistic:    NA on NA and NA DF,  p-value: NA



More information about the R-help mailing list