In AndySouth/resistanceGame: prototyping an insecticide resistance management computer game

This document demonstrates seasonal scenarios for vector populations and resistance to be used in the game. It follows on from IRM-prototype-game-scenarios4.

The game will modify input parameters to generate reasonable scenarios. The input parameters are simply a means to generate reasonable scenarios.

In the following plots time in days is represented on the x axis, the top panel shows insecticide use, the middle panel shows vector population and the lower panel shows frequency of resistance.

The code included is there merely to show us as developers how the scenarios were generated.

For an interactive version of the equations used to generate these plots see https://andysouth.shinyapps.io/shinyGame4.

#### load required packages
require(resistanceGame)

Remember that years are 336 days (7days * 4weeks * 12months), half=168, 2 years=672.

\newpage

1 constant emergence, no intervention, 2 years

emergence <- 0.3
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

2 constant emergence, intervention, no resistance

emergence <- 0.3
l_config <- read_config()
l_config2 <- config_plan(l_config, t_strt=c(168), t_stop=c(672),
                         control_id=c('irs_ops'))
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  l_config=l_config2,
                  insecticide_kill=0.5,
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

3 constant emergence, intervention, resistance metabolic

emergence <- 0.3
l_config <- read_config()
l_config2 <- config_plan(l_config, t_strt=c(168), t_stop=c(672),
                         control_id=c('irs_pyr'))
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  l_config=l_config2,
                  insecticide_kill=0.6, resist_freq_start = 0.1, resist_mech='metabolic',
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

4 constant emergence, intervention, resistance target

emergence <- 0.3
l_config <- read_config()
l_config2 <- config_plan(l_config, t_strt=c(168), t_stop=c(672),
                         control_id=c('irs_pyr'))
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  l_config=l_config2,
                  insecticide_kill=0.6, resist_freq_start = 0.1, resist_mech='target',
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

5 seasonal emergence, no intervention, 2 years

emergence <- expand_season(season_string="6:0.3;6:0.1;6:0.3;6:0.1")
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

6 seasonal emergence, intervention, no resistance

emergence <- expand_season(season_string="6:0.3;6:0.1;6:0.3;6:0.1")
l_config <- read_config()
l_config2 <- config_plan(l_config, t_strt=c(168), t_stop=c(672),
                         control_id=c('irs_ops'))
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  l_config=l_config2,
                  insecticide_kill=0.5,
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

7 seasonal emergence, intervention, resistance metabolic

emergence <- expand_season(season_string="6:0.3;6:0.1;6:0.3;6:0.1")
l_config <- read_config()
l_config2 <- config_plan(l_config, t_strt=c(168), t_stop=c(672),
                         control_id=c('irs_pyr'))
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  l_config=l_config2,
                  insecticide_kill=0.6, resist_freq_start = 0.1, resist_mech='metabolic',
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

8 seasonal emergence, intervention, resistance target

emergence <- expand_season(season_string="6:0.3;6:0.1;6:0.3;6:0.1")
l_config <- read_config()
l_config2 <- config_plan(l_config, t_strt=c(168), t_stop=c(672),
                         control_id=c('irs_pyr'))
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  l_config=l_config2,
                  insecticide_kill=0.6, resist_freq_start = 0.1, resist_mech='target',
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

9 seasonal emergence, intervention, resistance metabolic, change intervention to susceptible

emergence <- expand_season(season_string="6:0.3;6:0.1;6:0.3;6:0.1")
l_config <- read_config()
l_config2 <- config_plan(l_config, t_strt=c(1,169), t_stop=c(168,336),
                         control_id=c('irs_pyr','irs_ops'))
plot_sim( run_sim(num_tsteps=672, emergence=emergence, survival=0.7,
                  l_config=l_config2,
                  insecticide_kill=0.6, resist_freq_start = 0.1, resist_mech='metabolic',
                  resist_incr=0.02, resist_decr = 0.01), 
                  plot_emergence=TRUE )

\newpage

AndySouth/resistanceGame documentation built on May 5, 2019, 6:01 a.m.