An R script to perform an "in silico" ddRADseq