- data2<-read.table("q-gdpdef.txt",header=TRUE)
- X<-diff(data2$gdpdef,lag=1,diffences=1) #对原始数据进行差分
- ar1=ar(X,order.max=3,method='yw')#以AR(3)模型估计参数
- ar1$resid
q-gdpdef.txt
(4.55 KB)
对原始数据做了一次差分得到X
对X进行AR(3)的参数估计
ar1$resid应该是残差吧,但是得到的数据里前三个为NA缺失值,怎么会这样?求帮忙!!!!
还有可以帮忙解释下tsdiag这个函数吗?为什么ar1不能作为其参数,即tsdiag(ar1)非法


雷达卡


京公网安备 11010802022788号







