# NOT RUN {
# }
# NOT RUN {
ls <- landscape(population = egk_pop, suitability = egk_hab, carrying_capacity = egk_k)
pd <- population_dynamics(change = growth(egk_mat),
dispersal = kernel_dispersal(max_distance = 2000,
dispersal_kernel = exponential_dispersal_kernel(
distance_decay = 1000)),
density_dependence = ceiling_density())
sim <- simulation(landscape = ls,
population_dynamics = pd,
habitat_dynamics = NULL,
timesteps = 20)
# Extract the population raster for the second life-stage in the first
# replicate and ninth timestep
extract_spatial(sim, timestep = 9, stage = 2)
# }
Run the code above in your browser using DataLab