RT
我用的程序如下：

library(pdR)
data<-read.csv("VA.csv",header=TRUE)
x<-data[1:156,5]
y<-ts(x, start=c(1998, 1), frequency=12)
y
test<- HEGY.test(wts=y, itsd = c(1, 1, c(1:11)), regvar = 0, selectlags = list(mode = "bic", Pmax = 12))
test


输出的结果如下：

$stats
              Stat. p-value
tpi_1     -1.955926    0.10
tpi_2     -3.652335    0.01
Fpi_3:4   11.453916    0.10
Fpi_5:6    2.687860    0.01
Fpi_7:8    3.780324    0.01
Fpi_9:10  12.279330    0.10
Fpi_11:12 12.948845    0.10
Fpi_2:12  26.647945      NA
Fpi_1:12  24.491568      NA


感觉有点匪夷所思啊……求大神赐教。