#' calculate_r0
#'
#' @param th_in Binary indicator variable if model is age structured (1) or not (0) between children and adults. Defaults to NULL in which case no sampling happens
#' @param sus_c Suceptible child human hosts
#' @param sus_a Suceptible adult human hosts
#' @param sm_c Suceptible child moisquitoes
#' @param sm_a Suceptible adult moisquitoes
#' @param b_vary Defaults to 1
#' @export
calculate_r0_adam <- function(th_in,sus_c=1,sus_a=1,sm_c=1,sm_a=1,b_vary=1,control=1){
# Rate humans get infected -- FORMULATION WITH ONE MOSQUITO POP AND DIFFERENT BITING RATES
b_hv = b_vary * th_in$beta_h
b_hh = 0
# Rate vectors get infected
b_vh = th_in$beta_v * b_hv
b_vv = 0
rr_hh=rep(0,length(b_vary));
rr_vv = rr_hh;
delta_v <- th_in$MuV
exp_v <- th_in$Vex
exp_h <- th_in$Exp
inf_p <- th_in$Inf.
tau <- th_in$tau
m <- th_in$m
rr_hh=rep(0,length(b_vary));
rr_vv = rr_hh;
rr_hv = (sus_c*b_hv/delta_v)*(exp_v/(delta_v+exp_v))
rr_vc = (sm_c*b_vh/inf_p)
r0_hh=rep(0,length(b_vary));
r0_vv = rr_hh;
r0_hv = (b_hv/delta_v)*(exp_v/(delta_v+exp_v))
r0_vc = (b_vh/inf_p)
rr_post = NULL
r0_post = NULL
for(ii in 1:length(b_vary)){
rr_post=c(rr_post,( max(Re(eigen(matrix(c(rr_hh[ii],rr_hv[ii],rr_vc[ii],rr_vv[ii]),nrow=2))$values))) )
r0_post=c(r0_post,( max(Re(eigen(matrix(c(r0_hh[ii],r0_hv[ii],r0_vc[ii],r0_vv[ii]),nrow=2))$values))) )
}
return( list(r0_out=r0_post, rr_out=rr_post,rr_mat = matrix(c(rr_hh[1],rr_hv[1],rr_vc[1],rr_vv[1]),nrow=2,byrow=T)) )
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.