ラグモデル 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