Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Statsmodels arima model returns error

I have started to learn statsmodels package and can't implement basic forecasting with arima.

Error is

ValueError: Given a pandas object and the index does not contain dates

I am trying this as a version:

df = make_df(filename_data)

y = []
x = []

# here I am preparing day by day sequence as that I have inconsistent data and I set 0 to NAN values

start_date = df[date_col].min()
end_date = df[date_col].max()



while start_date <= end_date:

    x.append(start_date)

    try:
        y.append(
            df[df[date_col] == start_date][rev_col].values[0])
    except:
        y.append(0)

    start_date += datetime.timedelta(days=1)

y = np.array(y)
x = np.array(x)

y = pd.TimeSeries(y, index=x)
print(y)
arma_mod = sm.tsa.ARMA(y, order=(2,2))
arma_res = arma_mod.fit(trend='nc', disp=-1)

Before that I tried

df = make_df(filename_data)

y = np.array(df[rev_col])
x = np.array(df[date_col])

y = pd.TimeSeries(y, index=x)

Why is it happening?

The date - revenue data looks OK:

2014-08-04      59477
2014-08-05      29989
2014-08-06      29989
2014-08-07     116116
like image 754
paveltr Avatar asked Mar 13 '23 13:03

paveltr


1 Answers

You can simply transform your DataFrame with as_matrix().

Example working code:

from statsmodels.tsa.arima_model import ARIMA
import numpy as np

def plot_residuals(data, ord=(2, 0, 1)):
    model = ARIMA(endog=data, order=(ord[0], 0, ord[1])).fit()
    plt.plot(model.resid)
    plt.show()

data = np.log(data.values) - np.log(data.values.shift()).to_frame().dropna().as_matrix()
plot_residuals(data, (2, 0, 1))

As statsmodels has many unresolved issues, it will help you only temporarily.

like image 117
Tadas Talaikis Avatar answered Mar 19 '23 07:03

Tadas Talaikis