In ian-hastings/rotations: Simulating the Evolution of Insecticide Resistance under Rotations

library(knitr)
opts_chunk$set(fig.height=7, fig.width=7)
# opts_chunk$set(dev="tiff",
#                dev.args=list(compression="lzw"),
#                dpi=300,
#                cache=FALSE,
#                fig.path='pap2figs/')

  library(rotations)
  library(ggplot2)

Current state of play

1. Flexible number of insecticides between 1 & 5
2. Exposure can be set for each insecticide or same for all.
3. Starting frequency can be set for each insecticide or same for all.

Some example plots follow, these act as a visual test that the model is doing what we expect.

#expo_hi=c(1,0.1)
dfr <- run_rot(n_insecticides=2, rot_interval=100, expo_hi=c(1,0.4))

# Interval 50, migration 0.01

dfr <- run_rot(rot_interval=50, max_gen = 300, migration=0.01)

# Decrease migration to 0.001

dfr <- run_rot(rot_interval=50, max_gen = 300, migration=0.001)

### Rotate when resistant, migration 0.01

dfr <- run_rot(rot_interval=0, max_gen = 300, expo_hi=0.5, migration=0.01)

run_rot( n_insecticides = 3 , max_gen = 150 , start_freqs = 0.1 , rot_interval = 50 , eff = 0.6 , dom_sel = 0.5 , rr = 0.5 , expo_hi = 0.5 , coverage = 0.8 , migration = 0.01 , cost = 0.1, no_r_below_start =FALSE )

```r

run_rot( n_insecticides = 3 , max_gen = 150 , start_freqs = 0.1 , rot_interval = 50 , eff = 0.6 , dom_sel = 0.5 , rr = 0.5 , expo_hi = 0.5 , coverage = 0.8 , migration = 0.01 , cost = 0 )

ian-hastings/rotations documentation built on Dec. 14, 2020, 11:42 p.m.