library(ATNr)
set.seed(123)
masses <- runif(20, 10, 100) #body mass of species
L <- create_Lmatrix(masses, 10, Ropt = 10)
L[L > 0] <- 1
mod <- create_model_Unscaled(20, 10, masses, L)
mod <- initialise_default_Unscaled(mod)
Run the code above in your browser using DataLab