abkineticsFunction.Rd
This function is defined by the user and describes antibody kinetics given the model.
A function that returns estimate titre value
Add information here.
# Example usage: we define two functional form examples, one describing waning until infection, and then the other describing kinetics post infection.
noInfSerumKinetics <- function(titre_est, timeSince, pars) {
titre_est <- titre_est - pars[1] * (timeSince)
titre_est
}
infSerumKinetics <- function(titre_est, timeSince, pars) {
a <- pars[1]
b <- pars[2]
c <- pars[3]
if (timeSince < 14) {
titre_est <- titre_est + (log(exp(a) + exp(c)) * (timeSince) / 14) * 1;
} else {
titre_est <- titre_est + (log(exp(a) * exp(-b/10 * (timeSince - 14))) + exp(c)) * 1;
}
titre_est
}