readsdr (version 0.3.0)

sd_what_if_from_time: What if from time t we change the value of some parameters

Description

What if from time t we change the value of some parameters

Usage

sd_what_if_from_time(
  time,
  up_to_time = Inf,
  par_list,
  ds_inputs,
  start_time = NULL,
  stop_time = NULL,
  timestep = NULL,
  integ_method = "euler"
)

Value

A data frame

Arguments

time

Time at which the parameter values change

up_to_time

Time from which the original values are restored.

par_list

A list that indicates which parameters change from time t. For instance, if you wanted to change the value of parameter c to 4, you would provide the list(c = 4)

ds_inputs

A list of deSolve inputs generated by read_xmile

start_time

A number indicating the time at which the simulation begins.

stop_time

A number indicating the time at which the simulation ends.

timestep

A number indicating the time interval for the simulation. Also known as dt.

integ_method

A string indicating the integration method. It can be either "euler" or "rk4"

Examples

  filepath       <- system.file("models/", "SIR.stmx", package = "readsdr")
  mdl            <- read_xmile(filepath)
  ds_components  <- mdl$deSolve_components
  output         <- sd_what_if_from_time(3, Inf, list(c = 4), ds_components)