model.pasin {REMixed}R Documentation

Model from Pasin and al.,2019

Description

Generate trajectory of the Humoral Immune Response to a Prime-Boost Ebola Vaccine.

Usage

model.pasin(t, y, parms)

Arguments

t

vector of time ;

y

initial condition, named vector of form c(Ab=<...>,S=<...>,L=<...>) ;

parms

named vector of model parameter ; should contain "theta_S","theta_L","delta_Ab","delta_S","delta_L".

Details

The model correspond to the dynamics of the humoral response, from 7 days after the boost immunization with antibodies secreting cells -S and L, characterized by their half lives- that produces antibodies -AB- at rate \theta_S and \theta_L. All these biological entities decay at rate repectively \delta_S, \delta_L and \delta_{Ab}. Model is then defined as

\left\{\begin{matrix}\frac{d}{dt} Ab(t) &=& \theta_S S(t) + \theta_L L(t) - \delta_{Ab} Ab(t) \\ \frac{d}{dt}S(t) &=& -\delta_S S(t) \\ \frac{d}{dt} L(t) &=& -\delta_L L(t)\end{matrix}\right.

Value

Matrix of time and observation of antibody titer Ab, and ASCs S and L.

References

Pasin C, Balelli I, Van Effelterre T, Bockstal V, Solforosi L, Prague M, Douoguih M, ThiƩbaut R, for the EBOVAC1 Consortium. 2019. Dynamics of the humoral immune response to a prime-boost Ebola vaccine: quantification and sources of variation. J Virol 93: e00579-19. https://doi.org/10.1128/JVI.00579-19

See Also

indParm, model.clairon.

Examples

y = c(Ab=0,S=5,L=5)
parms = c(theta_S = 611,
          theta_L = 3.5,
          delta_Ab = 0.025,
          delta_S = 0.231,
          delta_L = 0.000152)

t = seq(0,100,5)
res <- model.pasin(t,y,parms)
plot(res)


[Package REMixed version 0.1.0 Index]