# NOT RUN {
## IASA model
## Parameters and intial conditions.
data(dogs)
dogs_iasa <- GetDataIASA(dogs,
destination.label = "Pinhais",
total.estimate = 50444)
# Solve for point estimates.
solve_iasa_pt <- SolveIASA(pars = dogs_iasa$pars,
init = dogs_iasa$init,
time = 0:15,
alpha.owned = TRUE,
method = 'rk4')
solve_iasa_rg <- SolveIASA(pars = dogs_iasa$pars,
init = dogs_iasa$init,
time = 0:10,
alpha.owned = TRUE,
s.range = seq(0, .4, l = 15),
a.range = c(0, .2),
alpha.range = c(0, .05),
v.range = c(0, .1),
method = 'rk4')
## Plot unowned population sizes using point estimates
# }
# NOT RUN {
PlotModels(solve_iasa_pt, variable = "ns2")
## Plot all scenarios and change the label for the scenarios.
## Not run
PlotModels(solve_iasa_rg, variable = "ns")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab