sim1 <- sim_poisson_community(s_pool = 200, n_sim = 20000, sad_type = "lnorm",
sad_coef = list("cv_abund" = 2))
sad1 <- community_to_sad(sim1)
plot(sad1, method = "rank")
plot(sad1, method = "octave")
Run the code above in your browser using DataLab