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
感觉有点匪夷所思啊……求大神赐教。


雷达卡




京公网安备 11010802022788号







