做数据分析为何要学统计学（11）——如何进行时间序列分析

时间序列是由随时间变化的值构成，如产品销量、气温数据等等。通过对时间序列展开分析，能够回答如下问题：

• （1）被研究对象的活动特征是否有周期性（也称季节性）
• （2）被研究对象的活动特征是否有趋势性（上升或下降）

时间序列分析有多种方法，常用方法为STL(Seasonal and Trend decomposition using Loess)算法。该方法是一种把时间序列分解为趋势项(trend component)、季节项(seasonal component)和残差项(remainder/residual component/random)的过滤过程。如下图：

STL算法在1990年由密歇根大学的R. B. Cleveland教授以及AT&T Bell实验室的W. S. Cleveland等人研发。其特点是：

• 稳健的估计趋势项和季节项，而不会被数据中的异常行为扭曲
• 可以指定季节项的周期为采样时间间隔任意大于一的整数倍
• 可以分解有缺失值的时间序列

 以下用纽约1946年1月到1959年12月的每月新生儿数作为时间序列分析人口增长的规律。代码如下：

#读入时间序列数据
import pandas as pd
X=pd.read_csv("https://robjhyndman.com/tsdldata/data/nybirths.dat",header=None,names=["birth"])
#绘制折线图
X.plot()

 

#进行时间序列分析
import  statsmodels.api as sm
import matplotlib.pyplot as plt
res = sm.tsa.seasonal_decompose(X,period=12)
res.plot()
plt.xlabel("Month")
plt.show()

结果如下

可以看人口出生在第38个月左右后出生有明显上升趋势，而且从每年3月-6期间，出生人口开始进入高峰期，从8月份左右开始出生人口急速下降。