ARTICLE AD BOX
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from statsmodels.tsa.seasonal import seasonal_decompose
from statsmodels.tsa.forecasting.theta import ThetaModel
units = pd.read_csv('EDC_units_IN.csv')
units['Date'] = pd.to_datetime(units['Date'], dayfirst=True)
units = units.set_index('Date')
units = units.asfreq('MS')
units.plot()
mod = ThetaModel(units, period=12, method='additive')
mod = mod.fit()
mod.summary()
mod.forecast(12)
mod.forecast_components(12)
