Implement Univariate SARIMA

[Ibrahim-Ola]

The current implementation of SARIMA uses exogenous variables (i.e., SARIMAX). We would like to forecast the $y$ variable using its own lagged values and the lagged values of its error term. In other words, let's implement SARIMA.

Ideally, we would implement this in a separate branch and merge once we're done.