Rで刈込み移動平均 trimed mean

Rで刈込み移動平均を計算する．
zooパッケージの移動関数適応関数rollapplyを利用する．

#刈込なし
library(zoo)
v<-rep(1,1000)+rnorm(1000)
v[c(10,200,500,700)]<-10^5;
x<-(rollapply(zoo(v),50,mean));
v2<-rep(NA,length(v))
v2[index(x)]<-x
plot(v2)

#刈込あり
v<-rep(1,1000)+rnorm(1000)
v[c(100,200,500,700)]<-10^5;
x<-(rollapply(zoo(v),50,function(x){mean(x,0.1)}));
v2<-rep(NA,length(v))
v2[index(x)]<-x
plot(v2)

現在はrcpprollを使ったほうがよいらしいです。
2014-03-10　window functionまとめ
R　http://d.hatena.ne.jp/dichika/20140310/p1 …