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.

[![enter image description here][1]][1]

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.

  [1]: https://i.sstatic.net/r8Ad8.png