ARIMA模型、随机游走模型RW模拟和预测时间序列趋势可视化

article/2025/9/11 18:28:26

原文链接: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

c9ea57e387da869c14fac81cce22dd7f.png

正如我们可以清楚地观察到的,这是一个非平稳序列,它的均值和标准偏差随时间变化不是恒定的。


点击标题查阅往期内容

4ae36fe60ed469ab87a88d136805e9b0.jpeg

R语言结合新冠疫情COVID-19股票价格预测:ARIMA,KNN和神经网络时间序列分析

outside_default.png

左右滑动查看更多

outside_default.png

01

0eb2ba440dd0061b5a5438e2067a865f.png

02

663c940e26482f7593fce92ec8a66663.png

03

90b6cf38909a87a9fae17f4220729380.png

04

fa5e7be4ed6d795b48711aa408ee751b.png

一阶差分序列

为了使序列平稳,我们取序列的一阶差分。

if <- diff

绘制时,您会注意到差分序列类似于白噪声。

2cc3582e8c4175484337123193dd2ed6.png

该 Rf 序列的统计数据计算如下:

> mean(Rf)
> sd(Wf)

cccd2a4e44fe1e363423315d5c594a3a.png

21c600bd468432c0ad0a473259f7bd2f.png

带偏移的随机游走序列

我们模拟的上述随机游走序列在均值附近上下徘徊。但是,我们可以让随机游走系列跟随上升或下降趋势,称为偏移。为此,我们为函数提供了一个额外的参数均值/截距  。这个截距是模型的斜率。我们还可以更改模拟序列的标准差。在下面的代码中,我们提供了 1 的平均值和 5 的标准差。

> Rt <- arima.sim
> plot.ts

fbbe85fac4c3d7c00999664e17fe2817.png

估计随机游走模型

为了拟合具有时间序列偏移的随机游走模型,我们将遵循以下步骤

  • 取数据的一阶差分。

  • arima使用阶数为 的函数 将白噪声模型拟合到差分数据 c(0,0,0)

  • 绘制原始时间序列图。

  • abline通过提供通过将白噪声模型拟合为斜率得到的截距,使用该函数添加估计趋势 。

1. 一阶差分

为了使这个数列平稳,我们将取数列的差值。

> plot.ts

feda72aac2891fcf4e2626a85aab4be4.png

2. 将白噪声模型拟合到差分数据

我们现在可以使用 arima将白噪声模型拟合到差分数据。

> whodl <- arima

84c25a45e925c4f44ddaeef072f76ed7.png

我们可以看到拟合的白噪声模型的截距为 0.67。

3. 绘制原始随机游走数据

这可以使用以下命令完成:

> plot.ts

80ed8570ea2659a92d45f47c55c65480.png

4.添加估计趋势

现在在同一个图上,我们要添加估计的趋势。在本课开始时,我们解释了随机游走序列如何是零均值白噪声序列的累积和(即积分)。因此,截距实际上是我们随机游走序列的斜率。

我们可以使用函数绘制趋势线 ,其中 a 是截距,b 是线的斜率。在我们的例子中,我们将指定白噪声模型的“a=0”和“b=intercept”。

> abline

估计的趋势线将添加到我们的图中。

20c73d81223d6ce114d01f21c162c0c5.png



0162078439057b2203ac652af6ea1d70.png

本文摘选《R语言模拟和预测ARIMA模型、随机游走模型RW时间序列趋势可视化》,点击“阅读原文”获取全文完整资料。


点击标题查阅往期内容

经济学:动态模型平均(DMA)、动态模型选择(DMS)预测原油时间序列价格

时变马尔可夫区制转换MRS自回归模型分析经济时间序列

R语言多维数据层次聚类散点图矩阵、配对图、平行坐标图、树状图可视化城市宏观经济指标数据

【数据分享】工人收入工资及其社会经济影响因素数据

R语言从经济时间序列中用HP滤波器,小波滤波和经验模态分解等提取周期性成分分析

R语言计量经济学:工具变量法(两阶段最小二乘法2SLS)线性模型分析人均食品消费时间序列数据和回归诊断

R语言计量经济学:虚拟变量(哑变量)在线性回归模型中的应用

PYTHON用时变马尔可夫区制转换(MRS)自回归模型分析经济时间序列

R语言使用ARIMAX预测失业率经济时间序列数据

【视频】Python和R语言使用指数加权平均(EWMA),ARIMA自回归移动平均模型预测时间序列

R语言时间序列GARCH模型分析股市波动率

R语言时变向量自回归(TV-VAR)模型分析时间序列和可视化

Python和R用EWMA,ARIMA模型预测时间序列

R语言用LASSO,adaptive LASSO预测通货膨胀时间序列

Python中的ARIMA模型、SARIMA模型和SARIMAX模型对时间序列预测

R语言arima,向量自回归(VAR),周期自回归(PAR)模型分析温度时间序列

Python用ARIMA和SARIMA模型预测销量时间序列数据

R语言线性回归和时间序列分析北京房价影响因素可视化案例

R语言使用ARIMA模型预测股票收益时间序列

R语言分布滞后线性和非线性模型(DLMs和DLNMs)分析时间序列数据

R语言广义相加模型 (GAMs)分析预测CO2时间序列数据

Python | ARIMA时间序列模型预测航空公司的乘客数量

R语言ARIMA,SARIMA预测道路交通流量时间序列分析:季节性、周期性

ARIMA模型预测CO2浓度时间序列-python实现

R语言用多元ARMA,GARCH ,EWMA, ETS,随机波动率SV模型对金融时间序列数据建模

欲获取全文文件,请点击左下角“阅读原文”。

f3d25a99e1bbe428907c3064f2951ea5.gif

64ae26007fa694de8ecbb48867ac6fb7.png

a5e31265abebd237233a41d93f19daa7.jpeg

934ce7b8d651daccf45ca43ba7b50cca.png

欲获取全文文件,请点击左下角“阅读原文”。

c460460969e06a233267e5636a9d8eea.gif


http://chatgpt.dhexx.cn/article/Z17UJ1i5.shtml

相关文章

随机游走(Random Walk)模型

Random Walk Model 1 模型及性质简介 给定一随机变量 u ( i ) { 1 , − 1 } u(i){\{1, -1\}} u(i){1,−1} 随机游走模型可表示为随时间 t t t变化的函数 y ( t ) ∑ i 1 t u ( i ) y(t)\sum_{i1}^{t} u(i) y(t)i1∑t​u(i) 几条随机游走可视化路线如下 性质一&#xff1a;…

读《PROSOSPEECH: ENHANCING PROSODY WITH QUANTIZED VECTOR PRE-TRAINING IN TEXT-TO-SPEECH》

当下韵律建模存在的问题&#xff1a; 1 提取的基音pitch信息存在误差&#xff0c;导致韵律合成出现问题 2 对韵律生成的相关要素 如基频 时长 能量等相互依存(dependent on each other) 共同产生了韵律相关的特征 3 韵律信息较高的可变性和高质量数据数目较少 导致不能完全学习…

UE4官方文档_Light Propagation Volumes_LPV方案

光线传播体积&#xff08;Light Propagation Volumes&#xff09;功能仍在开发中&#xff0c;不适用于生产。 本页面的内容&#xff1a; 启用光线传播体积基础场景设置光线传播体积设置 调整外观和性能 定向光源设置查看全局照明显示光线传播体积GI 替换材质切换其他注意事项 启…

Ue4 使用lpv快速增强间接光照效果

LPV缩写Light Propagation VolumesUe4自带&#xff0c;效果还可以&#xff0c;能快速在项目中实现不需要烘焙的间接光照效果主要原理使用光照生成点云进行对物体表面间接光进行计算测试版本4.16.3如何开启把r.LightPropagation1 加入到 consolevariables.ini 文件最后 &#…

实时GI方案概述

LPV CryTek原创的&#xff0c;但是貌似因为漏光的问题&#xff0c;没有广泛应用起来。 SVO VXGI Enlighten Enlighten的实时GI解决方案用的时预计算实时全局照明 (Precomputed Realtime GI)&#xff0c;这是一种允许交互式更新场景照明的技术&#xff0c;采用的是辐射度算…

IPVLAN

IPVLAN 一、拓扑图二、实验内容三、配置信息 一、拓扑图 二、实验内容 假设S1交换机由于某种原因无法配置&#xff0c;利用IP地址划分在S3做相应配置使得PC能供与服务器正常通信 三、配置信息 1、接口信息配置 S1的0/0/1和0/0/2接口无任何配置&#xff0c;0/0/3接口配置了a…

LPI

概述 LPI全称是Locality-specific Peripheral Interrupts&#xff08;LPIs&#xff09;&#xff0c;GICv3有两种方式支持LPIs&#xff1a; 1&#xff09;使用ITS把从设备发送的EventID转换成LPI INTID 2&#xff09;直接转发LPI INTID到Redistributors&#xff08;GICR_SETL…

系统辨识和自适应控制

系统辨识知识要点 1.为什么采用负反馈技术 2.什么是自适应控制&#xff0c;为什么采用自适应控制&#xff0c;指出自适应控制的使用场合 3.学习了什么辨识方法&#xff0c;这些方法之间的联系 4.最小二乘中的无偏性和一致性指的是什么 5.什么是白噪声 白噪声是一种具有…

【状态估计】用于描述符 LTI 和 LPV 系统的分析、状态估计和故障检测的算法(Matlab代码实现)

&#x1f4a5; &#x1f4a5; &#x1f49e; &#x1f49e; 欢迎来到本博客 ❤️ ❤️ &#x1f4a5; &#x1f4a5; &#x1f3c6; 博主优势&#xff1a; &#x1f31e; &#x1f31e; &#x1f31e;博客内容尽量做到思维缜密&#xff0c;逻辑清晰&#xff0c;为了方便读者。 …

Global Illumination_Light Propagation Volumes (LPV)

文章具体参照 https://ericpolman.com/ 本方法的思想就是把场景分成很多的小格子&#xff0c;然后计算每一个小格子里面的光照&#xff08;LPV&#xff09;。如果直接计算每个格子里面的光照那代价也是不可接受的&#xff0c;因此本算法用了一种很巧妙的方式来处理&#xff1…

LPV(Light Propagation Volumes)

lpv 测试了Light Propagation Volumes&#xff0c;全实时没有任何预处理的GI&#xff0c;而且可以适用任意场景。 文档很长&#xff0c;不过基本原理还是比较直白的&#xff1a; 生成reflect shadow map(rsm)。 将rsm信息用SH系数方式注入一个volumetexture中。 …

【GAMES-202实时渲染】4、3D空间全局光照(RSM、LPV、VXGI)

Lec7~8 1、Reflective Shadow Maps&#xff08;RSM&#xff09;2、Light Propagation Volumes&#xff08;LPV&#xff09;3、Voxel Global Illumination&#xff08;VXGI&#xff09; 1、Reflective Shadow Maps&#xff08;RSM&#xff09; RSM是一个特别经典的计算全局光照…

lpv

测试了Light Propagation Volumes&#xff0c;全实时没有任何预处理的GI&#xff0c;而且可以适用任意场景。 文档很长&#xff0c;不过基本原理还是比较直白的&#xff1a; 生成reflect shadow map(rsm)。 将rsm信息用SH系数方式注入一个volumetexture中。 在vol…

操作系统经典 pv过桥问题

Semophere bridge1; Semophere mutexNS1,mutexSN1;//用于保护countNS,countSN int countNS0,countSN0; Semophere s11,s20;//用于交替通过 StoN(){while(1){P(mutexSN);countSN;//来车了v(mutexSN);p(mutexSN);if(countNS0){//对面无车,则直接通过P(bridge);通过countSN--;V…

C语言解决四人/多人过桥问题

参加笔试的时候遇到一道经典的算法题&#xff0c;四人过桥问题。当时没写出来&#x1f605;。 四人过桥问题&#xff1a;在一个黑夜里&#xff0c;有四个人需要过桥&#xff0c;每次只能通过两人&#xff0c;其中一人必须拿着手电筒&#xff1b;但只有一个手电筒&#xff0c;所…

小明过桥问题

小明家必须要过一座桥。小明过桥最快要&#xff11;秒&#xff0c;小明的弟弟最快要&#xff13;秒&#xff0c;小明的爸爸最快要&#xff16;秒&#xff0c;小明的妈妈最快要&#xff18;秒&#xff0c;小明的爷爷最快要&#xff11;&#xff12;秒。每次此桥最多可过两人&…

过桥问题

在一个夜黑风高的晚上&#xff0c;有n&#xff08;n < 50&#xff09;个小朋友在桥的这边&#xff0c;现在他们需要过桥&#xff0c;但是由于桥很窄&#xff0c;每次只允许不大于两人通过&#xff0c;他们只有一个手电筒&#xff0c;所以每次过桥的两个人需要把手电筒带回来…

过桥问题的通解

问题一,一个典型过桥问题: 小明一家5口人在夜晚过一座桥,小明过桥要1分钟,小明的弟弟过桥要3分钟,小明的爸爸过桥要6分钟,小明的妈妈过桥要8分钟,小明的爷爷过桥要12分钟;这座桥每次只能过2个人,因是夜晚,过桥时必须提着灯,小明有一只灯,点燃后30分钟会熄灭,问怎…

如何打开tdms文件?

https://www.zhihu.com/question/305029962/answer/1203851780v 下载地址 http://www.ni.com/example/27944/en/ 还有一个综述&#xff0c;写的挺好&#xff01; https://wenku.baidu.com/view/c62700e4aa00b52acec7ca09.html

转载:TDM协议

转自http://www.wangdali.net/i2s/ 1. PCM简介 PCM (Pulse Code Modulation) 是通过等时间隔(即采样率时钟周期)采样将模拟信号数字化的方法。图11为4 bit 采样深度的PCM数据量化示意图。 图11. 4-bit PCM的采样量化 PCM数字音频接口,即说明接口上传输的音频数据通过PCM…