1- Change ref pop size min = 0 by min = 30
Thus modified
  N_ref[N_ref == 0] <- NaN
   by
  N_ref[N_ref < 30] <- NaN

2- added warning message (if too many iterations with N_ref < 30)

3- updated description

Updated plot_impact accordingly

Added imports (magrittr, popbio) to models M3 and M4 functions

Added piece "%>% sapply(min, 0.999)" in M3 and M4 to avoid survival parameter > 1.
Also added this in 'init_calib' function.

All models tested and work.

Made Model 3 (with DD) simulations faster by changing the calibrate_params function (based on optim) by an analytical approximation (using match_lam_delta function).

Tests reveals that the approx is largely good enough. Makes no differences in the results.

It works well but it's slow because of the optim procedure : lam_Nt > calibrate vr > Leslie matrix as f(Nt)

updated variables names in the run_analysis file

changed name survival_mean > survivals
changed name fecundities_mean > fecundities

devtools : added package on Git