[Rcode]时间序列(整理自R in action)【待更新】
来源:互联网 发布:mysql execute using 编辑:程序博客网 时间:2024/06/06 14:15
一 生成时序对象
sales<-c(18,33,41,7,34)tsales<-ts(sales,start=c(2003,1),frequency=12)###常规操作得到基础信息plot(tsales) #绘制时序图start(tsales);end(tsales) #起止点frequency(tsales) ###取子集tsales.subset<-window(tsales,start=c(2003,5),end=(2003,7));tsales.subset
二 时间序列的平滑化与季节性分解
#转化为对数数据可降低波动性同时可以使用相加模型方便做季节性分解 plot(AirPassengers) lAirPassengers<-log(AirPassengers) plot(lAirPassengers)#分解时间序列 fit<-stl(lAirPassengers,s.window = "period") plot(fit) fit$time.series #查看每个观测值各分解项的值 #注意这里s.window="period",从而每年的相同月份季节效应相同#把对数数据转化为原始尺度 exp(fit$time.series) #数据解读:7月的季节性效应为1.24,即乘客数增加了24%#季节的可视化 install.packages("forecast") library(forecast) opar<-par(no.readonly = T) par(mfrow=c(2,1)) monthplot(AirPassengers,xlab="",ylab="") seasonplot(AirPassengers,years.labels="TRUE",main="seasonplot")
三 指数预测模型
#指数模型预测(单,双,三指数模型) #单指数平滑 library(forecast) fit<-ets(nhtemp,model="ANN") fit #一步向前预测并预测 forecast(fit,1) plot(forecast(fit,1),xlab="Year",ylab=expression(paste("Temperature(",degree*F,")",)),main="New Haven Annual Mean Temperature") accuracy(fit) #得到准确性度量 #Holt指数平滑和Holt-Winters指数平滑 library(forecast) fit<-ets(log(AirPassengers),model="AAA") fit pred<-forecast(fit,5) pred #ets()函数与自动预测 library(forecast) fit<-ets(JohnsonJohnson) #我们没有指定模型时,软件自动搜索了一系列模型并且选择了最小化拟合标准的模型。 fit
四 ARIMA预测模型 (保证数据的平稳性)
#回归前的操作: #平稳性: #diff(ts,differences=d),forecast的ndiffs()函数可以帮我找到最优的d值 #保证方差为常数: #作对数变换或者Box-Cox变换 #example:Nile数据: library(forecast) library(tseries) plot(Nile) ndiffs(Nile) dNile<-diff(Nile) plot(dNile) adf.test(dNile) #tseries包的adf.test做ADF检验判断序列平稳性 #选择模型: Acf(dNile) Pacf(dNile) #拟合模型: library(forecast) fit<-arima(Nile,order=c(0,1,1)) fit accuracy(fit) #模型评价:(主要是残差的独立正态性) qqnorm(fit$residuals) qqline(fit$residuals) Box.test(fit$residuals,type="Ljung-Box") #预测: forecast(fit,3) plot(forecast(fit,3),xlab="Year",ylab="Annual Flow") #ARIMA的自动预测 library(forecast) fit<-auto.arima(sunspots) fit forecast(fit,3) accuracy(fit)
阅读全文
0 0
- [Rcode]时间序列(整理自R in action)【待更新】
- 时间序列分析笔记(待整理)
- R语言时间序列函数整理[不断更新]
- R语言时间序列函数整理
- R语言时间序列函数整理
- R语言时间序列函数整理
- R语言-时间序列函数整理
- R-时间序列自相关acf,偏自相关pacf
- jvisualvm监控使用(待更新整理)
- Mongodb命令整理(待更新)
- R in action 笔记(第一部分)
- R in action 笔记(第二部分)
- R in action
- R语言常见问题(持续更新整理)
- R语言做时间序列(未完)
- 时间序列完全教程(R)
- 时间序列完全教程(R)
- R----简述时间序列(ARIMA)
- 聊一聊 Spring 中的线程安全性
- android 给图片设置滤镜报错OutOfMemoryError
- YARN内存使用优化配置
- Hadoop的三种作业调度原理汇总
- FileUtils方法大全
- [Rcode]时间序列(整理自R in action)【待更新】
- Qt 隐藏鼠标指针
- Java解析excel操作
- PAT 1029. 旧键盘(20)
- 日历打印
- Happy Vertices
- [算法期中] 分组
- 8.11 应用举例
- 杭电ACM OJ 1031 Design T-Shirt 还是快速默写一波 快速排序 轻松搞定!