To account for the stochastic nature of the processes regulating both dynamics of mosquito and epidemic transmission, we used a discrete-time stochastic version of model M, with time-step where l is the force of infection (see for instance [38] [39] [40] All model parameters were calibrated to minimize the score function F~ffi