simulates an ode model with extra work — gsl_odeiv2

This function calls a C function which solves an initial value problem, calculates the sensitivity of the solution, log-likelihood value ll, gradient of ll amd Fisher-Information.

Usage

gsl_odeiv2_fi(
  name,
  experiments,
  p,
  abs.tol = 1e-06,
  rel.tol = 1e-05,
  initial.step.size = 0.001,
  method = 0
)

Arguments

name: the name of the ODE model to simulate (a shared library of the same name will be dynamically loaded and needs to be created first)
experiments: a list of N simulation experiments (time, parameters, initial value, events)
p: a matrix of parameters with M columns
abs.tol: absolute tolerance, real scalar
rel.tol: relative tolerance, real scalar
initial.step.size: initial value for the step size; the step size will adapt to a value that observes the tolerances, real scalar

Value

a list of the solution trajectories y(t;p) for all experiments (named like the experiments), as well as the output functions