3

このコードをベクトル表記で再現したい

getEMA2<-function(x,win){
k<-2/(win+1)
v<-vector()
for (i in 1:length(x)){
    if (i==1){
        v[i]<-x[i]
    }
    else{
        v[i]<-k*x[i]+(1-k)*v[i-1]
    }
}
return (v)
}
testOutput<-getEMA2(rnorm(100,0,1),5)

フィルター機能を使用してみましたが、再帰/畳み込みメソッドでこれを達成できるようには見えません

返信ありがとうございます。

4

1 に答える 1

7

filter計算するので

y[n] = x[n] + alpha * y[n-1]

結果を再スケーリングする必要があります。

f <- function(x,win) {
  alpha <- 2/(win+1)
  filter(x, 1-alpha, method="recursive", side=1, init=x[1]/alpha)*alpha
}
x <- 1:10
k <- 3
getEMA2(x,k)
f(x,k) # identical

これらのフィルターのほとんどは、TTRパッケージですでに定義されています。

于 2012-07-12T06:22:47.210 に答える