R语言模拟和预测ARIMA模型、随机游走模型RW时间序列趋势可视化
原文链接:http://tecdat.cn/?p=25122
当一个序列遵循随机游走模型时,就说它是非平稳的。我们可以通过对时间序列进行一阶差分来对其进行平稳化,这将产生一个平稳序列,即零均值白噪声序列。例如,股票的股价遵循随机游走模型,收益序列(价格序列的差分)将遵循白噪声模型。
让我们更详细地了解这种现象。
由于随机游走序列的差分是白噪声序列,我们可以说随机游走序列是零均值白噪声序列的累积和(即积分)。有了这些信息,我们可以以 ARIMA 模型的形式定义 Random Walk 系列,如下所示:
ARIMA(0,1,0)
其中
- 自回归部分,p = 0
- 积分部分,d = 1
- 移动平均部分,q = 0
模拟随机游走序列 我们现在可以通过为
arima.sim()
函数提供适当的参数来模拟 R 中的随机游走序列, 如下所示:R <- arima.sim
我们可以使用该
plot.ts()
函数绘制新生成的序列 。> plot.ts
文章图片
正如我们可以清楚地观察到的,这是一个非平稳序列,它的均值和标准偏差随时间变化不是恒定的。
一阶差分序列 为了使序列平稳,我们取序列的一阶差分。
if <- diff
绘制时,您会注意到差分序列类似于白噪声。
文章图片
该
RW_diff
序列的统计数据计算如下:> mean(Rf)
> sd(Wf)
文章图片
文章图片
带偏移的随机游走序列 我们模拟的上述随机游走序列在均值附近上下徘徊。但是,我们可以让随机游走系列跟随上升或下降趋势,称为偏移。为此,我们为函数提供了一个额外的参数均值/截距
arima.sim()
。这个截距是模型的斜率。我们还可以更改模拟序列的标准差。在下面的代码中,我们提供了 1 的平均值和 5 的标准差。> Rt <- arima.sim
> plot.ts
文章图片
估计随机游走模型 为了拟合具有时间序列偏移的随机游走模型,我们将遵循以下步骤
- 取数据的一阶差分。
arima()
使用阶数为 的函数 将白噪声模型拟合到差分数据c(0,0,0)
。- 绘制原始时间序列图。
abline()
通过提供通过将白噪声模型拟合为斜率得到的截距,使用该函数添加估计趋势 。
> plot.ts
文章图片
2. 将白噪声模型拟合到差分数据 我们现在可以使用该
arima()
函数将白噪声模型拟合到差分数据。> whodl <- arima
文章图片
我们可以看到拟合的白噪声模型的截距为 0.67。
3. 绘制原始随机游走数据 这可以使用以下命令完成:
> plot.ts
【R语言模拟和预测ARIMA模型、随机游走模型RW时间序列趋势可视化】
文章图片
4.添加估计趋势 现在在同一个图上,我们要添加估计的趋势。在本课开始时,我们解释了随机游走序列如何是零均值白噪声序列的累积和(即积分)。因此,截距实际上是我们随机游走序列的斜率。
我们可以使用函数绘制趋势线
abline(a,b)
,其中 a 是截距,b 是线的斜率。在我们的例子中,我们将指定白噪声模型的“a=0”和“b=intercept”。> abline
估计的趋势线将添加到我们的图中。
文章图片
文章图片
最受欢迎的见解
1.在python中使用lstm和pytorch进行时间序列预测
2.python中利用长短期记忆模型lstm进行时间序列预测分析
3.使用r语言进行时间序列(arima,指数平滑)分析
4.r语言多元copula-garch-模型时间序列预测
5.r语言copulas和金融时间序列案例
6.使用r语言随机波动模型sv处理时间序列中的随机波动
7.r语言时间序列tar阈值自回归模型
8.r语言k-shape时间序列聚类方法对股票价格时间序列聚类
9.python3用arima模型进行时间序列预测
推荐阅读
- 急于表达——往往欲速则不达
- 第三节|第三节 快乐和幸福(12)
- 20170612时间和注意力开销记录
- 2.6|2.6 Photoshop操作步骤的撤消和重做 [Ps教程]
- 对称加密和非对称加密的区别
- 眼光要放高远
- 樱花雨
- 前任
- 2020-04-07vue中Axios的封装和API接口的管理
- 烦恼和幸福