# exponential
model <- list(f=expon)
# single run
file<-"chp4/exp-pop-inp.csv"
out <- sim.comp(model,file)
# multiple runs
param <- list(plab="r", pval = seq(-0.02,0.02,0.01))
out.r <- sim.mruns(model, file, param)
# Monte Carlo run
model <- list(f=expon.rand)
file<-"chp4/exp-rnum-inp.csv"
out.rk <- sim.rnum(model,file)
# Harvest, multiple runs
exp.f <- list(f=expon.forced)
param <- list(plab="Ha",pval=c(0,-0.1,-0.2))
t.X <- sim(exp.f,"chp7/exp-harvest-inp.csv",param)
# Seasonality
exp.v <- list(f=expon.var)
param <- list(plab="rd",pval=c(0,-0.005,-0.01))
t.X <- sim(exp.v,"chp7/exp-var-inp.csv",param)
# sudden disturbance
exp.sud <- list(f=expon,z=expon.z,g=expon.g)
param <- list(plab="Ha",pval=c(-10,-20,-30))
t.X <- simd(exp.sud,file="chp7/exp-sud-inp.csv",param)
Run the code above in your browser using DataLab