http://stackoverflow.com/questions/8960517/how-to-calling-auto-arima-function-of-r-in-java-and-store-the-result-of-forecast

打开R,安装一个叫做"Forecast"的包。每次开启R后,使用之前用
library('forecast')
加载该包。

这里我用传说中的Airline Model数据。载入数据,转换成TS格式
airdata<-read.table('airline.dat')
airts<-ts(airdata,start=1949,frequency=12)

然后用forecast包中的auto.arima自动拟合Arima模型。
arima1<-auto.arima(airts,trace=T)
显示的结果如下:
ARIMA(2,0,2)(1,1,1)[12] with drift : 974.1468
ARIMA(0,0,0)(0,1,0)[12] with drift : 1077.823
ARIMA(1,0,0)(1,1,0)[12] with drift : 974.92
ARIMA(0,0,1)(0,1,1)[12] with drift : 1022.198
ARIMA(2,0,2)(0,1,1)[12] with drift : 967.1033
ARIMA(2,0,2)(0,1,0)[12] with drift : 966.755
ARIMA(1,0,2)(0,1,0)[12] with drift : 964.3004
ARIMA(1,0,1)(0,1,0)[12] with drift : 963.9208
ARIMA(1,0,1)(0,1,0)[12] : 971.225
ARIMA(1,0,1)(1,1,0)[12] with drift : 972.4003
ARIMA(1,0,1)(0,1,1)[12] with drift : 963.9781
ARIMA(1,0,1)(1,1,1)[12] with drift : 971.7862
ARIMA(0,0,1)(0,1,0)[12] with drift : 1022.291
ARIMA(2,0,1)(0,1,0)[12] with drift : 965.183
ARIMA(1,0,0)(0,1,0)[12] with drift : 966.9728

Best model: ARIMA(1,0,1)(0,1,0)[12] with drift

结果是一个AR(1),MA(1)和季节差分一次的Arima模型。Arima模型自动拟合的关键就是定阶,以前用的办法是EACF(extended (sample) autocorrelation function)来定阶,不过现在一般用AIC,AICc,BIC等统计量来定阶。例如上面的974.1468 等就是该模型的AIC

然后可以预测了。
airfore<-forecast(arima1,h=30,fan=T)
预测了30个月,另外计算了置信区间50到99的预测值。
plot(airfore)
最后绘图。如果需要获取其中的预测数据,可以用: airefore$mean 获取。

http://blog.macro2.org/?p=609

R forcast auto arima用法相关推荐

  1. 独家 | 利用Auto ARIMA构建高性能时间序列模型(附Python和R代码)

    作者:AISHWARYA SINGH 翻译:陈之炎 校对:丁楠雅 本文共3400字,建议阅读10+分钟. 本文介绍了ARIMA的概念,并带你用Python和R训练一个数据集实现它. 简介 想象你现在有 ...

  2. R语言拟合ARIMA模型:使用forecast包中的auto.arima函数自动搜索最佳参数组合、模型阶数(p,d,q)、设置seasonal参数指定在模型中是否包含季节信息

    R语言拟合ARIMA模型:使用forecast包中的auto.arima函数自动搜索最佳参数组合.模型阶数(p,d,q).设置seasonal参数指定在模型中是否包含季节信息 目录

  3. R语言 时间序列arima模型

    基本理论知识    ARMA模型称为自回归移动平均模型,是时间序列里常用的模型之一.ARMA模型是对不含季节变动的平稳序列进行建模.它将序列值表示为过去值和过去扰动项的加权和.模型形式如下: yt=c ...

  4. R语言时间序列ARIMA新手教程

    R语言时间序列ARIMA新手教程 首先说一下ARMA回归的底层逻辑,所谓的AR模型和MA模型都是ARMA模型的一种特殊情况,有点类似正方形和长方形都是矩形.ARMA模型的表达式为: p为自回归部分的滞 ...

  5. Python数据分析案例-分别使用时间序列ARIMA、SARIMAX模型与Auto ARIMA预测国内汽车月销量

    1. 前言 模型: ARIMA模型(英语:Autoregressive Integrated Moving Average model),差分整合移动平均自回归模型,又称整合移动平均自回归模型(移动也 ...

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

    "预测非常困难,特别是关于未来".丹麦物理学家尼尔斯·波尔(Neils Bohr),最近我们被要求撰写关于arima的研究报告,包括一些图形和统计输出. 很多人都会看到这句名言.预 ...

  7. R语言 - seasonal ARIMA与带傅里叶修正项的ARIMA预测及比较

    最近接触到了带傅里叶修正项的ARIMA模型(ARIMA with fourier modification, 以下简称FARIMA),学校里面虽然上过一些时间序列,知道time series deco ...

  8. 总结:min-height:100px; height:auto;的用法(新浪博客 )

    总结:min-height:100px; height:auto;的用法 1.IE6有个BUG,就是如果指定了一个盒子的高度,但里面内容又溢出了,那这个盒子的高度则会失去作用,盒子的高度会随着内容的增 ...

  9. T1275找出井字棋的获胜者+auto的用法

    题目:找出井字棋的获胜者 A 和 B 在一个 3 x 3 的网格上玩井字棋. 井字棋游戏的规则如下: 玩家轮流将棋子放在空方格 (" ") 上. 第一个玩家 A 总是用 " ...

最新文章

  1. NSHelper.showAlertTitle的两种用法 swift
  2. 再学 GDI+[11]: DrawCurve - 绘制曲线
  3. 2020-08-24绘制ROC   PR曲线 核心方法总结 ,计算AUC核心方法
  4. 也做SQL查询:班级总成绩 前三名,总成绩有相等的情况
  5. 22.PATH环境变量
  6. ORA-28000: the account is locked
  7. 以下表示中 不能用作c语言常量的是0UL,C语言笔试卷.doc
  8. 基于CSE的微服务工程实践-Native API先行
  9. 丰巢快递将强制收费引网友吐槽,全国第一驿站表示继续免费
  10. nyoj 471:好多的树(容斥原理)
  11. 硬盘无法识别怎么办 硬盘无法识别怎么恢复数据
  12. java定义vip顾客继承顾客_Java初级教频教程 - JavaSE - Java - 私塾在线 - 只做精品视频课程服务...
  13. [论文总结] 深度学习在农业领域应用论文笔记10
  14. k易语言html导入超级列表框,易语言超级列表框导入TXT内容的方法
  15. 视频教程-一小时搞定简单VBA编程 Excel宏编程快速入门-Office/WPS
  16. WEB前端学习day-6-盒子,浮动,学成在线案例
  17. PHP 判断一个数是否是质数
  18. C语言执行时进行窗口隐藏
  19. MOOS-ivp简介
  20. 努力前端【LeetCode-10】448. 找到所有数组中消失的数字 442. 数组中重复的数据(中等) 41. 缺失的第一个正数(困难) [鸽笼原理,数组,Map,类似No.645]

热门文章

  1. 安卓应用 - 公开市场上传投放
  2. 微信早安,利用uniCloud阿里云的云函数实现定时推送
  3. 部署外网网站(二)——宝塔搭建服务器
  4. 2006-10-30 18:37:00 著名Linux内核程序员大鹰 ox啊
  5. 联想一体机电源键不亮_联想力压华为拿下6·18 PC“六冠王”,背后有何秘诀?...
  6. 华夏ERP没有找到新增功能
  7. OpenGL画蜗型线、心形线、三叶曲线、四叶曲线、螺旋线
  8. sign-file: certs/signing_key.pem: 解决方法
  9. Oops是什么有什么用
  10. Hibernate 列映射 – 解决其他类型映射命名冲突