plt.figure(figsize=(16,8)) model_ARIMA = ARIMA(df_log, order=(7,1,7)) #Using p=7, d=1, q=7 results_ARIMA = model_ARIMA.fit() plt.plot(df_shift) plt.plot(results_ARIMA.fittedvalues, color='red') plt.title('ARIMA Model - RMSE: %.4f'% mean_squared_error(results_ARIMA.fittedvalues,df_shift['Close'], squared=False)) plt.show()