f) Making Runge and Mara Figure 28.6

Load the package

library(redstart)

Figure 28.6

F28.6_range_winter <- param_ranges(figure = 28.6,scenario = "winter.limited")
#> Setting breeding 'source' carrying capacity K.bc to 800
#> Setting winter   'good'   carrying capacity K.wg to 484
#> Setting carry over effect to vary from 1 to 2
F28.6_range_intermediate <- param_ranges(figure = 28.6,scenario = "intermediate")
#> Setting breeding 'source' carrying capacity K.bc to 224
#> Setting winter   'good'   carrying capacity K.wg to 580
#> Setting carry over effect to vary from 1 to 2
F28.6_range_summer <- param_ranges(figure = 28.6,scenario = "summer.limited")
#> Setting breeding 'source' carrying capacity K.bc to 205
#> Setting winter   'good'   carrying capacity K.wg to 900
#> Setting carry over effect to vary from 1 to 2

head(F28.6_range_winter)
#>            min     max
#> gamma      5.0     5.0
#> co.        1.0     2.0
#> K.bc     800.0   800.0
#> K.bk   10000.0 10000.0
#> K.wg     485.0   485.0
#> S.w.mg     0.8     0.8

F28.6_seq_winter <- param_seqs(F28.6_range_winter)
F28.6_seq_intermed <- param_seqs(F28.6_range_intermediate)
F28.6_seq_summer <- param_seqs(F28.6_range_summer)

F28.6_grid_winter <- param_grid(param.seqs = F28.6_seq_winter)
#> The dimension of the fully expanded dataframe is: 
#> 10 by 30
F28.6_grid_intermed <- param_grid(param.seqs = F28.6_seq_intermed)
#> The dimension of the fully expanded dataframe is: 
#> 10 by 30
F28.6_grid_summer <- param_grid(param.seqs = F28.6_seq_summer)
#> The dimension of the fully expanded dataframe is: 
#> 10 by 30

F28.6_winter <-runFAC_multi(param.grid = F28.6_grid_winter, 
                     verbose = F)
#> 
#> Model at equilibrium after 149 iterations
#> 
#> Model at equilibrium after 110 iterations
#> 
#> Model at equilibrium after 82 iterations
#> 
#> Model at equilibrium after 79 iterations
#> 
#> Model at equilibrium after 78 iterations
#> 
#> Model at equilibrium after 77 iterations
#> 
#> Model at equilibrium after 76 iterations
#> 
#> Model at equilibrium after 75 iterations
#> 
#> Model at equilibrium after 75 iterations
#> 
#> Model at equilibrium after 74 iterations
F28.6_intermed <-runFAC_multi(param.grid = F28.6_grid_intermed, 
                     verbose = F)
#> 
#> Model at equilibrium after 83 iterations
#> 
#> Model at equilibrium after 80 iterations
#> 
#> Model at equilibrium after 78 iterations
#> 
#> Model at equilibrium after 76 iterations
#> 
#> Model at equilibrium after 75 iterations
#> 
#> Model at equilibrium after 74 iterations
#> 
#> Model at equilibrium after 73 iterations
#> 
#> Model at equilibrium after 73 iterations
#> 
#> Model at equilibrium after 73 iterations
#> 
#> Model at equilibrium after 71 iterations
F28.6_summer <-runFAC_multi(param.grid = F28.6_grid_summer, 
                     verbose = F)
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations
#> 
#> Model at equilibrium after 94 iterations

plot_Fig28_6(winter.lim   = F28.6_winter$multiFAC.out.df.RM,
             summer.lim   = F28.6_intermed$multiFAC.out.df.RM,
             intermediate = F28.6_summer$multiFAC.out.df.RM)

Nathan L. Brouwer

2022-03-07

Figure 28.6