R时间序列TSA包eacf

1. library(TSA)
   
2. x=c(1944,1369,3707,3577,3460,2626,2846,2933,3163,2015,2359,2649,1162,1579,3817)
   
3. x1=ts(x,frequency=12,start=c(2013,1))
   
4. eacf(x1)

复制代码

eacf()最后会报错，Error in AA %*% t(X) : requires numeric/complex matrix/vector arguments。这是怎么回事呢？

x1是下面的内容，估计你已经看出来了：

           Jan    Feb    Mar    Apr    May   Jun   Jul    Aug   Sep   Oct   Nov   Dec
2013 1944 1369 3707 3577 3460 2626 2846 2933 3163 2015 2359 2649
2014 1162 1579 3817


报错提示这个是“需要数值/复数矩阵/矢量参数的错误”。
原始数据不能直接扔进来，R一般都会被处理成虚拟变量或者哑变量，但是ts不处理，就直接报错了。
而且你这里是8个singularity，函数最大的AR秩默认值是7。

请问楼主知道怎么解决了吗 求告知啊TT

数据量不够，导致AR和MA部分阶数上去后缺少足够的数据，你可以手动将AR和MA的最高阶数调整一下，比如ar.max=3, ma.max=3
计算出来的eacf矩阵就是：
AR/MA
  0 1 2 3
0 o o o o
1 o o o o
2 o o o o
3 o o o o

具体代码是这个：eacf(x1, ar.max = 3, ma.max = 3)
也许调整成其他的参数也可以，但是似乎两个阶数和大于等于7就不行了