R语言时间序列代码整理
时间序列:
data<-read.csv('/Users/ji_fanyang/Desktop/test2.csv')
myseries<- ts(data[1],start=c(2017,1),end=c(2020,7),frequency=12)
install.packages('forecast')
library(forecast)
1.时序的平滑化
平滑曲线,居中移动平均
opar<-par(no.readonly=TRUE)
par(mfrow=c(2,2))
ylim<-c(min(myseries),max(myseries))
plot(myseries,main="Raw time series")
plot(ma(myseries,3),main="simple moving average(k=3)",ylim=ylim)
plot(ma(myseries,7),main="simple moving average(k=7)",ylim=ylim)
plot(ma(myseries,15),main="simple moving average(k=15)",ylim=ylim)
2.季节性分解(趋势因子,季节因子,随机因子)
stl=(ts,s.window=,t.window=),为相加模型
S.window:控制季节效应变化的速度,s.window=“periodic”,使得季节效应在每个月都是一致的
t.window:控制趋势变化的速度
myseries2<-log(myse