右メニュー


※上記の広告は60日以上更新のないWIKIに表示されています。更新することで広告が下部へ移動します。

更新履歴

取得中です。



#library
par(mfrow=c(2,1))

set.seed(2)

ewma <- function (x, lambda, init)
{
 y <- filter(lambda*x, filter=1-lambda, method="recursive", init=init)
 return (y)
}

x = ts(diffinv(rnorm(100)))

ts.plot(x, ylim=c(-20,20))


for(i in 1:10){

a = (i-1)*0.1+0.1
y1 = ewma(x,lambda=a, init=x[1])
ts.plot(y1, ylim=c(-20,20), col=2)
par(new=T)
print(a)
}