Question: RSVSim cannot simulate duplications
gravatar for bdelolmo
4.4 years ago by
bdelolmo10 wrote:


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, size, replace, prob) : 
cannot take a sample larger than the population when 'replace = FALSE'

Does any know why is this happening?

ADD COMMENTlink modified 4.2 years ago by choosehappy10 • written 4.4 years ago by bdelolmo10
gravatar for choosehappy
4.2 years ago by
choosehappy10 wrote:

I had the same issue with RSVsim version 1.10 (it was the default installed by bioconductor for R3.2). I installed RSVsim version 1.12 from source, and the issue was resolved.

ADD COMMENTlink written 4.2 years ago by choosehappy10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1418 users visited in the last hour