I am trying to simulate some SVs using RSVSim. I can simulate both deletions and insertions at the same time, but the problem comes when I want to introduce duplications.
This is the code that I use:
sim = simulateSV(output='.', genome=genome, chrs = "chr12", ins = 150, sizeIns = sizeins, dels = 150, sizeDels = sizedels, dups = 150, sizeDups = sizedups)
The genome variable is the plain text ucsc hg 19 version. Size params are vectors of 150 random lengths between 20 to 200 bp.
The error that it throws is:
Loading genome from FASTA file: /home/bernat/reference/ucsc.hg19.fasta Breakpoints will be distributed uniformly across the genome Error in sample.int(length(x), size, replace, prob) : cannot take a sample larger than the population when 'replace = FALSE'
Does any know why is this happening?