ラグモデル fを任意の関数に。 lag_model<-function(tmax=100,lag=1,f=mean){ x<-0 tmp<-unif(lag+1) for(i in 1:tmax){ x[i+1]<-eval(f)(tmp) tmp<-c(tmp[2:length(tmp)],x[i+1]) } }

Rで関数を変数のように扱って実行 a<-sum; eval(a)(1:100); a<-sd; eval(a)(1:100) 参考 http://d.hatena.ne.jp/tsutatsutatsuta/20120114/1326542583

Rでbox平均。 corcor2<-0 sdsdsd2<-0 IQRIQR3<-0 lenlen3<-0 sdsdsd3<-0 meanmean3<-0 jjj<-flow_volume for(i in 1:length(filelist)){ try(ll<-mat_volume[,i]) try(corcor2[i]<-cor(ll,jjj,use="complete.obs")) try(sdsdsd2[i]<-sd(ll/jjj,na.rm=T)) ll2…

distribution_of_runs a123bvector(kkkb,c(1,-1)) b123bvector(kkkb,c(-1,1)) z diff(z) }