Awesome
Automated Time Series Models in Python (AtsPy)
Finance Quant Machine Learning
- ML-Quant.com - Automated Research Repository
Easily develop state of the art time series models to forecast univariate data series. Simply load your data and select which models you want to test. This is the largest repository of automated structural and machine learning time series models. Please get in contact if you want to contribute a model. This is a fledgling project, all advice appreciated.
Install
pip install atspy
Automated Models
ARIMA
- Automated ARIMA ModellingProphet
- Modeling Multiple Seasonality With Linear or Non-linear GrowthHWAAS
- Exponential Smoothing With Additive Trend and Additive SeasonalityHWAMS
- Exponential Smoothing with Additive Trend and Multiplicative SeasonalityNBEATS
- Neural basis expansion analysis (now fixed at 20 Epochs)Gluonts
- RNN-based Model (now fixed at 20 Epochs)TATS
- Seasonal and Trend no Box CoxTBAT
- Trend and Box CoxTBATS1
- Trend, Seasonal (one), and Box CoxTBATP1
- TBATS1 but Seasonal Inference is Hardcoded by PeriodicityTBATS2
- TBATS1 With Two Seasonal Periods
Why AtsPy?
- Implements all your favourite automated time series models in a unified manner by simply running
AutomatedModel(df)
. - Reduce structural model errors with 30%-50% by using LightGBM with TSFresh infused features.
- Automatically identify the seasonalities in your data using singular spectrum analysis, periodograms, and peak analysis.
- Identifies and makes accessible the best model for your time series using in-sample validation methods.
- Combines the predictions of all these models in a simple (average) and complex (GBM) ensembles for improved performance.
- Where appropriate models have been developed to use GPU resources to speed up the automation process.
- Easily access all the models by using
am.models_dict_in
for in-sample andam.models_dict_out
for out-of-sample prediction.
AtsPy Progress
- Univariate forecasting only (single column) and only monthly and daily data have been tested for suitability.
- More work ahead; all suggestions and criticisms appreciated, use the issues tab.
- Here is a Google Colab to run the package in the cloud and here you can run all the models.
Documentation by Example
Load Package
from atspy import AutomatedModel
Pandas DataFrame
The data requires strict preprocessing, no periods can be skipped and there cannot be any empty values.
import pandas as pd
df = pd.read_csv("https://raw.githubusercontent.com/firmai/random-assets-two/master/ts/monthly-beer-australia.csv")
df.Month = pd.to_datetime(df.Month)
df = df.set_index("Month"); df
<table class="dataframe">
<thead>
<tr style="text-align: right;">
<th></th>
<th>Megaliters</th>
</tr>
<tr>
<th>Month</th>
<th></th>
</tr>
</thead>
<tbody>
<tr>
<th>1956-01-01</th>
<td>93.2</td>
</tr>
<tr>
<th>1956-02-01</th>
<td>96.0</td>
</tr>
<tr>
<th>1956-03-01</th>
<td>95.2</td>
</tr>
<tr>
<th>1956-04-01</th>
<td>77.1</td>
</tr>
<tr>
<th>1956-05-01</th>
<td>70.9</td>
</tr>
</tbody>
</table>
AutomatedModel
AutomatedModel
- Returns a class instance.forecast_insample
- Returns an in-sample forcasted dataframe and performance.forecast_outsample
- Returns an out-of-sample forcasted dataframe.ensemble
- Returns the results of three different forms of ensembles.models_dict_in
- Returns a dictionary of the fully trained in-sample models.models_dict_out
- Returns a dictionary of the fully trained out-of-sample models.
from atspy import AutomatedModel
model_list = ["HWAMS","HWAAS","TBAT"]
am = AutomatedModel(df = df , model_list=model_list,forecast_len=20 )
Other models to try, add as many as you like; note ARIMA
is slow: ["ARIMA","Gluonts","Prophet","NBEATS", "TATS", "TBATS1", "TBATP1", "TBATS2"]
In-Sample Performance
forecast_in, performance = am.forecast_insample(); forecast_in
<table class="dataframe">
<thead>
<tr style="text-align: right;">
<th></th>
<th>Target</th>
<th>HWAMS</th>
<th>HWAAS</th>
<th>TBAT</th>
</tr>
<tr>
<th>Date</th>
<th></th>
<th></th>
<th></th>
<th></th>
</tr>
</thead>
<tbody>
<tr>
<th>1985-10-01</th>
<td>181.6</td>
<td>161.962148</td>
<td>162.391653</td>
<td>148.410071</td>
</tr>
<tr>
<th>1985-11-01</th>
<td>182.0</td>
<td>174.688055</td>
<td>173.191756</td>
<td>147.999237</td>
</tr>
<tr>
<th>1985-12-01</th>
<td>190.0</td>
<td>189.728744</td>
<td>187.649575</td>
<td>147.589541</td>
</tr>
<tr>
<th>1986-01-01</th>
<td>161.2</td>
<td>155.077205</td>
<td>154.817215</td>
<td>147.180980</td>
</tr>
<tr>
<th>1986-02-01</th>
<td>155.5</td>
<td>148.054292</td>
<td>147.477692</td>
<td>146.773549</td>
</tr>
</tbody>
</table>
performance
<table class="dataframe">
<thead>
<tr style="text-align: right;">
<th></th>
<th>Target</th>
<th>HWAMS</th>
<th>HWAAS</th>
<th>TBAT</th>
</tr>
</thead>
<tbody>
<tr>
<th>rmse</th>
<td>0.000000</td>
<td>17.599400</td>
<td>18.993827</td>
<td>36.538009</td>
</tr>
<tr>
<th>mse</th>
<td>0.000000</td>
<td>309.738878</td>
<td>360.765452</td>
<td>1335.026136</td>
</tr>
<tr>
<th>mean</th>
<td>155.293277</td>
<td>142.399639</td>
<td>140.577496</td>
<td>126.590412</td>
</tr>
</tbody>
</table>
Out-of-Sample Forecast
forecast_out = am.forecast_outsample(); forecast_out
<table class="dataframe">
<thead>
<tr style="text-align: right;">
<th></th>
<th>HWAMS</th>
<th>HWAAS</th>
<th>TBAT</th>
</tr>
<tr>
<th>Date</th>
<th></th>
<th></th>
<th></th>
</tr>
</thead>
<tbody>
<tr>
<th>1995-09-01</th>
<td>137.518755</td>
<td>137.133938</td>
<td>142.906275</td>
</tr>
<tr>
<th>1995-10-01</th>
<td>164.136220</td>
<td>165.079612</td>
<td>142.865575</td>
</tr>
<tr>
<th>1995-11-01</th>
<td>178.671684</td>
<td>180.009560</td>
<td>142.827110</td>
</tr>
<tr>
<th>1995-12-01</th>
<td>184.175954</td>
<td>185.715043</td>
<td>142.790757</td>
</tr>
<tr>
<th>1996-01-01</th>
<td>147.166448</td>
<td>147.440026</td>
<td>142.756399</td>
</tr>
</tbody>
</table>
Ensemble and Model Validation Performance
all_ensemble_in, all_ensemble_out, all_performance = am.ensemble(forecast_in, forecast_out)
all_performance
<table class="dataframe">
<thead>
<tr style="text-align: right;">
<th></th>
<th>rmse</th>
<th>mse</th>
<th>mean</th>
</tr>
</thead>
<tbody>
<tr>
<th>ensemble_lgb__X__HWAMS</th>
<td>9.697588</td>
<td>94.043213</td>
<td>146.719412</td>
</tr>
<tr>
<th>ensemble_lgb__X__HWAMS__X__HWAMS_HWAAS__X__ensemble_ts__X__HWAAS</th>
<td>9.875212</td>
<td>97.519817</td>
<td>145.250837</td>
</tr>
<tr>
<th>ensemble_lgb__X__HWAMS__X__HWAMS_HWAAS</th>
<td>11.127326</td>
<td>123.817378</td>
<td>142.994374</td>
</tr>
<tr>
<th>ensemble_lgb</th>
<td>12.748526</td>
<td>162.524907</td>
<td>156.487208</td>
</tr>
<tr>
<th>ensemble_lgb__X__HWAMS__X__HWAMS_HWAAS__X__ensemble_ts__X__HWAAS__X__HWAMS_HWAAS_TBAT__X__TBAT</th>
<td>14.589155</td>
<td>212.843442</td>
<td>138.615567</td>
</tr>
<tr>
<th>HWAMS</th>
<td>15.567905</td>
<td>242.359663</td>
<td>136.951615</td>
</tr>
<tr>
<th>HWAMS_HWAAS</th>
<td>16.651370</td>
<td>277.268110</td>
<td>135.544299</td>
</tr>
<tr>
<th>ensemble_ts</th>
<td>17.255107</td>
<td>297.738716</td>
<td>163.134079</td>
</tr>
<tr>
<th>HWAAS</th>
<td>17.804066</td>
<td>316.984751</td>
<td>134.136983</td>
</tr>
<tr>
<th>HWAMS_HWAAS_TBAT</th>
<td>23.358758</td>
<td>545.631579</td>
<td>128.785846</td>
</tr>
<tr>
<th>TBAT</th>
<td>39.003864</td>
<td>1521.301380</td>
<td>115.268940</td>
</tr>
</tbody>
</table>
Best Performing In-sample
all_ensemble_in[["Target","ensemble_lgb__X__HWAMS","HWAMS","HWAAS"]].plot()
Future Predictions All Models
all_ensemble_out[["ensemble_lgb__X__HWAMS","HWAMS","HWAAS"]].plot()
And Finally Grab the Models
am.models_dict_in
{'HWAAS': <statsmodels.tsa.holtwinters.HoltWintersResultsWrapper at 0x7f42f7822d30>,
'HWAMS': <statsmodels.tsa.holtwinters.HoltWintersResultsWrapper at 0x7f42f77fff60>,
'TBAT': <tbats.tbats.Model.Model at 0x7f42d3aab048>}
am.models_dict_out
{'HWAAS': <statsmodels.tsa.holtwinters.HoltWintersResultsWrapper at 0x7f9c01309278>,
'HWAMS': <statsmodels.tsa.holtwinters.HoltWintersResultsWrapper at 0x7f9c01309cf8>,
'TBAT': <tbats.tbats.Model.Model at 0x7f9c08f18ba8>}
Follow this link if you want to run the package in the cloud.
AtsPy Future Development
- Additional in-sample validation steps to stop deep learning models from over and underfitting.
- Extra performance metrics like MAPE and MAE.
- Improved methods to select the window length to use in training and calibrating the model.
- Add the ability to accept dirty data, and have the ability to clean it up, interpolation etc.
- Add a function to resample to a larger frequency for big datasets.
- Add the ability to algorithmically select a good enough chunk of a large dataset to balance performance and time to train.
- More internal model optimisation using AIC, BIC an AICC.
- Code annotations for other developers to follow and improve on the work being done.
- Force seasonality stability between in and out of sample training models.
- Make AtsPy less dependency heavy, currently it draws on tensorflow, pytorch and mxnet.
Citations
If you use AtsPy in your research, please consider citing it. I have also written a small report that can be found on SSRN.
BibTeX entry:
@software{atspy,
title = {{AtsPy}: Automated Time Series Models in Python.},
author = {Snow, Derek},
url = {https://github.com/firmai/atspy/},
version = {1.15},
date = {2020-02-17},
}
@misc{atspy,
author = {Snow, Derek},
title = {{AtsPy}: Automated Time Series Models in Python (1.15).},
year = {2020},
url = {https://github.com/firmai/atspy/},
}