2015-10-07から1日間の記事一覧

Rで刈込み移動平均 trimed mean

R

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)…