model_trajectory: Model Model trajectory calc

Description Usage Arguments Value See Also Examples

View source: R/model.R

Description

Calculates the ferret model trajectory for a single infection event. Uses an R implementation so easy to code and has named parameter vectors. An equivalent function using Rcpp is implemented at model_trajectory_cpp. Run parameter_descriptions for documentation of the model parameter vector

Usage

1
model_trajectory(pars, times, logSigma = FALSE)

Arguments

pars

the named vector of model parameters

times

the vector of time in days to solve over

logSigma

if TRUE, uses the exponential of pars["sigma"] (useful if MCMC is exploring over log space). Default FALSE.

Value

a vector of antibody titres

See Also

Other model functions: model_func_groups, model_func, model_trajectory_cpp

Examples

1
2
3
4
5
6
7
pars <- c("lower_bound"=0,"S"=1,"EA"=0,"MAX_TITRE"=13,
          "mu"=8,"tp"=12,"dp"=0.5,"ts"=10,"m"=0.003,"beta"=0.6, "c"=4,
          "sigma"=1,"y0_mod"=-20,"boost_limit"=-1,"tau"=0.05,
          "order"=1,"primed"=0,"mod"=1,
          "x"=0,"t_i"=10,"y0"=0,"eff_y0"=0)
times <- seq(0,100,by=10)
y <- model_trajectory(pars,times)

jameshay218/antibodyKinetics documentation built on Nov. 8, 2019, 11 a.m.