# logistic
logis <- list(f=logistic)
# single run
t.X <- sim(logis,file="chp6/logistic-inp.csv")
# multiple runs
param <- list(plab="r", pval = seq(0.2,0.6,0.2))
t.X <- sim(logis,"chp6/logistic-inp.csv",param)
# Harvest, multiple runs
logis.f <- list(f=logistic.forced)
param <- list(plab="Ha",pval=c(0,-0.1,-0.2,-0.3))
t.X <- sim(logis.f,"chp7/logis-harvest-inp.csv",param)
# Seasonality
logis.v <- list(f=logistic.var)
param <- list(plab="Kd",pval=c(0,-0.5,-0.6))
t.X <- sim(logis.v,"chp7/logis-var-inp.csv",param)
# sudden disturbance
logis.sud <- list(f=logistic,z=logistic.z,g=logistic.g)
t.X <- simd(logis.sud,file="chp7/logis-sud-inp.csv")
# M3 decay multiple runs
m3decay<- list(f=M3)
param <- list(plab="Kmax", pval = seq(-20,-60,-10))
t.X <- sim(m3decay,"chp6/m3decay-inp.csv", param)
# monod batch multiple runs
monod <- list(f=monod.batch)
param <- list(plab="Kmax", pval = seq(0.5,1.5,0.5))
t.X <- sim(monod,"chp6/monod-batch-inp.csv", param)
Run the code above in your browser using DataLab