compare_emp: Compare minimum expected populations

Description Usage Arguments Examples

Description

Compare minimum expected populations from two or more 'simulation_results' objects.

Usage

1
compare_emp(x, ..., interval = 95, all_points = FALSE, simulation_names = NULL)

Arguments

x

a simulation_results object

...

additional simulation results objects

interval

the desired confidence interval representing the uncertainty around the expected minimum population estimates from simulation comparisons; expressed as a whole integer between 0 and 100 (default value is 95).

all_points

should the expected minimum populations from all simulation replicates be shown on the plot?

simulation_names

an optional character vector of simulation names to override the defaults

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
## Not run: 
ls <- landscape(population = egk_pop, suitability = egk_hab, carrying_capacity = egk_k)

# Create populations dynamics with and without ceiling density dependence
pd1 <- population_dynamics(change = growth(egk_mat),
                           dispersal = kernel_dispersal(max_distance = 2000,
                           dispersal_kernel = exponential_dispersal_kernel(distance_decay = 1000)),
                           density_dependence = ceiling_density())
pd2 <- population_dynamics(change = growth(egk_mat),
                           dispersal = kernel_dispersal(max_distance = 2000,
                           dispersal_kernel = exponential_dispersal_kernel(distance_decay = 1000)))

# Run first simulation with ceiling density dependence and three replicates
sim1 <- simulation(landscape = ls,
                   population_dynamics = pd1,
                   habitat_dynamics = NULL,
                   timesteps = 20,
                   replicates = 3)
                   
# Run second simulation without ceiling density dependence and three replicates
sim2 <- simulation(landscape = ls,
                   population_dynamics = pd2,
                   habitat_dynamics = NULL,
                   timesteps = 20,
                   replicates = 3)

compare_emp(sim1, sim2)

## End(Not run)

steps documentation built on March 3, 2021, 5:06 p.m.