Map N(0,1) stochastic perturbations to an LNA path.

Map N(0,1) stochastic perturbations to an LNA path.

map_draws_2_lna(
  pathmat,
  draws,
  lna_times,
  lna_pars,
  lna_param_vec,
  lna_param_inds,
  lna_tcovar_inds,
  init_start,
  param_update_inds,
  stoich_matrix,
  forcing_inds,
  forcing_tcov_inds,
  forcings_out,
  forcing_transfers,
  svd_d,
  svd_U,
  svd_V,
  step_size,
  lna_pointer,
  set_pars_pointer
)

Arguments

pathmat

matrix where the LNA path should be stored

draws

matrix of N(0,1) draws to be mapped to an LNA path

lna_times

vector of interval endpoint times

lna_pars

numeric matrix of parameters, constants, and time-varying covariates at each of the lna_times

init_start

index in the parameter vector where the initial compartment volumes start

param_update_inds

logical vector indicating at which of the times the LNA parameters need to be updated.

stoich_matrix

stoichiometry matrix giving the changes to compartments from each reaction

forcing_inds

logical vector of indicating at which times in the time-varying covariance matrix a forcing is applied.

svd_d

vector in which to store SVD singular values

svd_U

matrix in which to store the U matrix of the SVD

svd_V

matrix in which to store the V matrix of the SVD

step_size

initial step size for the ODE solver (adapted internally, but too large of an initial step can lead to failure in stiff systems).

lna_pointer

external pointer to LNA integration function.

set_pars_pointer

external pointer to the function for setting the LNA parameters.

forcing_matrix

matrix containing the forcings.

Value

fill out pathmat with the LNA path corresponding to the stochastic perturbations.