This version of the code was used for the second revision of the manuscript submitted to Plos One in January 2021.
