R script of an MC simulation modeling the DGP of temporal frequency distributions