Question: Error message from 'scanForHeterozygotes' function in AllelicImbalance package
gravatar for JK
3.1 years ago by
United States
JK10 wrote:

I generated BAM files using the STAR RNA-Seq aligner and the work stream described in GATK:


My R command with the error I am running into with these BAM files:


> heterozygotePositions <- scanForHeterozygotes(reads, maxReads=1e+07, verbose=T)

Investigating chromosome chr4

Investigating sample 1 out of 33

Error in new[start[i]:(start[i] + length(toKeep[[i]]) - 1), i] <- charList[[i]][toKeep[[i]]] :

  number of items to replace is not a multiple of replacement length


I would appreciate some help with this error.

Thank you,


ADD COMMENTlink modified 2.8 years ago by Jesper Robert Gadin20 • written 3.1 years ago by JK10
gravatar for Jesper Robert Gadin
2.8 years ago by
European Union
Jesper Robert Gadin20 wrote:


Sorry to be late with my reply!

This error was solved some months ago, see


ADD COMMENTlink written 2.8 years ago by Jesper Robert Gadin20

Yes, this issue is resolved. Thank you.

ADD REPLYlink written 2.8 years ago by JK10
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.2.0
Traffic: 356 users visited in the last hour