Error message from 'scanForHeterozygotes' function in AllelicImbalance package
1
0
Entering edit mode
JK ▴ 10
@jk-6972
Last seen 6.2 years ago
United States

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

 

https://www.broadinstitute.org/gatk/guide/article?id=3891

 

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,

Joe

allelicimbalance • 791 views
ADD COMMENT
0
Entering edit mode
@jesper-robert-gadin-6326
Last seen 6.3 years ago
European Union

Hello,

Sorry to be late with my reply!

This error was solved some months ago, see https://github.com/pappewaio/AllelicImbalance/issues/3

Jesper

ADD COMMENT
0
Entering edit mode

Yes, this issue is resolved. Thank you.

ADD REPLY

Login before adding your answer.

Traffic: 418 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6