Temp.functions: Functions to define temperature response.

Description Usage Arguments Details Value Note Author(s) References See Also Examples

Description

Functions parab, kt.rate, Q10 are for single temperature response. Function Q10.mruns explores sensitivity to parameters of Q10.

Usage

1
2
3
4
parab(Temp,a,b)
kT.rate(Temp,k0,q10)
Q10(Temp, param)
Q10.mruns(Temp, param.nom, Topt.sens, Tmax.sens, q10.sens)

Arguments

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

Details

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.

Value

Temperature response

Note

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.

Author(s)

Miguel F. Acevedo Acevedo@unt.edu

References

Acevedo M.F. 2012. Simulation of Ecological and Environmental Models. CRC Press.

See Also

Simulation functions sim, simd

Examples

 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)

seem documentation built on April 14, 2017, 9:12 p.m.