Description Usage Arguments Details Value Note Author(s) References See Also Examples
Functions parab, kt.rate, Q10 are for single temperature response. Function Q10.mruns explores sensitivity to parameters of Q10.
1 2 3 4 |
Temp |
temperature |
a |
parameter |
b |
parameter |
k0 |
parameter |
q10 |
parameter |
param |
parameter array c(Topt,Tmax,q10) where Topt is temp for which optimum occurs, Tmax is maximum temperature, and q10 is coefficient |
param.nom |
nominal values of parameters c(Topt,Tmax,q10) |
Topt.sens |
values to vary parameter Topt |
Tmax.sens |
values to vary parameter Tmax |
q10.sens |
values to vary parameter q10 |
The three functions parab, kT.rate, and Q10, calculate the response to temperature using the parabolic, kT by q10, and the gamma-like q10.
Function Q10.mruns explores sensitivity to parameters of Q10.
Temperature response
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.
Miguel F. Acevedo Acevedo@unt.edu
Acevedo M.F. 2012. Simulation of Ecological and Environmental Models. CRC Press.
Simulation functions sim
, simd
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 38 | ## Not run:
# parabolic temp response
Temp <- seq(0,50); nT <- length(Temp)
Tu <- c(40,45,50);np <- length(Tu)
Tl <- 0
# dimension array
y <- matrix(nrow=nT,ncol=np)
# loop for the various optimal temperatures and then for temp range
for(j in 1:np){
for(i in 1:nT) y[i,j] <- parab(Temp[i],Tl,Tu[j])
}
# rate kT
Temp <- seq(20,30); nT <- length(Temp)
k0 <- 1
q10 <- c(1,1.5,2); np <- length(q10)
kT <- matrix(nrow=nT,ncol=np)
for(i in 1:np){
kT[,i] <- kT.rate(Temp,k0,q10[i])/kT.rate(Temp[1],k0,q10[i])
}
# Q10
Temp <- seq(0,50); nT <- length(Temp)
Topt= 25; Tmax <- 40; q10 <- 1.5
fT <- Q10(Temp, param=c(Topt,Tmax,q10))
Temp <- seq(0,50); nT <- length(Temp)
# param Topt, Tmax,q10
param.nom <- c(25,40,1.5)
Topt.sens <- c(20,25,30)
Tmax.sens <- c(40,45,50)
q10.sens <- c(1,1.5,2)
np <- length(q10)
test <- Q10.mruns(Temp,param.nom,Topt.sens,Tmax.sens,q10.sens)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.