I obtained the ARIMA forcast below from the following code:

ARIMAfit = auto.arima((data.ts), seasonal = TRUE) plot(residuals(ARIMAfit)) plot(forecast(ARIMAfit, h = 200))

where data.ts is a ts object with frequency 7 in R. The data set is the daily usage of bicycles in a bike-share program.

As you can imagine, I want a better forecast of what would happen. Do you have any suggestions on how to improve my forecast?

The data passes both the DF and DW tests after I apply the auto.arima.

Let me know if you need the data. I appreciate your help.