R/calculate_m.R

In lorecatta/DENVclimate:

calculate_m <- function(R0_values, k, beta_mh, beta_hm, mu, eip, parms) {

  ec <- parms$ec

  T_inf <- parms$T_inf

  R0_values * mu / ((k^2 * beta_hm * beta_mh * T_inf * exp(-mu*eip)) + ec)

}