AutoARIMA#

class AutoARIMA(start_p=2, d=None, start_q=2, max_p=5, max_d=2, max_q=5, start_P=1, D=None, start_Q=1, max_P=2, max_D=1, max_Q=2, max_order=5, sp=1, seasonal=True, stationary=False, information_criterion='aic', alpha=0.05, test='kpss', seasonal_test='ocsb', stepwise=True, n_jobs=1, start_params=None, trend=None, method='lbfgs', maxiter=50, offset_test_args=None, seasonal_test_args=None, suppress_warnings=False, error_action='warn', trace=False, random=False, random_state=None, n_fits=10, out_of_sample_size=0, scoring='mse', scoring_args=None, with_intercept=True, update_pdq=True, time_varying_regression=False, enforce_stationarity=True, enforce_invertibility=True, simple_differencing=False, measurement_error=False, mle_regression=True, hamilton_representation=False, concentrate_scale=False)[source]#

Auto-(S)ARIMA(X) forecaster, from pmdarima package.

Includes automated fitting of (S)ARIMA(X) hyper-parameters (p, d, q, P, D, Q).

Exposes pmdarima.arima.AutoARIMA [1] under the sktime interface. Seasonal ARIMA models and exogeneous input is supported, hence this estimator is capable of fitting auto-SARIMA, auto-ARIMAX, and auto-SARIMAX.

The auto-ARIMA algorithm seeks to identify the most optimal parameters for an ARIMA model, settling on a single fitted ARIMA model. This process is based on the commonly-used R function, forecast::auto.arima.

Auto-ARIMA works by conducting differencing tests (i.e., Kwiatkowski–Phillips–Schmidt–Shin, Augmented Dickey-Fuller or Phillips–Perron) to determine the order of differencing, d, and then fitting models within ranges of defined start_p, max_p, start_q, max_q ranges. If the seasonal optional is enabled, auto-ARIMA also seeks to identify the optimal P and Q hyper-parameters after conducting the Canova-Hansen to determine the optimal order of seasonal differencing, D.

In order to find the best model, auto-ARIMA optimizes for a given information_criterion, one of (‘aic’, ‘aicc’, ‘bic’, ‘hqic’, ‘oob’) (Akaike Information Criterion, Corrected Akaike Information Criterion, Bayesian Information Criterion, Hannan-Quinn Information Criterion, or “out of bag”–for validation scoring–respectively) and returns the ARIMA which minimizes the value.

Note that due to stationarity issues, auto-ARIMA might not find a suitable model that will converge. If this is the case, a ValueError will be thrown suggesting stationarity-inducing measures be taken prior to re-fitting or that a new range of order values be selected. Non- stepwise (i.e., essentially a grid search) selection can be slow, especially for seasonal data. Stepwise algorithm is outlined in Hyndman and Khandakar (2008).

Parameters:

start_pint, optional (default=2)

The starting value of p, the order (or number of time lags) of the auto-regressive (“AR”) model. Must be a positive integer.

dint, optional (default=None)

The order of first-differencing. If None (by default), the value will automatically be selected based on the results of the test (i.e., either the Kwiatkowski–Phillips–Schmidt–Shin, Augmented Dickey-Fuller or the Phillips–Perron test will be conducted to find the most probable value). Must be a positive integer or None. Note that if d is None, the runtime could be significantly longer.

start_qint, optional (default=2)

The starting value of q, the order of the moving-average (“MA”) model. Must be a positive integer.

max_pint, optional (default=5)

The maximum value of p, inclusive. Must be a positive integer greater than or equal to start_p.

max_dint, optional (default=2)

The maximum value of d, or the maximum number of non-seasonal differences. Must be a positive integer greater than or equal to d.

max_qint, optional (default=5)

he maximum value of q, inclusive. Must be a positive integer greater than start_q.

start_Pint, optional (default=1)

The starting value of P, the order of the auto-regressive portion of the seasonal model.

Dint, optional (default=None)

The order of the seasonal differencing. If None (by default, the value will automatically be selected based on the results of the seasonal_test. Must be a positive integer or None.

start_Qint, optional (default=1)

The starting value of Q, the order of the moving-average portion of the seasonal model.

max_Pint, optional (default=2)

The maximum value of P, inclusive. Must be a positive integer greater than start_P.

max_Dint, optional (default=1)

The maximum value of D. Must be a positive integer greater than D.

max_Qint, optional (default=2)

The maximum value of Q, inclusive. Must be a positive integer greater than start_Q.

max_orderint, optional (default=5)

Maximum value of p+q+P+Q if model selection is not stepwise. If the sum of p and q is >= max_order, a model will not be fit with those parameters, but will progress to the next combination. Default is 5. If max_order is None, it means there are no constraints on maximum order.

spint, optional (default=1)

The period for seasonal differencing, sp refers to the number of periods in each season. For example, sp is 4 for quarterly data, 12 for monthly data, or 1 for annual (non-seasonal) data. Default is 1. Note that if sp == 1 (i.e., is non-seasonal), seasonal will be set to False. For more information on setting this parameter, see Setting sp. (link to http://alkaline-ml.com/pmdarima/tips_and_tricks.html#period)

seasonalbool, optional (default=True)

Whether to fit a seasonal ARIMA. Default is True. Note that if seasonal is True and sp == 1, seasonal will be set to False.

stationarybool, optional (default=False)

Whether the time-series is stationary and d should be set to zero.

information_criterionstr, optional (default=’aic’)

The information criterion used to select the best ARIMA model. One of pmdarima.arima.auto_arima.VALID_CRITERIA, (‘aic’, ‘bic’, ‘hqic’, ‘oob’).

alphafloat, optional (default=0.05)

Level of the test for testing significance.

teststr, optional (default=’kpss’)

Type of unit root test to use in order to detect stationarity if stationary is False and d is None.

seasonal_teststr, optional (default=’ocsb’)

This determines which seasonal unit root test is used if seasonal is True and D is None.

stepwisebool, optional (default=True)

Whether to use the stepwise algorithm outlined in Hyndman and Khandakar (2008) to identify the optimal model parameters. The stepwise algorithm can be significantly faster than fitting all (or a random subset of) hyper-parameter combinations and is less likely to over-fit the model.

n_jobsint, optional (default=1)

The number of models to fit in parallel in the case of a grid search (stepwise=False). Default is 1, but -1 can be used to designate “as many as possible”.

start_paramsarray-like, optional (default=None)

Starting parameters for ARMA(p,q). If None, the default is given by ARMA._fit_start_params.

trendstr, optional (default=None)

The trend parameter. If with_intercept is True, trend will be used. If with_intercept is False, the trend will be set to a no- intercept value.

methodstr, optional (default=’lbfgs’)

The method determines which solver from scipy.optimize is used, and it can be chosen from among the following strings:

‘newton’ for Newton-Raphson
‘nm’ for Nelder-Mead
‘bfgs’ for Broyden-Fletcher-Goldfarb-Shanno (BFGS)
‘lbfgs’ for limited-memory BFGS with optional box constraints
‘powell’ for modified Powell’s method
‘cg’ for conjugate gradient
‘ncg’ for Newton-conjugate gradient
‘basinhopping’ for global basin-hopping solver

The explicit arguments in fit are passed to the solver, with the exception of the basin-hopping solver. Each solver has several optional arguments that are not the same across solvers. These can be passed as **fit_kwargs

maxiterint, optional (default=50)

The maximum number of function evaluations.

offset_test_argsdict, optional (default=None)

The args to pass to the constructor of the offset (d) test. See pmdarima.arima.stationarity for more details.

seasonal_test_argsdict, optional (default=None)

The args to pass to the constructor of the seasonal offset (D) test. See pmdarima.arima.seasonality for more details.

suppress_warningsbool, optional (default=False)

Many warnings might be thrown inside of statsmodels. If suppress_warnings is True, all of the warnings coming from ARIMA will be squelched.

error_actionstr, optional (default=’warn’)

If unable to fit an ARIMA due to stationarity issues, whether to warn (‘warn’), raise the ValueError (‘raise’) or ignore (‘ignore’). Note that the default behavior is to warn, and fits that fail will be returned as None. This is the recommended behavior, as statsmodels ARIMA and SARIMAX models hit bugs periodically that can cause an otherwise healthy parameter combination to fail for reasons not related to pmdarima.

tracebool, optional (default=False)

Whether to print status on the fits. A value of False will print no debugging information. A value of True will print some. Integer values exceeding 1 will print increasing amounts of debug information at each fit.

randombool, optional (default=’False’)

Similar to grid searches, auto_arima provides the capability to perform a “random search” over a hyper-parameter space. If random is True, rather than perform an exhaustive search or stepwise search, only n_fits ARIMA models will be fit (stepwise must be False for this option to do anything).

random_stateint, long or numpy RandomState, optional (default=None)

The PRNG for when random=True. Ensures replicable testing and results.

n_fitsint, optional (default=10)

If random is True and a “random search” is going to be performed, n_iter is the number of ARIMA models to be fit.

out_of_sample_sizeint, optional (default=0)

The number of examples from the tail of the time series to hold out and use as validation examples. The model will not be fit on these samples, but the observations will be added into the model’s endog and exog arrays so that future forecast values originate from the end of the endogenous vector. For instance:

y = [0, 1, 2, 3, 4, 5, 6]
out_of_sample_size = 2
> Fit on: [0, 1, 2, 3, 4]
> Score on: [5, 6]
> Append [5, 6] to end of self.arima_res_.data.endog values

scoringstr, optional (default=’mse’)

If performing validation (i.e., if out_of_sample_size > 0), the metric to use for scoring the out-of-sample data. One of (‘mse’, ‘mae’)

scoring_argsdict, optional (default=None)

A dictionary of key-word arguments to be passed to the scoring metric.

with_interceptbool, optional (default=True)

Whether to include an intercept term.

update_pdqbool, optional (default=True)

whether to update pdq parameters in update True: model is refit on all data seen so far, potentially updating p,d,q False: model updates only ARIMA coefficients via likelihood, as in pmdarima

Further arguments to pass to the SARIMAX constructor:

- time_varying_regressionboolean, optional (default=False)

Whether or not coefficients on the exogenous regressors are allowed to vary over time.

- enforce_stationarityboolean, optional (default=True)

Whether or not to transform the AR parameters to enforce stationarity in the auto-regressive component of the model. - enforce_invertibility : boolean, optional (default=True) Whether or not to transform the MA parameters to enforce invertibility in the moving average component of the model.

- simple_differencingboolean, optional (default=False)

Whether or not to use partially conditional maximum likelihood estimation for seasonal ARIMA models. If True, differencing is performed prior to estimation, which discards the first \(s D + d\) initial rows but results in a smaller state-space formulation. If False, the full SARIMAX model is put in state-space form so that all datapoints can be used in estimation. Default is False.

- measurement_error: boolean, optional (default=False)

Whether or not to assume the endogenous observations endog were measured with error. Default is False.

- mle_regressionboolean, optional (default=True)

Whether or not to use estimate the regression coefficients for the exogenous variables as part of maximum likelihood estimation or through the Kalman filter (i.e. recursive least squares). If time_varying_regression is True, this must be set to False. Default is True.

- hamilton_representationboolean, optional (default=False)

Whether or not to use the Hamilton representation of an ARMA process (if True) or the Harvey representation (if False). Default is False.

- concentrate_scaleboolean, optional (default=False)

Whether or not to concentrate the scale (variance of the error term) out of the likelihood. This reduces the number of parameters estimated by maximum likelihood by one, but standard errors will then not be available for the scale parameter.

Attributes:

cutoff: Cut-off = “present time” state of forecaster.
fh: Forecasting horizon that was passed.
is_fitted: Whether fit has been called.

See also

ARIMA
StatsForecastAutoARIMA

References

[1]

https://alkaline-ml.com/pmdarima/modules/generated/pmdarima.arima.AutoARIMA.html

Examples

>>> from sktime.datasets import load_airline
>>> from sktime.forecasting.arima import AutoARIMA
>>> y = load_airline()
>>> forecaster = AutoARIMA(
...     sp=12, d=0, max_p=2, max_q=2, suppress_warnings=True
... ) 
>>> forecaster.fit(y)  
AutoARIMA(...)
>>> y_pred = forecaster.predict(fh=[1,2,3])  

Methods

`check_is_fitted`()	Check if the estimator has been fitted.
`clone`()	Obtain a clone of the object with same hyper-parameters.
`clone_tags`(estimator[, tag_names])	Clone tags from another estimator as dynamic override.
`create_test_instance`([parameter_set])	Construct Estimator instance if possible.
`create_test_instances_and_names`([parameter_set])	Create list of all test instances and a list of names for them.
`fit`(y[, X, fh])	Fit forecaster to training data.
`fit_predict`(y[, X, fh, X_pred])	Fit and forecast time series at future horizon.
`get_class_tag`(tag_name[, tag_value_default])	Get a class tag's value.
`get_class_tags`()	Get class tags from the class and all its parent classes.
`get_config`()	Get config flags for self.
`get_fitted_params`([deep])	Get fitted parameters.
`get_param_defaults`()	Get object's parameter defaults.
`get_param_names`()	Get object's parameter names.
`get_params`([deep])	Get a dict of parameters values for this object.
`get_tag`(tag_name[, tag_value_default, ...])	Get tag value from estimator class and dynamic tag overrides.
`get_tags`()	Get tags from estimator class and dynamic tag overrides.
`get_test_params`([parameter_set])	Return testing parameter settings for the estimator.
`is_composite`()	Check if the object is composed of other BaseObjects.
`load_from_path`(serial)	Load object from file location.
`load_from_serial`(serial)	Load object from serialized memory container.
`predict`([fh, X])	Forecast time series at future horizon.
`predict_interval`([fh, X, coverage])	Compute/return prediction interval forecasts.
`predict_proba`([fh, X, marginal])	Compute/return fully probabilistic forecasts.
`predict_quantiles`([fh, X, alpha])	Compute/return quantile forecasts.
`predict_residuals`([y, X])	Return residuals of time series forecasts.
`predict_var`([fh, X, cov])	Compute/return variance forecasts.
`reset`()	Reset the object to a clean post-init state.
`save`([path, serialization_format])	Save serialized self to bytes-like object or to (.zip) file.
`score`(y[, X, fh])	Scores forecast against ground truth, using MAPE (non-symmetric).
`set_config`(**config_dict)	Set config flags to given values.
`set_params`(**params)	Set the parameters of this object.
`set_random_state`([random_state, deep, ...])	Set random_state pseudo-random seed parameters for self.
`set_tags`(**tag_dict)	Set dynamic tags to given values.
`summary`()	Summary of the fitted model.
`update`(y[, X, update_params])	Update cutoff value and, optionally, fitted parameters.
`update_predict`(y[, cv, X, update_params, ...])	Make predictions and update model iteratively over the test set.
`update_predict_single`([y, fh, X, update_params])	Update model with new data and make forecasts.

classmethod get_test_params(parameter_set='default')[source]#

Return testing parameter settings for the estimator.

Parameters:

parameter_setstr, default=”default”: Name of the set of test parameters to return, for use in tests. If no special parameters are defined for a value, will return "default" set.

Returns:

paramsdict or list of dict

check_is_fitted()[source]#

Check if the estimator has been fitted.

Raises:

NotFittedError: If the estimator has not been fitted yet.

clone()[source]#

Obtain a clone of the object with same hyper-parameters.

A clone is a different object without shared references, in post-init state. This function is equivalent to returning sklearn.clone of self.

Raises:

RuntimeError if the clone is non-conforming, due to faulty __init__.

Notes

If successful, equal in value to type(self)(**self.get_params(deep=False)).

clone_tags(estimator, tag_names=None)[source]#

Clone tags from another estimator as dynamic override.

Parameters:

estimatorestimator inheriting from :class:BaseEstimator
tag_namesstr or list of str, default = None: Names of tags to clone. If None then all tags in estimator are used as tag_names.

Returns:

Self: Reference to self.

Notes

Changes object state by setting tag values in tag_set from estimator as dynamic tags in self.

classmethod create_test_instance(parameter_set='default')[source]#

Construct Estimator instance if possible.

Parameters:

parameter_setstr, default=”default”: Name of the set of test parameters to return, for use in tests. If no special parameters are defined for a value, will return “default” set.

Returns:

instanceinstance of the class with default parameters

Notes

get_test_params can return dict or list of dict. This function takes first or single dict that get_test_params returns, and constructs the object with that.

classmethod create_test_instances_and_names(parameter_set='default')[source]#

Create list of all test instances and a list of names for them.

Parameters:

parameter_setstr, default=”default”: Name of the set of test parameters to return, for use in tests. If no special parameters are defined for a value, will return “default” set.

Returns:

objslist of instances of cls: i-th instance is cls(**cls.get_test_params()[i])
nameslist of str, same length as objs: i-th element is name of i-th instance of obj in tests convention is {cls.__name__}-{i} if more than one instance otherwise {cls.__name__}

property cutoff[source]#

Cut-off = “present time” state of forecaster.

Returns:

cutoffpandas compatible index element, or None: pandas compatible index element, if cutoff has been set; None otherwise

property fh[source]#: Forecasting horizon that was passed.

fit(y, X=None, fh=None)[source]#

Fit forecaster to training data.

State change:: Changes state to “fitted”.

Writes to self:

Sets fitted model attributes ending in “_”, fitted attributes are inspectable via get_fitted_params.

Sets self.is_fitted flag to True.

Sets self.cutoff to last index seen in y.

Stores fh to self.fh if fh is passed.

Parameters:

ytime series in sktime compatible data container format.

Time series to which to fit the forecaster.

Individual data formats in sktime are so-called mtype specifications, each mtype implements an abstract scitype.

Series scitype = individual time series, vanilla forecasting. pd.DataFrame, pd.Series, or np.ndarray (1D or 2D)
Panel scitype = collection of time series, global/panel forecasting. pd.DataFrame with 2-level row MultiIndex (instance, time), 3D np.ndarray (instance, variable, time), list of Series typed pd.DataFrame
Hierarchical scitype = hierarchical collection, for hierarchical forecasting. pd.DataFrame with 3 or more level row MultiIndex (hierarchy_1, ..., hierarchy_n, time)

For further details on data format, see glossary on mtype. For usage, see forecasting tutorial examples/01_forecasting.ipynb

fhint, list, np.array or ForecastingHorizon, optional (default=None)

The forecasting horizon encoding the time stamps to forecast at. If self.get_tag("requires-fh-in-fit") is True, must be passed in fit, not optional

Xtime series in sktime compatible format, optional (default=None).

Exogeneous time series to fit the model to. Should be of same scitype (Series, Panel, or Hierarchical) as y. If self.get_tag("X-y-must-have-same-index"), X.index must contain y.index.

Returns:

selfReference to self.

fit_predict(y, X=None, fh=None, X_pred=None)[source]#

Fit and forecast time series at future horizon.

Same as fit(y, X, fh).predict(X_pred). If X_pred is not passed, same as fit(y, fh, X).predict(X).

State change:: Changes state to “fitted”.

Writes to self:

Sets fitted model attributes ending in “_”, fitted attributes are inspectable via get_fitted_params.

Sets self.is_fitted flag to True.

Sets self.cutoff to last index seen in y.

Stores fh to self.fh.

Parameters:

ytime series in sktime compatible data container format

Time series to which to fit the forecaster.

Individual data formats in sktime are so-called mtype specifications, each mtype implements an abstract scitype.

Series scitype = individual time series, vanilla forecasting. pd.DataFrame, pd.Series, or np.ndarray (1D or 2D)
Panel scitype = collection of time series, global/panel forecasting. pd.DataFrame with 2-level row MultiIndex (instance, time), 3D np.ndarray (instance, variable, time), list of Series typed pd.DataFrame
Hierarchical scitype = hierarchical collection, for hierarchical forecasting. pd.DataFrame with 3 or more level row MultiIndex (hierarchy_1, ..., hierarchy_n, time)

For further details on data format, see glossary on mtype. For usage, see forecasting tutorial examples/01_forecasting.ipynb

fhint, list, np.array or ForecastingHorizon (not optional)

The forecasting horizon encoding the time stamps to forecast at.

Xtime series in sktime compatible format, optional (default=None).

Exogeneous time series to fit the model to. Should be of same scitype (Series, Panel, or Hierarchical) as y. If self.get_tag("X-y-must-have-same-index"), X.index must contain y.index.

X_predtime series in sktime compatible format, optional (default=None)

Exogeneous time series to use in prediction. If passed, will be used in predict instead of X. Should be of same scitype (Series, Panel, or Hierarchical) as y in fit. If self.get_tag("X-y-must-have-same-index"), X.index must contain fh index reference.

Returns:

y_predtime series in sktime compatible data container format: Point forecasts at fh, with same index as fh. y_pred has same type as the y that has been passed most recently: Series, Panel, Hierarchical scitype, same format (see above)

classmethod get_class_tag(tag_name, tag_value_default=None)[source]#

Get a class tag’s value.

Does not return information from dynamic tags (set via set_tags or clone_tags) that are defined on instances.

Parameters:

tag_namestr: Name of tag value.
tag_value_defaultany: Default/fallback value if tag is not found.

Returns:

tag_value: Value of the tag_name tag in self. If not found, returns tag_value_default.

classmethod get_class_tags()[source]#

Get class tags from the class and all its parent classes.

Retrieves tag: value pairs from _tags class attribute. Does not return information from dynamic tags (set via set_tags or clone_tags) that are defined on instances.

Returns:

collected_tagsdict: Dictionary of class tag name: tag value pairs. Collected from _tags class attribute via nested inheritance.

get_config()[source]#

Get config flags for self.

Returns:

config_dictdict: Dictionary of config name : config value pairs. Collected from _config class attribute via nested inheritance and then any overrides and new tags from _onfig_dynamic object attribute.

get_fitted_params(deep=True)[source]#

Get fitted parameters.

State required:: Requires state to be “fitted”.

Parameters:

deepbool, default=True

Whether to return fitted parameters of components.

If True, will return a dict of parameter name : value for this object, including fitted parameters of fittable components (= BaseEstimator-valued parameters).
If False, will return a dict of parameter name : value for this object, but not include fitted parameters of components.

Returns:

fitted_paramsdict with str-valued keys

Dictionary of fitted parameters, paramname : paramvalue keys-value pairs include:

always: all fitted parameters of this object, as via get_param_names values are fitted parameter value for that key, of this object
if deep=True, also contains keys/value pairs of component parameters parameters of components are indexed as [componentname]__[paramname] all parameters of componentname appear as paramname with its value
if deep=True, also contains arbitrary levels of component recursion, e.g., [componentname]__[componentcomponentname]__[paramname], etc

classmethod get_param_defaults()[source]#

Get object’s parameter defaults.

Returns:

default_dict: dict[str, Any]: Keys are all parameters of cls that have a default defined in __init__ values are the defaults, as defined in __init__.

classmethod get_param_names()[source]#

Get object’s parameter names.

Returns:

param_names: list[str]: Alphabetically sorted list of parameter names of cls.

get_params(deep=True)[source]#

Get a dict of parameters values for this object.

Parameters:

deepbool, default=True

Whether to return parameters of components.

If True, will return a dict of parameter name : value for this object, including parameters of components (= BaseObject-valued parameters).
If False, will return a dict of parameter name : value for this object, but not include parameters of components.

Returns:

paramsdict with str-valued keys

Dictionary of parameters, paramname : paramvalue keys-value pairs include:

always: all parameters of this object, as via get_param_names values are parameter value for that key, of this object values are always identical to values passed at construction
if deep=True, also contains keys/value pairs of component parameters parameters of components are indexed as [componentname]__[paramname] all parameters of componentname appear as paramname with its value
if deep=True, also contains arbitrary levels of component recursion, e.g., [componentname]__[componentcomponentname]__[paramname], etc

get_tag(tag_name, tag_value_default=None, raise_error=True)[source]#

Get tag value from estimator class and dynamic tag overrides.

Parameters:

tag_namestr: Name of tag to be retrieved
tag_value_defaultany type, optional; default=None: Default/fallback value if tag is not found
raise_errorbool: whether a ValueError is raised when the tag is not found

Returns:

tag_valueAny: Value of the tag_name tag in self. If not found, returns an error if raise_error is True, otherwise it returns tag_value_default.

Raises:

ValueError if raise_error is True i.e. if tag_name is not in
self.get_tags().keys()

get_tags()[source]#

Get tags from estimator class and dynamic tag overrides.

Returns:

collected_tagsdict: Dictionary of tag name : tag value pairs. Collected from _tags class attribute via nested inheritance and then any overrides and new tags from _tags_dynamic object attribute.

is_composite()[source]#

Check if the object is composed of other BaseObjects.

A composite object is an object which contains objects, as parameters. Called on an instance, since this may differ by instance.

Returns:

composite: bool: Whether an object has any parameters whose values are BaseObjects.

property is_fitted[source]#: Whether fit has been called.

classmethod load_from_path(serial)[source]#

Load object from file location.

Parameters:

serialresult of ZipFile(path).open(“object)

Returns:

deserialized self resulting in output at path, of cls.save(path)

classmethod load_from_serial(serial)[source]#

Load object from serialized memory container.

Parameters:

serial1st element of output of cls.save(None)

Returns:

deserialized self resulting in output serial, of cls.save(None)

predict(fh=None, X=None)[source]#

Forecast time series at future horizon.

State required:: Requires state to be “fitted”, i.e., self.is_fitted=True.

Accesses in self:

Fitted model attributes ending in “_”.

self.cutoff, self.is_fitted

Writes to self:: Stores fh to self.fh if fh is passed and has not been passed previously.

Parameters:

fhint, list, np.array or ForecastingHorizon, optional (default=None): The forecasting horizon encoding the time stamps to forecast at. Should not be passed if has already been passed in fit. If has not been passed in fit, must be passed, not optional
Xtime series in sktime compatible format, optional (default=None): Exogeneous time series to use in prediction. Should be of same scitype (Series, Panel, or Hierarchical) as y in fit. If self.get_tag("X-y-must-have-same-index"), X.index must contain fh index reference.

Returns:

y_predtime series in sktime compatible data container format: Point forecasts at fh, with same index as fh. y_pred has same type as the y that has been passed most recently: Series, Panel, Hierarchical scitype, same format (see above)

predict_interval(fh=None, X=None, coverage=0.9)[source]#

Compute/return prediction interval forecasts.

If coverage is iterable, multiple intervals will be calculated.

State required:: Requires state to be “fitted”, i.e., self.is_fitted=True.

Accesses in self:

Fitted model attributes ending in “_”.

self.cutoff, self.is_fitted

Writes to self:: Stores fh to self.fh if fh is passed and has not been passed previously.

Parameters:

fhint, list, np.array or ForecastingHorizon, optional (default=None): The forecasting horizon encoding the time stamps to forecast at. Should not be passed if has already been passed in fit. If has not been passed in fit, must be passed, not optional
Xtime series in sktime compatible format, optional (default=None): Exogeneous time series to use in prediction. Should be of same scitype (Series, Panel, or Hierarchical) as y in fit. If self.get_tag("X-y-must-have-same-index"), X.index must contain fh index reference.
coveragefloat or list of float of unique values, optional (default=0.90): nominal coverage(s) of predictive interval(s)

Returns:

pred_intpd.DataFrame

Column has multi-index: first level is variable name from y in fit,

second level coverage fractions for which intervals were computed.: in the same order as in input coverage.

Third level is string “lower” or “upper”, for lower/upper interval end.

Row index is fh, with additional (upper) levels equal to instance levels,

from y seen in fit, if y seen in fit was Panel or Hierarchical.

Entries are forecasts of lower/upper interval end,

for var in col index, at nominal coverage in second col index, lower/upper depending on third col index, for the row index. Upper/lower interval end forecasts are equivalent to quantile forecasts at alpha = 0.5 - c/2, 0.5 + c/2 for c in coverage.

predict_proba(fh=None, X=None, marginal=True)[source]#

Compute/return fully probabilistic forecasts.

Note: currently only implemented for Series (non-panel, non-hierarchical) y.

State required:: Requires state to be “fitted”, i.e., self.is_fitted=True.

Accesses in self:

Fitted model attributes ending in “_”.

self.cutoff, self.is_fitted

Writes to self:: Stores fh to self.fh if fh is passed and has not been passed previously.

Parameters:

fhint, list, np.array or ForecastingHorizon, optional (default=None): The forecasting horizon encoding the time stamps to forecast at. Should not be passed if has already been passed in fit. If has not been passed in fit, must be passed, not optional
Xtime series in sktime compatible format, optional (default=None): Exogeneous time series to use in prediction. Should be of same scitype (Series, Panel, or Hierarchical) as y in fit. If self.get_tag("X-y-must-have-same-index"), X.index must contain fh index reference.
marginalbool, optional (default=True): whether returned distribution is marginal by time index

Returns:

pred_distsktime BaseDistribution: predictive distribution if marginal=True, will be marginal distribution by time point if marginal=False and implemented by method, will be joint

predict_quantiles(fh=None, X=None, alpha=None)[source]#

Compute/return quantile forecasts.

If alpha is iterable, multiple quantiles will be calculated.

State required:: Requires state to be “fitted”, i.e., self.is_fitted=True.

Accesses in self:

Fitted model attributes ending in “_”.

self.cutoff, self.is_fitted

Writes to self:: Stores fh to self.fh if fh is passed and has not been passed previously.

Parameters:

fhint, list, np.array or ForecastingHorizon, optional (default=None): The forecasting horizon encoding the time stamps to forecast at. Should not be passed if has already been passed in fit. If has not been passed in fit, must be passed, not optional
Xtime series in sktime compatible format, optional (default=None): Exogeneous time series to use in prediction. Should be of same scitype (Series, Panel, or Hierarchical) as y in fit. If self.get_tag("X-y-must-have-same-index"), X.index must contain fh index reference.
alphafloat or list of float of unique values, optional (default=[0.05, 0.95]): A probability or list of, at which quantile forecasts are computed.

Returns:

quantilespd.DataFrame

Column has multi-index: first level is variable name from y in fit,: second level being the values of alpha passed to the function.
Row index is fh, with additional (upper) levels equal to instance levels,: from y seen in fit, if y seen in fit was Panel or Hierarchical.
Entries are quantile forecasts, for var in col index,: at quantile probability in second col index, for the row index.

predict_residuals(y=None, X=None)[source]#

Return residuals of time series forecasts.

Residuals will be computed for forecasts at y.index.

If fh must be passed in fit, must agree with y.index. If y is an np.ndarray, and no fh has been passed in fit, the residuals will be computed at a fh of range(len(y.shape[0]))

State required:: Requires state to be “fitted”. If fh has been set, must correspond to index of y (pandas or integer)
Accesses in self:: Fitted model attributes ending in “_”. self.cutoff, self._is_fitted
Writes to self:: Nothing.

Parameters:

ytime series in sktime compatible data container format

Time series with ground truth observations, to compute residuals to. Must have same type, dimension, and indices as expected return of predict.

If None, the y seen so far (self._y) are used, in particular:

if preceded by a single fit call, then in-sample residuals are produced
if fit requires fh, it must have pointed to index of y in fit

Xtime series in sktime compatible format, optional (default=None)

Exogeneous time series for updating and forecasting Should be of same scitype (Series, Panel, or Hierarchical) as y in fit. If self.get_tag("X-y-must-have-same-index"), X.index must contain both fh index reference and y.index.

Returns:

y_restime series in sktime compatible data container format: Forecast residuals at fh`, with same index as ``fh. y_res has same type as the y that has been passed most recently: Series, Panel, Hierarchical scitype, same format (see above)

predict_var(fh=None, X=None, cov=False)[source]#

Compute/return variance forecasts.

State required:: Requires state to be “fitted”, i.e., self.is_fitted=True.

Accesses in self:

Fitted model attributes ending in “_”.

self.cutoff, self.is_fitted

Writes to self:: Stores fh to self.fh if fh is passed and has not been passed previously.

Parameters:

fhint, list, np.array or ForecastingHorizon, optional (default=None): The forecasting horizon encoding the time stamps to forecast at. Should not be passed if has already been passed in fit. If has not been passed in fit, must be passed, not optional
Xtime series in sktime compatible format, optional (default=None): Exogeneous time series to use in prediction. Should be of same scitype (Series, Panel, or Hierarchical) as y in fit. If self.get_tag("X-y-must-have-same-index"), X.index must contain fh index reference.
covbool, optional (default=False): if True, computes covariance matrix forecast. if False, computes marginal variance forecasts.

Returns:

pred_varpd.DataFrame, format dependent on cov variable

If cov=False:

Column names are exactly those of y passed in fit/update.: For nameless formats, column index will be a RangeIndex.
Row index is fh, with additional levels equal to instance levels,: from y seen in fit, if y seen in fit was Panel or Hierarchical.

Entries are variance forecasts, for var in col index. A variance forecast for given variable and fh index is a predicted

variance for that variable and index, given observed data.

If cov=True:

Column index is a multiindex: 1st level is variable names (as above): 2nd level is fh.
Row index is fh, with additional levels equal to instance levels,: from y seen in fit, if y seen in fit was Panel or Hierarchical.
Entries are (co-)variance forecasts, for var in col index, and: covariance between time index in row and col.

Note: no covariance forecasts are returned between different variables.

reset()[source]#

Reset the object to a clean post-init state.

Using reset, runs __init__ with current values of hyper-parameters (result of get_params). This Removes any object attributes, except:

hyper-parameters = arguments of __init__

object attributes containing double-underscores, i.e., the string “__”

Class and object methods, and class attributes are also unaffected.

Returns:

self: Instance of class reset to a clean post-init state but retaining the current hyper-parameter values.

Notes

Equivalent to sklearn.clone but overwrites self. After self.reset() call, self is equal in value to type(self)(**self.get_params(deep=False))

save(path=None, serialization_format='pickle')[source]#

Save serialized self to bytes-like object or to (.zip) file.

Behaviour: if path is None, returns an in-memory serialized self if path is a file location, stores self at that location as a zip file

saved files are zip files with following contents: _metadata - contains class of self, i.e., type(self) _obj - serialized self. This class uses the default serialization (pickle).

Parameters:

pathNone or file location (str or Path): if None, self is saved to an in-memory object if file location, self is saved to that file location. If:

path=”estimator” then a zip file estimator.zip will be made at cwd. path=”/home/stored/estimator” then a zip file estimator.zip will be stored in /home/stored/.
serialization_format: str, default = “pickle”: Module to use for serialization. The available options are “pickle” and “cloudpickle”. Note that non-default formats might require installation of other soft dependencies.

Returns:

if path is None - in-memory serialized self
if path is file location - ZipFile with reference to the file

score(y, X=None, fh=None)[source]#

Scores forecast against ground truth, using MAPE (non-symmetric).

Parameters:

ypd.Series, pd.DataFrame, or np.ndarray (1D or 2D): Time series to score
fhint, list, array-like or ForecastingHorizon, optional (default=None): The forecasters horizon with the steps ahead to to predict.
Xpd.DataFrame, or 2D np.array, optional (default=None): Exogeneous time series to score if self.get_tag(“X-y-must-have-same-index”), X.index must contain y.index

Returns:

scorefloat: MAPE loss of self.predict(fh, X) with respect to y_test.