knitr::opts_chunk$set( collapse = TRUE, comment = "#>", eval = FALSE )
library(roleR) set.seed(1)
# set-up parameters p <- roleParams( individuals_local = 1000, individuals_meta = 10000, species_meta = 50, speciation_local = 0.5, speciation_meta = 1, extinction_meta = 0.8, trait_sigma = 2, env_sigma = 1, comp_sigma = 1, dispersal_prob = 0.5, mutation_rate = 0.01, equilib_escape = 1, num_basepairs = 250, init_type = 'oceanic_island', niter = 1000, niterTimestep = 100 ) # initialize the model m <- roleModel(p) # run the model m <- runRole(m) # get results getSumStats(m, list(rich = richness, hill = hillAbund))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.