Python 和 Scipy 中的季节调整

Question

我希望使用 Python 季节性调整每月数据。从这些系列中可以看出：www.emconfidential.com，季节性很高数据的组成部分。我想对此进行调整，以便我可以更好地判断系列趋势是上升还是下降。有人知道如何使用 scipy 或其他 Python 库轻松做到这一点吗？

Answer 1

统计模型可以做到这一点。它们具有基本的季节性分解以及 Census X13 调整的包装。您还可以使用 rpy2 访问 R 的一些优秀的 SA 库。这是 statsmodels 季节性分解：

import pandas as pd
import statsmodels.api as sm
import matplotlib.pyplot as plt
pd.options.display.mpl_style = 'default'
%matplotlib inline

dta = sm.datasets.co2.load_pandas().data.resample("M").fillna(method="ffill")

res = sm.tsa.seasonal_decompose(dta)

fig = res.plot()
fig.set_size_inches(10, 5)
plt.tight_layout()

http://statsmodels.sourceforge.net/0.6 .0/release/version0.6.html

Answer 2

现在有一个软件包似乎正是您正在寻找的！查看 seasonal 包，这里是链接 。我个人觉得非常有用，想知道其他人的想法。

Answer 3

没有神奇的 Python 库可以为你进行季节性调整。执行此类操作的应用程序往往相当大。

您需要自己计算出数学，然后使用 scipy 计算其余部分为你。

Answer 4

我建议使用 Facebook 数据科学团队开发的 Prophet 。它具有 Python+R API，用于时间序列预测，尽管您可以仅使用它将序列分解为其组件（趋势与季节性）。您可以轻松调整和可视化分解：

from fbprophet import Prophet
import numpy as np
import pandas as pd

# Create series
np.random.seed(0)
x = np.arange(0, 10, .285)
y_periodic = np.sin(x*np.pi)
y_random = np.random.normal(size=len(x))
y_trend = x / 10.
df = pd.DataFrame({'ds': pd.date_range('01-01-2017', periods=len(x)),
                    'y': y_periodic})
df.head() # has to be a DataFrame with columns "ds" and "y"
df.set_index('ds').plot(style='-*')

# Estimate the model
m = Prophet()
m.fit(df);
forecast = m.predict(df)
m.plot_components(forecast);

Answer 5

不确定这方面的编程，但我会认真考虑移动平均线来解决这个问题。