eode_sensitivity_proj {ecode}

Sensitivity Analysis

Description

Run a sensitivity analysis on an ODE system.

Usage

eode_sensitivity_proj(x, valueSpace, N, step = 0.01)

Arguments

x

object of class "pc" that represents the ODE system under consideration.

valueSpace

a list indicating initial conditions and parameters. Model variables must be included to specify initial values of each variable. Values can be a vector, indicating all the potential values to be considered in the sensitivity analysis.

N

number of iterations

step

interval of time for running simulations. Parameter of the function "eode_proj()".

Value

an object of "pcfamily" class, each component having three sub-components: $grid_var: variables or parameters whose values are changed throughout the sensitivity analysis. $fixed_var: variables whose values are not changed. $pc: phase curve. An object of "pc" class.

Examples

eq1 <- function(x, y, r1 = 4, a11 = 1, a12 = 2) (r1 - a11 * x - a12 * y) * x
eq2 <- function(x, y, r2 = 1, a21 = 2, a22 = 1) (r2 - a21 * x - a22 * y) * y
x <- eode(dxdt = eq1, dydt = eq2, constraint = c("x<100", "y<100"))
eode_sensitivity_proj(x, valueSpace = list(x = c(0.2, 0.3, 0.4), y = 0.1, a11 = 1:3), N = 100)