(data,lags=20,ax=ax1)
plt.title('原始数据的自相关图')
ax2 = fig.add_subplot(212)
fig = sm.graphics.tsa.plot_pacf...("二阶差分")
plt.show()
fig = plt.figure(figsize=(12,8))
ax1=fig.add_subplot(211)
fig = sm.graphics.tsa.plot_acf...(diff1,lags=20,ax=ax1)
plt.title('一阶差分的自相关图')
ax2 = fig.add_subplot(212)
fig = sm.graphics.tsa.plot_pacf...(diff1,lags=20,ax=ax2)
plt.title('一阶差分的偏相关图')
plt.show()
model1 = sm.tsa.ARMA(data,(7,1)).fit()
print...aic,model1.bic,model1.hqic)
pre = model1.predict('1740','1760',dynamic=True)
print(pre)
model2 = sm.tsa.ARIMA