SetRanges: Parameter ranges for global sensitivity analysis

Description

Define the minimum and maximum values for parameters whose global sensitivities are to be assesses with CalculateGlobalSens or sensRange functions.

Usage

SetRanges(pars = NULL, range = 0.1)

Arguments

pars

the same pars vector used in one of the following functions: SolveSI or SolveIASA.

range

scale factor to define the minimum and maximum for each parameter. The default is 0.1, which set the minimum and maximum as 10 percent lesser and greater than the pars values.

Value

data.frame with the complete set of parameter ranges.

References

Soetaert K and Petzoldt T (2010). Inverse modelling, sensitivity and monte carlo analysis in R using package FME. Journal of Statistical Software, 33(3), pp. 1-28.

Reichert P and Kfinsch HR (2001). Practical identifiability analysis of large environmental simulation models. Water Resources Research, 37(4), pp. 1015-1030.

http://oswaldosantos.github.io/capm

Examples

Run this code

data(psu.ssu)
data(survey.data)

#####################
## SolveIASA model ##
#####################

# Parameters and initial conditions.
pars.solve.iasa <- c(
   b1 = 21870.897, b2 = 4374.179,
   df1 = 0.104, dm1 = 0.098, df2 = 0.1248, dm2 = 0.1176,
   sf1 = 0.069, sf2 = 0.05, sm1 = 0.028, sm2 = 0.05,
   k1 = 98050.49, k2 = 8055.456, h1 = 1, h2 = .5,
   ab = 0.054, ad = 0.1, v = 0.1)

# Set ranges 10 \% greater and lesser than the
# point estimates.
rg.solve.iasa <- SetRanges(pars.solve.iasa)