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