首页 > 运维 > linux运维 > 正文

ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

蓮花仙者
发布: 2025-07-13 10:32:20
原创
793人浏览过

原文链接:http://tecdat.cn/?p=25180 时间序列分析

对于时间序列分析,有两种数据格式: ts (时间序列)和 xts (可扩展时间序列)。前者不需要时间戳,可以直接从向量转换。后者非常重视日期和时间,因此只能使用日期和/或时间列来定义。我们涵盖了基本的时间序列模型,即 arima、garch 和 var。

时间序列数据

函数 ts 将任何向量转换为时间序列数据。

代码语言:javascript代码运行次数:0运行复制
price
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

我们首先为估计定义一个时间序列(ts)对象。请注意, ts 与 xts类似, 但没有日期和时间。

代码语言:javascript代码运行次数:0运行复制
df <- ts(df)df
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
可扩展的时间序列数据xts

要处理高频数据(分秒),我们需要包 xts。该包定义可扩展时间序列 ( xts ) 对象。

以下代码安装并加载 xts 包。

代码语言:javascript代码运行次数:0运行复制
library(xts)
登录后复制

考虑我们的可扩展时间序列的以下数据

代码语言:javascript代码运行次数:0运行复制
date  time  price
登录后复制

现在我们准备定义 xts 对象。代码 as.POSIXct() 将字符串转换为带有分钟和秒的日期格式。

代码语言:javascript代码运行次数:0运行复制
df <-data.framedf$daime <-pastedf$dttime <-as.POSIXctdf <- xts
登录后复制

对于仅使用日期的转换,我们使用 POSIXlt() 而不是 POSIXct()。

代码语言:javascript代码运行次数:0运行复制
df$date <- as.POSIXctdf$price <-as.numericprice <-xts
登录后复制
自回归移动平均模型arima

可以使用 arima() 函数估计自回归移动平均模型。

以下代码估计了一个 AR(1) 模型:

代码语言:javascript代码运行次数:0运行复制
 AR1
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

以下代码估计了一个 AR(2) 模型:

代码语言:javascript代码运行次数:0运行复制
AR2 <- arimaAR2
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

以下代码估计一个 MA(1) 模型:

代码语言:javascript代码运行次数:0运行复制
MA1 <- arimaMA1
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

以下代码估计一个 MA(2) 模型:

代码语言:javascript代码运行次数:0运行复制
MA2 <- arima
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

以下代码估计了一个 ARMA(1,1) 模型:

代码语言:javascript代码运行次数:0运行复制
ARMA11 <- arima
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

有时,我们只想保留系数。

代码语言:javascript代码运行次数:0运行复制
coef #得到系数
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

以下代码显示了残差图。

代码语言:javascript代码运行次数:0运行复制
plot
登录后复制
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

R 有一个方便的函数来 autofit() 拟合ARIMA 模型的参数。

现在寻找最好的 ARIMA 模型了。

代码语言:javascript代码运行次数:0运行复制
autoarma
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

时间序列模型的一项重要功能是预测。以下代码给出了两步的预测:

代码语言:javascript代码运行次数:0运行复制
teFoast <-predict
登录后复制

下面显示了预测图。

代码语言:javascript代码运行次数:0运行复制
plot.ts#可视化预测
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARCH 和 GARCH模型

要估计 ARCH 和 GARCH 模型,我们需要安装garch。

我们将在生成随机数时使用 ARMA(1,1) 估计 GARCH(1,1)

代码语言:javascript代码运行次数:0运行复制
a <- runif #随机数Spec <-ugarchspec
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

为了获得 GARCH 模型的具体结果,我们使用以下代码:

代码语言:javascript代码运行次数:0运行复制
coffnt <-coefvoy <- sigma
登录后复制
VAR模型

以下数据将用于估计 VAR 模型。

以下代码估计 VAR(2) 模型。

代码语言:javascript代码运行次数:0运行复制
abr<-VAR #运行 VAR(2)coef       #VAR的系数公式
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
代码语言:javascript代码运行次数:0运行复制
summary   #VAR的摘要
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

生成系数图

以下代码为 VAR 模型生成系数图:

代码语言:javascript代码运行次数:0运行复制
plot
登录后复制
登录后复制
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列

以上就是ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列的详细内容,更多请关注php中文网其它相关文章!

最佳 Windows 性能的顶级免费优化软件
最佳 Windows 性能的顶级免费优化软件

每个人都需要一台速度更快、更稳定的 PC。随着时间的推移,垃圾文件、旧注册表数据和不必要的后台进程会占用资源并降低性能。幸运的是,许多工具可以让 Windows 保持平稳运行。

下载
相关标签:
来源:php中文网
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn
最新问题
开源免费商场系统广告
热门教程
更多>
最新下载
更多>
网站特效
网站源码
网站素材
前端模板
关于我们 免责申明 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送
PHP中文网APP
随时随地碎片化学习
PHP中文网抖音号
发现有趣的

Copyright 2014-2025 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号