时序平滑化处理

  • 经过简单移动平均进行平滑处理

时序数据中一般有很明显的随机或偏差成分,为了辨明数据中的规律,画出一条平滑曲线,该曲线消除或者减弱时序中不规则和其余变更,呈现现象的 长期趋势,最简单的方法就是 简单移动平均函数

其中 St 是时间点t的平滑值, k =2q +1 是每次用来平均观测值的个数,通常咱们会将其设为一个奇数,居中移动平均法的代价是,每一个时序集中咱们会损失最后的 (k-1)/2个观测值spa

  R中有几个函数均可以作简单的移动平均,包括 TTR包 中的 SMA()函数,zoo包中的 rollmean()函数,forecast() 包中的 ma() 函数code

#以下图1
> library(forecast)
> opar <- par(no.readonly=TRUE)
> par(mfrow=c(2,2))
> ylim <- c(min(Nile), max(Nile))
> plot(Nile, main="Raw time series")
> plot(ma(Nile, 3), main="Simple Moving Averages (k=3)", ylim=ylim)
> plot(ma(Nile, 7), main="Simple Moving Averages (k=7)", ylim=ylim)
> plot(ma(Nile, 15), main="Simple Moving Averages (k=15)", ylim=ylim)
> par(opar)

                      图1it

相关文章
相关标签/搜索