APPENDIX 1 We find the intrinsic growth for species 1 and species 2 simultaneously
intrinsic_growth <- pop_data_1_ for (col in c(“species_1”, “species_2”)) { intrinsic_growth[[col]] <- c( pop_data_1_[[col]][2:nrow(pop_data_1_)] / pop_data_1_[[col]][1:(nrow(pop_data_1_) - 1)] ) }