Description Details Author(s) References See Also Examples
Used in Acevedo, M.F. 2012. "Simulation of Ecological and Environmental Models", CRC press.
Package: | seem |
Type: | Package |
Version: | 1.0 |
Date: | 2012-06-05 |
License: | GPL-2 |
Package to perform simulations, including sensitivity analysis. It provides a suite of simulator functions: sim.comp, sim.mruns, sim.rnum, sim, simd, simr. Model functions are employed to define the ODE to be simulated by the simulation functions. Nominal parameter values are defined in input files. Variation of param values are defined in lists.
Input files are in 'datafiles.zip' in directory 'datafiles' and organized by chapters of Acevedo (2012). Input files are required to run the examples below.
In addition, seem includes wrappers to use several simulation programs written in Fortran and C. These include cerio.F, semi.F, zelig.F, river.C, and forsucc.C. The zelig simulator used here is based on the version 2.3 by D.L. Urban.
M.F. Acevedo <Acevedo@unt.edu>
Acevedo, M.F. 2012. "Simulation of Ecological and Environmental Models", CRC press.
Model functions expon
, expon.z
, expon.g
, logistic
, two.stage.cont
, two.stage.cont.delay
, two.stage.x0
, and many others; Methods RK4
, euler
, RK4D
, ramos
, Simulator output functions onerun.out
, mruns.out
, rnum.out
, vars.out
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 | ## Not run:
model <- list(f=expon);file<-"chp4/exp-pop-inp.csv"
out.eu <- sim.comp(model,file,method="euler",lab.out="Euler")
model <- list(f=expon);file<-"chp4/exp-pop-inp.csv"
param <- list(plab="r", pval = seq(-0.02,0.02,0.01))
out.r <- sim.mruns(model,file, param, pdfout=TRUE)
model<-list(f=expon.rand); file<-"chp4/exp-rnum-inp.csv"
out.rk <- sim.rnum(model,file)
logis <- list(f=logistic)
t.X <- sim(logis,"chp6/logistic-inp.csv")
param <- list(plab="r", pval = seq(0.2,0.6,0.2))
t.X <- sim(logis,"chp6/logistic-inp.csv",param)
exp.sud <- list(f=expon,z=expon.z,g=expon.g)
t.X <- simd(exp.sud,file="chp7/exp-sud-inp.csv")
logis.sud <- list(f=logistic,z=logistic.z,g=logistic.g)
param <- list(plab="r",pval=c(0.1,0.2))
out.f <- simd(logis.sud,file="chp7/logis-sud-inp.csv",param)
two.stage.delay <-list(f=two.stage.cont.delay,x0=two.stage.x0)
param <- list(plab="nd", pval = c(2,5,6,7))
y.s <- simt(two.stage.delay,"chp10/two-stage-delay-inp.csv", param)
ga <- list(f=green.ampt.ramos);param <- list(plab="Ks",pval=c(10,20,30,40))
t.X <- simr(ga,"chp15/ga-inp.csv", param, pdfout=TRUE)
spp <- c("Role1", "Role2", "Role3", "Role4")
fileprefix <- "rolesm"; label <- "Semi-Markov Order 1"
fileout <- semi.F("chp16",fileprefix)
x <- read.plot.semi.out(fileout, spp, label)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.