我正在尝试将基于SAS的ARIMA MOdel转换为基于Python的arima模型
SAS上的代码如下:
proc arima data=temp;
identify var=c3 nlag=12;run;
identify var=c3(1);run;
estimate p=1 q=1/*noint*/ method=ml;run;
estimate p=2 q=0 /*noint*/ method=ml;run;
outlier;run;
forecast lead=&lead id=tc_date interval=week printall out=forecast_ARIMA;
run;
退出;
数据临时包含2列:日期和c3评级(我需要通过ARIMA预测)
我为以下代码创建的python代码是:
从statsmodels导入tsa为tsa
p1 = 2
d1 = 0
q1 = 0
temp = tsa.arima_model.ARIMA(temp.c3.values, [p1, d1, q1])
model1 = temp.fit(trend='nc' , disp = False)
pred1 = model1.predict(start=0, end=len(ts_data)-1)
但是我从SA和python进程得到的预测值是完全不同的, 我的流程是否正确?如果是,为什么计算出的值不同,如果不是,python进程应该是什么?
谢谢