时间序列预测网站流量增长趋势(ARIMA)

在文章《浅析网站流量趋势预测》中，我们曾对网站的直接流量进行过分析和预测,直接流量整体较为平稳，但在工作日和周末存在较大差异。因此我们之前的方法是将工作日和周末的访问量数据分开处理，通过分布情况单独进行预估。这种方法得到的预测值范围较大，精确度也较差。本篇文章我们将使用时间序列的方法对网站直接流量的变化趋势进行预测。

下面是4月4日——5月1日的网站直接流量会话数据，我们在表格的周索引列中单独标明了工作日（weekday）和休息日（weekend）用以说明会话所属的类别。可以看出在这组数据中存在着明显的周末趋势，周末的会话数据要明显低于工作日的会话数据。

使用时间序列对数据进行分析和预测的第一步是观察数据的形态，因此我们将会话数据绘制到折线图中观察变化趋势。从直接流量变化趋势图上可以看出整体趋势较为平稳，没有明显的增长或下降。同时也能更明显的看出会话数据共分为4周，并且每一周都有明显的周末趋势。工作日每天的会话量在300—400之间，周末的会话量则为200以下，两者之间差距较大。如何对这两个差距较大的数据进行预测？我们使用一元线性回归的方法，一元线性回归可以反映一个因变量(Y)与一个自变量(X)之间的线性关系。对应到这里直接流量的会话就是因变量(Y)，我们只需要再找到一个自变量(X)就可以建立一元回归线性方程来预测直接流量的变化趋势了。

从前面的趋势图中可以看出，直接流量整体趋势比较平稳，只有在周末的时候才会出现较大的变化，造成周末数据与工作日数据的差异。因此造成直接流量变化的主要原因是所在日期的周属性，当所在日期为工作日(weekday)时，会话数据在300-400之间。而当所在日期为周末(weekend)时，会话数据下降到200以下。我们将日期的周属性设定为自变量(X)，取名“周变量”，这个周变量(X)共有两个取值，当日期为工作日时，取值为0，当日期为休息日时，取值为1。

设置好周变量后，我们在直接流量的源数据表中增加一列“周变量”，并按日期的周属性输入“周变量”的值。Weekday=0，weekend=1。现在我们有了因变量(Y)和自变量(X)，可以使用一元线性回归建立方程求两者之间的关系。

通过Excel中的回归分析功能，将会话作为因变量(Y)输入，将“周变量”(X)作为自变量输入，我们获得了以下分析结果。在这个结果中，我们主要关注“周变量”的系数和截距(Intercept)，通过这两个值可以建立“周变量”与会话的一元线性回归方程。

这是一元线性回归方程，我们将前面的分析结果代入到这个方程中，就可以计算出直接流量的值。其中y是直接流量的会话值，x是“周变量”的取值，a是“周变量”系数，b是截距。

前面说过，“周变量”有两个取值，工作日为0，休息日为1。因此在对直接流量进行预测时我们也分别按照日期的周属性取值来计算工作日和休息日的直接流量会话值，以下为两个取值的计算公式。

通过一元线性回归方程我们获得了直接流量在工作日和休息日的会话预测值，每天具体实际数据和预测数据展现在下面表格中。

将直接流量实际的会话数据与预测数据绘制成折线图进行对比，在较为平稳的趋势中预测数据较好的模拟并预测了直接流量中工作日和休息日的数据变化。

以上是一种较为简单且理想化的流量变化情况。实际工作中情况会更加复杂，直接流量除了自身的规律以外，还会受到外部因素的影响，呈现上升或下降的趋势。这时候，直接流量的整体趋势不再平稳，并且还同时掺杂自身的周末趋势。这种情况下如何对直接流量进行分析和预测呢？

下面是4月4日——5月1日的直接流量会话数据，由于有周末趋势存在我们同样增加了周索引列来标识会话所属日期的周属性。

第一步将直接流量的会话数据汇总折线图观察流量的形态。直接流量除了自身的周末趋势以外，还有一个明显的上涨趋势。我们在进行分析和预测时要兼顾直接流量的周末趋势和上涨趋势。这里一元线性回归已经无法满足我们的需求了，需要使用多元线性回归。多元线性回归是一元线性回归的升级版，可以反应一个因变量(Y)与多个自变量(X)之间的关系。这里直接流量的会话数据是我们要求的因变量(Y)，周末趋势是一个自变量(X)，另一个自变量(X)就是直接流量的增长趋势。

我们首先在直接流量的原始数据中增加两个字段作为两个自变量X。第一个字段是周变量，设置的方法和取值范围与前面例子中介绍的完全一致。第二个字段是趋势变量，用来描述直接流量的增长趋势。在前面的折线图中已经看到直接流量呈线性增长，因此我们使用从1到28作为趋势变量中每一天的取值。现在我们有了因变量(Y)和两个自变量(X),可以通过建立多元线性回归方程来计算并预测直接流量的会话值了。

使用Excel中的回归功能，将会话作为因变量(Y)输入，将趋势变量和周变量作为自变量(X)输入，获得以下分析结果。在这个结果中我们主要关注三个值，截距(Intercept)，趋势变量系数和周变量系数。

这是多元线性回归方程，我们将前面的分析结果代入到这个方程中，就可以获得直接流量的值。其中y是直接流量的会话值，x1是趋势变量，x2是周变量，b1是趋势变量系数，b2是周变量系数，a是截距。

通过多元线性回归方程我们获得了直接流量的预测值，每天的实际数据和预测数据展现在下面表格中。

将直接流量实际的会话数据与预测数据绘制成折线图进行对比，通过多元线性回归预测数据同时兼顾了直接流量中的周末趋势和整体的增长趋势。如果流量中还包括有其他影响因素只需通过绘制折线图观察影响因素的变化趋势和规律，然后增加相应的自变量(X)即可。

通过时间序列我们对直接流量的增长趋势进行了预测，很多时候除了知道趋势以外我们还希望对流量进行分解，了解影响流量变化的因素，这些因素在变化中所起的作用以及消除影响因素后流量真实的变化情况。我们使用时间序列对直接流量进行分解。下面是直接流量在最近4周时间的会话数据。

直接流量的变化可能受多个因素影响。例如前面的例子中，直接流量自身有明显的周末趋势，同时，受外部因素影响还有明显的增长趋势，此外可能还有一部分难以解释的随机波动因素。因此，我们将直接流量分为三个主要部分，分别为趋势值，周因素和随机波动。并对直接流量按这三部分进行分解。以下是时间序列乘法分解模型。这里需要说明的是，除了趋势值以外，周因素和随机波动将被分解为相对度量。数值大于1表示他们的影响在趋势值之上，数值小于1表示他们的影响在趋势值以下。

直接流量有明显的周末趋势，因此首先我们通过移动平均来消除直接流量中的周末趋势。由于周末趋势的变化是以周为单位的，因此我们选择7天的数据进行移动平均计算。下面是移动平均的计算公式。

我们按照移动平均公式对直接流量消除周末趋势，只保留直接流量的趋势值。以下是移动平均的计算结果。

移动平均通过对7天数据的平均对直接流量进行了平滑处理，消除了周末趋势及随机波动的影响，获得了直接流量的趋势值。下图是直接流量的实际值与趋势值的对比。

获得直接流量趋势值后，我们可以通过时间序列乘法分解模型获得周因素和随机波动的相对度量。将时间序列乘法分解模型两边同时除以趋势值，可以获得周因素和随机波动的组合度量。

分别将每一天的直接流量会话实际值除以移动平均后的趋势值，获得了周因素和随机波动的组合度量值。下表中汇总了整个时间段的周因素和随机波动的组合度量。

我们根据每一天的周因素和随机波动组合度量值计算出周指数，方法是将上表星期索引中同一天的值进行平均。下表是直接流量一周的指数值。其中大于1表示周因素的影响在趋势值之上，小于1表示周因素的影响在趋势值之下。从下表中可以明显看出，周一到周五周因素的影响都要大于趋势，而周末两天周因素的影响小于趋势。

将周指数还原到每一天中，获得直接流量每一天会话的星期指数。使用实际会话值除以星期指数获得直接流量消除星期因素的会话量数据。也就是直接流量排除周末因素影响后的变化趋势。

将消除周末趋势后的直接流量变化趋势绘制到折线图上，使用趋势线来拟合消除周因素影响的直接流量与时间周期。到目前为止我们分解出了直接流量中排除周因素和随机波动的会话值，排除周因素的会话值以及周因素对于趋势值的影响。下面我们将使用排除周因素的会话值对直接流量进行预测，并与周因素进行整合，预测出带有周因素的直接流量变化趋势。

我们将消除周因素后的直接流量会话值作为因变量(Y)，将时间作为自变量(X)，通过一元线性回归来描述会话值与时间的关系。将因变量(Y)与自变量(X)输入到回归分析中，获得以下结果。其中星期ID是时间X的系数，Intercept是截距。

将回归分析的结果导入到一元线性回归方程，计算出排除周因素的直接流量预测值。下表中显示了消除周末因素的直接预测值的计算结果。

将消除周末因素的预测值乘以星期指数，计算出直接流量实际值的预测趋势，下表中显示了计算结果和预测值。

将直接流量的实际值与增加周因素的预测值绘制到折线图中，由于我们是先剥离出了直接流量的趋势值，然后在趋势值和随机波动值上增加的周因素，因此预测值同时兼顾了直接流量的增长趋势和周末趋势。文章开始部分的预测相比，更加符合直接流量的真实变化趋势。

本篇文章的名字是使用时间序列预测网站流量增长趋势，但整篇文章只介绍了对直接流量的预测。其实这些方法适用于任何一种流量渠道，找到这些流量渠道的规律和影响因素就可以套用本文中的方法对流量趋势进行分解和预测。